Case Collection and Histopathology Information
The Institutional Review Board of Nanfang Hospital affiliated to Southern Medical University approved the study. 11 PXA, 2 A-PXA, and 1 E-GBM diagnosed by pathologists from 2016 to 2019 in Nanfang Hospital were collected. In this study, the clinical data of age, gender, pathological diagnosis (according to the WHO 2016 Central nervous system tumor classification), anatomic location of glioma, tumor size, and pathological molecular markers were collected. Three neuro-oncologist reviews medical records for neuroimaging features, treatment options, and survival data. All specimens’ slides were reviewed by pathologists to confirm the diagnosis of all cases. For the diagnosis of PXA, the criteria determined by Brat et al. [3] were used. All our cases satisfied the diagnostic criteria for PXA and A-PXA, respectively: mitotic activity of PXA was ≤ 5 mitoses per 10 high-power fields (HPF) [12], A-PXA based on either mitotic rate > 5/10 HPFs or presence of necrosis [2, 11]. GBM displays nuclear atypia, cellular pleomorphism, mitotic activity, a diffuse growth pattern, microvascular proliferation and/or necrosis [10]. Besides, we conducted glial fibrillary acidic protein (GFAP) immunotherapy in all cases. Immunocytochemical detection of GFAP can identify the astrological nature of tumor cells and thus support the exclusion of non-glial neoplasms resembling PXA or A-PXA [6].
Molecular Detection
The pathological molecular markers mainly included GFAP (OriGene, monoclonal, 1:100 dilution), epithelial membrane antigen (EMA, OriGene, monoclonal, 1:100 dilution), O6-methylguanine-DNA methyltransferase (MGMT, OriGene, monoclonal, 1:100 dilution), neuronal nuclei (Neu-N, OriGene, monoclonal, 1:100 dilution), oligodendrocyte transcription factor 2 (Oligo-2, OriGene, monoclonal, 1:100 dilution), epidermal growth factor receptor (EGFR, OriGene, monoclonal, 1:100 dilution), vascular endothelial growth factor (VEGF, OriGene, polyclonal, 1:100 dilution), isocitrate dehydrogenase 1 (IDH1, OriGene, monoclonal, 1:100 dilution), Ki-67 (OriGene, monoclonal, 1:100 dilution), ATRX chromatin remodeler (ATRX, OriGene, polyclonal, 1:100 dilution), CD34 molecule (CD34, OriGene, monoclonal, 1:100 dilution), Synaptophysin (Syn, OriGene, monoclonal, 1:100 dilution), BRAF V600E (Roche, monoclonal, 1:100 dilution), H3K27M (OriGene, polyclonal, 1:100 dilution), Neurofilament (NF, OriGene, monoclonal, 1:100 dilution), Chromogranin A (CgA, OriGene, monoclonal, 1:100 dilution), CD20 (OriGene, monoclonal, 1:100 dilution), Cytokeratin (CK, OriGene, monoclonal, 1:100 dilution), CD30 (OriGene, monoclonal, 1:100 dilution), Vimentin (OriGene, monoclonal, 1:100 dilution), Capicua (CIC, OriGene, monoclonal, 1:400 dilution), far upstream element binding protein 1 (FUBP1, OriGene, monoclonal, 1:400 dilution), Desmin (OriGene, monoclonal, 1:100 dilution) and Reticular Fibres (abcam, monoclonal, 1:500 dilution) .
Fluorescence in situ hybridization (FISH) Detection
We performed FISH on the formalin-fixed and paraffin-embedded (FFPE) tissue sections using the same commercial 1p/19q probes (Vysis paired probes 1p36/1q25 and 19q13/19p13, Abbott Laboratories). Briefly, 5-µm-thick formalin-fixed, paraffin-embedded sections were deparaffinized, pre-treated with Vysis Pretreatment Solution (1N, NaSCN, Vysis 32-801200, Abbott Molecular Inc.) and digested in protease solution (Vysis 32-801200, Abbott Molecular Inc.). The probe mix (5 to 15 µl) was added to each slide according to the manufacturer’s instructions. Target DNA and probes were codenatured at 74℃ for 5 minutes and incubated at 37℃ overnight in a humidified hybridization chamber (ThermoBrite, Abbott Molecular Inc.). Post-hybridization washes were performed in 2x SSC/0.3% NP-40 for 2 minutes at 75℃. Finally, the slides were air-dried and counterstained with 4,6-diamidino-2-phenylindole (DAPI, Cell Signaling Technology, #4083).
Tissue Processing and NGS Panel Genotyping
For the IDH1, IDH2, ATRX, telomerase reverse transcriptase (TERT), tumor protein p53 (TP53), MGMT, H3.3 histone A (H3F3A), H3 clustered histone 2 (HIST1H3B), B-Raf proto-oncogene (BRAF), EGFR, and phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA) of the special case of the E-GBM patient, the NGS Panel sequencing and analysis were performed at the Genomics Laboratory of GenomicCare Biotechnology (Shanghai, China). For thawed soft tissue or blood, DNA was extracted using the Maxwell RSC Blood DNA Kit (cat# AS1400, Promega, Madison, WI, USA) on a Maxwell RSC system (cat# AS4500, Promega). For the FFPE tissue, DNA was extracted using the MagMAX FFPE DNA/RNA Ultra Kit (cat# A31881, ThermoFisher, Waltham, MA, USA) on a KingFisher Flex system (ThermoFisher). The extracted DNA was sheared using a Covaris L220 sonicator, then the library preparation and capture was done using Tecan EVO 150 (Thermofisher), and sequenced on an Ion S5 sequencer (Thermofisher) to generate paired-end reads. After removing adapters and low-quality reads, the reads were aligned to the National Center for Biotechnology Information (NCBI) human genome reference assembly hg19 using the Burrows-Wheeler Aligner alignment algorithm and further processed using the Genome Analysis Toolkit (GATK, version 3.5), including the GATK Realigner Target Creator to identify regions that needed to be realigned. Somatic single-nucleotide variants (SNV), insertion and deletion (Indel), and copy number variation (CNV) were determined using the MuTect/ANNOVAR/dbNSFP31, VarscanIndel, and CNV nator software respectively as reported in Zang et al. [24].
10 ng DNA from tumor tissue was extracted with Ion AmpliSeq™ Direct FFPE DNA Kit according to the manufacturer’s instruction for sequencing. The Ion AmpliSeq™ Library Kit was used for the amplification of ten genes (IDH1, IDH2, ATRX, TERT, TP53, H3F3A, HIST1H3B, BRAF, EGFR, and PIK3CA) and the preparation of amplicon libraries. Amplicon libraries were processed using the Ion 540™ Kit-Chef and then sequenced by the Thermo Fisher Scientific Ion S5™ System panel platform. Sequencing reads were mapped to hg19 with BWA using the default parameters. SNV annotation was performed by ANNOVAR and indels detection was analyzed with Varscan 2.4.2. Candidate somatic SNV and indels were identified if they satisfied the following conditions: the sequencing depth was larger than 300, then they were supported by at least 30 reads; the mutation frequency was larger than 1%.
MGMT promoter methylation status detection was carried out on FFPE tissue by methylation-specific polymerase chain reaction (MS-PCR). 200 ng purified DNA was subjected to bisulfite conversion using the EZ-96 DNA Methylation-Lightning MagPrep kit (cat# D5046, Zymo Research, Irvine, CA USA) and amplified using EpiMark Hot Start Taq DNA polymerase (cat#M0490, NEB) in PCR with sense (5’-CCATCTCATCCCTGCGTGTCTCCGACTCAG-Index-GATATGTTGGGATAGTT-3’) and anti-sense (5’- CCTCTCTATGGGCAGTCGGTGAT -CCCAAACACTC
ACCAAAT-3’) primer oligos. The sequences matched to the MGMT promoter region were underlined. The PCR products were purified using Agencourt AMPure XP magnetic beads (cat# A63881, Beckman Coulter, Brea, CA, USA) and subjected to library preparation and sequencing similar to the above in the NGS panel section. Sequenced reads were aligned to NCBI human genome reference assembly hg19 using Bowtie2 (version 2.2.8), Samtools (version 1.3.1), and Bismark version 0.16.3 to identify methylated cytosine.
Multiplex PCR based on the NGS was used to confirm true chromosome (chr) 1p / 19q status. First of all, 10 single nucleotide polymorphisms (SNPs) are selected on 1p and 10 SNP is located in 19q.13.3 section on chr 19q. There are 5 SNP sections located in 19q.13.3 section on chr 19q. These two regions were used to judge whether loss of heterozygosity occurs in common deletion regions, and combined with other SNP. We can judge whether loss of heterozygosity occurs in the whole long arm or severed arm. The heterozygosity distribution of each SNP in the East Asian population refers to the gnomAD_exome_EAS database, and the selection criteria were that the probability that the SNP was heterozygous in the population is close to 50%. The distance of SNP in the genome is more than 300 kb. Then primers and amplification fragments were designed. According to the selected sites, primers suitable for multiple amplification were designed, and the T at the beginning of the second base at the 3' end of each primer was modified to U. The samples were amplified by Platinum Multiplex PCR Master Mix of Thermofisher and treated with uracil DNA glycosylase (UDG) to make the single strand of U form an Apyrimidine (AP) site. Then the library is constructed and sequenced, the Ion torrent of Thermofisher is used for sequencing, and the IonXpress series connector of Ion torrent is used to construct the sequencing library. Finally, the sequencing results are analyzed, and the sequence data were compared with the hg19 sequence to get the allele frequency (AF) percentage of each SNP. In this study, the homozygous type is defined as 95% − 100% or 0–5%, and the heterozygous type is defined as 45% -55%. If it is outside the above range, it is considered that loss of heterozygosity (LOH) has occurred in the SNP. If both 1p and 19q are LOH, it is considered to be chr 1p/19q co-LOH.
Statistical analysis
The SPSS statistical software package (version 25, IBM Corp.) was used for all analyses. Categorical variables are shown as numbers and percentages, while continuous variables are shown as the mean and standard deviation (SD).