Material and Methods
This is a prospective cross-sectional study done at the National Center of Neurological Sciences (NCNS) Khartoum, from 2018 to 2021. The study included all cerebral tumors that were radio- logically diagnosed preoperatively as meningioma. Ethical approval was taken from Institute of Endemic disease, Khartoum University, and written consent was obtained from each patient.
The tumor specimens were obtained immediately from the operating room, before surgery was completed, and sent to the research laboratory in the center within 2-5 minutes. At arrival to the laboratory, small portion from 20 tissue samples were processed for DNA extraction, three of which were selected for whole exome analysis and each was processed further for detection of gene of centrality by using PCR and thus sequencing analysis. Three non meningioma brain tumor were used as control
Whole Exome Sequencing (WES)
Sample Preparation
After evaluating the detailed histopathology of our tumor specimens, 3 samples were selected for DNA extraction for whole exome analysis and each was proceeded further for detection genes of centrality.DNA quality and quantity were assessed using Nano-Drop ND-1000 (Nano-Drop Products, Wilmington, DE, USA) spectrophotometer at wavelength spectrum of 220-750nm and standard 1% agarose gel electrophoresis. Each sample was labeled with a unique accession number, and DNA was sent to Beijing Genomics Institute (BGI) for whole exome sequencing.
Exploratory Analysis
Genes set were prepared by filtering for all genes harboring exonic and splice site variations stratified into SNVs and INDELs set genes for each sample. In addition to that, we used 1000G online database to explore novel mutations. And also we explored the excess variant of allele frequency per 3 samples.
To predict possible impact of an amino acid substitution on the structure and function of a human protein, SIFT and PolyPhen-2 programs were used.
In addition to this, shared genes filtering for all genes which stratified into two classes (SNVs and INDELs) per 3 sample, were constructed by utilizing R software (version 4.1.0, 2021-05-18)
Network Analysis Using Network Analyst 2021
In this study, Network analyst 2021 was used to put up individual functional interaction network of genes per each sample, and for gene set which harboring damaging effect. Although there was number of mutated genes, they were not all included in network analysis, we managed shared gene analysis in both class (INDEL and SNVs) per 3 sample, when needed , high and moderate, low and modifier impact were defined in SNVs. On topography of network analysis, genes of centrality can be conveyed.
To address meticulous genes interaction, and whether these interactions with particular centrality, all shared genes within the class INDEL and SNVs were used as input on Network Analyst (GenMania, Reactome and STRING) to predict protein- protein interaction and functional gene onto
Pathway and GO Enrichment Analysis
Based on enrichment analysis online tools (Network Analyst 2021), which was used to examine pathways enrichment for sets of genes in various data base, and application to forecast diseases associated.
Additionally, shared mutations within genes of centrality were detected using polymerase chain reaction by means of specific forward and reverse primers.
The primers were designed using Primer 3 online program. The table below is containing the list of primers used in this study.
AKT1 |
F |
CCCCTCAGATGATCTCTCCA |
R |
TGAAGAATTTGGAGGGAAGG |
CD44 |
F |
AGGAACAGTGGTTTGGCAAC |
R |
AACTGGCTTGTATCCATTCCT |
Discussion:
Meningioma is the most common brain tumor in Sudan and mostly reported in female, [
3] similarly, in this current study, there was clear cut towards female predominance. Most of these tumors are slowly growing and associated with symptoms of gradual increase in intracranial pressure, in this context, headache was the main cardinal feature among our patients. Understanding of the normal development of the nervous system has noticeably improved in recent years, it has a complex differentiation hierarchy ranging from a neural stem cell that can give rise to all of the major lineages in the brain parenchyma to lineage committed progenitors that have a more restricted differentiation potential to terminally differentiated cells. Stem cell concepts can influence the understanding of brain cancer. The cancer stem cell hypothesis proposes that established tumors consist of a cellular hierarchy with a subpopulation of tumor cells able to maintain and propagate the tumor. In this current study, thousands of mutations were predicted by whole exome sequencing analysis, of these high C>T of transversions mutations followed by G>A were observed per fibrous, meningiothelial and fibrous recurrence, our results goes with international results, that different cancer type have different mutational signature similar to gastric cancer in which the prevalence of C>T and C>A transversions were high.[
19]What is more, among our samples, numerous mutations that harboring deleterious poly phen 2 effects were found, which highlight the genetic complexity on this type of brain tumor. Are these thousands mutation are being claimed with development of tumorigenesis process on this type of brain tumor, a question is not merely answered. Despite the similarity, further investigations are considered necessary in our study to investigate prospective functions.
In recent years, the availability of large-scale network and genomic data sets allow studying the association between the position of proteins within molecular networks and their patterns of molecular evolution,[
20] the consequence of these studies have shown the individual genes are affected by the position that their encoded products occupy in molecular networks. Genes acting at the center of protein–protein interaction networks and metabolic networks evolve under higher levels of purifying selection than those acting at fundamentals of a network periphery, and those genes are so named genes/ proteins of Centrality. Hence Protein interactions play important roles in vital biological processes such as cell cycle control, metabolic and signaling pathways and disease pathways. These interactions can be represented as complex networks, where the nodes are the proteins and the edges represent the interactions between the pairs of proteins they connect.[
21]
With regard to the potential of our datasets to guide precision analysis, genes/ protein of centrality are
Akt1,
UBC and numerous stemness genes like
CD44, which emphasize the significance of network.The significance of
Akt1 gene encodes one of the three members of the human
Akt serine-threonine protein kinase family which are often referred to as protein kinase B alpha, beta, and gamma. These proteins are phosphorylated by phosphoinositide 3-kinase (PI3K), AKT/PI3K forms a key in component of many signaling pathways that involve the binding of membrane-bound ligands resembling receptor tyrosine kinases, G-protein coupled receptors, and integrin-linked kinase. Therefore, these
Akt proteins regulate wide variety of cellular functions of which cell proliferation, survival, metabolism, and angiogenesis together in normal and malignant cells, this is along with,
Akt proteins are recruited to the cell membrane by phosphatidylinositol 3, 4, 5-trisphosphate (PIP3, 4, 5).Mutations in this gene are associated with multiple types of cancer like breast, colorectal, and ovarian cancers.[
22]
UBC genewhich is located on 12q24.3, consisting of 2 exons,the promoter of it contains heat shock element which renders and mediates
UBC stimulation upon stress. It also plays key role in maintaining cellular ubiquitin level under stress condition.[
23]
Ubiquitin (
Ub), was first identified by Gideon Goldstein et al. in 1975 and further confirmed over the next several decades in the human genome. It is highly conserved regulatory protein containing 76 amino acids, can be covalently tagged to target proteins via cascade of enzymatic reactions and regulates numerous physiological and pathological functions.[
24,
25]
Moreover, this gene is associated with various signaling pathways, of which important nutrient and environmental stimulus. The amino acid plays critical role in the mechanistic target of rapamycin complex 1 (mTORC1) signaling pathway. [
26,
27] In cancerous stem cell, it regulates and maintains stemness of genes expression likewise
Oct4,
Sox2 and
Nanog.[
28]It also plays fundamental role in CSC characteristics, such as self-renewal, maintenance, differentiation and tumor-genesis.[
29]
Notably, this gene is up regulated in several cancer, meanwhile the role of UBC in meningioma should be addressed, in this current study, UBC per 2 classes is mutated and represented gene of centrality within each class individually, however, it is not the case, when all shared being the case.
Tumor is the disease which is characterized by multiple phenotypes of which is abnormal cell growth and miscellaneous biological progression. The interaction of biological entities of genes, proteins and metabolites is of great interest in cancer development, as shown in our results, signaling pathway that yelled protein of Akt family. Recently, our results on using next generation sequencing, identified multiple genetic signaling pathways of which is Pathways in cancer, Proteoglycans in cancer and PI3K-Akt signaling pathway.
It is known that signaling pathways control cell cycle, however, the mechanism and interactions among these pathways vary between individual tumor and tumor types. For instance, our results among meningioma histological subtypes, taking into account pathways, more than 10 pathways are significantly co expressed and have been compared in a number of pathway data base. Interestingly, these pathways are of considerable similarity between our 3 meningioma sample, and this might denote genetic resemblance, since most of input genes did not show variation. Furthermore, our exome analysis highlight number of pathways that related to cancer development; of these, proteoglycan pathway which provided connection between cell membrane and surrounding extra cellular membrane, thus playing vital role in regulating cell cancer, adhesion and migration. Another detected pathway in our sample was miRNA signaling pathway, miRNAs are family of small non coding RNAs which regulate wide range of biological processes that encoding carcinogenesis, and thus have been known dysregulated in cancer. Generously, abnormal expression of miRNA in tumor, it believed that this tumor dysregulation might affect several cancers hallmarks on other words, initiation, and progression. Besides, functionally these miRNAs incorporate into cell cycle proliferation pathway, and the dysregulated miRNAs are accountable for evading growth suppresser and hence sustaining proliferating signaling pathway in cancer.
Nevertheless, Akt pathway is intracellular pathway promoted metabolism, proliferation, cell survival growth and angiogenesis in response to extracellular signals. Yet, this pathwayregulates different cell functions and play conservative role through mTORcomplex.
Remarkably, stemness is regarded as the key factor of carcinogenesis and resistance to chemotherapy. [
30]Although the cancer stem cells were yet to be isolated from meningioma, cancer stem cell and embryonic stem cell markers have been widely identified from meningioma and strongly associated withpatients’ outcomes. [
31]
The majority of exome analysis focused on analyzing cancer etiology depending on the mutations, however with no or little look over nature of these genes, few of these studies questioned stem cell genes.
Outstandingly, from our results, stem cell genes were mutated equally per our 3 meningioma samples which might suggest that, cancerous stem cell were probably drive the force of meningioma development, and taken together the distinct histological variant of meningioma ( fibrous, meningiothelial and fibrous recurrence) which are the most common variant, suggesting the biological behavior other than recurrence. In absolute numbers of meningioma, most recurrence meningiomaare corresponded to histological benign variant (WHO grade1), however, their potential aggressive behavior still need to be evaluated.
Various cancers are associated with mutations on
Akt1 gene, and some mutations have been associated with the meningioma recurrence. Interestingly, in this current study,
Akt1 has shown number of mutations among 20 meningioma variants. While
Akt1 plays role in modulating synaptic plasticity and signal transduction pathways,[
32] deficiency in
Akt1 may lead to abnormal prefrontal cortical structure and deficits in cognitive functions and risk for schizophrenia.[
33]This might be claimed to behavioral changes among our patients with this kind of brain tumor.
Overall, our finding implicated rs17846829 which was seen in fibrous and atypical meningioma variant and rs58565216 that was seen within the fibrous variant might explain cognitive decline in meningioma patients. A recent robust study using array comparative genomic hybridization (aCGH), exposed that, identification of
SMO and
Akt1 in meningioma increased the possibility of targeting therapy among meningioma patients.[
34]Another German a study, concluded that, mutation of
Akt1 in skull base meningioma designate shorter time of recurrence and activate mTOR and ERK1/2 signaling and targeting mTOR pathway might be of benefit to patients with recurrent skull base meningioma.[
35] Furthermore, case study from Switzerland exhibited presence of activated mutations of
Akt1in sub-population of meningioma by using sequencing, and culturing of meningioma discovered sensitivity to AZD5363, which is a selective inhibitor to kinase activity. Thus clinical therapy of AZD5363 was given to the patients for tumor control. [
36] Our study promotes the view, that mutations in
Akt1 denote recurrence and might be future challenges in clinical therapy; however, big data might be needed to consolidate our findings. Moreover,
CD44 A>G (rs9666607 mutation was detected in all our tumor samples, this mutation is located within stretch of residues annotated in UniProt as a distinct region, and the differences in amino acid properties can disturb this region and eventually disturb its function. Furthermore, this mutation is located in a region with known splice variants, and is described K > R (in dbSNP: rs9666607). Despite the mutant residue among the observed residue types at this position in other homologous sequences, and sometimes suggests that, this variant was not damaging for the protein’s structure and function, however, this mutation was already established as having harmful effect.[
37]In one study which analyzed SNPs within stemness-related genes for association with overall survival in the prostate, lung, colorectal and ovarian cancer screening and thus identified SNPs of
CD44 (rs9666607)were evaluated for association with gene expression.[
38]In light of the above this rs9666607 in
CD44among our meningioma samples was detected in
CD44 involves with
Akt1 gene in a proteoglycan pathway in cancer.[
39]Nevertheless,
CD44 frame shift mutation (D419Tfs*61) in our samples, which is located within stretch of residues annotated in UniProt as special region consequently differences in amino acid properties can disturb this region and disturb its function, and as well this mutation showed difference in charge between the wild-type and mutant amino acid. Additionally, the charge of the wild-type residue was lost, this be able to cause loss of interactions with other molecules or residues.[
40]