Identification and comparison of m6A modifications in glioblastoma non-coding RNAs with MeRIP-seq and Nanopore dRNA-seq

The most prominent RNA modification - N6-methyladenosine (m6A) - affects gene regulation and cancer progression. The extent and effect of m6A on long non-coding RNAs (lncRNAs) is, however, still not clear. The most established method for m6A detection is methylated RNA immunoprecipitation and sequencing (MeRIP-seq). However, Oxford Nanopore Technologies recently developed direct RNA-seq (dRNA-seq) method, allowing m6A identification at higher resolution and in its native form. We performed whole transcriptome sequencing of the glioblastoma cell line U87-MG with both MeRIP-seq and dRNA-seq. For MeRIP-seq, m6A peaks were identified using nf-core/chipseq, and for dRNA-seq - EpiNano pipeline. MeRIP-seq analysis revealed 5086 lncRNAs transcripts, while dRNA-seq identified 336 lncRNAs transcripts from which 556 and 198 were found to be m6A modified, respectively. While 24 lncRNAs with m6A overlapped between two methods. Gliovis database analysis revealed that the expression of the major part of identified overlapping lncRNAs was associated with glioma grade or patient survival prognosis. We found that the frequency of m6A occurrence in lncRNAs varied more than 9-fold throughout the provided list of 24 modified lncRNAs. The highest m6A frequency was detected in MIR1915HG, THAP9-AS1, MALAT1, NORAD1, and NEAT1 (49-88nt), while MIR99AHG, SNHG3, LOXL1-AS1, ILF3-DT showed the lowest m6A frequency (445-261nt). Taken together, (1) we provide a high accuracy list of 24 m6A modified lncRNAs of U87-MG cells; (2) we conclude that MeRIP-seq is more suitable for an initial m6A screening study, due to its higher lncRNA coverage, whereas dRNA-seq is most useful when more in-depth analysis of m6A quantity and precise location is of interest.Abbreviations: (dRNA-seq) direct RNA-seq, (GBM) glioblastoma, (LGG) low-grade glioma, (lncRNAs) long non-coding RNAs, (m6A) N6-methyladenosine, (MeRIP-seq) methylated RNA immunoprecipitation and sequencing, (ncRNA) non-coding RNA, (ONT) Oxford Nanopore Technologi; Lietuvos Mokslo Taryba.

Effect of Organic Selenium on the Homeostasis of Trace Elements, Lipid Peroxidation, and mRNA Expression of Antioxidant Proteins in Mouse Organs

(1) In this study we determined the effect of long-term selenomethionine administration on the oxidative stress level and changes in antioxidant protein/enzyme activity; mRNA expression; and the levels of iron, zinc, and copper. (2) Experiments were performed on 4-6-week-old BALB/c mice, which were given selenomethionine (0.4 mg Se/kg b.w.) solution for 8 weeks. The element concentration was determined via inductively coupled plasma mass spectrometry. mRNA expression of SelenoP, Cat, and Sod1 was quantified using real-time quantitative reverse transcription. Malondialdehyde content and catalase activity were determined spectrophotometrically. (3) After long-term SeMet administration, the amount of Se increased by 12-fold in mouse blood, 15-fold in the liver, and 42-fold in the brain, as compared to that in the control. Exposure to SeMet decreased amounts of Fe and Cu in blood, but increased Fe and Zn levels in the liver and increased the levels of all examined elements in the brain. Se increased malondialdehyde content in the blood and brain but decreased it in liver. SeMet administration increased the mRNA expression of selenoprotein P, dismutase, and catalase, but decreased catalase activity in brain and liver. (4) Eight-week-long selenomethionine consumption elevated Se levels in the blood, liver, and especially in the brain and disturbed the homeostasis of Fe, Zn, and Cu. Moreover, Se induced lipid peroxidation in the blood and brain, but not in the liver. In response to SeMet exposure, significant up-regulation of the mRNA expression of catalase, superoxide dismutase 1, and selenoprotein P in the brain, and especially in the liver, was determined.

Radiomic features of amygdala nuclei and hippocampus subfields help to predict subthalamic deep brain stimulation motor outcomes for Parkinson‘s disease patients

Background and purpose

The aim of the study is to predict the subthalamic nucleus (STN) deep brain stimulation (DBS) outcomes for Parkinson’s disease (PD) patients using the radiomic features extracted from pre-operative magnetic resonance images (MRI).

Methods

The study included 34 PD patients who underwent DBS implantation in the STN. Five patients (15%) showed poor DBS motor outcome. All together 9 amygdalar nuclei and 12 hippocampus subfields were segmented using Freesurfer 7.0 pipeline from pre-operative MRI images. Furthermore, PyRadiomics platform was used to extract 120 radiomic features for each nuclei and subfield resulting in 5,040 features. Minimum Redundancy Maximum Relevance (mRMR) feature selection method was employed to reduce the number of features to 20, and 8 machine learning methods (regularized binary logistic regression (LR), decision tree classifier (DT), linear discriminant analysis (LDA), naive Bayes classifier (NB), kernel support vector machine (SVM), deep feed-forward neural network (DNN), one-class support vector machine (OC-SVM), feed-forward neural network-based autoencoder for anomaly detection (DNN-A)) were applied to build the models for poor vs. good and very good STN-DBS motor outcome prediction.

Results

The highest mean prediction accuracy was obtained using regularized LR (96.65 ± 7.24%, AUC 0.98 ± 0.06) and DNN (87.25 ± 14.80%, AUC 0.87 ± 0.18).

Conclusion

The results show the potential power of the radiomic features extracted from hippocampus and amygdala MRI in the prediction of STN-DBS motor outcomes for PD patients.

The Role of CASC2 and miR-21 Interplay in Glioma Malignancy and Patient Outcome

Recently long non-coding RNAs (lncRNAs) were highlighted for their regulatory role in tumor biology. The novel human lncRNA cancer susceptibility candidate 2 (CASC2) has been characterized as a potential tumor suppressor in several tumor types. However, the roles of CASC2 and its interplay with miR-21 in different malignancy grade patient gliomas remain unexplored. Here we screened 99 different malignancy grade astrocytomas for CASC2, and miR-21 gene expression by real-time quantitative polymerase chain reaction (RT-qPCR) in isocitrate dehydrogenase 1 (IDH1) and O-6-methylguanine methyltransferase (MGMT) assessed gliomas. CASC2 expression was significantly downregulated in glioblastomas (p = 0.0003). Gliomas with low CASC2 expression exhibited a high level of miR-21, which was highly associated with the higher glioma grade (p = 0.0001), IDH1 wild type gliomas (p < 0.0001), and poor patient survival (p < 0.001). Taken together, these observations suggest that CASC2 acts as a tumor suppressor and potentially as a competing endogenous RNA (ceRNA) for miR-21, plays important role in IDH1 wild type glioma pathogenesis and patients' outcomes.

Significance of Amphiregulin (AREG) for the Outcome of Low and High Grade Astrocytoma Patients

Background: Amphiregulin (AREG) is one of the ligands of the epidermal growth factor receptor which levels was shown to have a tight coherence with various types of cancer. AREG was also designated to be a promising marker for several types of cancer however precious little data about AREG role in the most frequent and generally lethal human brain tumours - astrocytomas reported up to date. The aim of the study was to investigate how AREG changes at epigenetic and expression levels reflect on astrocytoma malignancy and patient outcome. Methods: In total 205 low and high grade astrocytoma samples (15 pilocytic astrocytomas, 56 diffuse astrocytomas, 32 anaplastic astrocytomas and 102 glioblastomas) were used for target mRNA, protein expression and DNA methylation analysis applying qRT-PCR, Western-Blot and MS-PCR assays, respectively. Results: Present research revealed that AREG expression level and methylation in cancer tissue is dependent on the grade of astrocytoma. GBM tissue disclosed elevated AREG mRNA expression but reduced AREG protein level as compared to grade II and grade III astrocytomas (p<0.001). Increased methylation frequency was also more abundant in GBM (74%) than grade I, II and III astrocytomas (25%, 34%, and 36%, respectively). The survival analysis revealed relevant differences in patient overall survival between AREG methylation, mRNA and protein expression groups. Kaplan-Meier analysis encompassing only malignant tumours showed similar results indicating that AREG is associated with astrocytoma patient survival independently from astrocytoma grade. Conclusions: Current findings demonstrate that AREG appearance is associated with patient survival as well as astrocytomas malignancy indicating its influence on tumour progression and suggest its applicability as a promising marker.

GFAP expression is influenced by astrocytoma grade and rs2070935 polymorphism

Glial fibrillary acidic protein (GFAP) is an intermediate filament that provides mechanical support to astrocytes. Rs2070935 is a single nucleotide polymorphism (SNP) located in the promoter region of the GFAP gene. The aim of this pilot study is to investigate GFAP expression at mRNA, protein levels and rs2070935 polymorphism in 50 different grade human astrocytoma samples. GFAP expression at mRNA level was measured using quantitative reverse transcription polymerase chain reaction (qRT-PCR) with SYBR Green dye, whereas the translational activity of the following gene was detected using western blot assay. Furthermore, genotypes of rs2070935 were identified using qPCR with TaqMan probes. As a result, GFAP mRNA and protein expression was found to be declining with increasing astrocytoma grade (p < 0.05). A tendency was observed between increased GFAP mRNA expression and shorter grade IV astrocytoma patient survival (p = 0.2117). The rs2070935 CC genotype was found to be associated with increased GFAP translational activity in grade II astrocytoma (p = 0.0238). Possible links between rs2070935 genotypes and alternative splicing of GFAP were also observed. The rs2070935 AA genotype was found to be associated with poor clinical outcome for grade IV astrocytoma patients (p = 0.0007), although the following data should be checked in a larger sample size of astrocytoma patients.

N-myc downstream-regulated gene 2 (NDRG2) promoter methylation and expression in pituitary adenoma

BACKGROUND

Pituitary adenoma (PA) is a benign primary tumor that arises from the pituitary gland and is associated with ophthalmological, neurological and endocrinological abnormalities. However, causes that increase tumor progressing recurrence and invasiveness are still undetermined. Several studies have shown N-myc downstream regulated gene 2 (NDRG2) as a tumor suppressor gene, but the role of NDRG2 gene in pituitary adenoma pathogenesis has not been elucidated. The aim of our research has been to examine NDRG2 mRNA expression in PA and to determine the associations between the NDRG2 gene epigenetic changes and the development of recurrence or invasiveness of PA and patient clinical data.

METHODS

The MS-PCR was used for NDRG2 promoter methylation analysis and gene mRNA expression levels were evaluated by qRT-PCR in 68 non-functioning and 73 functioning adenomas. Invasiveness was evaluated using magnetic resonance imaging with Hardy's modified criteria. Statistical analysis was performed to find correlations between NDRG2 gene mRNA expression, promoter methylation and patient clinical characteristics and PA activity.

RESULTS

The NDRG2 mRNA expression was significantly lower in the case of acromegaly (GH and IGF-1 hypersecretion) than in other diagnoses of PAs (p < 0.05). Also, the NDRG2 expression was significantly higher in prolactinoma (PRL hypersecretion) than in in other diagnoses of PAs (p < 0.05). The promoter of NDRG2 was methylated in 22.69% (12/58 functioning and 15/61 non-functioning) of patients with PA. However, the NDRG2 gene mRNA expression was not significantly related to its methylation status. Clinical factors, such as: age, gender, relapse and diagnoses of Cushing syndrome were of no significance for NDRG2 promoter methylation and mRNA expression levels, as well as secreting or non-secreting PAs and the invasiveness of PAs.

CONCLUSION

The different NDRG2 promoter methylation and expression levels in PA samples showed tumor heterogeneity and indicates a potential role of this gene in pituitary adenoma pathogenesis, but the corresponding details require intensive research.

MMP-14 and TGFβ-1 methylation in pituitary adenomas

Pituitary adenoma (PA) is one of the most common abnormalities in the sellar region. Despite the fact that PA is a benign monoclonal neoplasm, it can cause serious complications, including ophthalmological, neurological and endocrinological abnormalities. Currently, the causes that increase the progression of tumors are unknown. Epigenetic silencing of the matrix metalloproteinase-14 (MMP-14) and transforming growth factor beta-1 (TGFβ-1) genes may be associated with the development of PA, since these genes are important in the processes of tumor metastasis and angiogenesis. The purpose of the present study was to determine if the methylation status of the MMP-14 and TGFβ-1 promoters is associated with PA development. In the present study, 120 tissue samples of PA were used. The methylation status of the MMP-14 and TGFβ-1 promoters was investigated by methylation specific-polymerase chain reaction. Statistical analysis was conducted to investigate the associations between the methylation status, age and gender of PA patients, PA tumoral activity, recurrence and invasiveness. The MMP-14 gene was methylated in 30.00% (17/56 functioning and 19/64 non-functioning) of patients with PA, while the TGFβ-1 gene was methylated in 13.33% (9/56 functioning and 7/64 non-functioning) of patients with PA. It was also observed that promoter methylation of MMP-14 correlated with the male gender (58.8 vs. 35.7%, P=0.022), while unmethylated (non-silenced) MMP-14 correlated with the female gender (64.3 vs. 41.7%, P=0.027). Associations between the promoter methylation status of the MMP-14 and TGFβ-1 genes and PA functioning or recurrence were not identified. The present study reveals that silencing of the MMP-14 gene correlates with patients' gender. However, MMP-14 and TGFβ-1 promoter methylation cannot be considered as a prognostic marker in PAs.

Testin (TES) as a candidate tumour suppressor and prognostic marker in human astrocytoma

Astrocytomas are one of the most common brain tumours; however, the current methods used to characterize these tumours are inadequate. The establishment of molecular markers may identify variables required to improve tumour characterization and subtyping, and may aid to specify targets for improved treatment with essential prognostic value for patient survival. One such candidate is testin (TES), which was reported to have prognostic value for glioblastoma patients. However, the role of TES protein in gliomagenesis is currently unknown. In the present study, the methylation status of the TES promoter was investigated in post-operative astrocytoma tumours of different malignancy grade, and its association with the survival of astrocytoma patients was evaluated. In addition, the expression of TES protein was investigated in the same set of astrocytoma tumours tissue, and the association of protein expression with glioma patients survival was evaluated. The methylation status of TES was assessed by methylation-specific polymerase chain reaction in 138 different grade astrocytoma samples. Western blot analysis was used to characterize the expression pattern of TES in 86 different grade astrocytoma specimens: 13 of pathological grade I, 31 of pathological grade II, 17 of pathological grade III and 25 of pathological grade IV (glioblastoma). Statistical analyses were conducted to investigate the association between tumour molecular pattern, patient clinical variables and overall survival. The methylation analysis of the TES promoter exhibited a distinct profile between astrocytomas of different malignancy grade (P<0.001). Furthermore, gene promoter methylation was significantly associated with patients' age, survival and pathological grade (P<0.001). The protein expression level of TES was significantly lower in glioblastoma (grade IV astrocytoma) than in lower grade (II–III) astrocytoma tissue (P=0.028 and P=0.04, respectively). Additionally, short overall survival of patients was markedly associated with low TES protein expression (P=0.007). However, no association between TES methylation and TES protein expression was noticed. The present study demonstrated that decreased expression of TES may be important in tumour progression and prognosis in human astrocytomas. TES may be a useful marker for predicting the clinical outcome of astrocytoma patients.

MGMT, GATA6, CD81, DR4, and CASP8 gene promoter methylation in glioblastoma

BACKGROUND

Methylation of promoter region is the major mechanism affecting gene expression in tumors. Recent methylome studies of brain tumors revealed a list of new epigenetically modified genes. Our aim was to study promoter methylation of newly identified epigenetically silenced genes together with already known epigenetic markers and evaluate its separate and concomitant role in glioblastoma genesis and patient outcome.

METHODS

The methylation status of MGMT, CD81, GATA6, DR4, and CASP8 in 76 patients with primary glioblastomas was investigated. Methylation-specific PCR reaction was performed using bisulfite treated DNA. Evaluating glioblastoma patient survival time after operation, patient data and gene methylation effect on survival was estimated using survival analysis.

RESULTS

The overwhelming majority (97.3%) of tumors were methylated in at least one of five genes tested. In glioblastoma specimens gene methylation was observed as follows: MGMT in 51.3%, GATA6 in 68.4%, CD81 in 46.1%, DR4 in 41.3% and CASP8 in 56.8% of tumors. Methylation of MGMT was associated with younger patient age (p < 0.05), while CASP8 with older (p < 0.01). MGMT methylation was significantly more frequent event in patient group who survived longer than 36 months after operation (p < 0.05), while methylation of CASP8 was more frequent in patients who survived shorter than 36 months (p < 0.05). Cox regression analysis showed patient age, treatment, MGMT, GATA6 and CASP8 as independent predictors for glioblastoma patient outcome (p < 0.05). MGMT and GATA6 were independent predictors for patient survival in younger patients' group, while there were no significant associations observed in older patients' group when adjusted for therapy.

CONCLUSIONS

High methylation frequency of tested genes shows heterogeneity of glioblastoma epigenome and the importance of MGMT, GATA6 and CASP8 genes methylation in glioblastoma patient outcome.