In Cerebellar Atrophy of 12-Month-Old ATM-Null Mice, Transcriptome Upregulations Concern Most Neurotransmission and Neuropeptide Pathways, While Downregulations Affect Prominently Itpr1, Usp2 and Non-Coding RNA.

The autosomal recessive disorder Ataxia-Telangiectasia is caused by a dysfunction of the stress response protein, ATM. In the nucleus of proliferating cells, ATM senses DNA double-strand breaks and coordinates their repair. This role explains T-cell dysfunction and tumour risk. However, it remains unclear whether this function is relevant for postmitotic neurons and underlies cerebellar atrophy, since ATM is cytoplasmic in postmitotic neurons. Here, we used ATM-null mice that survived early immune deficits via bone-marrow transplantation, and that reached initial neurodegeneration stages at 12 months of age. Global cerebellar transcriptomics demonstrated that ATM depletion triggered upregulations in most neurotransmission and neuropeptide systems. Downregulated transcripts were found for the ATM interactome component Usp2, many non-coding RNAs, ataxia genes Itpr1, Grid2, immediate early genes and immunity factors. Allelic splice changes affected prominently the neuropeptide machinery, e.g., Oprm1. Validation experiments with stressors were performed in human neuroblastoma cells, where ATM was localised only to cytoplasm, similar to the brain. Effect confirmation in SH-SY5Y cells occurred after ATM depletion and osmotic stress better than nutrient/oxidative stress, but not after ATM kinase inhibition or DNA stressor bleomycin. Overall, we provide pioneer observations from a faithful A-T mouse model, which suggest general changes in synaptic and dense-core vesicle stress adaptation.

PhosY-secretome profiling combined with kinase-substrate interaction screening defines active c-Src-driven extracellular signaling.

c-Src tyrosine kinase is a renowned key intracellular signaling molecule and a potential target for cancer therapy. Secreted c-Src is a recent observation, but how it contributes to extracellular phosphorylation remains elusive. Using a series of domain deletion mutants, we show that the N-proximal region of c-Src is essential for its secretion. The tissue inhibitor of metalloproteinases 2 (TIMP2) is an extracellular substrate of c-Src. Limited proteolysis-coupled mass spectrometry and mutagenesis studies verify that the Src homology 3 (SH3) domain of c-Src and the P31VHP34 motif of TIMP2 are critical for their interaction. Comparative phosphoproteomic analyses identify an enrichment of PxxP motifs in phosY-containing secretomes from c-Src-expressing cells with cancer-promoting roles. Inhibition of extracellular c-Src using custom SH3-targeting antibodies disrupt kinase-substrate complexes and inhibit cancer cell proliferation. These findings point toward an intricate role for c-Src in generating phosphosecretomes, which will likely influence cell-cell communication, particularly in c-Src-overexpressing cancers.

Development and Characterization of Selective FAK Inhibitors and PROTACs with In Vivo Activity.

Focal adhesion kinase (FAK) is an attractive drug target due to its overexpression in cancer. FAK functions as a non-receptor tyrosine kinase and scaffolding protein, coordinating several downstream signaling effectors and cellular processes. While drug discovery efforts have largely focused on targeting FAK kinase activity, FAK inhibitors have failed to show efficacy as single agents in clinical trials. Here, using structure-guided design, we report the development of a selective FAK inhibitor (BSJ-04-175) and degrader (BSJ-04-146) to evaluate the consequences and advantages of abolishing all FAK activity in cancer models. BSJ-04-146 achieves rapid and potent FAK degradation with high proteome-wide specificity in cancer cells and induces durable degradation in mice. Compared to kinase inhibition, targeted degradation of FAK exhibits pronounced improved activity on downstream signaling and cancer cell viability and migration. Together, BSJ-04-175 and BSJ-04-146 are valuable chemical tools to dissect the specific consequences of targeting FAK through small-molecule inhibition or degradation.

Global Mass Spectrometry-Based Analysis of Protein Ubiquitination Using K-ε-GG Remnant Antibody Enrichment.

Ubiquitination is a post-translational modification that affects protein degradation as well as a variety of cellular processes. Methods that globally profile ubiquitination are powerful tools to better understand these processes. Here we describe an updated method for identification and quantification of thousands of sites of ubiquitination from cells, tissues, or other biological materials. The method involves cell lysis and digestion to peptides, immunoaffinity enrichment with an antibody recognizing di-glycine remnants left behind at ubiquitinated lysines, and liquid chromatography-tandem mass spectrometry analysis of the enriched peptides.

Global Proteome Profiling to Assess Changes in Protein Abundance Using Isobaric Labeling and Liquid Chromatography-Tandem Mass Spectrometry.

Protein degradation is a critical component of all facets of cell biology, and recently methods have been developed to make use of targeted protein degradation as both an investigative tool and a potential therapeutic avenue. Mass spectrometry-based proteomic studies have allowed detailed characterization of changes in protein level and the biology underlying growth, development, and disease. Current methods and instrumentation allow identification and quantitative analysis of thousands of proteins in a single assay. The method described here involves cell lysis and digestion to peptides, labeling peptides with isobaric tagging TMT reagents, basic reversed phase fractionation, and liquid chromatography-tandem mass spectrometry analysis of the enriched peptides.

Acetylation-mediated remodeling of the nucleolus regulates cellular acetyl-CoA responses.

The metabolite acetyl-coenzyme A (acetyl-CoA) serves as an essential element for a wide range of cellular functions including adenosine triphosphate (ATP) production, lipid synthesis, and protein acetylation. Intracellular acetyl-CoA concentrations are associated with nutrient availability, but the mechanisms by which a cell responds to fluctuations in acetyl-CoA levels remain elusive. Here, we generate a cell system to selectively manipulate the nucleo-cytoplasmic levels of acetyl-CoA using clustered regularly interspaced short palindromic repeat (CRISPR)-mediated gene editing and acetate supplementation of the culture media. Using this system and quantitative omics analyses, we demonstrate that acetyl-CoA depletion alters the integrity of the nucleolus, impairing ribosomal RNA synthesis and evoking the ribosomal protein-dependent activation of p53. This nucleolar remodeling appears to be mediated through the class IIa histone deacetylases (HDACs). Our findings highlight acetylation-mediated control of the nucleolus as an important hub linking acetyl-CoA fluctuations to cellular stress responses.

Widespread protein lysine acetylation in gut microbiome and its alterations in patients with Crohn's disease.

Lysine acetylation (Kac), an abundant post-translational modification (PTM) in prokaryotes, regulates various microbial metabolic pathways. However, no studies have examined protein Kac at the microbiome level, and it remains unknown whether Kac level is altered in patient microbiomes. Herein, we use a peptide immuno-affinity enrichment strategy coupled with mass spectrometry to characterize protein Kac in the microbiome, which successfully identifies 35,200 Kac peptides from microbial or human proteins in gut microbiome samples. We demonstrate that Kac is widely distributed in gut microbial metabolic pathways, including anaerobic fermentation to generate short-chain fatty acids. Applying to the analyses of microbiomes of patients with Crohn's disease identifies 52 host and 136 microbial protein Kac sites that are differentially abundant in disease versus controls. This microbiome-wide acetylomic approach aids in advancing functional microbiome research.

mTOR and S6K1 drive polycystic kidney by the control of Afadin-dependent oriented cell division.

mTOR activation is essential and sufficient to cause polycystic kidneys in Tuberous Sclerosis Complex (TSC) and other genetic disorders. In disease models, a sharp increase of proliferation and cyst formation correlates with a dramatic loss of oriented cell division (OCD). We find that OCD distortion is intrinsically due to S6 kinase 1 (S6K1) activation. The concomitant loss of S6K1 in Tsc1-mutant mice restores OCD but does not decrease hyperproliferation, leading to non-cystic harmonious hyper growth of kidneys. Mass spectrometry-based phosphoproteomics for S6K1 substrates revealed Afadin, a known component of cell-cell junctions required to couple intercellular adhesions and cortical cues to spindle orientation. Afadin is directly phosphorylated by S6K1 and abnormally decorates the apical surface of Tsc1-mutant cells with E-cadherin and α-catenin. Our data reveal that S6K1 hyperactivity alters centrosome positioning in mitotic cells, affecting oriented cell division and promoting kidney cysts in conditions of mTOR hyperactivity.

Phospho-PTM proteomic discovery of novel EPO- modulated kinases and phosphatases, including PTPN18 as a positive regulator of EPOR/JAK2 Signaling.

The formation of erythroid progenitor cells depends sharply upon erythropoietin (EPO), its cell surface receptor (erythropoietin receptor, EPOR), and Janus kinase 2 (JAK2). Clinically, recombinant human EPO (rhEPO) additionally is an important anti-anemia agent for chronic kidney disease (CKD), myelodysplastic syndrome (MDS) and chemotherapy, but induces hypertension, and can exert certain pro-tumorigenic effects. Cellular signals transduced by EPOR/JAK2 complexes, and the nature of EPO-modulated signal transduction factors, therefore are of significant interest. By employing phospho-tyrosine post-translational modification (p-Y PTM) proteomics and human EPO- dependent UT7epo cells, we have identified 22 novel kinases and phosphatases as novel EPO targets, together with their specific sites of p-Y modification. New kinases modified due to EPO include membrane palmitoylated protein 1 (MPP1) and guanylate kinase 1 (GUK1) guanylate kinases, together with the cytoskeleton remodeling kinases, pseudopodium enriched atypical kinase 1 (PEAK1) and AP2 associated kinase 1 (AAK1). Novel EPO- modified phosphatases include protein tyrosine phosphatase receptor type A (PTPRA), phosphohistidine phosphatase 1 (PHPT1), tensin 2 (TENC1), ubiquitin associated and SH3 domain containing B (UBASH3B) and protein tyrosine phosphatase non-receptor type 18 (PTPN18). Based on PTPN18's high expression in hematopoietic progenitors, its novel connection to JAK kinase signaling, and a unique EPO- regulated PTPN18-pY389 motif which is modulated by JAK2 inhibitors, PTPN18's actions in UT7epo cells were investigated. Upon ectopic expression, wt-PTPN18 promoted EPO dose-dependent cell proliferation, and survival. Mechanistically, PTPN18 sustained the EPO- induced activation of not only mitogen-activated protein kinases 1 and 3 (ERK1/2), AKT serine/threonine kinase 1-3 (AKT), and signal transducer and activator of transcription 5A and 5B (STAT5), but also JAK2. Each effect further proved to depend upon PTPN18's EPO- modulated (p)Y389 site. In analyses of the EPOR and the associated adaptor protein RHEX (regulator of hemoglobinization and erythroid cell expansion), wt-PTPN18 increased high molecular weight EPOR forms, while sharply inhibiting the EPO-induced phosphorylation of RHEX-pY141. Each effect likewise depended upon PTPN18-Y389. PTPN18 thus promotes signals for EPO-dependent hematopoietic cell growth, and may represent a new druggable target for myeloproliferative neoplasms.

CBP/p300 Drives the Differentiation of Regulatory T Cells through Transcriptional and Non-Transcriptional Mechanisms.

Regulatory T cells (Treg) are immunosuppressive and negatively impact response to cancer immunotherapies. CREB-binding protein (CBP) and p300 are closely related acetyltransferases and transcriptional coactivators. Here, we evaluate the mechanisms by which CBP/p300 regulate Treg differentiation and the consequences of CBP/p300 loss-of-function mutations in follicular lymphoma. Transcriptional and epigenetic profiling identified a cascade of transcription factors essential for Treg differentiation. Mass spectrometry analysis showed that CBP/p300 acetylates prostacyclin synthase, which regulates Treg differentiation by altering proinflammatory cytokine secretion by T and B cells. Reduced Treg presence in tissues harboring CBP/p300 loss-of-function mutations was observed in follicular lymphoma. Our findings provide novel insights into the regulation of Treg differentiation by CBP/p300, with potential clinical implications on alteration of the immune landscape. SIGNIFICANCE: This study provides insights into the dynamic role of CBP/p300 in the differentiation of Tregs, with potential clinical implications in the alteration of the immune landscape in follicular lymphoma.

Heat stress induced, ligand-independent MET and EGFR signalling in hepatocellular carcinoma.

PURPOSE: The aims of the present study were 2-fold: first, to test the hypothesis that heat stress induces MET and EGFR signalling in hepatocellular carcinoma (HCC) cells and inhibition of this signalling decreases HCC clonogenic survival; and second, to identify signalling pathways associated with heat stress induced MET signalling. MATERIALS AND METHODS: MET+ and EGFR+ HCC cells were pre-treated with inhibitors to MET, EGFR, PI3K/mTOR or vehicle and subjected to heat stress or control ± HGF or EGF growth factors and assessed by colony formation assay, Western blotting and/or quantitative mass spectrometry. IACUC approved partial laser thermal or sham ablation was performed on orthotopic N1S1 and AS30D HCC tumours and liver/tumour assessed for phospho-MET and phospho-EGFR immunostaining. RESULTS: Heat-stress induced rapid MET and EGFR phosphorylation that is distinct from HGF or EGF in HCC cells and thermal ablation induced MET but not EGFR phosphorylation at the HCC tumour ablation margin. Inhibition of the MET and EGFR blocked both heat stress and growth factor induced MET and EGFR phosphorylation and inhibition of MET decreased HCC clonogenic survival following heat stress. Pathway analysis of quantitative phosphoproteomic data identified downstream pathways associated with heat stress induced MET signalling including AKT, ERK, Stat3 and JNK. However, inhibition of heat stress induced MET signalling did not block AKT signalling. CONCLUSIONS: Heat-stress induced MET and EGFR signalling is distinct from growth factor mediated signalling in HCC cells and MET inhibition enhances heat stress induced HCC cell killing via a PI3K/AKT/mTOR-independent mechanism.

Identification of Post-Translational Modifications from Serum/Plasma by Immunoaffinity Enrichment and LC-MS/MS Analysis Without Depletion of Abundant Proteins.

Immunoaffinity enrichment combined with LC-MS/MS enables identification of Post-translational modifications (PTMs) from serum/plasma samples without abundant protein depletion. Here we described the workflow in details in identifying various types of PTMs such as lysine acetylation and arginine methylation from cancer serum. The method described is compatible with all common proteomic analysis platforms and quantitative methods.

Role of the E3 ubiquitin ligase RNF157 as a novel downstream effector linking PI3K and MAPK signaling pathways to the cell cycle.

The interconnected PI3K and MAPK signaling pathways are commonly perturbed in cancer. Dual inhibition of these pathways by the small-molecule PI3K inhibitor pictilisib (GDC-0941) and the MEK inhibitor cobimetinib (GDC-0973) suppresses cell proliferation and induces cell death better than either single agent in several preclinical models. Using mass spectrometry-based phosphoproteomics, we have identified the RING finger E3 ubiquitin ligase RNF157 as a target at the intersection of PI3K and MAPK signaling. We demonstrate that RNF157 phosphorylation downstream of the PI3K and MAPK pathways influences the ubiquitination and stability of RNF157 during the cell cycle in an anaphase-promoting complex/cyclosome-CDH1-dependent manner. Deletion of these phosphorylation-targeted residues on RNF157 disrupts binding to CDH1 and protects RNF157 from ubiquitination and degradation. Expression of the cyclin-dependent kinase 2 (CDK2), itself a downstream target of PI3K/MAPK signaling, leads to increased phosphorylation of RNF157 on the same residues modulated by PI3K and MAPK signaling. Inhibition of PI3K and MEK in combination or of CDK2 by their respective small-molecule inhibitors reduces RNF157 phosphorylation at these residues and attenuates RNF157 interaction with CDH1 and its subsequent degradation. Knockdown of endogenous RNF157 in melanoma cells leads to late S phase and G2/M arrest and induces apoptosis, the latter further potentiated by concurrent PI3K/MEK inhibition, consistent with a role for RNF157 in the cell cycle. We propose that RNF157 serves as a novel node integrating oncogenic signaling pathways with the cell cycle machinery and promoting optimal cell cycle progression in transformed cells.

Phosphoproteome analysis of the MAPK pathway reveals previously undetected feedback mechanisms.

The RAS-RAF-MEK-ERK (MAPK) pathway is prevalently perturbed in cancer. Recent large-scale sequencing initiatives profiled thousands of tumors providing insight into alterations at the DNA and RNA levels. These efforts confirmed that key nodes of the MAPK pathway, in particular KRAS and BRAF, are among the most frequently altered proteins in cancer. The establishment of targeted therapies, however, has proven difficult. To decipher the underlying challenges, it is essential to decrypt the phosphorylation network spanned by the MAPK core axis. Using mass spectrometry we identified 2241 phosphorylation sites on 1020 proteins, and measured their responses to inhibition of MEK or ERK. Multiple phosphorylation patterns revealed previously undetected feedback, as upstream signaling nodes, including receptor kinases, showed changes at the phosphorylation level. We provide a dataset rich in potential therapeutic targets downstream of the MAPK cascade. By integrating TCGA (The Cancer Genome Atlas) data, we highlight some downstream phosphoproteins that are frequently altered in cancer. All MS data have been deposited in the ProteomeXchange with identifier PXD003908 (http://proteomecentral.proteomexchange.org/dataset/PXD003908).

Quantitative phosphoproteomic analysis of the PI3K-regulated signaling network.

The PI3K pathway is commonly activated in cancer. Only a few studies have attempted to explore the spectrum of phosphorylation signaling downstream of the PI3K cascade. Such insight, however, is imperative to understand the mechanisms responsible for oncogenic phenotypes. By applying MS-based phosphoproteomics, we mapped 2509 phosphorylation sites on 1096 proteins, and quantified their responses to activation or inhibition of PIK3CA using isogenic knock-in derivatives and a series of targeted inhibitors. We uncovered phosphorylation changes in a wide variety of proteins involved in cell growth and proliferation, many of which have not been previously associated with PI3K signaling. A significant update of the posttranslational modification database PHOSIDA (http://www.phosida.com) allows efficient use of the data. All MS data have been deposited in the ProteomeXchange with identifier PXD003899 (http://proteomecentral.proteomexchange.org/dataset/PXD003899).

Quantitative Profiling of Post-translational Modifications by Immunoaffinity Enrichment and LC-MS/MS in Cancer Serum without Immunodepletion.

A robust method was developed and optimized for enrichment and quantitative analysis of posttranslational modifications (PTMs) in serum/plasma samples by combining immunoaffinity purification and LC-MS/MS without depletion of abundant proteins. The method was used to survey serum samples of patients with acute myeloid leukemia (AML), breast cancer (BC), and nonsmall cell lung cancer (NSCLC). Peptides were identified from serum samples containing phosphorylation, acetylation, lysine methylation, and arginine methylation. Of the PTMs identified, lysine acetylation (AcK) and arginine mono-methylation (Rme) were more prevalent than other PTMs. Label-free quantitative analysis of AcK and Rme peptides was performed for sera from AML, BC, and NSCLC patients. Several AcK and Rme sites showed distinct abundance distribution patterns across the three cancer types. The identification and quantification of posttranslationally modified peptides in serum samples reported here can be used for patient profiling and biomarker discovery research.

A respiratory chain controlled signal transduction cascade in the mitochondrial intermembrane space mediates hydrogen peroxide signaling.

Reactive oxygen species (ROS) such as hydrogen peroxide (H2O2) govern cellular homeostasis by inducing signaling. H2O2 modulates the activity of phosphatases and many other signaling molecules through oxidation of critical cysteine residues, which led to the notion that initiation of ROS signaling is broad and nonspecific, and thus fundamentally distinct from other signaling pathways. Here, we report that H2O2 signaling bears hallmarks of a regular signal transduction cascade. It is controlled by hierarchical signaling events resulting in a focused response as the results place the mitochondrial respiratory chain upstream of tyrosine-protein kinase Lyn, Lyn upstream of tyrosine-protein kinase SYK (Syk), and Syk upstream of numerous targets involved in signaling, transcription, translation, metabolism, and cell cycle regulation. The active mediators of H2O2 signaling colocalize as H2O2 induces mitochondria-associated Lyn and Syk phosphorylation, and a pool of Lyn and Syk reside in the mitochondrial intermembrane space. Finally, the same intermediaries control the signaling response in tissues and species responsive to H2O2 as the respiratory chain, Lyn, and Syk were similarly required for H2O2 signaling in mouse B cells, fibroblasts, and chicken DT40 B cells. Consistent with a broad role, the Syk pathway is coexpressed across tissues, is of early metazoan origin, and displays evidence of evolutionary constraint in the human. These results suggest that H2O2 signaling is under control of a signal transduction pathway that links the respiratory chain to the mitochondrial intermembrane space-localized, ubiquitous, and ancient Syk pathway in hematopoietic and nonhematopoietic cells.

Neuroblastoma tyrosine kinase signaling networks involve FYN and LYN in endosomes and lipid rafts.

Protein phosphorylation plays a central role in creating a highly dynamic network of interacting proteins that reads and responds to signals from growth factors in the cellular microenvironment. Cells of the neural crest employ multiple signaling mechanisms to control migration and differentiation during development. It is known that defects in these mechanisms cause neuroblastoma, but how multiple signaling pathways interact to govern cell behavior is unknown. In a phosphoproteomic study of neuroblastoma cell lines and cell fractions, including endosomes and detergent-resistant membranes, 1622 phosphorylated proteins were detected, including more than half of the receptor tyrosine kinases in the human genome. Data were analyzed using a combination of graph theory and pattern recognition techniques that resolve data structure into networks that incorporate statistical relationships and protein-protein interaction data. Clusters of proteins in these networks are indicative of functional signaling pathways. The analysis indicates that receptor tyrosine kinases are functionally compartmentalized into distinct collaborative groups distinguished by activation and intracellular localization of SRC-family kinases, especially FYN and LYN. Changes in intracellular localization of activated FYN and LYN were observed in response to stimulation of the receptor tyrosine kinases, ALK and KIT. The results suggest a mechanism to distinguish signaling responses to activation of different receptors, or combinations of receptors, that govern the behavior of the neural crest, which gives rise to neuroblastoma.

Complementary PTM Profiling of Drug Response in Human Gastric Carcinoma by Immunoaffinity and IMAC Methods with Total Proteome Analysis.

Gaining insight into normal cellular signaling and disease biology is a critical goal of proteomic analyses. The ability to perform these studies successfully to extract the maximum value and discovery of biologically relevant candidate biomarkers is therefore of primary importance. Many successful studies in the past have focused on total proteome analysis (changes at the protein level) combined with phosphorylation analysis by metal affinity enrichment (changes at the PTM level). Here, we use the gastric carcinoma cell line MKN-45 treated with the c-Met inhibitor SU11274 and PKC inhibitor staurosporine to investigate the most efficient and most comprehensive strategies for both total protein and PTM analysis. Under the conditions used, total protein analysis yielded few changes in response to either compound, while analysis of phosphorylation identified thousands of sites that changed differentially between the two treatments. Both metal affinity and antibody-based enrichments were used to assess phosphopeptide changes, and the data generated by the two methods was largely complementary (non-overlapping). Label-free quantitation of peptide peak abundances was used to accurately determine fold-changes between control and treated samples. Protein interaction network analysis allowed the data to be placed in a biologically relevant context, and follow-up validation of selected findings confirmed the accuracy of the proteomic data. Together, this study provides a framework for start-to-finish proteomic analysis of any experimental system under investigation to maximize the value of the proteomic study and yield the best chance for uncovering actionable target candidates.

Interrogating cAMP-dependent kinase signaling in Jurkat T cells via a protein kinase A targeted immune-precipitation phosphoproteomics approach.

In the past decade, mass-spectrometry-based methods have emerged for the quantitative profiling of dynamic changes in protein phosphorylation, allowing the behavior of thousands of phosphorylation sites to be monitored in a single experiment. However, when one is interested in specific signaling pathways, such shotgun methodologies are not ideal because they lack selectivity and are not cost and time efficient with respect to instrument and data analysis time. Here we evaluate and explore a peptide-centric antibody generated to selectively enrich peptides containing the cAMP-dependent protein kinase (PKA) consensus motif. This targeted phosphoproteomic strategy is used to profile temporal quantitative changes of potential PKA substrates in Jurkat T lymphocytes upon prostaglandin E2 (PGE2) stimulation, which increases intracellular cAMP, activating PKA. Our method combines ultra-high-specificity motif-based immunoaffinity purification with cost-efficient stable isotope dimethyl labeling. We identified 655 phosphopeptides, of which 642 (i.e. 98%) contained the consensus motif [R/K][R/K/X]X[pS/pT]. When our data were compared with a large-scale Jurkat T-lymphocyte phosphoproteomics dataset containing more than 10,500 phosphosites, a minimal overlap of 0.2% was observed. This stresses the need for such targeted analyses when the interest is in a particular kinase. Our data provide a resource of likely substrates of PKA, and potentially some substrates of closely related kinases. Network analysis revealed that about half of the observed substrates have been implicated in cAMP-induced signaling. Still, the other half of the here-identified substrates have been less well characterized, representing a valuable resource for future research.