Abstract
Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is a hereditary cancer syndrome in which affected individuals are at risk for development of cutaneous and uterine leiomyomas and an aggressive form of type II papillary kidney cancer. HLRCC is characterized by germline mutation of the tricarboxylic acid (TCA) cycle enzyme, fumarate hydratase (FH). FH-deficient kidney cancer is characterized by impaired oxidative phosphorylation and a metabolic shift to aerobic glycolysis, a form of metabolic reprogramming referred to as the Warburg effect. Increased glycolysis generates ATP needed for increased cell proliferation. In FH-deficient kidney cancer, levels of AMP-activated protein kinase (AMPK), a cellular energy sensor, are decreased resulting in diminished p53 levels, decreased expression of the iron importer, DMT1, leading to low cellular iron levels, and to enhanced fatty acid synthesis by diminishing phosphorylation of acetyl CoA carboxylase, a rate-limiting step for fatty acid synthesis. Increased fumarate and decreased iron levels in FH-deficient kidney cancer cells inactivate prolyl hydroxylases, leading to stabilization of hypoxia-inducible factor (HIF)-1α and increased expression of genes such as VEGF and glucose transporter 1 (GLUT1) to provide fuel needed for rapid growth demands. Several therapeutic approaches for targeting the metabolic basis of FH-deficient kidney cancer are under development or are being evaluated in clinical trials, including the use of agents such as metformin, which would reverse the inactivation of AMPK, approaches to inhibit glucose transport, lactate dehydrogenase A (LDHA), the antioxidant response pathway, the heme oxygenase pathway, and approaches to target the tumor vasculature and glucose transport with agents such as bevacizumab and erlotinib. These same types of metabolic shifts, to aerobic glycolysis with decreased oxidative phosphorylation, have been found in a wide variety of other cancer types. Targeting the metabolic basis of a rare cancer such as FH-deficient kidney cancer will hopefully provide insights into the development of effective forms of therapies for other, more common forms of cancer. Clin Cancer Res; 19(13); 3345–52. ©2013 AACR.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
CME Staff Planners' Disclosures
The members of the planning committee have no real or apparent conflict of interest to disclose.
Learning Objective(s)
Upon completion of this activity, the participant should have a better understanding of the clinical manifestations of hereditary leiomyomatosis and renal cell carcinoma, the biochemical pathway of fumarate hydratase–deficient kidney cancer, and the potential therapeutic approaches for patients with advanced cases of this disease.
Acknowledgment of Financial or Other Support
This activity does not receive commercial support.
Background
Hereditary leiomyomatosis and renal cell carcinoma
Hereditary leiomyomatosis and renal cell carcinoma (HLRCC) is an autosomal dominant hereditary cancer syndrome in which affected individuals are at risk for the development of cutaneous and uterine leiomyomas and kidney cancer (1, 2). Leiomyoma, which is a hallmark of HLRCC, is a benign smooth muscle neoplasm. We have detected HLRCC-associated cutaneous leiomyomas in more than 75% of our affected patients. The number of cutaneous leiomyomas per patient ranged from zero to more than 100 lesions (3). Cutaneous leiomyomas most often occur in the arrectores pilorum muscles, which are attached to the hair follicles (4), The arrector pilori can be considered an energy sensing organ; when it is cold, the arrectores pilorum contract to trap air and provide insulation. These raised lesions can be painful and are often sensitive to touch or cold. No effective systemic or topical treatment currently exists for HLRCC-associated cutaneous leiomyomas. The treatment is most often symptomatic, sometimes including surgical resection. Females affected with HLRCC are at risk for development of multiple early-onset uterine fibroids. Uterine fibroids are detected in more than 90% of affected HLRCC females, 70% to 80% of whom will undergo either a myomectomy (surgical removal of the uterine leiomyoma) or hysterectomy (removal of the uterus; refs. 3, 5).
HLRCC-associated kidney cancer
Patients affected with HLRCC are at risk for the development of early-onset, bilateral, and/or multifocal renal cysts and papillary kidney tumors. Papillary kidney cancer has been reported to occur in 60% of HLRCC families (3, 5). Papillary kidney cancer can be classified as type I or type II. Type I papillary kidney cancer is characterized by an indolent growth pattern and is the histologic subtype found in patients with hereditary papillary renal carcinoma, which is associated with germline MET gene mutation (6, 7). Type II papillary kidney cancer, however, is a lethal disease associated with rapidly aggressive growth and early metastasis. HLRCC tumors are classified as type II papillary kidney cancer and are characterized histologically by orangophilic nucleoli and prominent perinucleolar halo formation (8). Because HLRCC renal tumors may occur in young individuals and have a propensity to spread when the tumors are very small, affected individuals have a life-long need to undergo annual abdominal imaging starting at the age of 10 years (9). Patients affected with other familial kidney cancer syndromes, such as von Hippel–Lindau (VHL; germline VHL gene mutation), hereditary papillary renal cell carcinoma (germline MET gene mutation), or Birt–Hogg–Dubé syndrome (germline FLCN gene mutation) are managed with active surveillance until the largest renal tumor reaches the 3-cm threshold, at which time renal parenchyma-sparing surgical resection is recommended (10, 11). However, due to the aggressive nature of HLRCC renal tumors, active surveillance is not recommended for the management of even small HLRCC-associated renal tumors; wide surgical excision is recommended when an HLRCC-associated renal tumor is detected (9).
Fumarate hydratase: the HLRCC gene
HLRCC is characterized by germline mutation of the tricarboxylic acid (TCA; Krebs) cycle enzyme gene, fumarate hydratase (FH; ref. 12; Fig. 1). FH, which exists in both a mitochondrial and a cytosolic form, is a homotetramer that catalyzes the hydration of fumarate to malate. Mutations of the FH gene, which may be missense, frameshift, or complete or partial deletion, are detected in 90% of HLRCC families (3, 5). FH, localized to chromosome 1q42.2, is subject to a 2-hit, loss of function in HLRCC tumors; loss of the somatic allele usually represents the second hit that is found in HLRCC-associated kidney tumors and uterine leiomyomas (13–15).
Shift to aerobic glycolysis in FH-deficient kidney cancer: the Warburg effect
FH-deficient kidney cancer undergoes a metabolic shift to aerobic glycolysis. The FH-deficient kidney cancer cell line UOK262, which is tumorigenic in a nude mouse model, is characterized by impaired oxidative phosphorylation and increased levels of glycolysis, as assessed by decreased oxygen consumption rate and increased extracellular acidification rate. Both the glycolytic and tumorigenic features of the FH-deficient cell line can be reversed by restoration of FH activity. In contrast with the cell line model derived from other forms of kidney cancer, such as VHL-deficient clear cell kidney cancer, FH-deficient kidney cancer cells are uniformly dependent on glucose and glycolysis for ATP production needed for rapid proliferation (16–18). This is shown in patients by the observation that, in contrast with other genetically defined types of kidney cancer, metastatic HLRCC-associated kidney cancer exhibits consistently high fluorodeoxyglucose uptake on positron emission tomography (PET) scan imaging (ref. 14; Fig. 1).
Inhibition of HIF prolyl hydroxylase results in stabilization of HIF-1α
Hypoxia-inducible factor (HIF)-1α and HIF-2α are transcription factors whose stability is regulated by HIF prolyl hydroxylase (PHD). In normoxic conditions, PHD (in the presence of α-ketoglutarate and iron) hydroxylate 2 HIF prolines, allowing the VHL complex to bind and facilitate ubiquitin-mediated degradation. In hypoxia (or with an inactivating mutation in VHL), HIF degradation is impaired and HIF-1α and HIF-2α accumulate (19). HIF-1α and HIF-2α are transcription factors that increase the transcription of a number of genes critical to a cancer dependent on aerobic glycolysis, such as the glucose transporter (GLUT1) and VEGF, which would increase glucose transport and tumor vascularity. FH-deficient kidney cancer has an increased need for both vascularity and glucose transport to provide increased nutrients for rapid growth and proliferation and for ATP production to compensate for a decrease in oxidative phosphorylation associated with a TCA cycle gene mutation. In FH-deficient cells, increased levels of fumarate inhibit the HIF PHDs and thus the activity of the VHL ubiquitination complex, leading to stabilization of HIF-1α and increased transcription of its targets such as VEGF and GLUT1 (ref. 20; Fig. 1).
Diminished AMPK, a negative regulator of the Warburg effect in FH-deficient kidney cancer, facilitates increased fatty acid and protein biosynthesis needed for rapid tumor growth
The AMP-activated protein kinase (AMPK) is a cellular energy sensor that reflects cellular energy status by undergoing phosphorylation and increasing activity when AMP levels increase and ATP levels decrease, indicative of energy deficiency (21). Although HLRCC tumor cells might be expected to show AMPK activation because of their dependence on less efficient glycolysis and their rapid consumption of ATP during growth, unexpectedly, the reverse was observed, with markedly reduced phopho-AMPK levels, but also reduced levels of AMPK subunits, AMPKα and AMPKβ1, at both transcriptional and protein levels. The reason for diminished AMPK levels is not known, but it was shown to result in reduced phosphorylation of acetyl CoA carboxylase (pACC), a state that normally renders the enzyme inactive, that would be expected to increase the synthesis of a key fatty acid biosynthetic intermediate, malonyl-CoA. In addition, the phospho-S6 ribosomal protein, an mTOR downstream effector, was also found to be activated. These results suggest that enhanced fatty acid and protein biosynthesis in FH-deficient kidney cancer results from decreased AMPK activity that promotes anabolic growth associated with cell growth and proliferation in tumors (17).
Decreased AMPK leads to decreased iron and increased HIF-1α levels.
The glycolytic shift in FH-deficient kidney cancer is associated with low AMPK levels, which decrease p53 (14) and the iron transporter, DMT1, which seems to require p53 for expression (14). In the FH-deficient kidney cancer cell line, UOK262, the iron responsive proteins, IRP1 and IRP2, as well as the IRP target, transferrin receptor protein 1 (TFRC), are elevated, indicating that cytosolic iron concentrations decrease secondary to decreased DMT1 activity (14). Changes in cytosolic iron levels and corresponding IRP activities potentially have opposite effects on HIF-1α and HIF-2α. Because PHDs require iron for their catalytic activity, decreased cellular iron would be expected to stabilize both HIF-1α and HIF-2α (22). However, although both FH-deficient kidney cancer cell line models, UOK262 and UOK268, have increased HIF-1α levels, both also have decreased HIF-2α levels. Decreased HIF-2α levels are consistent with increased IRP1 levels, as IRP1 is known to target the iron responsive element in the HIF-2α 5′-untranslated region and to inhibit HIF-2α translation when cytosolic iron is low, in vitro (23) and in vivo (24).
Activation of AMPK reduces the invasive potential of FH-deficient kidney cancer.
As noted above, AMPK is decreased in UOK262 cells, resulting in activation of phospho-S6 ribosomal protein and a decrease in pACC (which would lead to an increase in fatty acid synthesis). To evaluate the potential role of activation of AMPK in FH-deficient kidney cancer, UOK262 cells were treated with either metformin or 5-aminoimidazole-4-carboxamide-1-β-D-ribofuranoside (AICAR). Both metformin and AICAR were found to significantly impair the invasiveness of FH-deficient kidney cancer cells (14).
Silencing of HIF-1α reduces the tumorigenic potential of FH-deficient kidney cancer.
Some studies have suggested that increased HIF-2α levels are essential for tumorigenesis in VHL-deficient clear cell kidney cancer cells (25, 26), and a recent genetic and functional study has concluded that HIF-1α actually functions as a tumor suppressor gene in clear cell kidney cancer (27). However, HLRCC cancer does not conform to this model. In UOK262 FH-deficient kidney cancer cells, HIF-1α is increased, whereas HIF-2α is decreased. It could be that the role of HIF-1α in suppressing or enabling cancer development is disease specific. HIF-1α promotes glucose uptake, increases glycolysis, and activates pyruvate dehydrogenase kinase 1 (PDK1), which inactivates the TCA cycle enzyme, pyruvate dehydrogenase (PDH), inhibiting the flow of glucose into the TCA cycle via the conversion of pyruvate to acetyl-CoA (28). Blocking the conversion of pyruvate to acetyl-CoA would increase lactate production and could block usage of glucose to provide carbon backbones in the anabolic synthesis of lipids (29, 30). In a cancer that is dependent on oxidative metabolism of glucose for lipid production, such as a VHL-deficient kidney cancer that has an intact TCA cycle, HIF-1α could potentially have a suppressive effect on cell growth by reducing energy production. However, in a cancer cell with a remodeled metabolic pathway that is (i) extremely dependent on efficient use of glucose for glycolysis for generation of ATP, and (ii) is dependent on carbons from glutamine, rather than from glucose, for lipid production and for carbon backbones of amino acid precursors such as oxaloacetate, increased HIF-1α seems to be supportive of tumorigenesis. To evaluate the antitumor effectiveness of targeting the HIF-1α pathway in FH-deficient kidney cancer, silencing of HIF-1α by siRNA was conducted in UOK262 cells. Silencing of HIF-1α was found to significantly impair UOK262 invasiveness, suggesting that HIF-1α is a critical oncoprotein in FH-deficient kidney cancer (17).
FH-deficient kidney cancer is characterized by a metabolic shift to a reductive carboxylation, glutamine-dependent pathway
13C metabolite analysis was used to further characterize the metabolic shift in FH-deficient kidney cancer. Even thought FH-deficient kidney cancer has a significant decrease in the flux of glucose-derived pyruvate to the mitochondria, it is still critically dependent on mitochondrial metabolic flux for critical biosynthetic reactions. Unlike other cancers that conduct oxidative metabolism of both glucose and glutamine, FH-deficient kidney cancer is characterized by glutamine-dependent reductive carboxylation (see red arrows; Fig. 1), rather than the normal procession of the TCA cycle, for citrate formation for the citrate/acetyl-CoA/lipid synthesis pathway. FH-deficient kidney cancer cells require glutamine to proliferate, and glutamine was found to be the major carbon source for the increased fatty acid synthesis required for rapid cellular proliferation of this most aggressive form of kidney cancer (31).
Clinical–Translational Advances
KEAP1 succination and Nrf2 signaling: antioxidant response in FH-deficient kidney cancer
In complementary studies, using different approaches, Adam and colleagues and Ooi and colleagues delineated the role of the Kelch-like ECH-associated protein 1 (KEAP1) succination and nuclear factor E2-related factor 2 (Nrf2) signaling in FH-deficient kidney cancer and the role of fumarate as an oncometabolite (32, 33). Nrf2 (is a transcription factor that induces expression of a number of genes involved in antioxidant response, such as to counteract the effect of reactive oxygen species (ROS) and that facilitate heme synthesis. KEAP1 is an electrophile that combines with cullin 3 (CUL3) to facilitate degradation of Nrf2 (34). KEAP1 mutations have been found in a wide variety of cancers, including 34% of squamous lung cancers, and Nrf2 has been found to be elevated in a number of cancers, including head and neck, gallbladder, lung, and pancreatic cancer (35, 36). In a study of HLRCC (as well as sporadic) type II papillary kidney cancer, Ooi and colleagues found a gene expression pattern consistent with upregulation of antioxidant response elements (ARE), including the aldo-keto reductase family member B10 gene (AKR1B10) as well as other ARE-driven genes, NAD(P)H dehydrogenase (encoded by quinone 1 also known as NQO1) and thioredoxin reductase 1 (TXNRD1; ref. 33). Nrf2, which is known to bind the ARE enhancer and the electrophile sensor, KEAP1, forms a complex with the CUL3 ubiquitin ligase to initiate ubiquitin-imediated degradation (37). An exposed KEAP1 residue, Cys-151, reacts with fumarate, which has been shown to be elevated in FH-deficient HLRCC uterine fibroids (38) as well in FH-deficient kidney cancer cells (33), resulting in a conformational change that inhibits KEAP1 from binding to NRF2, resulting in NRF2 accumulation, and increases expression of antioxidant response genes (33, 39, 40). The findings of Ooi and colleagues show that increased levels of fumarate induce succination of KEAP1, resulting in stabilization of NRF2 and increased expression of AKR1B10. Alderson and colleagues showed that S-(2-Succinyl) cysteine (2SC) is formed by increased levels of fumarate, which inactivates the sulfhydryl enzyme, glycetaldehyde-3-phosphate dehydrogenase (41). These findings identify fumarate as an endogenous electrophile and indicate a role for fumarate as an oncoprotein in metabolic regulation (41). Nagai and colleagues subsequently showed that increase in 2SC is the result of mitochondrial stress and that increased levels of fumarate modify cysteine residues in many proteins, a process called succination (42). Bardella and colleagues showed that 2SC accumulates to high levels in FH-deficient cells and tumors and that 2SC staining could be a useful marker for detection of this disease (43).
In a parallel study, Adam and colleagues showed that FH-deficient murine renal cyst formation is Hif and Phd independent and is associated with a pronounced upregulation of antioxidant pathways (32). In the Fh1-deficient mouse, activation of the NRF2 antioxidant pathway was found to arise as a direct consequence of FH inactivation. Critical residues of KEAP1 were found to be succinated in this model, indicating that NRF2 activation results from succination of KEAP1. Impairment of KEAP1 function by succination diminishes NRF2 proteosomal degradation and thereby increases NRF2 stability and activity (ref. 32; Fig. 1).
Dysregulation of NRF2 degradation, a critical adaptive response that seems to occur early in hyperplastic cystic and tumor formation, could provide the foundation for the development of novel approaches for therapy of FH-deficient kidney cancer. Decreased KEAP1 function increases NRF2 activity, which would reduce oxidative stress and could provide a growth advantage to aggressive type II papillary kidney cancer cells. Thus, potential therapeutic approaches to target the oxidative response pathway could include inhibition of AKR1B10 or TXNRD1 (33).
FH deficiency induces HIF-1α stabilization by generation of reactive oxygen species
Another potential therapeutic approach to target FH-deficient kidney cancer involves targeting the activated antioxidant pathway by increasing reactive oxygen. Sudarshan and colleagues showed that inactivating mutations of FH-deficient kidney cancer cells results in glucose-mediated generation of cellular ROS and ROS-dependent HIF-1α stabilization and that the metabolic shift to aerobic glycolysis is critical to HIF stabilization (44). To evaluate a potential therapeutic approach for HLRCC-associated kidney cancer, Sourbier and colleagues investigated the cytotoxic effects of increasing ROS levels using the proteosome inhibitor, bortezamib, which inhibits NF-κB, in combination with cisplatin in FH-deficient cancer cells (45). Bortezamib was found to induce apoptosis in vitro and inhibited growth in vivo. Cellular ROS levels correlated with bortezomib-associated cytotoxicity. Combining other ROS inducers with bortezomib enchanced cytotoxicity, whereas combining a ROS scavenger with bortezomib inhibited its cytotoxic effect. When HLRCC murine xenografts were treated with cisplatin, a ROS inducer, and bortezomib, dependent tumor regression resulted (45).
Heme oxygenase is synthetically lethal in FH-deficient cells
Frezza and colleagues described a glutamine-to-bilirubin pathway in FH-deficient cells that involves the biosynthesis and degradation of heme, which enables partial mitochondrial NADH production and ends with bilirubin excretion from the FH-deficient cells (46). They showed that inhibition of heme oxygenase is synthetically lethal with FH-deficient cells, meaning that the FH cells cannot survive without intact heme oxygenase. Thus, inhibition of heme oxygenase may cause specific death of tumor cells in HLRCC-associated kidney cancer while sparing cells that are not FH deficient (46).
Targeting glycolysis and LDHA as a therapeutic strategy for FH-deficient kidney cancer
FH-deficient kidney cancer has undergone a metabolic shift to aerobic glycolysis and is dependent on high levels of glycolysis for ATP production needed for the rapid proliferation of this aggressive form of cancer (17, 18). A number of different approaches have been tried for targeting glucose metabolism in a Warburg cancer, including targeting critical glycolytic enzymes such as HKII, PKM2, or LDHA (reviewed in Hamanaka and Chandel; ref. 47). To evaluate whether FH-deficient cells would be sensitive to LDHA blockade, Xie and colleagues showed that LDHA is overexpressed in HLRCC kidney cancer, LDHA inhibition increased apoptosis in FH-deficient cells, this effect is ROS mediated, and LDHA knockdown in FH-deficient cells resulted in a significant decrease in xenograft tumor growth (48). Efforts are currently under way to identify an effective LDHA inhibitor for evaluation in preclinical and clinical trials.
Englerin A: limiting the cell's access to glucose
Englerin A, a sesquiterpen isolated from the root bark and stem bark of Phyllanthus engleri Pax (49), has been shown to have selective inhibition against RCC growth (50, 51). Sourbierand colleagues showed that englerin A binds to and activates protein kinase C-θ, inducing an insulin-resistant phenotype, whereas limiting the cell's access to glucose. Simultaneously, englerin A induces PKCθ activation of the transcription factor, HSF1, which induces glucose dependence. By starving the cells of glucose while simultaneously inducing glucose dependence, englerin A is lethal to a number of types of kidney cancer, including FH-deficient RCC. FH-deficient RCC cells are very sensitive to treatment with englerin A, whereas cells with FH replaced are resistant (52). Studies to further assess the effect of englerin A in FH deficient as well as other forms of kidney cancer are in progress.
Bevacizumab and erlotinib: targeting FH-deficient kidney cancer tumor vasculature with anti-VEGF reagents in patients with metastatic HLRCC-associated kidney cancer
FH-deficient kidney cancer is characterized by impaired oxidative phosphorylation, a shift to aerobic glycolysis, and a dependence on glucose for ATP generation and oxidative stress. Increased oxidative stress and/or increased levels of fumarate inhibit HIF PHD, resulting in HIF-1α stabilization, leading to increased transcription of a number of hypoxia-induced pathway genes, such as VEGF and GLUT1. HLRCC-associated kidney cancer, which is readily detected by [18] fluoro-deoxyglucose-based PET scanning, is a fast growing tumor that has a propensity to metastasize when the primary tumors are as small as .5 cm in size (9). In our experience managing patients with HLRCC since 1989, we have not found this disease to be responsive to chemotherapy, immunotherapy, or the conventional targeted therapeutic approaches in use for patients with other forms of advanced kidney cancer. On the basis of the knowledge that these cancer cells are critically dependent on tumor vasculature to ensure a high level of glucose transport needed for the increased glycolytic demand of this rapidly growing cancer, a clinical trial evaluating the effect of bevacizumab and erlotinib in patients with advanced HLRCC-kidney cancer is currently under way (53).
Conclusions
FH-deficient HLRCC-associated type II papillary kidney cancer, the most aggressive form of kidney cancer, may have an Achilles' heel, which can be successfully exploited for therapy. This early-onset cancer, which has a high propensity for metastasis, has “precisely the properties expected of a Warburg tumor” (54). In FH-deficient kidney cancer, oxidative phosphorylation is impaired because the TCA cycle is blocked by FH deficiency; FH− cells undergo a metabolic shift to aerobic glycolysis to provide the ATP needed for maintenance of highly proliferating cells. Increased levels of fumarate and decreased levels of iron in FH-deficient RCC (14) inhibit HIF PHD, resulting in stabilization of HIF-1α levels (20, 55), which would drive vascularization and glucose transport. FH-deficient kidney cancer is characterized by a shift to a reductive glutamine-dependent pathway that provides citrate, which generates the acetyl CoA needed for lipid production in rapidly growing cells (31). A number of potential therapeutic approaches to target this most aggressive form of kidney cancer, such as activation of AMPK, inhibition of LDHA, repressing the antioxidant response and the heme oxygenase pathway, have been developed. One hopes that these approaches, as well as the clinical trial under way targeting the vasculature and glucose transport in HLRCC-associated kidney cancer with bevacizumab and erlotinib, will provide the foundation for the development of an effective form of therapy for patients affected with this most aggressive form of Krebs cycle mutation cancer and could provide insight into the management of patients with other forms of cancer characterized by a metabolic shift to aerobic glycolysis and impaired oxidative phosphorylation.
Authors' Contributions
Conception and design: W.M. Linehan
Development of methodology: W.M. Linehan, T.A. Rouault
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): W.M. Linehan
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): W.M. Linehan, T.A. Rouault
Writing, review, and/or revision of the manuscript: W.M. Linehan, T.A. Rouault
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): W.M. Linehan
Acknowledgments
The authors thank Georgia Flaten Shaw for the outstanding editorial and graphics support and Christopher J. Ricketts for critical review of the manuscript.
Grant Support
This work was supported by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research.