Skip to main content

Tumor microenvironment in a minipig model of spinal cord glioma



Spinal cord glioma (SCG) is considered an orphan disease that lacks effective treatment options with margins that are surgically inaccessible and an overall paucity of literature on the topic. The tumor microenvironment is a critical factor to consider in treatment and modeling design, especially with respect to the unresectable tumor edge. Recently, our group developed a high-grade spinal cord glioma (SCG) model in Göttingen minipigs.


Immunofluorescence and ELISA were performed to explore the microenvironmental features and inflammation cytokines in this minipig SCG model. Protein carbonyl assay and GSH/GSSG assay were analyzed in the core and edge lesions in the minipig SCG model. The primary core and edge cells proliferation rate were shown in vitro, and the xenograft model in vivo.


We identified an elevated Ki-67 proliferative index, vascular and pericyte markers, CD31 and desmin in the tumor edge as compared to the tumor core. In addition, we found that the tumor edge demonstrated increased pro-inflammatory and gliomagenic cytokines including TNF-α, IL-1β, and IL-6. Furthermore, the mediation of oxidative stress is upregulated in the tumor edge. Hypoxic markers had statistically significant increased staining in the tumor core, but were notably still present in the tumor edge. The edge cells cultures derived from SCG biopsy also demonstrated an increased proliferative rate compared to core cell cultures in a xenotransplantation model.


Our study demonstrates heterogeneity in microenvironmental features in our minipig model of high-grade SCG, with a phenotype at the edge showing increased oxidative stress, proliferation, inflammatory cytokines, neovascularization, and decreased but present staining for hypoxic markers. These findings support the utility of this model as a means for investigating therapeutic approaches targeting the more aggressive and surgically unresectable tumor border.


Spinal Cord Gliomas (SCGs) account for 4–8% of central nervous system (CNS) tumors, where high-grade lesions have a poor prognosis [1]. Most commonly, gross total resection is not an option in such lesions due to a lack of resectable margins and the risk of worsening morbidity [2]. Importantly, SCG is considered an orphan disease that is overlooked in comparison to its intracranial counterpart whereby human specimens are limited owing to this relative rarity. Consequently, there is a paucity of literature on SCG pathogenesis, intratumoral heterogeneity, and tumor microenvironment despite some recent small retrospective studies [3,4,5,6].

While rodent models have been utilized in the preclinical setting for the study of SCG, there are known limitations especially with regard to anatomic constraints, drug delivery, and surgical translatability [7, 8]. In the past, efforts to investigate surgical techniques, drug delivery, and de novo tumor development have depended on spontaneous canine gliomas [9] which are rare and present with ethical and logistical issues. Moreover, SCG is extremely rare in dogs. As such, development of new therapeutic strategies for SCG has been hindered by the lack of a relevant preclinical animal model for understanding the delivery, therapeutic response, efficacy, toxicity, and tumor microenvironment of SCG.

Our group has developed a minipig model of high-grade SCG using lentiviral gene transfer of oncogenic transgenes including platelet derived growth factor beta (PDGF-B), HRAS-G12V, and small hairpin-loop RNA targeting P53 (shRNA-P53) [10] as well a vetted methodology for inducing glioma in prior rodent studies[11,12,13] and validated this approach in the rat spinal cord [14]. We employed this approach given involvement of the RTK/RAS/PI3K and p14/CDKN2A and P53 pathways in a large majority of high-grade gliomas [10, 11]. This model has histopathologic and radiologic features consistent with human disease processes in a highly penetrant manner [10]. The advantage of utilizing a large animal model is to facilitate translation given the similarity of human's metabolic, oxidative, genetic, and immunologic features to that of the swine [15]. However, the tumor microenvironment and heterogeneity of our model is unknown. In an effort to further characterize and advance preclinical modeling, the present study reports key microenvironmental features of the model, including oxidative stress, vascular markers, and intratumoral heterogeneity between the tumor core and edge.

Among these features is oxidative stress, a well-known contributor to glioma pathogenesis. It is exerted by reactive oxygen species (ROS) that accumulate due to an imbalance between ROS generation and elimination, which can consequently promote glioma growth and progression [16]. In the tumor microenvironment, an interplay between tumor cells and surrounding parenchyma can enable tumor progression by sustaining malignant cells and providing adaptation to restricted oxygen supply [17]. In addition to adaptation to a hypoxic environment, ROS generation can lead to lipid, protein, and DNA damage that can further contribute to glioma pathogenesis. Furthermore, it is well documented that the tumor core in high grade supratentorial gliomas have a hypoxic core and over-expression and stabilization of hypoxia inducible factor-1 (HIF-1α) [18, 19]. Downstream of HIF-1α are hypoxia responsive genes, including glucose transporter 1 (GLUT-1) and CA-IX (Carbonic Anhydrase 9) which are expressed in cells undergoing anaerobic metabolism and demonstrated to be overexpressed in glioblastoma cells [20]. In addition, this hypoxic microenvironment leads to a variety of downstream effects including a pro-angiogenic environment, metabolic reprogramming, and local immunosuppression [21]. Consequently, understanding whether this phenomenon is present in this model is crucial to its translational relevance.

The differences in oxidative stress, region specific factors and subsequent molecular features may be found depending on the geographic location within the tumor including well-reported differences in proliferation [22], inflammatory cytokines [23], and microvascular structure [24, 25]. Certainly, molecular subtypes have been identified from a transcriptomic standpoint in supratentorial high-grade glioma [26, 27] and correlated with characteristic immunohistochemical profiles in regions of hypoxia and hyper-proliferation amongst others [28]. Given these studies in clinical and preclinical models in the brain, it stands to reason that modeling SCG requires robust characterization of factors. This is especially crucial from a clinical standpoint in the spinal cord, whereby the unresectable edge of the tumor is a major cause for recurrence.

In our present study, we compare and contrast features of the tumor core and edge to evaluate differences in redox homeostasis, hypoxic markers, and validate the aggressive proliferation of tumor cells. We further described recapitulation of these features upon evaluation in porcine SCG-derived culture models both in vitro and on xenotransplantation. These data provide the basis for future work employing this large animal model as a platform for translational surgical studies, pharmacotherapeutic studies targeting redox homeostasis, and comparative studies with patient derived samples.

Materials and methods

Gӧttingen minipig spinal cord glioma model

Four Gӧttingen minipigs (± 15–20 kg, 8–9 months old) were employed in the present study. The surgical procedures were performed as previously described for administration of oncogenic vectors. Third-generation and VSV-G pseudotyped lentiviral vectors were employed as previously described [10]. Briefly, three individual vector constructs (1 – EF1a-PDGFB-IRES-eGFP, 2 – EF1a-HRAS-G12V-mPlum, and 3 – H6 & U6 shRNA-P53 with PGK-mCherry) were produced for viral packaging at > 10^9 IU/ml in a standard fashion. Animals were monitored using physical examinations at baseline and post-operatively three times per week until endpoint. All experiments were approved by the Institutional Animal Care and Use Committee (IACUC) in coordination with the Division of Animal Resources (DAR) and veterinary staff [10].

Magnetic resonance imaging (MRI)

Pigs were scanned used a clinical 3 T Trio MRI with an integrated spine coil (Siemens Medical Solutions, Malvern, PA, USA) using standard clinical sequences including T1-weighted, T2-weighted, and T1-weighted post-gadolinium (0.1 mmol/kg intravenous, Multihance, Bracco Diagnostics, Italy). Scans were reviewed and processed using RadiAnt DICOM Viewer (Medixant, Poznań, Poland) [10].

Nude mice xenograft model

Athymic nude mice were housed and maintained with the protocols approved by the Emory Institutional Animal Care and Use Committee (IACUC). The nude mice are 5–6 weeks old, half male and half female, with eight mice per group, total 16 mice. Pig SCG core and edge cells (3 × 106) in 100 μl of PBS were inoculated subcutaneously into the mice on the dorsal flank. Body weight and tumor growth were recorded once every three days. The total tumor volume (TV) was calculated by the following formula: TV (mm3) = a * b2 /2, where "a" is the minimum diameter, and "b" denotes the maximum diameter. The mice were euthanized after three weeks [29].

Hematoxylin–eosin (H&E) staining and immunohistochemistry (IHC) staining

The minipig SCG tumors and xenograft tumors were fixed overnight with 4% PFA followed by paraffin embedding. H&E staining was performed to evaluate histopathologic features in pig SCG tumors and xenografts. Immunohistochemistry was performed using primary antibodies for Ki67 (ab15580, ABCAM, USA), SOX2 (ab92494, ABCAM, USA), Olig2 (P21954, Invitrogen, USA), NG2 (ab129051, ABCAM, USA), NQO1 (#3187, Cell Signaling Technology, USA), 4-HNE (MA5-27,570, Invitrogen, USA), HIF-1α (PA1-16,601, Invitrogen, USA), GLUT-1 (ab14683, ABCAM, USA), CA-IX (66,243, Proteintech, USA) with Vectastain Kit (PK-8200, Vector Laboratories, Inc., USA). All images were acquired using whole-slide scanning at 40 × magnification (Leica Aperio AT2 Slide Scanner) [14].

Immunofluorescence staining

Paraffin embedded sections were also stained with primary antibodies for PDGFB (ab23914, ABCAM, USA), HRAS (sc-29, SANTA CRUZ, USA), p53 (#2524, Cell Signaling Technology, USA), Ki67 (ab15580, ABCAM, USA), NQO1 (#3187, Cell Signaling Technology, USA), 4-HNE (MA5-27,570, Invitrogen, USA), Desmin (MA106401, Invitrogen, USA), CD31 (MCA1746GA, Invitrogen, USA), Collagen IV (ab6586, ABCAM, USA) with Opal Detection Kit (AKOYA BIOSCIENCES, USA) and the images were taken by Vectra Polaris Multispectral Imaging System (AKOYA BIOSCIENCES, USA) [29, 30].

Enzyme-linked immunosorbent (ELISA) assays

Tissue lysates of the core and edge of minipig SCG samples were harvested and quantified for cytokines levels using ELISA Kits, including IL-1β (ab100754, ABCAM, USA), IL-6 (ab100755, ABCAM USA) and TNF-α (88 ab100756, ABCAM, USA) in a standard fashion according to manufacturer specifications [31].

8-hydroxy 2 deoxyguanosine (8-OH-dG) Level, Protein Carbonyl (PC) Assay and GSH/GSSG Assay

A8-OH-dG level and ratio of GSH/GSSG allows for evaluation of the status of oxidative stress in the biological systems [32], while protein carbonyl content can measure Oxidative damage to proteins. The tissue lysates of the core and edge of minipig spinal cord glioma samples were harvested and quantified for 8-hydroxy 2 deoxyguanosine Kit (ab201734, ABCAM, USA), GSH/GSSG-Glo™ (Promega Corporation, USA) and Protein Carbonyl Assay Kit (#ab126287, ABCAM, USA) following with manufacturer's assays procedure.

ATP Assay and Lactate dehydrogenase (LDH) release

Tissue lysates of minipig SCG samples from core and edge locations were evaluated for total ATP and LDH release using ATP Assay Kit (ab83355, ABCAM, USA) and LDH assay kits (J2380, Promega Corporation, USA) according to the protocols and read using a microplate reader (BioTek Instruments Inc., USA) [31, 33].

Primary cell cultures of minipig SCG Core and edge samples

After sacrificed the minipigs, the tumors were collected from the spinal cord of the animal. Core and edge cell samples were cultured in Dulbecco's Modified Eagle Medium (DMEM), with addition of media components including 10 ng/mL PDGFAA (# 78,095, STEMCELL Technologies, USA),10 ng/mL bFGF (# 78,003.1, STEMCELL Technologies, USA), N2 Supplement (07152, STEMCELL Technologies, USA), 20% FBS and penicillin (100units/mL)-streptomycin (100 μg/mL) (Hyclone, USA) in a humidified incubator with 5% CO2 at 37 °C according to a previously described protocol [34].

Cell proliferation assay

In order to evaluate cell proliferation, core and edge cells were seeded (2 × 105 cells per mL) in growth medium in a 96-well plate and analyzed by MTT Cell Proliferation Assay Kit (ab211091, ABCAM, USA). Briefly, 50 μL of MTT Reagent was added to each well, and the plates were incubated for 3 h at 37 °C. After incubation, 150 µL of MTT Solvent were added to each well. The absorbance at OD 590 nm was recorded using a microplate reader (BioTek Instruments Inc., USA) at specified timepoints [35].

Statistical analysis

Data visualizations were performed with GraphPad Prism 6 (GraphPad Software Inc., La Jolla, CA, USA) and BioRender (BioRender, Toronto, Canada). The analyses were performed using either Student’s t-test or one-way ANOVA. Significant Difference among groups was assessed as *p < 0.05; **p < 0.01; ***p < 0.001.


Establishing minipig SCG model and comparison of core and edge

The development of our minipig spinal cord glioma model was reported previously [10]. In short, 14 days post-operatively, animals (N = 4/4) developed motor deficits in accordance with oncogene injection (Fig. 1A and B). 21 days post-operatively, MRI was used to validate the presence of an enhancement mass lesion (Fig. 1C) and animals underwent necropsy for tissue harvesting. H&E was used to confirm histopathologic features including hypercellularity, nuclear atypia, mitotic indices, necrosis, hemorrhage, increased vascularity, and astrocytic morphology. An outline of the definition of core and edge regions of the tumor is presented in Fig. 1D and grossly defined as 1 mm border circumferential border [36].

Fig. 1
figure 1

Establishment a Cohort of Minipig SCG and Defining Core and Edge Regions. a Lentiviral vectors expressing PDGFB, HRASG12V, and shRNA targeting P53 (sh787, sh944) were designed based on prior literature. b Schematic representation of the timeline of tumor progression in the spinal cord minipig model. c Representative endpoint MRI of T1w post-contrast demonstrating a contrast enhancing lesion in the oncogene injection site. d Representative H&E staining showing edge and core tumor regions in the spinal cord glioma. Scale bar: 2 mm and 100 μm

Enrichment of SCG edge in pericytes and endothelial cells

One crucial component to microenvironmental differences in high-grade glioma is differences in vascular structure, including endothelial cells, pericytes, and basement membrane. We co-stained pig SCG tumors with anti-CD31 (PECAM-1), anti-Desmin antibodies, anti–collagen IV antibodies, and DAPI (Fig. 2A). We found a statistically significant elevation in CD-31 and Desmin positivity (p < 0.01) in the edge compare to the core cells. No significant difference was observed in the overall cellularity or Collagen IV (Fig. 2B). In the tumor edge, we also observed a high-degree of pericyte coverage of blood vessels and vascular endothelium in a disorganized fashion.

Fig. 2
figure 2

Vascular Microenvironment of SCG edge Demonstrates Enrichment in Endothelial Cells and Pericytes. a Immunofluorescence imaging of core and edge tumor sections co-stained for anti-CD31(endothelial cells, red), anti-Desmin (pericytes, green), and anti-collagen IV (basement membrane, grey). Scale bars are 50 μm. b Quantification of DAPI signal area, CD31, Desmin, and collagen IV-positive areas. Data are shown as means ± SEM. (**p < 0.01, one-way ANOVA, n = 4)

Increased proliferative glial progenitors and cytokines in the edge

In comparing the tumor core and edge regions, we found statistically significant elevations in proliferative marker Ki-67 as well as glial progenitor markers SOX2, Olig2, and NG2, (p < 0.05, Fig. 3A). In order to further characterize gliomagenic cytokines in this model, we evaluated TNF-a, IL-1B, and IL-6. On ELISA of tumor lysates of core and edge samples, edge lysates demonstrated statistically significant elevations in cytokines TNF-α, IL-1β and IL-6 compared to core lysates (Fig. 3 B, C, D).

Fig. 3
figure 3

Tumor Edge Demonstrates Elevated Proliferation and Inflammatory Cytokines. a Glioma markers and proliferative index were analyzed by immunohistochemistry staining. Scale bar: 100 μm. b IL-1β, c IL-6, and d TNF-α levels in the core and edge lysate cells of the minipig SCG model. The spinal cord glioma lysates were subjected to TNF-α, IL-1β, and IL-6 ELISA analysis. Data are shown as means ± SEM. (*p < 0.05, one-way ANOVA, n = 4)

Elevated markers of oxidative stress in the SCG edge cells

Given significant elevations in proliferative index, we investigated markers of oxidative stress in the core and edge regions in the tissues. 8-OH-dg, a biomarker of DNA damage in response to oxidative stress, was elevated two-fold on average in edge versus core lysates (p < 0.05, Fig. 4A). Protein Carbonyl expression levels, a marker of oxidative stress and protein oxidation, were also elevated in the edge samples (p < 0.01, Fig. 4B). On examination of GSH to GSSG ratios, core samples were significantly elevated (p < 0.05, Fig. 4C) compared to edge samples, indicating increased oxidative stress in the edge samples. Overall, this suggests elevated ROS and oxidative stress in the edge samples compared to core samples. ATP and LDH assay revealed enhanced cell death in edge cells compared to core cells (Fig. 4 D and E). Moreover, 4-HNE and NQO1 mirrored the Ki67 expression patterns on immunofluorescent co-staining (Fig. 4F). Hence, these data support elevated ROS in the edge region of the minipig SCG model. Examining microenvironmental features of redox homeostasis would be incomplete without also characterizing markers of hypoxia in the present model, given its well documented role in glioma progression and impact on glioma cell metabolism [17, 20, 21]. As such, we further examined hypoxic markers through immunohistochemical staining and found that HIF-1α, CA-IX, and GLUT-1 staining was increased in the core compared to the tumor edge (Fig. 4G). This suggests that the present model recapitulates a largely hypoxic tumor core, and that the tumor edge has a degree of hypoxia, albeit less than the tumor core.

Fig. 4
figure 4

SCG Tumor Edge Exhibits Increased Markers of Oxidative Stress. a Expression of 8-hydroxy-2-deoxyguanosine. b Carbonyl expression levels. c GSH/GSSG ratios. d ATP levels. e LDH ratios. Data are means ± SEM (*p < 0.05; **p < 0.01; ***p < 0.001, n = 4) f Representative immunofluorescent staining of edge and core minipig SCG tissues for 4-HNE (green), Ki67 (red), NQO1 (cyan) and DAPI (blue). Scale bars: 100 μm. n = 4. g HIF-1α, CA-IX, GLUT-1 staining in tumor core and edge of minipig SCG tissues. Scale bars: 100 μm. n = 4

Proliferative phenotype of edge cells is recapitulated in culture and pig-derived xenograft

In order to explore whether edge cells recapitulate their proliferative phenotype in vitro, we established the core and edge tissue culture models derived from fresh biopsies of the minipig SCGs collected at the time of euthanasia. We confirmed presence of vector reporters by immunofluorescent staining (Fig. 5A). Next, we compared the rate of proliferation and found that edge cells exhibited higher cell proliferation rates than core cells (Fig. 5B).

Fig. 5
figure 5

Cell proliferation of core and edge minipig spinal cord glioma cells in a nude mice xenograft model a Core and edge cells were staining with GFP and mCherry. Scale bar: 50 μm; b Cell proliferation rates of two cell lines were determined by MTT assay for 5 days. c Body weight of nude mice with flank injection of core and edge cells. d Tumor volume curve in nude mice by core and edge cells. e Representative macroscopic images of tumor growth in core and edge nude mice models. f H&E, Ki67, 4-HNE and NQO1 staining of tumor slices in various groups. Scale Bars: 50 μm. Significant differences are considered when *p < 0.05; **p < 0.01. n = 8

We then looked to examine this phenotype in vivo in a xenograft model. Both cell lines were inoculated in nude mice subcutaneously, and tumor volumes were quantitatively measured every three days. The tumor growth did not significantly affect the nude mice's body weight (Fig. 5C). Notably, the edge xenografts exhibited significantly increased tumor growth compared to the core cells (Fig. 5D and E). After termination, H&E staining showed a high-grade glioma with features consistent with prior analyses in both core and edge xenografts. Ki-67 staining once more showed that cell proliferation was strongly elevated in the edge xenografts. ROS staining assay with 4-HNE (4-hydroxynonenal), a marker for ROS level, showed that edge cells evidently ameliorated ROS, coupled with NQO1 immunohistochemical staining (Fig. 5F). We then examined if the xenograft model would reproduce the findings on hypoxic markers, and found similar results on immunohistochemical staining with increased HIF-1α, CA-IX, and GLUT-1 staining in the tumor core vs the tumor edge (Additional file 1: Figure S1).


Our present study corroborated several key features of our previously described minipig model of high-grade SCG. From a practical standpoint, we continued to report a highly penetrant model system with tumor induction following vector administration. We further described and confirmed the expression of relevant transgenes on immunofluorescence with no major differences in the tumor core compared to the tumor edge. Of note, we observed significant elevation in proliferation in the SCG tumor edge. This finding was recapitulated on establishing tissue culture from biopsies of the core and edge, as well as on xenotransplantation into nude mice. This highlights a point of clinical interest, - it is often tumor cells at the periphery which have higher proliferation rates that are surgically inaccessible given the infiltrative nature of eloquent CNS regions like SCG. The inability to perform a gross total resection given the infiltrative nature of SCGs represents one of main barriers to development of an effective treatment strategies and is the key reason for our interest in the heterogeneity between cells in the tumor core and edge.

The findings in our model were in agreement with the literature in describing increased proliferation, increased vascularity, and oxidative stress at the edge in supratentorial glioma as well as increased hypoxic markers in the tumor core [20, 37,38,39]. It has been widely reported that the tumor edge in supratentorial high-grade glioma exhibits a more aggressive, proliferative, and invasive phenotype [22, 39]. We examined vascular markers and found that edge regions were enriched in endothelial cells, pericytes, and overall vasculature consistent with prior reports [24]. In non-disease states, pericytes have a multitude of roles, including maintenance of the blood–brain-barrier through interaction with astrocytes, development and maturation of vasculature, paracrine signaling, as well as exhibiting stem cell properties [40,41,42,43]. It has also been reported that pericytes directly promote glioma progression in in vitro and in vivo models [44]. Moreover, microvascular proliferation is a well-known feature of high-grade gliomas. In order to explain this finding, the increased vascularity at the tumor edge in this model might simply be due to the innate nature of microvascular proliferation secondary to rising metabolic needs locally in the more mitotically active edge region. This proliferation is conventionally attributed to the fact that hypoxia was not necessarily restricted solely to the tumor core [37]. This concept is supported by the fact that while the tumor edge had, as expected, a significantly lower degree of staining for hypoxic markers (HIF-1a, CA-IX, GLUT-1), up to 30% of the cells in the tumor edge still stained strongly for these markers (Fig. 4G, Additional file 1: Figure S1). As such, this region may represent a site of transitional hypoxia in the setting of rapid tumor proliferation that should be examined in future studies including identifying specific cell populations and highlighting potential translational targets [37]. This is of particular interest given that since pericytes and endothelial cells both express PDGF receptors, it is also possible that our minipig model with the PDGFB transgene have been driven to elevated levels of pericytes and microvascular proliferation [24].

In examining select cytokines in the core and edge regions, we observed an elevation in TNF-α, IL-1β, and IL-6 in the edge compared to the tumor core. IL-6 is a cytokine that is known to trigger JAK-STAT3 signaling in glioblastoma, gliomagenesis, MGMT methylation, and is correlated with tumor grade and overall patient survival [45,46,47]. TNF-α and IL-1β are pro-inflammatory cytokines that are found in supratentorial high-grade gliomas and are involved in neuroinflammation, gliomagenesis, and development of radio-resistance [23, 48]. Overall, the elevated presence of these cytokines throughout the tumor was consistent with reports on supratentorial high-grade gliomas. It was of interest with utilizing this model as a platform for radio-resistance, drug resistance, or aggressive phenotype in this particularly challenging region.

For evaluation of oxidative stress, we examined markers of oxidative DNA damage (8-OH-dg), protein oxidative damage (protein carbonyl content), as well as non-enzymatic antioxidant defense through ROS scavenging (GSH/GSSG ratio) [49]. We found that the tumor edge has a significantly elevated level of markers of ROS, including a high 8-OH-dG and low GSH/GSSG ratio, as well as oxidative damage with an elevated protein carbonyl content in Edge tumor lysates. This supports the notion of elevated oxidative stress in the setting of a highly proliferative tumor edge. On histologic evaluation, we further observed an increase in NQO1 and 4-HNE staining patterns at the tumor edge. Next steps toward further examining molecular characteristics of core and edge regions of this model will include rigorous evaluation of hypoxic markers including , PDK1 and others. In addition to redox mechanisms contribution to glioma pathophysiology [17], the difference in oxidative stress between the core and edge in this model is of particular interest given that there are a variety of druggable targets in the realm of oxidative stress [16, 29, 50,51,52,53].

The development of advanced pre-clinical disease modeling has increasingly employed pigs as biomedical models given their increased anatomic, genetic, and immunologic similarity to the human [54]. Porcine models, including minipigs, have successfully been utilized to model numerous pathologies including colorectal cancer, osteosarcoma, cystic fibrosis, muscular dystrophy, as well as numerous others [55, 56]. Notably, the use of TALEN mediated gene editing was utilized to generate a minipig model of neurofibromatosis type 1 that better recapitulates features and natural history of the human disease [10, 57]. Specific to glioma models, as of 2020 there are only three published models including two orthotopic xenografts [58, 59] using commercially available cell lines, as well as our presently reported spinal cord model using vector driven gene transfer [60]. Further work to develop enhanced porcine models of glioma may create a platform better suited for pre-clinical translational efforts.


Overall, the data presented in this study demonstrated multiple distinct microenvironmental features within the pig model of high-grade SCG. Future work on the development of this model system could explore each of these individual components, including but not limited to profiling the inflammatory microenvironment, transcriptomic characterization, or further elucidating the cause for findings of intra-tumoral heterogeneity between the core, edge, perivascular niche, or other regions. However, at present this model represents a step forward for facilitating surgical translation, as well as a potential space for evaluating proof-of-concept for surgical resection followed by select druggable targets that could prove more efficient in treating the infiltrative tumor edge which is surgically inaccessible. It would be advantageous to also initially employ culture and xenograft models derived from the pig SCG with high-throughput screening methodologies to identify existing FDA-approved or novel candidates that may be able to better target the tumor edge cells.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



Spinal cord glioma


Central nervous system


Reactive oxygen species








8-Hydroxy 2 deoxyguanosine


Protein Carbonyl


Lactate dehydrogenase




NAD(P)H quinone dehydrogenase 1


  1. Hsu S, Quattrone M, Ostrom Q, Ryken TC, Sloan AE, Barnholtz-Sloan JS. Incidence patterns for primary malignant spinal cord gliomas: a Surveillance, Epidemiology, and End Results study. J Neurosurg Spine. 2011;14(6):742–7.

    PubMed  PubMed Central  Google Scholar 

  2. Abd-El-Barr MM, Huang KT, Chi JH. Infiltrating spinal cord astrocytomas: Epidemiology, diagnosis, treatments and future directions. J Clin Neurosci. 2016;29:15–20.

    PubMed  Google Scholar 

  3. Teng YD, Abd-El-Barr M, Wang L, Hajiali H, Wu L, Zafonte RD. Spinal cord astrocytomas: progresses in experimental and clinical investigations for developing recovery neurobiology-based novel therapies. Exp Neurol. 2019;311:135–47.

    PubMed  Google Scholar 

  4. Misove A, Vicha A, Broz P, Vanova K, Sumerauer D, Stolova L, Sramkova L, Koblizek M, Zamecnik J, Kyncl M, et al. Integrated genomic analysis reveals actionable targets in pediatric spinal cord low-grade gliomas. Acta Neuropathol Commun. 2022;10(1):143.

    PubMed  Google Scholar 

  5. Konovalov N, Timonin S, Asyutin D, Raevskiy M, Sorokin M, Buzdin A, Kaprovoy S. Transcriptomic Portraits and Molecular Pathway Activation Features of Adult Spinal Intramedullary Astrocytomas. Front Oncol. 2022;12: 837570.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Zhang M, Iyer RR, Azad TD, Wang Q, Garzon-Muvdi T, Wang J, Liu A, Burger P, Eberhart C, Rodriguez FJ, et al. Genomic landscape of intramedullary spinal cord gliomas. Sci Rep. 2019;9(1):18722.

    CAS  PubMed  PubMed Central  Google Scholar 

  7. Ellis JA, Castelli M, Assanah M, Bruce JN, Canoll P, Ogden AT. Unique microenvironmental responses to PDGF stimulation in brain and spinal cord gliomas determine tumor phenotype. J Neurooncol. 2015;123:27–33.

    CAS  PubMed  Google Scholar 

  8. Flisikowska T, Kind A, Schnieke A. The new pig on the block: modelling cancer in pigs. Transgenic Res. 2013;22:673–80.

    CAS  PubMed  Google Scholar 

  9. Dickinson PJ, LeCouteur RA, Higgins RJ, Bringas JR, Larson RF, Yamashita Y, Krauze MT, Forsayeth J, Noble CO, Drummond DC. Canine spontaneous glioma: a translational model system for convection-enhanced delivery. Neuro Oncol. 2010;12(9):928–40.

    PubMed  PubMed Central  Google Scholar 

  10. Tora MS, Texakalidis P, Neill S, Wetzel J, Rindler RS, Hardcastle N, Nagarajan PP, Krasnopeyev A, Roach C, James R. Lentiviral vector induced modeling of high-grade spinal cord glioma in minipigs. Sci Rep. 2020;10(1):5291.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Lynes J, Wibowo M, Koschmann C, Baker GJ, Saxena V, Muhammad AK, Bondale N, Klein J, Assi H, Lieberman AP, et al. Lentiviral-induced high-grade gliomas in rats: the effects of PDGFB, HRAS-G12V, AKT, and IDH1-R132H. Neurotherapeutics. 2014;11(3):623–35.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Assanah M, Lochhead R, Ogden A, Bruce J, Goldman J, Canoll P. Glial progenitors in adult white matter are driven to form malignant gliomas by platelet-derived growth factor-expressing retroviruses. J Neurosci. 2006;26(25):6781–90.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Rahme GJ, Luikart BW, Cheng C, Israel MA. A recombinant lentiviral PDGF-driven mouse model of proneural glioblastoma. Neuro Oncol. 2018;20(3):332–42.

    CAS  PubMed  Google Scholar 

  14. Nagarajan PP, Tora MS, Neill SG, Federici T, Texakalidis P, Donsante A, Canoll P, Lei K, Boulis NM. Lentiviral-induced spinal cord gliomas in rat model. Int J Mol Sci. 2021;22(23):12943.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Ribitsch I, Baptista PM, Lange-Consiglio A, Melotti L, Patruno M, Jenner F, Schnabl-Feichter E, Dutton LC, Connolly DJ, van Steenbeek FG. Large animal models in regenerative medicine and tissue engineering: to do or not to do. Front Bioeng Biotechnol. 2020;8:972.

    PubMed  PubMed Central  Google Scholar 

  16. Liu S, Dong L, Shi W, Zheng Z, Liu Z, Meng L, Xin Y, Jiang X. Potential targets and treatments affect oxidative stress in gliomas: an overview of molecular mechanisms. Front Pharmacol. 2022;13: 921070.

    CAS  PubMed  PubMed Central  Google Scholar 

  17. Krawczynski K, Godlewski J, Bronisz A. Oxidative stress-part of the solution or part of the problem in the hypoxic environment of a brain tumor. Antioxidants (Basel). 2020;9(8):747.

    CAS  PubMed  Google Scholar 

  18. Rong Y, Durden DL, Van Meir EG, Brat DJ. “Pseudopalisading” necrosis in glioblastoma: a familiar morphologic feature that links vascular pathology, hypoxia, and angiogenesis. J Neuropathol Exp Neurol. 2006;65(6):529–39.

    PubMed  Google Scholar 

  19. Renfrow JJ, Soike MH, West JL, Ramkissoon SH, Metheny-Barlow L, Mott RT, Kittel CA, D’Agostino RB Jr, Tatter SB, Laxton AW, et al. Attenuating hypoxia driven malignant behavior in glioblastoma with a novel hypoxia-inducible factor 2 alpha inhibitor. Sci Rep. 2020;10(1):15195.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. Domenech M, Hernandez A, Plaja A, Martinez-Balibrea E, Balana C. Hypoxia: The Cornerstone of Glioblastoma. Int J Mol Sci. 2021;22(22):12608.

    CAS  PubMed  PubMed Central  Google Scholar 

  21. Colwell N, Larion M, Giles AJ, Seldomridge AN, Sizdahkhani S, Gilbert MR, Park DM. Hypoxia in the glioblastoma microenvironment: shaping the phenotype of cancer stem-like cells. Neuro Oncol. 2017;19(7):887–96.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Bastola S, Pavlyukov MS, Yamashita D, Ghosh S, Cho H, Kagaya N, Zhang Z, Minata M, Lee Y, Sadahiro H, et al. Glioma-initiating cells at tumor edge gain signals from tumor core cells to promote their malignancy. Nat Commun. 2020;11(1):4660.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Kore RA, Abraham EC. Inflammatory cytokines, interleukin-1 beta and tumor necrosis factor-alpha, upregulated in glioblastoma multiforme, raise the levels of CRYAB in exosomes secreted by U373 glioma cells. Biochem Biophys Res Commun. 2014;453(3):326–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. Behling K, Maguire WF, Di Gialleonardo V, Heeb LE, Hassan IF, Veach DR, Keshari KR, Gutin PH, Scheinberg DA, McDevitt MR. Remodeling the vascular microenvironment of glioblastoma with alpha-particles. J Nucl Med. 2016;57(11):1771–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Chen J, Mao S, Li H, Zheng M, Yi L, Lin JM, Lin ZX. The pathological structure of the perivascular niche in different microvascular patterns of glioblastoma. PLoS ONE. 2017;12(8): e0182183.

    PubMed  PubMed Central  Google Scholar 

  26. Behnan J, Finocchiaro G, Hanna G. The landscape of the mesenchymal signature in brain tumours. Brain. 2019;142(4):847–66.

    PubMed  PubMed Central  Google Scholar 

  27. Darmanis S, Sloan SA, Croote D, Mignardi M, Chernikova S, Samghababi P, Zhang Y, Neff N, Kowarsky M, Caneda C, et al. Single-cell RNA-Seq analysis of infiltrating neoplastic cells at the migrating front of human glioblastoma. Cell Rep. 2017;21(5):1399–410.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Bergmann N, Delbridge C, Gempt J, Feuchtinger A, Walch A, Schirmer L, Bunk W, Aschenbrenner T, Liesche-Starnecker F, Schlegel J. The intratumoral heterogeneity reflects the intertumoral subtypes of glioblastoma multiforme: a regional immunohistochemistry analysis. Front Oncol. 2020;10:494.

    PubMed  PubMed Central  Google Scholar 

  29. Lei K, Xia Y, Wang X-C, Ahn EH, Jin L, Ye K. C/EBPβ mediates NQO1 and GSTP1 anti-oxidative reductases expression in glioblastoma, promoting brain tumor proliferation. Redox Biol. 2020;34: 101578.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Behling K, Maguire WF, Di Gialleonardo V, Heeb LE, Hassan IF, Veach DR, Keshari KR, Gutin PH, Scheinberg DA, McDevitt MR. Remodeling the vascular microenvironment of glioblastoma with α-particles. J Nucl Med. 2016;57(11):1771–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  31. Ahn EH, Lei K, Kang SS, Wang Z-H, Liu X, Hong W, Wang YT, Edgington-Mitchell LE, Jin L, Ye K. Mitochondrial dysfunction triggers the pathogenesis of Parkinson’s disease in neuronal C/EBPβ transgenic mice. Mol Psychiatry. 2021;26(12):7838–50.

    CAS  PubMed  Google Scholar 

  32. Valavanidis A, Vlachogianni T, Fiotakis C. 8-hydroxy-2′-deoxyguanosine (8-OHdG): a critical biomarker of oxidative stress and carcinogenesis. J Environ Sci Health C. 2009;27(2):120–39.

    CAS  Google Scholar 

  33. Lei K, Kang SS, Ahn EH, Chen C, Liao J, Liu X, Li H, Edgington-Mitchell LE, Jin L, Ye K. C/EBPβ/AEP Signaling regulates the oxidative stress in malignant cancers, stimulating the metastasisC/EBPβ/AEP signaling regulates the oxidative stress in TNBC. Mol Cancer Ther. 2021;20(9):1640–52.

    CAS  PubMed  Google Scholar 

  34. Mayer-Pröschel M. Isolation and generation of oligodendrocytes by immunopanning. Curr Protocol Neurosci. 1998. 

  35. Lei K, Gu X, Alvarado AG, Du Y, Luo S, Ahn EH, Kang SS, Ji B, Liu X, Mao H. Discovery of a dual inhibitor of NQO1 and GSTP1 for treating glioblastoma. J Hematol Oncol. 2020;13:1–21.

    Google Scholar 

  36. Chowdhury FN, Reisinger J, Gomez KE, Chimed T-S, Thomas CM, Le PN, Miller B, Morton JJ, Nieto CM, Somerset HL. Leading edge or tumor core: intratumor cancer stem cell niches in oral cavity squamous cell carcinoma and their association with stem cell function. Oral Oncol. 2019;98:118–24.

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Lin CM, Yu CF, Huang HY, Chen FH, Hong JH, Chiang CS. Distinct tumor microenvironment at tumor edge as a result of astrocyte activation is associated with therapeutic resistance for brain tumor. Front Oncol. 2019;9:307.

    PubMed  PubMed Central  Google Scholar 

  38. Minata M, Audia A, Shi J, Lu S, Bernstock J, Pavlyukov MS, Das A, Kim SH, Shin YJ, Lee Y, et al. Phenotypic plasticity of invasive edge glioma stem-like cells in response to ionizing radiation. Cell Rep. 2019;26(7):1893–905.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. Huang Y, Qi L, Kogiso M, Du Y, Braun FK, Zhang H, Huang LF, Xiao S, Teo WY, Lindsay H, et al. Spatial dissection of invasive front from tumor mass enables discovery of novel microRNA drivers of glioblastoma invasion. Adv Sci (Weinh). 2021;8(23): e2101923.

    PubMed  Google Scholar 

  40. Cheng L, Huang Z, Zhou W, Wu Q, Donnola S, Liu JK, Fang X, Sloan AE, Mao Y, Lathia JD, et al. Glioblastoma stem cells generate vascular pericytes to support vessel function and tumor growth. Cell. 2013;153(1):139–52.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Zhang XN, Yang KD, Chen C, He ZC, Wang QH, Feng H, Lv SQ, Wang Y, Mao M, Liu Q, et al. Pericytes augment glioblastoma cell resistance to temozolomide through CCL5-CCR5 paracrine signaling. Cell Res. 2021;31(10):1072–87.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. Sun R, Kong X, Qiu X, Huang C, Wong PP. The emerging roles of pericytes in modulating tumor microenvironment. Front Cell Dev Biol. 2021;9: 676342.

    PubMed  PubMed Central  Google Scholar 

  43. Bergers G, Song S. The role of pericytes in blood-vessel formation and maintenance. Neuro Oncol. 2005;7(4):452–64.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Sena IFG, Paiva AE, Prazeres P, Azevedo PO, Lousado L, Bhutia SK, Salmina AB, Mintz A, Birbrair A. Glioblastoma-activated pericytes support tumor growth via immunosuppression. Cancer Med. 2018;7(4):1232–9.

    PubMed  PubMed Central  Google Scholar 

  45. West AJ, Tsui V, Stylli SS, Nguyen HPT, Morokoff AP, Kaye AH, Luwor RB. The role of interleukin-6-STAT3 signalling in glioblastoma. Oncol Lett. 2018;16(4):4095–104.

    PubMed  PubMed Central  Google Scholar 

  46. Wang H, Lathia JD, Wu Q, Wang J, Li Z, Heddleston JM, Eyler CE, Elderbroom J, Gallagher J, Schuschu J, et al. Targeting interleukin 6 signaling suppresses glioma stem cell survival and tumor growth. Stem Cells. 2009;27(10):2393–404.

    CAS  PubMed  Google Scholar 

  47. Bowman T, Broome MA, Sinibaldi D, Wharton W, Pledger WJ, Sedivy JM, Irby R, Yeatman T, Courtneidge SA, Jove R. Stat3-mediated Myc expression is required for Src transformation and PDGF-induced mitogenesis. Proc Natl Acad Sci U S A. 2001;98(13):7319–24.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Zhou W, Jiang Z, Li X, Xu Y, Shao Z. Cytokines: shifting the balance between glioma cells and tumor microenvironment after irradiation. J Cancer Res Clin Oncol. 2015;141(4):575–89.

    CAS  PubMed  Google Scholar 

  49. Ramirez-Exposito MJ, Martinez-Martos JM. The delicate equilibrium between oxidants and antioxidants in brain glioma. Curr Neuropharmacol. 2019;17(4):342–51.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Lei K, Gu X, Alvarado AG, Du Y, Luo S, Ahn EH, Kang SS, Ji B, Liu X, Mao H, et al. Discovery of a dual inhibitor of NQO1 and GSTP1 for treating glioblastoma. J Hematol Oncol. 2020;13(1):141.

    PubMed  PubMed Central  Google Scholar 

  51. Lei K, Xia Y, Wang XC, Ahn EH, Jin L, Ye K. C/EBPbeta mediates NQO1 and GSTP1 anti-oxidative reductases expression in glioblastoma, promoting brain tumor proliferation. Redox Biol. 2020;34: 101578.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Salazar-Ramiro A, Ramirez-Ortega D, Perez de la Cruz V, Hernandez-Pedro NY, Gonzalez-Esquivel DF, Sotelo J, Pineda B. Role of redox status in development of glioblastoma. Front Immunol. 2016;7:156.

    PubMed  PubMed Central  Google Scholar 

  53. Policastro LL, Ibanez IL, Notcovich C, Duran HA, Podhajcer OL. The tumor microenvironment: characterization, redox considerations, and novel approaches for reactive oxygen species-targeted gene therapy. Antioxid Redox Signal. 2013;19(8):854–95.

    CAS  PubMed  Google Scholar 

  54. Lunney JK, Van Goor A, Walker KE, Hailstock T, Franklin J, Dai C. Importance of the pig as a human biomedical model. Sci Transl Med. 2021;13(621):eabd5758.

    CAS  PubMed  Google Scholar 

  55. Perleberg C, Kind A, Schnieke A. Genetically engineered pigs as models for human disease. Dis Model Mech. 2018. 

  56. Ardan T, Baxa M, Levinska B, Sedlackova M, Nguyen TD, Klima J, Juhas S, Juhasova J, Smatlikova P, Vochozkova P. Transgenic minipig model of Huntington’s disease exhibiting gradually progressing neurodegeneration. Dis Model Mech. 2019.

  57. Isakson SH, Rizzardi AE, Coutts AW, Carlson DF, Kirstein MN, Fisher J, Vitte J, Williams KB, Pluhar GE, Dahiya S, et al. Genetically engineered minipigs model the major clinical features of human neurofibromatosis type 1. Commun Biol. 2018;1:158.

    PubMed  PubMed Central  Google Scholar 

  58. Selek L, Seigneuret E, Nugue G, Wion D, Nissou MF, Salon C, Seurin MJ, Carozzo C, Ponce F, Roger T, et al. Imaging and histological characterization of a human brain xenograft in pig: the first induced glioma model in a large animal. J Neurosci Methods. 2014;221:159–65.

    PubMed  Google Scholar 

  59. Khoshnevis M, Carozzo C, Bonnefont-Rebeix C, Belluco S, Leveneur O, Chuzel T, Pillet-Michelland E, Dreyfus M, Roger T, Berger F, et al. Development of induced glioblastoma by implantation of a human xenograft in Yucatan minipig as a large animal model. J Neurosci Methods. 2017;282:61–8.

    PubMed  Google Scholar 

  60. Hicks WH, Bird CE, Pernik MN, Haider AS, Dobariya A, Abdullah KG, Aoun SG, Bentley RT, Cohen-Gadol AA, Bachoo RM, et al. Large animal models of glioma: current status and future prospects. Anticancer Res. 2021;41(11):5343–53.

    CAS  PubMed  Google Scholar 

Download references


This research project was supported in part by the Viral Vector Core of the Emory Center for Neurodegeneration Disease, Emory Division of Animal Resources (DAR), Emory Center for Systems Imaging (CSI), Cancer Tissue and Pathology shared resource of Winship Cancer Institute of Emory University (NIH/NCI under award number P30CA138292) for their guidance and expertise in the conduct of this study. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.


This work is funded with Federal funds from the National Cancer Institute (NCI), National Institutes of Health (NIH) grant R01CA251393. Research reported in this publication was supported in part by the Emory Initiative, Biological Discovery through Chemical Innovation (BDCI)(00114072).

Author information

Authors and Affiliations



NB, KL and MT conceived the study. KL, MT, YL, HA, MZ and PN performed the experiments. SN analyzed the pathology. KL and MT analyzed the data. KL and MT prepared the figures. NB, KL and MT wrote the original manuscript draft. TF, JG, KH and YD reviewed and edited the paper. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Kecheng Lei or Nicholas M. Boulis.

Ethics declarations

Ethics approval and consent to participate

This study has been conducted in accordance with ethical standards and according to the Declaration of Helsinki and the national and international guidelines, and has been approved by Emory University.

Consent for publication

Not applicable.

Competing interests

No potential conflicts of interest were disclosed for the authors.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1

. IHC staining of edge and core groups for HIF-1α, GLUT1 and CAIX in a nude mice xenograft model. Scale bars: 100 μm. Significant differences are considered when *p < 0.05. n = 8.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tora, M.S., Neill, S.G., Lakhina, Y. et al. Tumor microenvironment in a minipig model of spinal cord glioma. J Transl Med 21, 667 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: