Impact Factor 4.848 | CiteScore 3.5
More on impact ›

Original Research ARTICLE

Front. Oncol., 14 August 2020 | https://doi.org/10.3389/fonc.2020.01382

High Copy-Number Variation Burdens in Cranial Meningiomas From Patients With Diverse Clinical Phenotypes Characterized by Hot Genomic Structure Changes

Junpeng Ma1, Yaqiang Hong2, Wei Chen3,4, Da Li1, Kaibing Tian1, Ke Wang1, Yang Yang1, Yuan Zhang1, Yujia Chen1, Lairong Song1, Liangpeng Chen1, Liwei Zhang1,5,6,7, Jiang Du8, Junting Zhang1,6,7, Zhen Wu1,5,6,7*, Dake Zhang3,4* and Liang Wang1,5,6*
  • 1Department of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Beijing, China
  • 2Tsinghua-Peking Center for Life Sciences, School of Life Sciences, Tsinghua University, Beijing, China
  • 3Beijing Advanced Innovation Center for Biomedical Engineering, School of Biological Science and Medical Engineering, Beihang University, Beijing, China
  • 4Key Laboratory of Genomic and Precision Medicine, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, China
  • 5China National Clinical Research Center for Neurological Diseases, Beijing, China
  • 6Center of Brain Tumor, Beijing Institute for Brain Disorders, Beijing, China
  • 7Beijing Key Laboratory of Brain Tumor, Beijing, China
  • 8Department of Neuropathology, Beijing Neurological Institute, Capital Medical University, Beijing, China

Meningiomas, as the most common primary tumor of the central nervous system, are known to harbor genomic aberrations that associate with clinical phenotypes. Here we performed genome-wide genotyping for cranial meningiomas in 383 Chinese patients and identified 9,821 copy-number variations (CNVs). Particularly, patients with diverse clinical features had distinct tumor CNV profiles. CNV burdens were greater in high-grade (WHO grade II and III) samples, recurrent lesions, large tumors (diameter >4.3 cm), and those collected from male patients. Nevertheless, the level of CNV burden did not relate to tumor locations, peritumoral brain edema, bone invasion, or multiple lesions. Overall, the most common tumor CNVs were the copy-number gain (CNG) at 22q11.1 and the copy-number losses (CNLs) at 22q13.2, 14q11.2, 1p34.3, and 1p31.3. Recurrent lesions were featured by the CNLs at 1p31.3, 6q22.31, 9p21.3, and 11p12, and high-grade samples had more CNVs at 4q13.3 and 6q22.31. Meanwhile, large tumors were more likely to have the CNVs at 1p31.3 and 1p34.3. Additionally, recurrence prediction indicated the CNLs at 4p16.3 (p = 0.009, hazard ratio = 5.69) and 10p11.22 (p = 0.037, hazard ratio = 4.53) were candidate independent risk factors.

Introduction

Meningiomas represent the most common primary intracranial tumor type, accounting for 37% of central nervous system neoplasms (1). They are believed to arise from progenitor cells of both the arachnoid cap cells of the arachnoid layer and fibroblasts that reside in the inner dura mater (2). Despite the identification of NF2 mutations or loss of function, recent sequencing studies also revealed mutations involving TRAF7, KLF4, AKT1, SMO, POLR2A, and the ARID1A and TERT promoters of in meningiomas (36). Cytogenetic changes, such as losses of chromosomes 22q, 14q, 1p, 1q, 10q, and 9q, are also commonly reported, some of which are related to tumor progression (712). Copy-number variations (CNVs) of cytobands located at 22q, 1p, and 14q were most common (1214). In tumor development, losses of 6q and 4q have been reported to be significantly associated with high-grade lesions (13). Furthermore, meningiomas in specific locations may have featured CNVs; for instance, those at anterior skull base are likely to have intact chromosome 22q, which loses tumor suppressor gene NF2 (15). However, as a relative benign tumor, meningioma had few data from a relatively large cohort to characterize genome-wide CNV changes, which limits efforts on applying them in tumor progression evaluation, prognosis, and the development of new treatments.

Although maximal but safe resection can cure the majority of meningiomas (16), tumor recurrence still occurs even after gross total resection (GTR) (17). The recurrence status cannot be completely predicted by histopathologic grade alone, as it is mainly based on histopathological characterizations of mitotic rate, cellular features of atypia, and local invasion (18). Meningiomas are well-known for their female-biased predominance (19), but tumors in male patients demonstrate not only a higher annual growth rate (20) but also a higher probability of recurrence (2123). Previous studies have proposed molecular markers for prognostic scoring systems in recent years (14, 21, 2426), and a better WHO classification of meningiomas integrated with independent molecular markers may help to predict the recurrence risk and adjust treatment plans for patients with meningiomas. Although genomic structure changes in neurologic tumors are common, extensive efforts are still required to evaluate roles of diverse recurrent CNVs in the models for tumor classification, prognosis scoring, and recurrence prediction.

To our knowledge, we here collected cranial meningiomas by far at the largest sample size in the Chinese population. We performed genome-wide genotyping for all these samples and identified diverse common CNVs. Along with detailed clinical information, we investigated their relations with gender difference, tumor location, grade classification, and recurrence, and we further proposed candidate predictors for tumor recurrence.

Materials and Methods

Sample Collection

This study was approved by the Institutional Review Board of Beijing Tiantan Hospital affiliated with Capital Medical University. Three hundred and eighty-three frozen meningioma samples were collected at Beijing Tiantan Hospital, Capital Medical University, between August 2008 and August 2017. Signed informed consent forms were acquired from all patients or their guardians before surgery. Tumor specimens from meningioma samples were stored in liquid nitrogen immediately following collection. Genomic DNA was purified from tumor samples using a Biomek 3000 automated workstation with an E.Z.N.A Mag-Bind Tissue DNA Kit (Omega Bio-Tek, Norcross, GA, USA). DNA quality and quantity were determined using a NanoDrop 1000 instrument (Thermo Scientific, Wilmington, DE, USA).

Clinical Data Collection and Follow-Up

Clinical information for 383 patients, including gender, age, primary or recurrent, degree of resection, tumor location, tumor diameter, bone invasion, peritumoral brain edema, pathological subtype, WHO grade, and follow-up results (recurrence and survival), was collected and summarized in Table 1. Pathological diagnosis was reviewed according to the 2016 WHO classification for meningiomas. Tumor recurrence was defined as tumor reemergence after GTR (gross total resection), or tumor regrowth with a minimum change of 25% increase of any tumor diameter after non-GTR based on contrast-enhanced MRIs (27). The degree of resection was decided according to the criteria of Simpson grading and classified as GTR (Simpson grade I to III) or STR (subtotal resection), verified by postoperative magnetic resonance images (MRIs) (28). Recurrence-free survival was defined as the period from the time of present surgery in our hospital to tumor recurrence (or last follow-up visit).

TABLE 1
www.frontiersin.org

Table 1. The clinicopathological features of the entire cohort and subcohort in the prognostic analysis for recurrence.

Whole-Genome Single-Nucleotide Polymorphism (SNP) Genotyping and Statistical Analysis

Whole-genome SNP array analysis for 383 meningioma samples was performed on Illumina Human Infinium CoreExome BeadChips (Illumina, San Diego, CA, USA). Raw intensity values were processed to obtain a normalized B allele frequency (BAF) and a log R ratio (LRR) for each probe using the GenomeStudio Software v2.0.4 (Illumina, San Diego, USA). LRR values were segmented with Genome Alteration Detection Analysis (GADA, Juan R. González, Barcelona, Spain) using parameters of T > 10 and segment lengths containing ≥50 continuous probes. For loss of heterozygosity (LOH) analysis, the sliding window approach was used with a window size of 100 informative SNPs. A window was considered to represent LOH if more than 80% of the SNPs had a minor allele frequency ≤0.9. A segment was defined either as normal or as having one of 3 types of alteration status based on the following criteria: (1) normal, |LRR| < 0.075 and retaining heterozygosity; (2) gain, LRR ≥ 0.075; (3) loss, LRR ≤ −0.075; and (4) copy-number neutral loss of heterozygosity (CNNLOH), |LRR| < 0.075.

To assess genome instability, the genomic fractions of CNVs, and CNNLOH were estimated by dividing the number of SNPs undergoing a specific alteration by the total number of SNPs present in the respective chromosome or in the respective sample. To identify minimal common regions (MCRs) of copy-number gains and losses, the Genomic Identification of Significant Targets in Cancer (GISTIC, Broad Institute, Boston, USA) algorithm was utilized. Thresholds of LRR were set at 0.1 and −0.1 to allow GISTIC to identify amplifications and deletions, respectively. Q-values of minimal common regions <0.01 were defined as significant, and 0.99 was used as the confidence level to determine regions that contained potential driver genes. For genes within candidate CNV markers, the differential gene expression analysis was performed using NCBI GEO2R for the dataset GSE74385 in Gene Expression Omnibus (GEO) (29); the survival analysis for diverse tumors using the gene expression data and clinical information in The Cancer Genome Atlas (TCGA) project via the portal UALCAN (30).

Statistics

Chi-squared tests and Wilcoxon rank-sum tests were performed using R scripts. The two-sided significance level was set at p ≤ 0.01, two-sided. Gene set enrichment analysis (GSEA) was performed using online tools (31), and the two-sided significance level was set at q ≤ 0.05. Prognostic analyses were conducted using Kaplan–Meier analysis and a Cox proportional regression model. Before conducting the prognostic analysis, patients with history of surgery, or with STR, or with postoperative radiotherapy were excluded. The two-sided significance level was set at p ≤ 0.05.

Results

Patient Characteristics

A total of 383 meningioma patients were enrolled in this study, with the female-to-male ratio at 2.24 (265:118), consistent with previous observation. Their average age was 49 years old (18 to 81 years old). According to the 2016 WHO meningioma grading classification, meningiomas of WHO grade I accounted for 86% (331/383) of tumors, among which 89 were from male patients. Meningiomas of WHO grade II constituted 12% (19 in females compared with 27 in males) of tumors, while meningiomas of grade III were only 1.6%, with 4 from females and 2 from males. The median of tumor diameter was 4.30 cm, and 239 meningiomas located in the skull base. In all, 338 patients had primary meningiomas and 45 patients suffered from tumor recurrence. GTR was achieved in 302 patients and STR in 81 patients; meanwhile, 118 patients had peritumoral brain edema, and 57 had bone invasion. For the pathological subtypes of these meningiomas, there were 140 transitional, 117 meningothelial, 66 fibrous, 18 atypical, 14 angiomatous, 6 anaplastic, 5 chondroid, 5 clear cell, and 12 other types. The summarized clinicopathological features are shown in Table 1.

Landscape of CNVs in Meningiomas

A total of 9,821 high-confidence CNVs were identified, and each sample had 26 CNVs on average, including 6,416 gains and 3,405 losses (detailed information for each CNV, Table S1). Their sizes ranged from 298 bp to 198 Mb, with 6,722 over 500 kb and 2,869 over 5 Mb. According to MCRs covered by diverse CNVs (Materials and Methods), we identified 36 common losses (39 kb to 79 Mb) in 27 chromosomes and 28 common gains (2 kb to 20 Mb) in 23 chromosomes in these 383 meningiomas (Figure 1 and Table 2). Copy-number losses (CNLs) were most likely to occur in 22q13.2 (31%), followed by 14q11.2 (29%). Moreover, chromosome 1 was also likely to lose fragments of its short arm, with three common losses of 1p34.3 (21%), 1p31.3 (19%), and 1p22.1 (16%). In comparison, copy-number gains (CNGs) were frequently detected at 22q11.1 (35%), 15q22.2 (16%), 14q11.2 (14%), 10q23.31 (14%), 8p11.22 (14%), and 7p12.3 (14%).

FIGURE 1
www.frontiersin.org

Figure 1. MCR Profiling of CNVs in meningiomas. The peaks in red and blue represent consistent regions of gains and losses in meningiomas based on Genomic Identification of Significant Targets in Cancer (GISTIC), respectively. The GISTIC Q-value is shown at the bottom. The green lines indicate the Q-value (0.01) considered significant in the analysis. Chromosomes are shown in the middle with odd-numbered chromosomes as white and even-numbered chromosomes as gray.

TABLE 2
www.frontiersin.org

Table 2. List of MCRs of CNVs in meningiomas.

Particularly, four CNLs at 1p22.1, 8p11.22, 14q32.2, and 22q13.2 were more common in non-skull-base meningiomas, and a CNG at 22q11.1 more frequently occurred in skull-base meningiomas (Figure 2). These CNVs commonly led to the deletion of 737 genes and the amplification of 146 genes. Pathways over-represented by deleted genes were G2M checkpoint, IL6 JAK STAT3 signaling, TNFA signaling via NFKB, epithelial–mesenchymal transition, inflammatory response, and KRAS signaling (GSEA, Materials and Methods, Table S1).

FIGURE 2
www.frontiersin.org

Figure 2. Landscape of CNVs in meningiomas. The solid bar in the same column represents a single MCR of a CNV with colors of red for gain, blue for loss, and white for normal. Different clinical information is represented at the top, which includes gender, age, WHO grade, death state, recurrent state, primary/recurrent lesions, pathology subtype, tumor location (skull base and non-skull base), bone invasion, peritumoral brain edema, multiple lesions, tumor size (>4.3 cm), and Simpson grade with colors as shown in the figure legends.

High CNV Burdens in Either WHO Grade II and III Meningiomas or Recurrent Lesions Featured by Large CNVs Over 500 kb

In the grade I meningiomas, 331 samples had 7,416 CNVs (22 per sample); in the 46 grade II meningiomas, there were 2,048 CNVs (45 per sample); and 6 grade III meningiomas had 357 CNVs (60 per sample). Therefore, the number of CNVs was similar in grade II and grade III meningiomas (II and III together defined as high grade, p = 0.374, Wilcoxon rank-sum test), which was significantly higher than that in grade I (p = 8.61 × 10−4 and 6.74 × 10−3, respectively, Wilcoxon rank-sum tests, Figure 3A). Moreover, the size of CNVs in high-grade meningioma was larger. The CNVs at 500–1,000 kb (5 vs. 6 vs. 2), 1–5 Mb (14 vs. 21 vs. 6), and >5 Mb (14 vs. 21 vs. 6) were more common in high-grade meningiomas than in grade I (P < 0.01, Wilcoxon rank-sum tests). Meanwhile, meningiomas at three grades had similar number of small CNVs (<500 kb) (grade I, 8; grade II, 11; and grade III, 11). Overall, large CNVs contributed to the higher CNV burden in WHO grade II and III meningiomas. In all, we found 15 common CNVs with differential incidences among different grades of meningiomas (Table 3). Interestingly, two CNLs at 4q13.3 and 6q22.31, both larger than 19 Mb, were most commonly observed in high-grade samples (P = 4.01 × 10−7, 4.05 × 10−10, respectively, Chi-square tests). A CNL at 4q13.3 covered 104 genes, which were over-represented in pathways including inflammatory response, IL6 JAK STAT3 signaling, TNFA signaling via NFKB, epithelial–mesenchymal transition, KRAS signaling up, and angiogenesis (Table S1). For 397 genes affected by the CNL at 6q22.31, the enriched pathways were Hypoxia, IL2 STAT5 signaling, and androgen response (Table S1).

FIGURE 3
www.frontiersin.org

Figure 3. CNV burdens in meningiomas with different clinical features. P-values were calculated with Wilcoxon rank-sum test at each size range. *P < 0.01; **P < 0.001. (A) Average number of CNVs in meningiomas of different WHO grades. (B) Average number of CNVs in meningiomas of primary and recurrent lesions. (C) Average number of CNVs in meningiomas of different tumor sizes. (D) Average number of CNVs in meningiomas of different gender. (E) Average number of CNVs in meningiomas of different positions. (F) Average number of CNVs in meningiomas with or without peritumoral brain edema. (G) Average number of CNVs in meningiomas with or without brain edema. (H) Average number of CNVs in single or multiple lesions.

TABLE 3
www.frontiersin.org

Table 3. MCRs of CNVs differently distributed in meningiomas of different WHO grade, history of surgery, and tumor size.

Recurrent meningiomas had significantly more CNVs (39 per sample) than primary ones (24 per sample, p = 2.08 × 10−3, Wilcoxon rank-sum test, Figure 3B). Meanwhile, recurrent lesions also had more CNVs over 500 kb (31 vs. 16, p = 3.03 × 10−4, Wilcoxon rank-sum test). Among seven common CNVs differentially distributed between recurrent and primary meningiomas (P < 0.01, chi-squared tests, Table 3), four CNVs of 1p31.3, 6q22.31, 9p21.3, and 11p12 were over 500 kb. Chromosome fragment losses at these sites were more common in recurrent meningiomas. For disrupted genes in each site (Sheet 3, Table S1), 6q22.31 covers nearly 400 genes. Other sites have much less genes affected, including four genes in 1p31.3 (USP1, ANGPTL3, ATG4C, and DOCK7), four in 9p21.3 (CDKN2A, CDKN2B, C9orf53, and CDKN2B-AS1), and seven in 11p12 (API5, TTC17, LRRC4C, HNRNPKP3, MIR129-2, MIR670, and LOC100507205).

High CNV Burdens in Meningiomas of Large Diameter and Male Patients

According to the median of tumor diameter (4.30 cm), we grouped these tumors into large (>4.3 cm) and small groups (≤4.3 cm). The large group had 191 samples with 33 CNVs on average, significantly higher than the observation in the small group (192 samples with 18 CNVs on average, P = 7.51 × 10−4, Wilcoxon rank-sum test, Figure 3C). Besides, large meningiomas also had more CNVs over 500 kb (24 vs. 11, p = 4.84 × 10−4, Wilcoxon rank-sum test). Large lesions had more CNVs of 1p31.3 and 1p34.3, and only one CNV, either loss or gain, was over 500 kb locating at 1p31.3 (four genes affected, Table S1).

Meningiomas from male patients had significantly more CNVs (118 samples, 39 CNVs on average) than those in female patients (265 samples, 20 CNVs on average; p = 4.11 × 10−6, Wilcoxon rank-sum test, Figure 3D). Moreover, these CNVs in male samples were larger, more of which were over 200 kb (male: 34 vs. female: 16, p = 3 × 10−6, Wilcoxon rank-sum test). Seven common CNVs showed significant gender difference (P < 0.01, chi-squared tests, Table 4). Five out of them were over 200 kb including the CNVs of 1p22.1, 1p31.3, 1p34.3, 14q23.1, and 19p12. The CNLs at these sites were more common in meningiomas from male patients. For genes affected by these CNVs (Table S1), there were four located in 1p22.1 (DR1, FNBP1L, CCDC18, and LOC100131564), four in 1p31.3 (USP1, ANGPTL3, ATG4C, and DOCK7), four in 14q23.1 (ARID4A, KIAA0586, TIMM9, and TOMM20L), five in 1p34.3 (SFPQ, ZMYM4, ZMYM6, ZMYM1, and ZMYM6NB), and 42 in 19p12.

TABLE 4
www.frontiersin.org

Table 4. MCRs of CNVs differently distributed in meningiomas of different gender, tumor location, and with or without peritumoral brain edema.

The Number of CNVs in Meningiomas Was Independent of Tumor Locations, Peritumoral Brain Edema, Bone Invasion, and Single or Multiple Lesions

Skull-base meningiomas (239 samples, 23 CNVs per sample) had a similar number of CNVs to with non-skull-base lesions (144 samples, 31 CNVs per sample; p = 0.013, Wilcoxon rank-sum test, Figure 3E). Overall, extremely large CNVs (>5 Mb) were more likely to present in non-skull-base meningiomas (p = 7.18 × 10−5, Wilcoxon rank-sum test). Nevertheless, a CNG at 22q11.1 of this type was preponderant in skull-base meningiomas (p = 3.45 × 10−5, chi-squared test, Table 4), affecting 27 genes in the region (Table S1).

Patients present with peritumoral brain edema (118 samples, 33 CNVs on average) showed no significantly difference in number of CNVs with those without (265 samples, 22 CNVs on average; P > 0.01, Wilcoxon rank-sum test, Figure 3F). However, more large CNVs (500 kb−1 Mb and >5 Mb) were observed in meningiomas with peritumoral brain edema (P = 3.48 × 10−3, 5.29 × 10−3, respectively, Wilcoxon rank-sum tests). The featured one with most significance located at 9p21.3, which was a CNL covering four genes (CDKN2A, CDKN2B, C9orf53, and CDKN2B-AS1) (p = 1.96 × 10−3, chi-square test; Table 4 and Table S1).

In tumors with or without bone invasions (P = 0.597) or single or multiple lesions (P = 0.869), the CNV burdens were similar (Wilcoxon rank-sum tests; Figures 3G,H). No CNVs were more prevalent in meningiomas with bone invasions. Nevertheless, a CNG at 10q23.31 had a higher incidence in multiple lesions (multiple, 33% vs. single, 12%; p = 3.49 × 10−3, chi-squared test, Table 5), which only covers one gene KIF20B.

TABLE 5
www.frontiersin.org

Table 5. The MCRs of CNVs differently distributed in meningiomas of single or multiple lesions.

Identification of Independently Significant Prognostic CNVs in Predicting Tumor Recurrence

Based on common CNVs, we tried to predict the tumor recurrence. After excluding patients with recurrent lesions, with subtotal resection, or having postoperative radiotherapy, 267 patients were included for further prognostic analysis, and the detailed clinicopathological features of this subcohort are shown in Table 1. In the follow-up (mean period, 60 months), 12 patients suffered from tumor recurrence. All common CNV regions and clinical features were included in univariate Cox analysis of tumor recurrence. As shown in Table 6, skull-base lesions (p = 0.040), loss of 1p22.1 (p = 0.039), 1p34.3 (p = 0.024), 4q13.3 (p = 0.029), 4p16.3 (p = 0.001), 7q11.21 (p = 0.015), 10p11.22 (p = 0.003), 14q23.1 (p = 0.032), 19q13.12 (p = 0.013), and 19p12 (p = 0.01) were significant risk factors for tumor recurrence. In particular, most significant independent risk factors for recurrence were loss of 4p16.3 (p = 0.009, HR = 5.69, multivariate Cox analysis) and 10p11.22 (p = 0.037, HR = 4.53). As shown in Figure 4, patients with losses of both 4p16.3 and 10p11.22 were more likely to suffer from tumor recurrence than patients with loss of either one, or patients with neither of these CNV changes. Calculated by Cox analysis, the hazard ratio (HR) increased by 5.10 (95% CI: 2.35–11.08, p = 3.7 × 10−5) for each additional prognostic CNV.

TABLE 6
www.frontiersin.org

Table 6. Significant factors for tumor recurrence of meningiomas in subcohort for prognostic analysis.

FIGURE 4
www.frontiersin.org

Figure 4. Kaplan–Meier plots of combined losses involving the CNLs at 4p16.3 and 10p11.22 as risk factor of tumor recurrence. Log-rank test was used. Both, with CNLs identified at both sites; Either, with CNLs identified at only one of two sites; None, no CNLs at both sites.

Eight genes located within these two CNLs: ZNF141, ABCA11P, ZNF595, ZNF721, ZNF718, ZNF876P, ZNF732 in 4p16.3, and CCDC7 in 10p11.22 (Table S1). In differential gene expression analysis between non-recurrent (grade I, 13 and grade II, 6) and recurrent (grade I, 7 and grade II, 8) lesions in a public gene expression dataset of meningioma (Method), ZNF141 and ZNF595 showed the tendency to have lower expression levels in recurrent samples (unadjusted p < 0.05, Table S1). We further examined the effects of expression levels of these eight genes on survival time for 31 tumor types in TCGA (Method) and identified 28 associations with significance (p < 0.05, Figure S1). All of these genes had lower expression levels in certain tumor types from patients with shorter survival, which indicated their decreased functions related to malignant phenotypes. Particularly, most of these genes (five out of eight) had the same effects on patient survival in head and neck squamous cell carcinoma (HNSC), and the low expression of ZNF718 (p = 0.0027), CCDC7 (p = 0.01), ZNF141 (p = 0.012), ZNF721 (p = 0.029), and ZNF732 (p = 0.045) all demonstrated significant associations with the shorter survival of patients.

Discussion

In the present study, we clarified the CNV characteristics of cranial meningiomas in 383 Chinese patients. Particularly, we compared the CNV burdens of meningiomas in diverse phenotypes. We found more CNVs in the samples of high-grades, recurrent lesions, tumor diameter over 4.3 cm, and samples from male patients. Meanwhile, CNV burden may not relate to tumor locations, peritumoral brain edema, bone invasion, and multiple lesions. Moreover, we also identified featured CNVs in each clinical group. Besides, we found two candidates as independent prognostic CNVs in predicting tumor recurrence.

Based on a relatively large cohort of cranial meningiomas, we observed that CNVs of 22q (61%), 14q (54%), and 1p (38%) were the most prevalent, followed by 15q (32%), 6p (30%), 8p (29%), 10q (29%), and 1q (25%). In previous studies, CNVs of 22q, 14q, and 1p are always among the most frequent CNVs of meningiomas (7, 10, 12, 14, 15, 24, 32), which is in accordance with our observation. For instance, a recent study reported the top three CNVs in their samples located at 1p (71%), 22q (64%), and 14q (42%) (14). By far, no differences were noted among different ethnic groups.

The CNVs frequently identified in patients with distinct clinical features hold clues for further functional studies. For instance, two CNLs of 1p31.3 and 1p34.3, commonly seen in meningiomas of high-grade and recurrent or large lesions, contain a lot of genes with functional importance. At 1p34.3, the SFPQ gene participates in transcriptional regulation, DNA double-strand break repairs, and suppression of RNA:DNA-hybrid-related telomere instability (33, 34). At 1p31.3, the USP1 gene, involved in multiple DNA repair pathways, can function as a key senescence regulator controlling genomic integrity (35); autophagy protein ATG4C participates in controlling the unregulated cell growth (36). Reduced levels of autophagy have been described as being linked to malignant tumors (37). Functional changes related to these genes may also contribute to the progression of meningiomas, which needs further studies for validation.

The CNG at 10q23.31 was the only CNV more commonly seen in multiple meningiomas rather than in single lesions. It covers only one gene, KIF20B, an oncogene involved in cytokinesis. A recent study suggested to target the KIF20B gene in the treatment for hepatocellular carcinoma (38). Inhibition of KIF20B can block mitosis at both metaphase and telophase, which enhance the cytotoxicity of two chemotherapeutic drugs, hydroxycamptothecin, and mitomycin C (39). The role of KIF20B in tumorigenesis of meningiomas, especially multiple lesions, suggests that its suppression might be a novel strategy in the treatment for multiple meningiomas in the future. Moreover, the CNG at 6p21.33, more frequently found in lesions from male patients, is where HLA-B and HLA-C are located, indicating the existence of immune factors underlying gender difference of meningioma occurrence; the CNG at 20q13.33, more frequently identified in patients with peritumoral brain edema, covers the SYCP2 gene, which is related to the depth of cervical invasion in squamous cell carcinoma (40).

Tumor recurrence is an important issue for patients with meningiomas, and patients with meningiomas prone to recurrence need adjuvant radiation therapy or close follow-up. Meanwhile, patients with low risk of tumor recurrence could be spared from the toxicity of radiation therapy. Nevertheless, these patients are not accurately identified by WHO grading (41). Here, we demonstrated the potential of CNV profiling in recurrence prediction. Loss of 1, 4, 9, and 10p and gain of 1q or other chromosomal regions have been revealed to be risk factors for tumor recurrence in previous studies (711, 14, 21, 24). In our observation, the CNLs of 4p16.3 and 10p11.22 were independent risk factors for cranial meningioma recurrence. The CNL of 4p16.3 covers MiR-571, ABCA11P, ZNF141, ZNF595, ZNF718, ZNF721, ZNF732, and ZNF876P. A recent study identified miR-571 as the first miRNA that prevents aberrant DNA replication, and the Cdk2-c-Myc-miR-571 axis was identified as a new pathway for regulating DNA replication, cell cycle, and genomic stability in cancer cells (42). As a result, loss of miR-571 may lead to genomic instability. Although some studies have reported differential expression or mutation occurrence of ZNF595 (43), ZNF721 (44), ZNF718 (45), and ZNF141 (46), their functions remain unclear. Besides, potential roles of ABCA11P, ZNF732, and ZNF876P are novel in meningioma recurrence. In the 10p11.22, CCDC7, also known as Biot2, highly expressed in CD133-positive stem cells, functions as a risk factor for poor prognosis in colorectal cancer (47, 48). In our study, the CNL at 10p11.22 (CCDC7) was an independent risk factor of tumor recurrence, and the underlying mechanisms need further investigation.

The cross-sectional analysis in the entire cohort compared primary and recurrent lesions from different groups of patients, and some primary tumors may also harbor CNVs contributing to tumor recurrence. It may undermine the ability to identify CNVs related to recurrence, which may explain the missing of the CNLs of 4p16.3 and 10p11.22 in the comparison. Meanwhile, the comparison results may be cofounded by differential CNVs present in the early stage of tumor development between two groups of patients. Therefore, the follow-up study provides us an opportunity to identify those CNVs related to recurrence. The recurrence rates of patients with these two CNVs were over 20% (loss of 4p16.3, 21%, 4/19; loss of 10p11.22, 27%, 3/11), significantly higher than the recurrent rate (about 4%) in patients without them. Nevertheless, only 12 patients (4.5%, 12/267) had tumor recurrence during a mean follow-up period of 5 years in our subcohort for prognostic analysis. Although it is similar to previous observation, which is 3% for WHO grade I meningiomas and 30% for WHO grade II meningiomas in GTR patients (28), the prediction power of these two candidate markers requires further evaluation in a larger group of patients with tumor recurrence in the future follow-up. Besides, the recurrence factors may have heterogeneity, and 4p16.3 and 10p11.22 together accounted for the 13% (loss of 4p16.3, 4/45; loss of 10p11.22, 4/45, losses of both, 2/45) recurrent lesions in the cross-sectional analysis. It needs further efforts to dissect other CNVs related to tumor recurrence.

Conclusions

Based on a large number of patients with cranial meningiomas, we identified that the CNVs of 22q, 1p, and 14q were the most prevalent. Meningiomas of high WHO grades, recurrent tumors, large size, and male gender were likely to have more CNVs, especially of large size (>500 kB). Additionally, the CNLs at 4p16.3 and 10p11.22 were promising candidates as independent risk factors for tumor recurrence prediction.

Data Availability Statement

The datasets in this study have been publicly deposited. They can be accessed at: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE147673.

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Beijing Tiantan Hospital Affiliated with Capital Medical University. Signed informed consent forms were acquired from all patients or their guardians.

Author Contributions

LW, DZ, JM, YH, ZW, JZ, LZ, and WC designed this study. JM, YH, KT, DL, KW, YY, and JD conducted the experiments and acquired data. JM, YZ, YC, LS, and LC conducted the follow-up and acquired data. JM and YH analyzed and interpreted data. JM, YH, KT, DL, KW, YY, JD, YZ, YC, LS, and LC drafted the manuscript. LW, DZ, ZW, JZ, WC, and LZ revised the manuscript. All authors have seen and approved the manuscript.

Funding

This work was funded by the Beijing Municipal Natural Science Foundation (7192056 and J180005, ZW), the Innovation Promotion Association CAS (2016098, DZ), and the Outstanding Talent Training Projects Foundation of Beijing (2017000021469G261, DL and 2017000021469G259, YY).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to the efforts from Junhua Li on the specimen storage.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2020.01382/full#supplementary-material

References

1. Kruchko C, Ostrom QT, Gittleman H, Barnholtz-Sloan JS. The CBTRUS story: providing accurate population-based statistics on brain and other central nervous system tumors for everyone. Neuro Oncol. (2018) 20:295–8. doi: 10.1093/neuonc/noy006

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Bi WL, Abedalthagafi M, Horowitz P, Agarwalla PK, Mei Y, Aizer AA, et al. Genomic landscape of intracranial meningiomas. J Neurosurg. (2016) 125:525–35. doi: 10.3171/2015.6.JNS15591

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Clark VE, Erson-Omay EZ, Serin A, Yin J, Cotney J, Ozduman K, et al. Genomic analysis of non-NF2 meningiomas reveals mutations in TRAF7, KLF4, AKT1, and SMO. Science. (2013) 339:1077–80. doi: 10.1126/science.1233009

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Abedalthagafi MS, Bi WL, Merrill PH, Gibson WJ, Rose MF, Du Z, et al. ARID1A and TERT promoter mutations in dedifferentiated meningioma. Cancer Genet. (2015) 208:345–50. doi: 10.1016/j.cancergen.2015.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Clark VE, Harmanci AS, Bai H, Youngblood MW, Lee TI, Baranoski JF, et al. Recurrent somatic mutations in POLR2A define a distinct subset of meningiomas. Nat Genet. (2016) 48:1253–9. doi: 10.1038/ng.3651

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sahm F, Schrimpf D, Olar A, Koelsche C, Reuss D, Bissel J, et al. TERT promoter mutations and risk of recurrence in meningioma. J Natl Cancer Inst. (2016) 108:djv377. doi: 10.1093/jnci/djv377

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Bello MJ, De Campos JM, Kusak ME, Vaquero J, Sarasa JL, Pestana A, et al. Allelic loss at 1p is associated with tumor progression of meningiomas. Genes Chromosomes Cancer. (1994) 9:296–8. doi: 10.1002/gcc.2870090411

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Lamszus K, Kluwe L, Matschke J, Meissner H, Laas R, Westphal M. Allelic losses at 1p, 9q, 10q, 14q, and 22q in the progression of aggressive meningiomas and undifferentiated meningeal sarcomas. Cancer Genet Cytogenet. (1999) 110:103–10. doi: 10.1016/S0165-4608(98)00209-X

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Perry A, Banerjee R, Lohse CM, Kleinschmidt-Demasters BK, Scheithauer BW. A role for chromosome 9p21 deletions in the malignant progression of meningiomas and the prognosis of anaplastic meningiomas. Brain Pathol. (2002) 12:183–90. doi: 10.1111/j.1750-3639.2002.tb00433.x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Maillo A, Orfao A, Espinosa AB, Sayagues JM, Merino M, Sousa P, et al. Early recurrences in histologically benign/grade I meningiomas are associated with large tumors and coexistence of monosomy 14 and del(1p36) in the ancestral tumor cell clone. Neuro Oncol. (2007) 9:438–46. doi: 10.1215/15228517-2007-026

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Jansen M, Mohapatra G, Betensky RA, Keohane C, Louis DN. Gain of chromosome arm 1q in atypical meningioma correlates with shorter progression-free survival. Neuropathol Appl Neurobiol. (2012) 38:213–9. doi: 10.1111/j.1365-2990.2011.01222.x

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Mcnulty SN, Schwetye K, Goldstein M, Carter J, Schmidt RE, Ansstas G, et al. Analysis of point mutations and copy number variation in grade II and III meningioma. Exp Mol Pathol. (2018) 105:328–33. doi: 10.1016/j.yexmp.2018.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Lee Y, Liu J, Patel S, Cloughesy T, Lai A, Farooqi H, et al. Genomic landscape of meningiomas. Brain Pathol. (2010) 20:751–62. doi: 10.1111/j.1750-3639.2009.00356.x

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Olar A, Wani KM, Wilson CD, Zadeh G, Demonte F, Jones DT, et al. Global epigenetic profiling identifies methylation subgroups associated with recurrence-free survival in meningioma. Acta Neuropathol. (2017) 133:431–44. doi: 10.1007/s00401-017-1678-x

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Kros J, De Greve K, Van Tilborg A, Hop W, Pieterman H, Avezaat C, et al. NF2 status of meningiomas is associated with tumour localization and histology. J Pathol. (2001) 194:367–72. doi: 10.1002/path.909

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Aghi MK, Carter BS, Cosgrove GR, Ojemann RG, Amin-Hanjani S, Martuza RL, et al. Long-term recurrence rates of atypical meningiomas after gross total resection with or without postoperative adjuvant radiation. Neurosurgery. (2009) 64:56–60; discussion 60. doi: 10.1227/01.NEU.0000330399.55586.63

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Marciscano AE, Stemmer-Rachamimov AO, Niemierko A, Larvie M, Curry WT, Barker FG II, et al. Benign meningiomas (WHO grade I) with atypical histological features: correlation of histopathological features with clinical outcomes. J Neurosurg. (2016) 124:106–14. doi: 10.3171/2015.1.JNS142228

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Viaene AN, Zhang B, Martinez-Lage M, Xiang C, Tosi U, Thawani JP, et al. Transcriptome signatures associated with meningioma progression. Acta Neuropathol Commun. (2019) 7:67. doi: 10.1186/s40478-019-0690-x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ostrom QT, Gittleman H, Fulop J, Liu M, Blanda R, Kromer C, et al. CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2008-2012. Neuro Oncol. (2015) 17(Suppl. 4):iv1–62. doi: 10.1093/neuonc/nov189

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Oya S, Kim SH, Sade B, Lee JH. The natural history of intracranial meningiomas. J Neurosurg. (2011) 114:1250–6. doi: 10.3171/2010.12.JNS101623

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Domingues PH, Sousa P, Otero A, Goncalves JM, Ruiz L, De Oliveira C, et al. Proposal for a new risk stratification classification for meningioma based on patient age, WHO tumor grade, size, localization, and karyotype. Neuro Oncol. (2014) 16:735–47. doi: 10.1093/neuonc/not325

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Hua L, Zhu H, Deng J, Tian M, Jiang X, Tang H, et al. Clinical and prognostic features of spinal meningioma: a thorough analysis from a single neurosurgical center. J Neurooncol. (2018) 140:639–47. doi: 10.1007/s11060-018-2993-3

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Juratli TA, Mccabe D, Nayyar N, Williams EA, Silverman IM, Tummala SS, et al. DMD genomic deletions characterize a subset of progressive/higher-grade meningiomas with poor outcome. Acta Neuropathol. (2018) 136:779–92. doi: 10.1007/s00401-018-1899-7

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ketter R, Henn W, Niedermayer I, Steilen-Gimbel H, Konig J, Zang KD, et al. Predictive value of progression-associated chromosomal aberrations for the prognosis of meningiomas: a retrospective study of 198 cases. J Neurosurg. (2001) 95:601–7. doi: 10.3171/jns.2001.95.4.0601

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kishida Y, Natsume A, Kondo Y, Takeuchi I, An B, Okamoto Y, et al. Epigenetic subclassification of meningiomas based on genome-wide DNA methylation analyses. Carcinogenesis. (2012) 33:436–41. doi: 10.1093/carcin/bgr260

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Sahm F, Schrimpf D, Stichel D, Jones DTW, Hielscher T, Schefzyk S, et al. DNA methylation-based classification and grading system for meningioma: a multicentre, retrospective analysis. Lancet Oncol. (2017) 18:682–94. doi: 10.1016/S1470-2045(17)30155-9

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Huang RY, Bi WL, Weller M, Kaley T, Blakeley J, Dunn I, et al. Proposed response assessment and endpoints for meningioma clinical trials: report from the response assessment in neuro-oncology working group. Neuro Oncol. (2019) 21:26–36. doi: 10.1093/neuonc/noy137

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Goldbrunner R, Minniti G, Preusser M, Jenkinson MD, Sallabanda K, Houdart E, et al. EANO guidelines for the diagnosis and treatment of meningiomas. Lancet Oncol. (2016) 17:e383–91. doi: 10.1016/S1470-2045(16)30321-7

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Schmidt M, Mock A, Jungk C, Sahm F, Ull AT, Warta R, et al. Transcriptomic analysis of aggressive meningiomas identifies PTTG1 and LEPR as prognostic biomarkers independent of WHO grade. Oncotarget. (2016) 7:14551–68. doi: 10.18632/oncotarget.7396

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravathi BVSK, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia. (2017) 19:649–58. doi: 10.1016/j.neo.2017.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Bi WL, Zhang M, Wu WW, Mei Y, Dunn IF. Meningioma genomics: diagnostic, prognostic, and therapeutic applications. Front Surg. (2016) 3:40. doi: 10.3389/fsurg.2016.00040

PubMed Abstract | CrossRef Full Text | Google Scholar

33. De Silva HC, Lin MZ, Phillips L, Martin JL, Baxter RC. IGFBP-3 interacts with NONO and SFPQ in PARP-dependent DNA damage repair in triple-negative breast cancer. Cell Mol Life Sci. (2019) 76:2015–30. doi: 10.1007/s00018-019-03033-4

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Petti E, Buemi V, Zappone A, Schillaci O, Broccia PV, Dinami R, et al. SFPQ and NONO suppress RNA:DNA-hybrid-related telomere instability. Nat Commun. (2019) 10:1001. doi: 10.1038/s41467-019-08863-1

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Ogrunc M, Martinez-Zamudio RI, Sadoun PB, Dore G, Schwerer H, Pasero P, et al. USP1 regulates cellular senescence by controlling genomic integrity. Cell Rep. (2016) 15:1401–11. doi: 10.1016/j.celrep.2016.04.033

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Korkmaz G, Le Sage C, Tekirdag KA, Agami R, Gozuacik D. miR-376b controls starvation and mTOR inhibition-related autophagy by targeting ATG4C and BECN1. Autophagy. (2012) 8:165–76. doi: 10.4161/auto.8.2.18351

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Long M, Mcwilliams TG. Monitoring autophagy in cancer: from bench to bedside. Semin Cancer Biol. (2019). doi: 10.1016/j.semcancer.2019.05.016. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Liu X, Zhou Y, Liu X, Peng A, Gong H, Huang L, et al. MPHOSPH1: a potential therapeutic target for hepatocellular carcinoma. Cancer Res. (2014) 74:6623–34. doi: 10.1158/0008-5472.CAN-14-1279

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liu X, Li Y, Zhang X, Liu XY, Peng A, Chen Y, et al. Inhibition of kinesin family member 20B sensitizes hepatocellular carcinoma cell to microtubule-targeting agents by blocking cytokinesis. Cancer Sci. (2018) 109:3450–60. doi: 10.1111/cas.13794

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Guo P, Wang D, Wu J, Yang J, Ren T, Zhu B, et al. The landscape of alternative splicing in cervical squamous cell carcinoma. Onco Targets Ther. (2015) 8:73–9. doi: 10.2147/OTT.S72832

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Suppiah S, Nassiri F, Bi WL, Dunn IF, Hanemann CO, Horbinski CM, et al. Molecular and translational advances in meningiomas. Neuro Oncol. (2019) 21:i4–17. doi: 10.1093/neuonc/noy178

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Zhang Y, Li Z, Hao Q, Tan W, Sun J, Li J, et al. The Cdk2-c-Myc-miR-571 axis regulates DNA replication and genomic stability by targeting geminin. Cancer Res. (2019) 79:4896–910. doi: 10.1158/0008-5472.CAN-19-0020

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Cui J, Yin Y, Ma Q, Wang G, Olman V, Zhang Y, et al. Comprehensive characterization of the genomic alterations in human gastric cancer. Int J Cancer. (2015) 137:86–95. doi: 10.1002/ijc.29352

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Shi X, Li Y, Yan P, Shi Y, Lai J. Weighted gene co-expression network analysis to explore the mechanism of heroin addiction in human nucleus accumbens. J Cell Biochem. (2020) 121:1870–9. doi: 10.1002/jcb.29422

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Arathimos R, Sharp GC, Granell R, Tilling K, Relton CL. Associations of sex hormone-binding globulin and testosterone with genome-wide DNA methylation. BMC Genet. (2018) 19:113. doi: 10.1186/s12863-018-0703-y

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Kalsoom UE, Klopocki E, Wasif N, Tariq M, Khan S, Hecht J, et al. Whole exome sequencing identified a novel zinc-finger gene ZNF141 associated with autosomal recessive postaxial polydactyly type A. J Med Genet. (2013) 50:47–53. doi: 10.1136/jmedgenet-2012-101219

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Shen YM, Arbman G, Sandstrom P, Gullstrand P, Wei YQ, Zhang H, et al. Novel gene hBiot2 is an independent prognostic factor in colorectal cancer patients. Oncol Rep. (2012) 27:376–82. doi: 10.3892/or.2011.1521

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Kim ST, Sohn I, Do IG, Jang J, Kim SH, Jung IH, et al. Transcriptome analysis of CD133-positive stem cells and prognostic value of survivin in colorectal cancer. Cancer Genomics Proteomics. (2014) 11:259–66. doi: 10.1016/S0959-8049(14)70327-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: copy number variation, female prominence, multiple meningiomas, oncogenic driver, recurrence, tumor location

Citation: Ma J, Hong Y, Chen W, Li D, Tian K, Wang K, Yang Y, Zhang Y, Chen Y, Song L, Chen L, Zhang L, Du J, Zhang J, Wu Z, Zhang D and Wang L (2020) High Copy-Number Variation Burdens in Cranial Meningiomas From Patients With Diverse Clinical Phenotypes Characterized by Hot Genomic Structure Changes. Front. Oncol. 10:1382. doi: 10.3389/fonc.2020.01382

Received: 06 January 2020; Accepted: 30 June 2020;
Published: 14 August 2020.

Edited by:

Hailiang Tang, Huashan Hospital Affiliated to Fudan University, China

Reviewed by:

Hua He, Shanghai Changzheng Hospital, China
Shi Jinlon. G, Affiliated Hospital of Nantong University, China

Copyright © 2020 Ma, Hong, Chen, Li, Tian, Wang, Yang, Zhang, Chen, Song, Chen, Zhang, Du, Zhang, Wu, Zhang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhen Wu, wuzhen1966@aliyun.com; Dake Zhang, dakezhang@gmail.com; Liang Wang, saintage7@126.com

These authors have contributed equally to this work