3. Discussion
This comprehensive bioinformatics assessment of the GABA system transcriptome in TCGA ACC patient tumors reveals that these tumors express transcripts encoding GABA shunt enzymes, multiple GABA receptor subunits, and GABA transporters. Of particular note, transcripts encoding GABA’s metabolizing enzyme, ABAT, are upregulated in more than 40% of ACC tumors, and its high expression is associated with favorable clinical attributes, including a three-fold increase in patient survival time. Tumors with upregulated ABAT transcript have enriched expression of genes in the gene ontology pathways of “mesenchymal-to-epithelial transition” and “interferon gamma signaling,” two pathways associated with limiting the spread of cancer. Our assessment of all GABA shunt transcripts, and those encoding other enzymes proximal to the shunt (such as those in the TCA cycle), implicates a GABA shunt metabolic phenotype that favors mitochondrial metabolism in patients with positive prognoses. Beyond these metabolic ramifications, ACC tumors express transcripts encoding multiple GABA transporters, suggesting that intracellular and intercellular tumoral GABA levels are highly regulated. Additionally, tumors express transcripts for both GABAA and GABAB receptor subunits. Regarding the former, upregulation of two transcripts encoding GABAA receptor subunits, GABRB2 and GABRD, are prognostic. This extensive GABA system transcriptome in ACC tumors is largely mirrored by the ACC cell line NCI-H295R, suggesting that this cell line is an appropriate in vitro model for studying the role of the GABA system in growth phenotypes for this cancer.
Approximately one-third of ACC patients have elevated serum cortisol levels [
6], with symptomatic hypercortisolism negatively associated with ACC patient survival even following tumor resection and/or mitotane treatment [
29]. For these patients, those few with upregulated tumoral
ABAT transcripts trend toward positive clinical outcomes vs those with unaltered
ABAT expression, while in patients with normal cortisol levels, upregulation of tumoral
ABAT levels is strongly correlated with positive clinical outcome. Enrichment analysis of ACC tumors with upregulated
ABAT identified two altered cortisol-related genes. Underexpression of
HSD3B2 (encoding 3β-hydroxysteroid dehydrogenase, 3β-HSD) was observed in tumors with upregulated
ABAT transcript expression, i.e., those with favorable patient prognoses. Consistent with this, high prostatic 3β-HSD protein levels are correlated with several unfavorable clinical attributes in prostate cancer patients [
30], while
HSD3B2 transcript expression has been proposed as a potential prognostic biomarker in lung squamous cell carcinoma [
31]. The involvement of steroidogenesis pathways, especially those leading to cortisol production and glucocorticoid receptor activation, may be relevant to cancer as the glucocorticoid receptor can regulate the expression of hundreds of genes, including those affecting cell proliferation, apoptosis, and gene transcription [
32,
33]. In this regard, overexpression of
NR3C1 transcripts (encoding the glucocorticoid receptor) was noted in tumors with upregulated
ABAT transcripts. Intriguingly, several glucocorticoid receptor binding sites are present within the
ABAT gene (our mining of the Gene Transcription Regulation Database [
34]), implying that the glucocorticoid receptor may regulate
ABAT transcription.
ABAT’s correlation with ACC patient hormone status and the expression of cortisol-related transcripts may represent an intersection of the GABA and steroid systems that could bear on ACC progression.
In primary tumors of other cancers, high
ABAT transcript expression is associated with favorable patient outcomes. Estrogen-receptor positive (ER
+) breast cancers express highly upregulated levels of
ABAT transcripts and its encoded protein, and have a more favorable patient prognosis, than low-
ABAT-expressing hormone-insensitive breast cancers [
35]. Consistent with these findings, low
ABAT transcript levels are associated with resistance to endocrine therapy in ER
+ breast cancer [
36], as well as a decrease in relapse-free survival [
35,
37,
38] and resistance to chemotherapy [
39] in both ER
+ and estrogen receptor negative (ER
−) tumors. Notably, triple-negative breast cancer has relatively low levels of
ABAT transcripts and protein, and a very poor patient prognosis [
35], with overexpression of
ABAT reducing migration and invasiveness in triple-negative cell lines [
39]. Additionally, low levels of
ABAT transcripts are associated with poor patient survival in mesenchymal glioblastoma [
40] as well as progression, recurrence, and reduced patient survival time in hepatocellular carcinoma (HCC) [
41]. In fact, ABAT protein is a highly sensitive biomarker for the identification of HCC and hepatoid adenocarcinomas [
42]. Recently, high
ABAT transcript levels were found to correlate with positive patient outcome in clear cell renal carcinoma, and to decrease cell proliferation and migration while increasing cell death in a kidney cancer cell line model [
43]. High
ABAT transcript expression is also favorably prognostic in renal (papillary), liver, and lung cancers (our mining of the Human Protein Atlas, v19.proteinatlas.org [
44]), and prevalently expressed in other cancers, such as thyroid carcinoma [
42], but is negatively prognostic for patients with esophageal squamous cell carcinoma [
45].
From an epigenetic standpoint,
ABAT gene-body methylation is negatively correlated with
ABAT transcript expression in ACC, implying methylation-driven transcriptional silencing of this gene. Accordingly, as increased
ABAT transcript expression is correlated with favorable ACC clinical outcome, increased levels of
ABAT gene body methylation are associated with poor patient prognosis. It is possible that
ABAT may be methylated by
DNMT1, as
DNMT1 upregulation is associated with increased levels of
ABAT methylation and is mutually exclusive with
ABAT transcript upregulation. Although
DNMT3B regulates
ABAT methylation in noncancerous cells [
46,
47], we found no correlation between
DNMT3B upregulation and
ABAT methylation level in ACC tumors. In other cancers, reversal of transcriptional silencing by demethylating agents in a glioblastoma cell line increases
ABAT transcript expression [
48]. High levels of
ABAT methylation are correlated with poor patient survival in myelodysplastic syndrome [
49]. While the complexities of
ABAT methylation remain to be understood, such regulation may impact cancer progression.
In general, the upregulation of
ABAT transcripts suggests an increase in the expression of its encoded protein, the aminotransferase ABAT. If so, tumors of ACC patients with the most favorable prognoses would likely display an increase in GABA catabolism, the consequence of which affects two distinct processes, as depicted in
Figure 11. Firstly, GABA activates GABA
A and GABA
B receptors. Therefore, GABA catabolism by ABAT lowers the amount of GABA available to activate GABA receptors. Secondly, the GABA shunt is a TCA cycle bypass mechanism beginning with the metabolism of glutamate to GABA by the activity of glutamic acid decarboxylase (GAD). Subsequently, the catabolic portion of the GABA shunt metabolizes GABA via ABAT activity into succinic semialdehyde (SSA), which is then metabolized by succinic semialdehyde dehydrogenase (SSADH) into succinate, a TCA cycle intermediate and substrate of the electron transport chain. The use of the GABA shunt to generate succinate has both anaplerotic and bioenergetic consequences.
The upregulation of tumoral
ABAT levels in a subset of primary ACC tumors reveals a prognostically favorable inter-patient metabolic heterogeneity among the ACC cohort that also includes the upregulation of
GAD1 and high expression of
ALDH5A1. This “GABA shunt metabolic phenotype” is accompanied by transcript variations among genes encoding enzymes proximal to the GABA shunt that can affect the availability of glutamate for GABA synthesis. Focusing on this synthesis, co-expression of
ABAT with
GLS, which encodes the glutamine-to-glutamate aminohydrolase glutaminase, suggests that increased glutamate may be available for GABA synthesis in tumors with high levels of
ABAT transcripts. In this context,
NAGS, encoding
N-acetylglutamate synthase (NAGS), is underexpressed in tumors with upregulated
ABAT. As NAGS catalyzes the generation of
N-acetylglutamate (NAG) from glutamate and acetyl-coA, decreased expression of
NAGS transcripts may result in more glutamate available as a substrate for GAD67 (encoded by
GAD1)—thus increasing GABA synthesis (
Figure 11). These data support the existence of a GABA shunt-centric inter-tumoral metabolic heterogeneity among the ACC cohort. Interestingly, an intra-tumoral heterogeneity of GABA levels has recently been reported in clonal cell lines of distinct metabolic phenotypes (e.g., glycolytic vs mitochondrial) derived from single tumor ex vivo mouse pancreatic cancer cells [
50]. Evaluation of both intra- and inter-tumoral heterogeneity of the GABA shunt may be important for understanding its role in mitochondrial metabolic heterogeneity, a feature of tumors that is increasingly recognized to influence cancer growth, progression, and drug treatment sensitivity [
51].
Focusing on the catabolic arm of the GABA shunt, high
ALDH5A1 transcript expression is present in nearly all primary ACC tumors. This gene encodes SSADH, which converts SSA to succinate, ensuring that the GABA shunt is driven to completion. Succinate is a substrate of succinate dehydrogenase (SDH), a TCA cycle enzyme and the electron transport chain’s complex II. The catalytically active subunit of SDH (SDHA) is encoded by the tumor-suppressor gene
SDHA [
52], and converts succinate into fumarate while generating a high-energy reducing equivalent. In ACC tumors, the co-occurrence of
SDHA and
ABAT upregulation supports the notion that GABA shunt activation contributes to succinate production. Additionally, the mutual exclusivity of
ABAT upregulation with that of
AKR7A2 (which encodes succinic semialdehyde reductase (SSR) and may lead to the diversion of GABA carbons out of the GABA shunt), and the under-expression of
AKR7A2 in ACC tumors with upregulated
ABAT transcripts may further ensure that the GABA shunt in these tumors culminates with the production of succinate. Unsurprisingly,
AKR7A2 upregulation is negatively prognostic for disease-free and progression-free survival, compared to patients with downregulated
AKR7A2 (unpublished). As the GABA shunt is commonly thought to act as a diversionary bypass of two steps of the TCA cycle, we were surprised to discover that expression of both
OGDH and
SUCLG2 transcripts, which encode the active subunits of the enzymes in these bypassed steps, were positively correlated with that of
ABAT, potentially indicating an increase in succinate production by both pathways. As both the TCA cycle and the GABA shunt may serve anaplerotic and bioenergetic purposes, these data imply an enhanced mitochondrial metabolism in a subset of primary ACC tumors in patients with favorable prognoses.
Metastases contribute to most cancer-related deaths [
53]. Our analysis of ACC patient clinical attributes revealed that higher levels of
ABAT transcript expression are correlated with lower-stage, organ-confined cases of ACC. In keeping with this finding, upregulated
ABAT expression is correlated with enriched expression of genes in the “mesenchymal to epithelial transition” (MET) pathway, which opposes the epithelial to mesenchymal transition (EMT) required for metastatic spread [
23]. Indeed, induction of the MET pathway in cancer cell lines reduces proliferation, migration, and/or invasion in vitro, as well as tumor growth and/or metastatic capability in mouse models [
54,
55,
56]. Of the MET pathway genes identified here, the majority are experimentally associated with maintaining an epithelial phenotype in cancer cells [
57,
58,
59,
60] and/or cell adhesion complexes that prevent migration and invasion [
61,
62,
63], or are known to be downregulated during EMT [
61,
64,
65]. Strikingly, when we evaluated the twenty enriched MET genes for their prognostic value in ACC patient outcome, expression of two genes—
PPL (a structural component of desmosomes) and
GPRC5A (a G protein-coupled receptor with an unidentified endogenous ligand)—were found to be associated with increased patient survival. In fact, no patients whose tumors demonstrated upregulated
PPL died during the 154-month observation period, whereas those with lower
PPL expression had a median survival time of 69 months.
The enrichment of MET pathway genes in tumors with upregulated
ABAT transcripts indicates a decreased likelihood of invasive cells exiting the primary tumor, suggesting that low
ABAT levels are associated with a predisposition for metastasis. Consistent with this notion, as shown here, primary ACC tumors in patients with metastases have significantly lower
ABAT transcript levels than primary tumors without metastases. This is in keeping with the correlation of low
ABAT transcript expression with both an increased risk of metastasis in triple-negative breast cancer tumors [
39] and an increased expression of a “metastasis gene signature” in hepatocellular carcinoma [
41]. While low
ABAT transcript levels may be important for EMT and the initiation of metastasis, the colonization of the metastatic niche requires activation of the MET pathway [
23]. Therefore, we predict that, in ACC metastatic spread, secondary-site colonization is associated with an increase in
ABAT transcript levels. Consistent with this prediction, both human breast-to-brain metastases [
66], and brain metastases of various cancers in xenograft mouse models [
67] show elevated
ABAT transcript expression relative to primary tumors. In fact, the ABAT inhibitor vigabatrin suppressed the proliferation and the number of brain metastases, and reduced the size of existing brain metastases, in xenograft mouse models [
67]. In humans, breast-to-brain metastases also express increased levels of
GAD1 transcripts [
66], while brain metastasis of ccRCC contain high levels of GAD67 protein [
68]. Mouse models of brain metastases similarly show elevated levels of
GAD1 transcripts [
67].
ABAT encodes a dual-function mitochondrial enzyme responsible for both GABA catabolism and maintenance of mitochondrial DNA copy number [
69]. Other genes encoding mitochondrial proteins that are prognostic for survival in ACC include
PINK1, encoding a mitochondrial serine/threonine kinase [
18,
70], and
FATE1, encoding a mitochondrial fission factor-family protein [
71,
72]. Although not tested in ACC clinical trials in the United States (
ClinicalTrials.gov), drugs targeting TCA cycle-associated enzymes are now in clinical trials for various other cancers [
73]. We determined that many transcripts encoding these enzymes are prognostic for ACC outcome (unpublished). In support of mitochondrial metabolism as an ACC therapeutic target, a drug blocking such metabolism in a xenograft mouse model of ACC inhibits the growth of tumors without overt toxicity [
74].
Our enrichment analysis also reveals that ACC tumors with upregulated
ABAT transcripts are enriched for transcripts of the “interferon gamma (IFN-γ) signaling” biological pathway [
75]. These enriched IFN-γ pathway genes are associated with decreased tumor progression and improved patient prognosis [
76,
77,
78,
79,
80,
81]. Of particular note are transcripts encoding HLA complex proteins (
Table 5), which are commonly associated with tumor regression [
82], as the expression of HLA proteins on the surface of cancer cells allows T-cells to recognize and destroy cancer cells [
83]. Interestingly, T-cells synthesize [
84] and secrete [
85] GABA, and express multiple components of the GABA system—including
GAD1 and
ABAT transcripts, as well as functional GABA transporters, which may be upregulated upon T-cell activation [
84]. Additionally, GABA
A receptor activation suppresses T-cell proliferation [
84,
86,
87].
ACC tumors express transcripts encoding three distinct plasma membrane GABA transporters—
SLC6A1,
SLC6A12, and
SLC6A13. Importantly, the expression of
SLC6A1 transcripts in most ACC tumors is both high and positively correlated with
ABAT transcript expression. Plasma membrane transporters in tumors may permit the exchange of GABA among tumoral cells, controlling GABA levels through bidirectional GABA transport as has been noted in both neurons [
88] and cultured T-cells [
85,
89]. In this context, extracellular GABA can be taken up either for use as a metabolic substrate and/or to terminate receptor activation, as in the central nervous system. Conversely, intracellular GABA may be released into the extracellular tumor microenvironment (TME) to activate GABA receptors, or for cooperative metabolism between cells of the tumor. Such exchange of metabolites between tumor cells has been shown for other metabolic substrates [
90,
91].
Despite GABA being primarily recognized as a neurotransmitter, its broader importance is evidenced by the use of the GABA shunt by bacteria, yeast, and plants—which utilize the GABA shunt to mitigate stressors including oxidative, pH, and nutrient deprivation stress [
92,
93,
94,
95,
96,
97,
98,
99,
100,
101,
102]. In particular, flux of carbons through the GABA shunt may affect ROS production, as supported by the demonstration that altered
ABAT expression and ABAT activity affects mitochondrial oxygen consumption in non-cancerous cells [
47,
103], and that the deletion of
ABAT orthologs in fungi causes increased sensitivity to oxidative stress [
96]. Modulation of ROS production by GABA shunt activity, to our knowledge, is unexplored with respect to cancer.
The presence of both GABA
A and GABA
B receptor-encoding transcripts in ACC tumors suggests that both types of GABA receptors may be operable in ACC. With respect to the GABA
A receptor,
GABRB2 and
GABRD, the genes encoding the receptor’s β2 and δ subunits, are both highly expressed in ACC tumors, and are positively and negatively prognostic for ACC patient survival, respectively. The most commonly contemporaneously-expressed transcripts within individual tumors are
GABRA3,
GABRA5,
GABRB2,
GABRB3,
GABRG1,
GABRD, and
GABRE—encoding the α3, α5, β2, β3, γ1, δ, and ε subunits, respectively—with the vast majority of all ACC tumors expressing transcripts encoding both δ and ε subunits. Given that the prototypic GABA
A receptor composition requires two α subunits, two β subunits, and a fifth position subunit (γ, δ, ε, θ, π), we propose that the most common putative pentameric receptor schemas found in ACC tumors are permutations of α
3α
5β
2β
3γ
1/δ/ε receptors (
Figure 7). Receptors composed of these subunits have been identified as being tentatively extant in brain [
104], while similar receptor subunit compositions have been proposed in lung [
105].
Correlation of GABA
A receptor subunit gene expression with cancer growth phenotypes in cancer cell lines, animal models, and patients strongly implicates GABA
A receptors in cancer progression. Expression of
GABRG3,
GABRQ,
GABRE, and
GABRR2 is correlated with inhibition of growth phenotypes in cell lines and with favorable patient outcomes [
106,
107,
108] (Broad Institute of MIT & Harvard, firebrowse.org), while expression of other transcripts (
GABRB2,
GABRP) is associated with cancer progression [
109,
110,
111]. Of particular relevance to our work, high expression of
GABRA3 transcripts, or its encoded α3 protein, is correlated with poor patient outcome in neuroblastoma [
112] and lung adenocarcinoma [
113,
114,
115]. In breast cancer cell culture and xenograft models,
GABRA3 transcripts undergo adenosine to inosine (A-to-I) RNA editing that decreases cell invasiveness vs unedited
GABRA3, by inducing intracellular retention of α3 subunit-containing receptors [
116]. Although
GABRA3 transcripts are not prognostic for ACC patients, their A-to-I editing status in ACC tumors has not been determined.
Pharmacological experiments with GABA
A receptor ligands demonstrate that receptor activation induces cell death in ex vivo medulloblastoma cells [
117] and a transgenic mouse model of neuroblastoma [
112], and inhibits the invasiveness and tumor formation of colon carcinoma cells in vitro and in xenograft mice [
118,
119]; while blocking GABA
A receptors promotes tumor growth of glioblastoma [
120]. Not all studies, however, have found GABA
A receptor activation to inhibit cancer cell growth, as receptor activation increases cell proliferation in prostate adenocarcinoma and oral squamous cell carcinoma cell lines [
121,
122,
123]. Such discrepant findings between studies of different cancers may be due to the varying composition of heteropentameric GABA
A receptors, as the 19 distinct subunits impart different functional characteristics (reviewed in [
104]).
In addition to mediating effects via GABA
A receptors, GABA also activates the GABA
B receptor, a heterodimeric G protein-coupled receptor. While two-thirds of ACC patient tumors express transcripts encoding both GABA
B receptor subunits, neither are correlated with patient prognosis. However, it is likely that GABA
B receptors are functional in ACC tumors, as shown in the ACC NCI-H295R cell line—albeit with undetermined effects on cell growth phenotype [
20]. In other cancers, activation of GABA
B receptors enhances the invasiveness of prostate adenocarcinoma cell lines [
124,
125,
126], while inhibiting migration and proliferation in cell lines of hepatocellular carcinoma [
127,
128], pancreatic cancer [
129], and lung adenocarcinoma [
115,
130], as well as decreasing tumor volume in hepatocellular carcinoma mouse xenografts [
127]. Given the ability of GABA
B receptors to affect growth phenotypes in other cancers, the presence of their transcripts in ACC tumors, and their functionality in the NCI-H295R ACC cell line, the impact of GABA
B receptors on cancer growth phenotypes warrants investigation.
The NCI-H295R cell line is a widely-used [
28,
131], and by far the most widely-cited (>800 citations in PubMed) ACC cell line model, that we show has a strikingly similar GABA system transcriptome to that of ACC tumors. Such similarity indicates that GABA system components likely reside, at a minimum, within the cancerous cells of the tumor. Like ACC tumors, the NCI-H295R cell line expresses transcripts encoding GABA shunt enzymes as well as a variety of GABA receptor subunits and multiple GABA transporters. In fact, these cells express
GAD1 transcripts that encode functional GAD67 protein, as determined by the production of
14C-CO
2 from
14C-glutamate [
20]. Of potential relevance to cancer, treatment of these cells with endocrine disrupting agents induces a four-fold upregulation of
GAD1 transcript levels [
132] (see the supplementary microarray data in Ref. 132) and a significant decrease in
ABAT transcript levels (unpublished, determined using Song, et al.’s methodology). Due to the similar GABA system transcriptomes between ACC tumors and the NCI-H295R cell line, and the ability of endocrine disruptors to affect this transcriptome, this cell line is a particularly appropriate in vitro ACC model for assessing the role of the GABA system in cancer cell growth phenotypes such as proliferation, migration, and invasion.
One notable difference between ACC tumors and the NCI-H295R cell line is that ACC tumors display high levels of
GABRD transcript expression, while we observed only relatively low levels expressed in the cell line. One explanation for this disparity is that
GABRD transcripts may be predominantly expressed in ACC tumor host cells, such as T-cells. This notion is consistent with an unsupervised bioinformatics study of TCGA tumors (not including ACC) that identified
GABRD expression as a pan-cancer marker that is overexpressed in tumor vs adjacent normal tissue
[27]. Our data mining effort extends this Gross et al., finding to include ACC tumors. Additionally, our results indicate that within the ACC cohort,
GABRD expression is upregulated in a subset of patient tumors with lower
ABAT transcript levels, and that
GABRD upregulation is not only unfavorably prognostic, but also may be associated with a decreased capacity to degrade GABA.
Our work identifies the expression levels and prognostic implications for transcripts encoding components of the GABA system—most notably, the GABA shunt and certain GABA
A receptor subunits, which largely appear together in the same subset of ACC tumors. Additionally, transcripts encoding GABA
B receptors and GABA transporters are ubiquitously present in ACC tumors. The use of human tumor data lends particular relevance to our findings. Beyond this strength, we present a full assessment of these genes as a system with respect to cancer, giving a context to the many unsupervised studies in which these genes have been individually identified. A limitation of this study is that bioinformatics data are correlational, and thus functional studies are required to determine whether any differences in transcript expression are biologically or phenotypically meaningful [
133].
Future functional studies will need to be conducted in model systems to determine whether GABA shunt function and/or GABA receptor activation affects cell growth phenotypes or cell metabolism. To this end, we have determined that the NCI-H295R cell line GABA system transcriptome mirrors that of ACC tumors. Future studies of this cell line should confirm the presence of GABA system proteins and evaluate whether pharmacological and/or genetic manipulation of this system affects cell proliferation, migration, or invasion. Additionally, the metabolic phenotype associated with GABA’s role as a potential anaplerotic substrate should be investigated with respect to mitochondrial metabolism and bioenergetics. If functional studies in cell lines determine that manipulation of the GABA system impacts cell growth and metabolism, in vivo ACC models can be used to characterize whether the arsenal of GABA system therapeutics—currently used to treat neurological and psychiatric disorders [
134,
135]—affects tumor formation, growth and metastasis. Such therapeutics include those that inhibit ABAT (Sabril
®), allosterically potentiate GABA
A receptors (Klonipin
®, Solfoton
®, and Brexanolone
®, among many others), activate GABA
B receptors (Lioresal
®), and inhibit GABA transporters (Gabitril
®). If functional studies prove fruitful, the clinical availability of these pharmacological agents may offer an accelerated path towards drug discovery for ACC treatment.