- Open Access
Differential gene methylation and expression of HOX transcription factor family in orbitofacial neurofibroma
Acta Neuropathologica Communications volume 8, Article number: 62 (2020)
Although most commonly benign, neurofibromas (NFs) can have devastating functional and cosmetic effects in addition to the possibility of malignant transformation. In orbitofacial neurofibromatosis type 1, NFs may cause progressive, disfiguring tumors of the lid, brow, temple, face and orbit. The purpose of this study was to identify biological differences between orbitofacial NFs and those occurring at other anatomic sites. We used Illumina Methylation EPIC BeadChip to study DNA methylation differences between orbitofacial NFs (N = 20) and NFs at other sites (N = 4). Global methylation differences were detected between the two groups and the top differentially methylated genes were part of the HOX (Homebox) family of transcription factors (HOXC8, HOXC4, HOXC6, HOXA6 and HOXD4), which were hypomethylated in orbitofacial NFs compared to the non-orbital NFs. Conversely, LTF (lactoferrin) was relatively hypermethylated in orbitofacial NF compared to non-orbitofacial NF. HOXC8 protein levels were higher in orbitofacial plexiform NFs (p = 0.04). We found no significant differences in the expression of HOXC4, HOXA6, or HOXD4 between the two groups. HOXC8 mRNA levels were also higher in orbitofacial NFs and HOXC8 overexpression in a non-neoplastic human Schwann cell line resulted in increased growth. In summary, we identified gene methylation and expression differences between orbitofacial NF and NFs occurring at other locations. Further investigation may be warranted, given that the HOX family of genes play an important role during development, are dysregulated in a variety of cancers, and may provide novel insights into therapeutic approaches.
Neurofibromatosis type 1 (NF1), also known as von Recklinghausen disease, is a relatively frequent (1 in 3000 live births), autosomal dominant, neurocutaneous disorder with heterogeneous clinicopathologic manifestations . The most prevalent tumor in NF1 patients is neurofibroma (NF), which is typically slow growing and benign, but can have devastating functional and cosmetic effects. Malignant transformation develops in 5–10% of NF1 patients, typically in a plexiform neurofibromas subtype .
Orbitofacial neurofibromatosis type 1 (OFNF) has been recognized as a unique clinical variant of NF1 for many years . In OFNF, which occurs in 1–22% of patients, orbitofacial NFs may cause progressive, disfiguring tumors of the lid, brow, temple, face and orbit [11, 33, 48]. Large NFs involving the orbit tend to be more aggressive and infiltrative compared to NFs located elsewhere in the body (Fig. 1) . The recurrence rate after excision is also high. In addition, this aggressiveness is most striking in infancy and early childhood and tends to improve somewhat as the individual ages [18, 23, 24].
These tumors not only affect facial appearance , but they can also disturb the growth patterns of the developing skull [19, 49], orbit [20, 43], and globe [15, 30]. Sphenoid and orbital dysplasia are prominent. Skull and facial features ipsilateral to neurofibromas include enlargement of the middle cranial fossa and defects in the greater sphenoid wing and lateral orbital wall. In addition, globe enlargement occurs on the affected side of patients with OFNF under certain circumstances [15, 30]. Both neuroimaging [19, 20, 49] and clinical reports  have now documented that NFs are frequently contiguous to globes that are abnormally large and to bony changes in the skull and orbit . However, the reasons that orbitofacial NFs are more aggressive than NFs at other locations, and the mechanisms by which they disturb developing skull and facial structures are currently unclear.
Orbitofacial NFs designation is usually applied to patients with often complex NFs that, given the infiltrative nature, involve the orbit and ocular adnexa, but also frequently adjacent facial structures. We previously reviewed our pathology experience with nerve sheath tumors involving the eye and ocular adnexa. Neurofibroma was the predominant nerve sheath tumor type involving this region. Of note, 63 neurofibromas at this location were studied, and they were predominantly plexiform, diffuse or with mixed patterns (68%) and 89% develop in patients with known NF1 . Although orbitofacial NFs often behave in an aggressive fashion, malignant transformation appears to be paradoxically rare. We only encountered a single case of malignant peripheral nerve sheath tumor (MPNST) in our review of 90 nerve sheath tumors from 67 patients involving the eye and ocular adnexa . This patient was a 63-year-old woman who underwent clinical progression of a childhood orbitofacial neurofibroma over several decades, with eventual malignant transformation, dissemination and death .
Recent advances in our understanding of human neoplasia have uncovered an important role for epigenetic changes in addition to DNA-level mutations and copy number changes. Furthermore, methylation profiling has emerged as a robust tool to aid in the classification of tumors of the nervous system , including those of peripheral nerve origin . DNA methylation profiles vary according to tumor subtypes. For example, MPNST, an aggressive malignant neoplasm that can develop in NF1 patients, is characterized by relative global hypomethylation . Another rationale to study epigenetic changes in subsets of NFs is that genetic alterations in these slow-growing tumors are essentially limited to NF1 mutations . We hypothesized that global methylation profiling would be helpful in studying biological differences between orbitofacial NFs and those developing in other body sites.
Materials and methods
Patients and samples
Formalin fixed, paraffin embedded tissues of 85 neurofibromas from 63 patients were studied by methylation array, qPCR and/or immunohistochemistry as described below. Only larger tumors with adequate material for DNA, RNA, and/or IHC analysis were selected, which excluded small localized tumors and enriched for larger more clinically significant tumors (more properly designated orbitofacial NF). Demographics and basic location characteristics were abstracted from retrospective review of clinical and pathologic records. Orbitofacial (n = 34) and non-orbitofacial (n = 51) tumors were studied. Pathologic subtypes included plexiform (n = 43), diffuse (n = 15), localized intraneural (n = 4) and localized cutaneous (n = 23). Clinicopathologic characteristics and tests performed per tumor are summarized in Supplementary table 1.
DNA methylation profiling
Formalin-fixed paraffin-embedded (FFPE) tissue from 20 orbitofacial and 4 non-orbital NFs were included for methylation profiling studies. In brief, FFPE sections were used to isolate genomic DNA with the QIAamp DNA FFPE Tissue Kit (Qiagen) and stored at 4 °C. Genomic DNA quality was assessed by low concentration agarose gel (0.6%) electrophoresis and fluorescent spectrometry with the PicoGreen DNA Kit (Life Technologies).
DNA bisulfite conversion was carried out using the EZ DNA Methylation Kit (Zymo Research) per manufacturer instructions with modifications for the Illumina Infinium Methylation Assay. Briefly, 400 ng of genomic DNA was first mixed with 5 μl of M-Dilution Buffer and incubated at 37 °C for 15 min, then mixed with 100 μl of CT Conversion Reagent prepared per the instruction manual. Mixtures were incubated in a thermocycler with 16 thermal cycles at 95 °C for 30 s and 50 °C for 1 h. Bisulfite-converted DNA samples were loaded onto 96-column plates provided in the kit for desulphonation and purification. Concentration of eluted DNA was measured using the Nanodrop-1000 spectrometer.
Bisulfite-converted DNA was analyzed using Illumina’s Infinium Human MethylationEPIC Beadchip Kit (WG-317-1002) per manufacturer instructions. The Beadchip contains over 850,000 CpG loci in human genome. Briefly, 4 μl of bisulfite-converted DNA was added to a 0.8 ml 96-well storage plate (Thermo Scientific), denatured in 0.014 N sodium hydroxide, and neutralized and amplified with kit-provided reagents and buffer at 37 °C for 20–24 h. Samples were fragmented using kit-provided reagents and buffer at 37 °C for 1 h and precipitated by adding 2-propanol. Re-suspended samples were denatured in a 96-well plate heat block at 95 °C for 20 min. 26 μl of each sample was loaded onto an 8-sample chip and the chips were assembled into hybridization chamber per manufacturer instructions. After incubation at 48 °C for 16–20 h, chips were briefly washed and then assembled and placed in a fluid flow-through station for primer-extension and staining procedures. Polymer-coated chips were image-processed in Illumina’s iScan scanner.
Array data acquisition and bioinformatic analysis
The quality of the raw data was assessed in an R environment. Briefly, the raw data was inputted and processed using the minfi R/Bioconductor  package to obtain the Beta- and M-values (log ratio of beta). Standard quality controls were performed (e.g. density plots and control probes expression profiles). Probes within 300 base pairs of each other were analyzed as a single cluster. Differentially methylated regions (DMRs) were detected by fitting a linear model with age as a covariate to adjust for differences between pediatric and adult patients. For each cluster using the bumphunter R/Bioconductor package  to compute the family-wide error rates (FWER) and p-values based on 500 bootstraps to estimate a null-profile. The top 5% regions with FWER ≤ 0.1 were selected as candidates for validation. Each DMR was annotated to its nearest gene based on the hg19 assembly.
RNA isolation and quantitative real time PCR analysis
The RNA was extracted using the RNeasy FFPE Kit (Qiagen). RNA concentration and purity were assessed using NanoDrop equipment (NanoDrop Technologies Inc.). One microgram of total RNA was used to produce cDNA with the super script III kit (INVITROGEN) and quantitative PCR was performed using SYBR Green reagents (Applied Biosystems, California) on an I-Cycler IQ Real-Time detection system (Bio-Rad, California) per manufacturer instructions. Primer sequences were as followed: huHOXC8 primers 5′-CTGTAAATCCTCCGCCAACAC-3′ (forward), 5′-CTTCAATCCGACGTTTTCGTG-3′ (reverse); huHOXC4 primers 5′-CCAGCAAGCAACCCATAGTC-3′ (forward), 5′-ATCCTTCTCCTTCGGGTCAG-3′ (reverse); huHOXA6 5′- CAAAGCACTCCATGACGAAGG-3′ (forward), 5′- CTCCTTCTCCAGCTCCAGTGT-3′ (reverse); GAPDH primers 5′- GCA GGG GGG AGC CAA AAG GGT-3′ (forward), 5′- TGG GTG GCA GTG ATG GCA TGG-3′ (reverse). The ΔΔCt method for quantitation of relative gene expression was used to determine the mean expression of each target gene normalized to the geometric mean of β-actin and GAPDH.
Tissue microarray construction and immunohistochemistry
A total of 77 NFs, including 44 orbitofacial and 31 non-periorbital, were examined by immunohistochemistry. Clinicopathologic features of the orbitofacial tumors were reported in a prior study . Tissue microarrays were constructed using four cores (0.6 mm diameter) from FFPE blocks. FFPE sections were deparaffinized in xylene, and antigen retrieval was performed thrice using sodium citrate buffer (pH 6.0) in a microwave oven at 5 min each. Sections were stained with rabbit polyclonal anti-HOXC4 (Abcam; 1:50 dilution, overnight at 4 °C), rabbit polyclonal anti-HOXC8 (Abcam; 1:50 dilution, overnight at 4 °C), rabbit polyclonal anti-HOXD4 (Novus Biologicals; 1:200, overnight at 4 °C); rabbit polyclonal anti-HOXA6 (Novus Biologicals; 1:200, overnight at 4 °C), followed by an anti-rabbit IgG secondary antibody staining with the VECTASTAIN® Elite® ABC HRP Kit (vectorlabs). After a diaminobenzidine reaction (Liquid DAB + Substrate Chromogen System, DakoCytomation), sections were counterstained with hematoxylin. All H&E slides were reviewed by one author (FJR), who scored the blinded samples. Immunoreactivity was scored semi-quantitatively based on the intensity of staining (0, no; 1, weak; 2, moderate; 3, intense). An overall immunohistochemical H-score (intensity x percentage of positive tumor cells) was assigned as low (H-score 0–100), medium (100–200), or high (200–300), as previously described . Proliferation index was evaluated as a visual estimate of the percentage of positive neoplastic nuclei using a Ki67 antibody (Ventana, clone 30–9).
Cell culture and transfection
A human Schwann cell line was obtained from ScienCell (Sciencell research laboratories) and maintained in complete Schwann Cell Medium (Sciencell research laboratories). All cell cultures were maintained at 37 °C in a humidified 5% CO2 atmosphere.
DNA constructs and viral infection
The construct encoding HOXC8 was purchased from Addgene (Addgene plasmid 21,001). The HOXC8 cDNA was subcloned into the retroviral vector pBABE_puro. Retroviral particles were produced by transfecting 293GT cells with VSVG envelope plasmid and the plasmid of interest, using Lipofectamine2000 (ThermoFisher) per manufacturer instructions. After 24 h, the medium was changed to human Schwann cell media, and supernatants were collected at 48 and 72 h. These supernatants were pooled, passed through a 0.45-mm filter, and then stored at − 80 °C until further needed. Human Schwann cells were incubated with retroviral particles in combination with polybrene for 24 h. After approximately 1 week in culture, puromycin (ThermoFisher) at 2 μg/ml was added to the cell culture medium to select for stable HOXC8-overexpressing cells. Control cells were infected with empty vector viral particles.
CellTiter-blue cell viability assay
Human Schwann cell line and human Schwann cell line transduced with HOXC8 were seeded in a 96-well plate (1000 cell/well). On the following day (day 0), day 4, and day 6, 40 μl of Cell titer Blue Reagent (Promega) was added to 200 μl culture medium for 4 h at 37 °C; fluorescence with excitation 560 nm and emission 590 nm was measured using the TECAN reader (TECAN).
Western blot analysis
Protein extraction, separation, and immunoblotting were performed as previously described . All Western blots were representative of three independent experiments. Antibodies used were as follows: rabbit polyclonal anti-HOXC8 (Abcam; 1:1000), β-ACTIN (Santa Cruz; 1:5000), and GAPDH (Santa Cruz; 1:5000). Anti-mouse IgG HRP-linked (#7076S) and anti-rabbit IgG HRP-linked (#7074S) (Cell Signaling Technology) were used as secondary antibodies. Overexpression lysates for recombinant human HOXC8 were used as positive controls (Origene; NM_022658).
HOX family of transcription factors are hypomethylated in orbitofacial neurofibromas
Differential methylation analysis identified 7 differential methylated regions between the orbitofacial (n = 20) and non-orbitofacial (n = 4) neurofibromas (FWER ≤0.1; P ≤ 0.0001) (Fig. 2 and Table 1). These DMRs were found predominantly within or next to the HOX family of transcription factors genes. HOXC8, HOXC4, HOXC6, HOXA6 and HOXD4 were found hypomethylated in the orbitofacial NF compared with the non-ocular (Fig. 3). Conversely, the transferrin family member LTF was relatively hypermethylated in orbitofacial NF compared to the non-orbital neurofibromas (Fig. 2). Gene methylation differences were further adjusted by age (Supplementary table 1).
HOXC8 protein levels are higher in orbitofacial compared to non-orbital neurofibromas
Immunohistochemistry of the 46 orbitofacial NFs and 31 non-orbital NFs showed no significant differences in protein expression between the two groups, although HOXC8 protein levels were higher in orbitofacial plexiform NFs (average H-score 126, range 0–180) compared to non-orbital plexiform NFs (average H score 100, range 10–160) (P = 0.04) (Fig. 4). We found no significant differences in protein expression between the orbitofacial NFs vs. non-orbital groups for HOXC4 (average H-score 170, range 80–190, vs. average H-score 172, range 160–190), HOXA6 (average H-score 155, range 90–190 vs. average H-score 158, range 90–190), or HOXD4 (average H-score 146, range 20–180 vs average H-score 155, range 50–190) (p > 0.05). Clinicopathologic features of the cohort are summarized in Supplementary Table 2. Next, we studied the ki67 proliferation index using immunohistochemistry. Proliferation index was slightly higher on average in the orbitofacial NFs compared with the non-orbital NFs (Supplementary Figure 1).
HOXC8 mRNA expression is increased in orbitofacial compared to non-orbital neurofibromas
Next we studied mRNA expression levels in orbitofacial NF (n = 13) samples compared to non-orbital neurofibromas (n = 7) using qRT-PCR. HOXC8 mRNA levels were ~ 1.3 fold higher in orbitofacial neurofibromas compared to non-orbitofacial neurofibromas (P = 0.0254) which was in line with the higher protein expression (Fig. 5). No significant differences were detected for mRNA levels in HOXC4 and HOXA6. Clinicopathologic features of the cases studied are summarized in Supplementary table 2.
HOXC8 increases human Schwann cell growth
We wanted to study the effect of overexpression of HOXC8 in a human Schwann cell line. For this we overexpressed HOXC8 through viral transduction. After selection, both HOXC8 human Schwann cell clones showed increased protein expression compared to the donor human Schwann cell line (Fig. 6a). To investigate if higher HOXC8 protein will have an effect on cell growth, we performed cell titer blue experiments over 6 days. Both clones of the HOXC8 overexpressing cell lines showed significant increase in cell viability compared with the donor human Schwann cell line (clone#1 P = 0.0029; clone#2 P = 0.0007) (Fig. 6b).
The importance of epigenetic changes has been increasingly recognized in a variety of physiologic and pathologic states, particularly neoplasms. In the current study, we identify important epigenetic differences between orbitofacial NF and those occurring at other sites, which may partially explain the increased aggressiveness of orbitofacial NFs. In particular, HOX family genes appear to be altered by differential methylation in NFs.
The HOX family encompass a group of genes arranged into four distinct clusters (HOXA, HOXB, HOXC, HOXD) in various chromosomal locations . They are known to play a critical role during development, encoding proteins that provide critical positional information during morphogenesis through variable protein composition, leading to regional patterning . The HOX genes are sequentially expressed during development, underlying the embryologic process known as the ‘temporal collinearity’ of the ‘HOX clock’ [9, 17].
Recent evidence has supported a role for HOX genes in physiologic and pathologic states involving stem cell biology, including cell hierarchy, positional identity, self-renewal and differentiation . Altered expression and mutations leading to genetic predisposition have been studied in a variety of cancer types, including gynecologic and hematologic malignancies . In our study, we found decreased DNA methylation in several HOX family genes in orbitofacial NFs. In particular, HOXC8 demonstrated increased mRNA and protein expression in orbitofacial NFs compared to NFs arising at other sites.
Altered HOXC8 levels have been documented in a variety of tumor types, and have been shown to be of biological and prognostic significance. For instance, high expression levels of HOXC8 and HOXC6 by immunohistochemistry were associated with poorer prognosis in a cohort of patients with esophageal squamous cell carcinoma , and increased gene expression was associated with worse outcome in cervical carcinoma . In another study, HOXC8 levels were increased in non-small cell lung cancer and was correlated with lymph node metastasis, differences that were thought to be mediated by TGFβ1 . HOXC8 has also been found to regulate proliferation, migration, and epithelial mesenchymal transition in triple negative breast cancer  and invasion in prostate cancer . However, a possible tumor suppressor role has also been reported in other studies, where overexpression of HOXC9 leads to decreased growth of nasopharyngeal carcinoma and stem cell properties in breast cancer cells . Also, HOXC8 expression has been shown to be inversely correlated with metastasis in pancreatic cancer .
The regulation of HOX genes during development and disease is complex, but epigenetic changes such as DNA methylation appear to be important mediators . For example, HOXC8 methylation seems to regulate cashmere fiber length in the cashmere goat . Variable HOX gene methylation has also been described in oral squamous cell carcinoma cell lines .
The results from our study suggest that HOX genes may play a role in the development of various NF subtypes. Of relevance, in a recent study Chen et al. demonstrated that inactivation of Nf1 in mice in a HOXB7 expressing neural crest lineage cell leads to the development of cutaneous and plexiform NFs in restricted anatomic distributions .
LTF (lactoferrin) is a member of the transferrin gene family. It is an iron binding protein which participates in the reduction of oxidative stress and inflammation . In our study it was the only significant gene to be differentially hypermethylated in orbitofacial NF. LTF has been found to be down-regulated in papillary thyroid cancer , gastric cancer , and underexpressed in association with hypermethylation in prostate cancer [36, 42] and nasopharyngeal carcinoma , supporting a role as a tumor suppressor gene in a variety of cancer types.
One of the caveats of our study is that clinical data, including precise age of onset, time of progression, extent of invasion as well as the precise NF1 mutations associated with their development were not studied. In addition, current efforts to more precisely characterize cutaneous neurofibromas, integrating clinical and pathological data are ongoing , and will have management implications. Orbitofacial neurofibromas are relatively rare, and larger cohorts will require multiinstitutional studies, with prospective clinical data that includes age of onset, precise clinicopathologic subtypes and NF1 mutations at a minimum, allowing to adjust for all these variables. Future studies into the molecular characterization of orbitofacial NF, and the relevance of HOX genes, will certainly include these variables, given their relevance to the study of NF1-associated neoplasia.
In summary, we report that there is differential methylation and expression of HOX genes in orbitofacial NFs when compared to NFs developing at other sites. Overexpression of HOXC8 leads to increased cell growth in non-neoplastic Schwann cells, suggesting a possible contribution to the increased aggressiveness in certain subtypes of NF. These findings warrant further investigation in clinical cohorts and animal models to continue the characterization of HOX genes and their role in the pathobiology of NFs, and develop rational therapies for tumors that although benign, are associated with disproportionate morbidity with current treatment regimens.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the GEO repository.
Adwan H, Zhivkova-Galunska M, Georges R, Eyol E, Kleeff J, Giese NA, Friess H, Bergmann F, Berger MR (2011) Expression of HOXC8 is inversely related to the progression and metastasis of pancreatic ductal adenocarcinoma. Br J Cancer 105:288–295. https://doi.org/10.1038/bjc.2011.217
Aryee MJ, Jaffe AE, Corrada-Bravo H, Ladd-Acosta C, Feinberg AP, Hansen KD, Irizarry RA (2014) Minfi: a flexible and comprehensive bioconductor package for the analysis of Infinium DNA methylation microarrays. Bioinformatics 30:1363–1369. https://doi.org/10.1093/bioinformatics/btu049
Axlund SD, Lambert JR, Nordeen SK (2010) HOXC8 inhibits androgen receptor signaling in human prostate cancer cells by inhibiting SRC-3 recruitment to direct androgen target genes. Mol Cancer Res 8:1643–1655. https://doi.org/10.1158/1541-7786.MCR-10-0111
Bai WL, Wang JJ, Yin RH, Dang YL, Wang ZY, Zhu YB, Cong YY, Deng L, Guo D, Wang SQet al (2017) Molecular characterization of HOXC8 gene and methylation status analysis of its exon 1 associated with the length of cashmere fiber in Liaoning cashmere goat. Genetica 145: 115–126 Doi https://doi.org/10.1007/s10709-017-9950-5
Capper D, Jones DTW, Sill M, Hovestadt V, Schrimpf D, Sturm D, Koelsche C, Sahm F, Chavez L, Reuss DE et al (2018) DNA methylation-based classification of central nervous system tumours. Nature 555:469–474. https://doi.org/10.1038/nature26000
Chaudhry IA, Morales J, Shamsi FA, Al-Rashed W, Elzaridi E, Arat YO, Jacquemin C, Oystreck DT, Bosley TM (2012) Orbitofacial neurofibromatosis: clinical characteristics and treatment outcome. Eye (Lond) 26:583–592. https://doi.org/10.1038/eye.2011.336
Chen Z, Mo J, Brosseau JP, Shipman T, Wang Y, Liao CP, Cooper JM, Allaway RJ, Gosline SJC, Guinney J et al (2019) Spatiotemporal loss of NF1 in Schwann cell lineage leads to different types of cutaneous Neurofibroma susceptible to modification by the hippo pathway. Cancer Discov 9:114–129. https://doi.org/10.1158/2159-8290.CD-18-0151
Deschamps J, Duboule D (2017) Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock. Genes Dev 31:1406–1416. https://doi.org/10.1101/gad.303123.117
Dolle P, Izpisua-Belmonte JC, Falkenstein H, Renucci A, Duboule D (1989) Coordinate expression of the murine Hox-5 complex homoeobox-containing genes during limb pattern formation. Nature 342:767–772. https://doi.org/10.1038/342767a0
Du YB, Dong B, Shen LY, Yan WP, Dai L, Xiong HC, Liang Z, Kang XZ, Qin B, Chen KN (2014) The survival predictive significance of HOXC6 and HOXC8 in esophageal squamous cell carcinoma. J Surg Res 188:442–450. https://doi.org/10.1016/j.jss.2014.01.017
Erb MH, Uzcategui N, See RF, Burnstine MA (2007) Orbitotemporal neurofibromatosis: classification and treatment. Orbit 26:223–228 Doi 788720091 [pii]. https://doi.org/10.1080/01676830600987227
Gong C, Zou J, Zhang M, Zhang J, Xu S, Zhu S, Yang M, Li D, Wang Y, Shi J et al (2019) Upregulation of MGP by HOXC8 promotes the proliferation, migration, and EMT processes of triple-negative breast cancer. Mol Carcinog 58:1863–1875. https://doi.org/10.1002/mc.23079
Grenier N, Guibert-Tranier F, Nicolau A, Caille JM (1984) Contribution of computerized tomography to the study of spheno-orbital dysplasia in neurofibromatosis. JNeuroradiol 11:201–211
Hosseinkhan N, Honardoost M, Blighe K, Moore CBT, Khamseh ME (2020) Comprehensive transcriptomic analysis of papillary thyroid cancer: potential biomarkers associated with tumor progression. J Endocrinol Invest. https://doi.org/10.1007/s40618-019-01175-7
Hoyt CS, Billson FA (1977) Buphthalmos in neurofibromatosis: is it an expression of regional giantism? JPediatrOphthalmol 14:228–234
Huang Y, Chen L, Guo A (2018) Upregulated expression of HOXC8 is associated with poor prognosis of cervical cancer. Oncol Lett 15:7291–7296. https://doi.org/10.3892/ol.2018.8200
Izpisua-Belmonte JC, Falkenstein H, Dolle P, Renucci A, Duboule D (1991) Murine genes related to the Drosophila AbdB homeotic genes are sequentially expressed during development of the posterior part of the body. EMBO J 10:2279–2289
Jackson IT, Carbonnel A, Potparic Z, Shaw K (1993) Orbitotemporal neurofibromatosis: classification and treatment. Plast Reconstr Surg 92:1–11
Jacquemin C, Bosley TM, Liu D, Svedberg H, Buhaliqa A (2002) Reassessment of sphenoid dysplasia associated with neurofibromatosis type 1. AJNR Am J Neuroradiol 23:644–648
Jacquemin C, Bosley TM, Svedberg H (2003) Orbit deformities in craniofacial neurofibromatosis type 1. AJNR Am J Neuroradiol 24:1678–1682
Jaffe AE, Murakami P, Lee H, Leek JT, Fallin MD, Feinberg AP, Irizarry RA (2012) Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. Int J Epidemiol 41:200–209. https://doi.org/10.1093/ije/dyr238
Kamkar F, Xaymardan M, Asli NS (2016) Hox-mediated spatial and temporal coding of stem cells in homeostasis and Neoplasia. Stem Cells Dev 25:1282–1289. https://doi.org/10.1089/scd.2015.0352
Lee V, Ragge NK, Collin JR (2004) Orbitotemporal neurofibromatosis. Clinical features and surgical management. Ophthalmology 111:382–388. https://doi.org/10.1016/j.ophtha.2003.05.029 S0161-6420(03)01193-X [pii]
Lee V, Ragge NK, Collin JR (2003) The surgical management of childhood orbito-temporal neurofibromatosis. Br J Plast Surg 56:380–387 Doi S0007122603001723 [pii]
Legrand D (2016) Overview of Lactoferrin as a natural immune modulator. J Pediatr 173(Suppl):S10–S15. https://doi.org/10.1016/j.jpeds.2016.02.071
Li B, Huang Q, Wei GH (2019) The role of HOX transcription factors in Cancer predisposition and progression. Cancers (Basel) 11. https://doi.org/10.3390/cancers11040528
Liu H, Zhang M, Xu S, Zhang J, Zou J, Yang C, Zhang Y, Gong C, Kai Y, Li Y (2018) HOXC8 promotes proliferation and migration through transcriptional up-regulation of TGFbeta1 in non-small cell lung cancer. Oncogenesis 7:1. https://doi.org/10.1038/s41389-017-0016-4
Luo G, Zhou Y, Yi W, Yi H (2015) Lactotransferrin expression is downregulated and affects the mitogen-activated protein kinase pathway in gastric cancer. Oncol Lett 9:2409–2413. https://doi.org/10.3892/ol.2015.3011
Macfarlane R, Levin AV, Weksberg R, Blaser S, Rutka JT (1995) Absence of the greater sphenoid wing in neurofibromatosis type 1: congenital or acquired: case report. Neurosurgery 37:129–133
Morales J, Chaudhry IA, Bosley TM (2009) Glaucoma and globe enlargement associated with neurofibromatosis type 1. Ophthalmology 116:1725–1730. https://doi.org/10.1016/j.ophtha.2009.06.019 S0161-6420(09)00647-2 [pii]
Natsumeda M, Maitani K, Liu Y, Miyahara H, Kaur H, Chu Q, Zhang H, Kahlert UD, Eberhart CG (2016) Targeting notch signaling and autophagy increases cytotoxicity in Glioblastoma Neurospheres. Brain Pathol 26:713–723. https://doi.org/10.1111/bpa.12343
Nix JS, Haffner MC, Ahsan S, Hicks J, De Marzo AM, Blakeley J, Raabe EH, Rodriguez FJ (2018) Malignant peripheral nerve sheath tumors show decreased global DNA methylation. J Neuropathol Exp Neurol 77:958–963. https://doi.org/10.1093/jnen/nly076
North K (1997) Neurofibromatosis type 1 in childhood. MacKeith Press, City
Ortonne N, Carroll S, Rodriguez FJ, Miller DC, Nazarian RM, Blakeley JO, Madaj ZB, Verma SK, Stemmer-Rachamimov A (2019) Assessing interobserver variability and accuracy in the histological diagnosis and classification of cutaneous neurofibromass. Neuro-Oncology Adv. https://doi.org/10.1093/noajnl/vdz050
Pemov A, Li H, Patidar R, Hansen NF, Sindiri S, Hartley SW, Wei JS, Elkahloun A, Chandrasekharappa SC, Program NCS et al (2017) The primacy of NF1 loss as the driver of tumorigenesis in neurofibromatosis type 1-associated plexiform neurofibromas. Oncogene 36:3168–3177. https://doi.org/10.1038/onc.2016.464
Porter CM, Haffner MC, Kulac I, Maynard JP, Baena-Del Valle JA, Isaacs WB, Yegnasubramanian S, De Marzo AM, Sfanos KS (2019) Lactoferrin CpG Island Hypermethylation and decoupling of mRNA and protein expression in the early stages of prostate carcinogenesis. Am J Pathol 189:2311–2322. https://doi.org/10.1016/j.ajpath.2019.07.016
Rasmussen SA, Friedman JM (2000) NF1 gene and neurofibromatosis 1. Am J Epidemiol 151:33–40
Riccardi VM (1992) Neurofibromatosis: phenotype, natural history and pathogenesis. Johns Hopkins University Press, City
Rodriguez EF, Blakeley J, Langmead S, Olivi A, Tufaro A, Tabbarah A, Berkenblit G, Sacks JM, Newsome SD, Montgomery E et al (2017) Low-grade Schwann cell neoplasms with leptomeningeal dissemination: clinicopathologic and autopsy findings. Hum Pathol 60:121–128. https://doi.org/10.1016/j.humpath.2016.08.010
Rohrich M, Koelsche C, Schrimpf D, Capper D, Sahm F, Kratz A, Reuss J, Hovestadt V, Jones DT, Bewerunge-Hudler M et al (2016) Methylation-based classification of benign and malignant peripheral nerve sheath tumors. Acta Neuropathol 131:877–887. https://doi.org/10.1007/s00401-016-1540-6
Shah M, Cardenas R, Wang B, Persson J, Mongan NP, Grabowska A, Allegrucci C (2017) HOXC8 regulates self-renewal, differentiation and transformation of breast cancer stem cells. Mol Cancer 16:38. https://doi.org/10.1186/s12943-017-0605-z
Shaheduzzaman S, Vishwanath A, Furusato B, Cullen J, Chen Y, Banez L, Nau M, Ravindranath L, Kim KH, Mohammed A et al (2007) Silencing of Lactotransferrin expression by methylation in prostate cancer progression. Cancer Biol Ther 6:1088–1095. https://doi.org/10.4161/cbt.6.7.4327
Tada M, Sawamura Y, Ishii N, Chin S, Abe H (1998) Massive plexiform neurofibroma in the orbit in a child with von Recklinghausen's disease. Childs Nerv Syst 14:210–212
Tao Y, Xi S, Briones V, Muegge K (2010) Lsh mediated RNA polymerase II stalling at HoxC6 and HoxC8 involves DNA methylation. PLoS One 5:e9163. https://doi.org/10.1371/journal.pone.0009163
Wong WW, Hirose T, Scheithauer BW, Schild SE, Gunderson LL (1998) Malignant peripheral nerve sheath tumor: analysis of treatment outcome. Int J Radiat Oncol Biol Phys 42:351–360 S0360-3016(98)00223-5 [pii]
Xavier FC, Destro MF, Duarte CM, Nunes FD (2014) Epigenetic repression of HOXB cluster in oral cancer cell lines. Arch Oral Biol 59:783–789. https://doi.org/10.1016/j.archoralbio.2014.05.001
Yi HM, Li H, Peng D, Zhang HJ, Wang L, Zhao M, Yao KT, Ren CP (2006) Genetic and epigenetic alterations of LTF at 3p21.3 in nasopharyngeal carcinoma. Oncol Res 16:261–272. https://doi.org/10.3727/000000006783981008
Zhang ML, Suarez MJ, Bosley TM, Rodriguez FJ (2017) Clinicopathological features of peripheral nerve sheath tumors involving the eye and ocular adnexa. Hum Pathol 63:70–78. https://doi.org/10.1016/j.humpath.2017.02.006
Zimmerman RA, Bilaniuk LT, Metzger RA, Grossman RI, Schut L, Bruce DA (1983) Computed tomography of orbitalfacial neurofibromatosis. Radiology 146:113–116
The authors thank Wayne Yu at the Microarray Core. Samples quality assessment and microarray analysis were conducted at The Sidney Kimmel Cancer Center Microarray Core Facility at Johns Hopkins University, supported by NIH grant P30 CA006973 entitled Regional Oncology Research Center.
This work was funded in part by a collaborative grant between the King Khaled Eye Specialist Hospital, Riyadh, Saudi Arabia and Wilmer Eye Institute, Baltimore, MD (to FJR) and NIH grant P30 CA006973 to the Sidney Kimmel Comprehensive Cancer Center (PI: W. Nelson).
Ethics approval and consent to participate
All studies were performed under Institutional Review Board (IRB) approval at Johns Hopkins (IRB#00079727). Appropriate consent or waiver of consent when applicable was obtained for all cases.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Differentiatially methylated probes and corresponding genes in orbitofacial vs. non-orbitofacial neurofibromas adjusted by age.
Clinicopathologic data of orbitofacial and non-orbitofacial neurofibromas.
Additional file 3.
About this article
Cite this article
Arnold, A., Imada, E.L., Zhang, M.L. et al. Differential gene methylation and expression of HOX transcription factor family in orbitofacial neurofibroma. acta neuropathol commun 8, 62 (2020). https://doi.org/10.1186/s40478-020-00940-7
- Plexiform neurofibromas
- HOX genes