2. Materials and Methods
Experimental Groups
Forty adult female dogs were distributed into the following groups – (G1) 8 female dogs with malignant mammary neoplasia (Simple Carcinomas); (G2) 8 female dogs with benign mammary neoplasia (Carcinoma in mixed tumor); (G3) 8 female dogs with tumor relapse after mastectomy (metastasis); and (G4) 8 female dogs in clinical remission (over a year following mastectomy of malignant tumors); and (G5) 8 healthy female controls. For G1 and G2 groups, blood was collected at the time of diagnosis of mammary tumor (G1.1 and G2.1) and at the time of surgical suture removal (G1.2 and G2.2), approximately 10 days after mastectomy. G3, G4, and G5 groups had only one collection. A summary of our methodological approach is shown in
Figure 1. The average age of the female dogs in all groups was 8.5 years, with mixed breed dogs being the most prevalent among all the breeds.
The inclusion criteria for the neoplastic group were adult dogs with a clinical and histological diagnosis of CMT, with or without metastasis. The control group consisted of healthy adult spayed patients without comorbidities or treatments. Histopathological grading was performed using the Goldschmidt classification (2011) [
11], which is based on Nottingham’s classification.
Liquid Biopsy
Three to five milliliters of blood were collected into tubes with EDTA anticoagulant and immediately centrifuged twice at 2,500g at 23˚C for 15 min to obtain platelet-poor plasma. Subsequently, the plasma underwent to three consecutive centrifugations at 4˚C, first at 300g for 10 min to remove cells, then at 2,500g for 10 min to remove cellular debris, and finally at 16,500g for 30 min to remove small extracellular vesicles (SEVs) larger than 150 nm. After five centrifugation cycles, the plasma was aliquoted and stored at -80˚C for subsequent isolation and characterization of the SEVs.
Tissue samples (tumor fragment) and histopathology
The tumor fragment obtained during surgical excision was fixed in 10% formaldehyde for 24 h and then transferred to 70% alcohol for subsequent histopathological diagnosis. Histopathological analysis of the canine patients was conducted at the PatAnimal Laboratory in São José do Rio Preto, SP. The samples were processed following the standard protocol, which included sectioning in a microtome and staining with hematoxylin and eosin (HE). The slides were examined by two independent veterinary pathologists using a common Nikon E20 optical microscope.
Exosomes Isolation
Exosome isolation was performed using ultracentrifugation (UC) (Optima XE-90 Ultracentrifuge; rotor 70 Ti; Beckman Coulter) at the Laboratory of Molecular Morphophysiology and Development (LMMD/ZMV), University of São Paulo, Pirassununga/SP/Brazil. To obtain exosome-enriched content, the pellets were filtered through a 0.22 μm pore filter and ultracentrifuged twice at 119,700 g for 70 minutes at 4 ºC. Following the second ultracentrifugation, the pellet was diluted in phosphate-saline (PBS) and stored at -80 ºC.
Nanoparticle Tracking Analysis (NTA)
The pellets containing SEVs were characterized for size, morphology, and quantity. Initially, they were diluted in 50 µL of magnesium-calcium-free PBS at a dilution factor of 1:500 in PBS. A laser beam was passed through a sampling chamber containing this liquid suspension. The suspended particles scattered the light, allowing visualization under a microscope at 20x magnification. A video camera (sCMOS in Camera Level 15 at a temperature of 37 ºC) captured 5 images of 30 seconds of these particles in Brownian motion. Using the Nano-Sight software (NS300; NTA 3.1 Build 3.1.45; Malvern), the particles were individually tracked, and the hydrodynamic diameter of the particles was calculated using the Stokes-Einstein equation.
Transmission Electron Microscopy (TEM)
Following ultracentrifugation, the pellet was diluted in a fixative solution (0.1M sodium cacodylate, 2.5% Glutaraldehyde, 4% paraformaldehyde, pH between 7.2 and 7.4) for 2 hours at room temperature (RT). Then, 2 mL of ultrapure water were added followed by centrifugation at 119,700 g for 70 minutes at 4 ºC to remove the fixative solution. The pellet was diluted in 20 uL of ultrapure water and stored at 4ºC. The contents were placed on a copper grid coated with pioloform for 5 min, and the excess was removed with moist filter paper. The grid was then inserted into a drop of 2% aqueous uranyl acetate for 3 min, removing the excess again with moist filter paper. Transmission electron microscopy (FEI 200kV, model Tecnai 20, emitter LAB6) was performed at the Multiuser Electron Microscopy Laboratory, USP/Ribeirão Preto/Brazil. Images were captured at 140,000 X magnification.
Proteomics analysis
Mass spectrometry analysis was conducted using the nanoElute nanoflow chromatographic system from Bruker Daltonics, Bremen, Germany, coupled online to a hybrid trapped ion mobility spectrometry-quadrupole time-of-flight mass spectrometer-timsTof Pro mass spectrometer, also from Bruker Daltonics at the Max Feffer Genetics of Plants Laboratory, Departament of Genetics, USP/Piracicaba/SP/Brazil. An aliquot (1µL) of the sample, equivalent to 200 ng of digested peptides, was injected into an Aurora 2 C18 trap column (1.6µm, 250mm x 75µm), from ionOpticks, Australia. A typical reverse-phase gradient (Solvent A: 0.1% AF, 99.9% H2O MilliQ; Solvent B: 0.1% AF, 99.9% CH3CN) was established in a liquid chromatography nanoflow system and separated at a flow rate of 400 nL.min. 1. The column temperature was maintained at 50°C. The chromatographic run was 120 min (2% to 15% of Solvent B for 60 min; raised to 25% at 90 min; raised to 37% at 100 min; raised to 95% at 110 min and finally 95% for 10 min for column washing). The column was coupled online with a timsTOF-Pro with a CaptiveSpray ion source, both from Bruker Daltonik GmbH. The temperature of the capillary ion transfer line was set at 180°C. Ion accumulation for 123ms and mobility separation were obtained with an input potential ramp from -160V to -20V within 123s. During acquisition, to enable the parallel accumulation-serial fragmentation (PASEF) method, precursor m/z and mobility information were first derived from a Tims-MS full scan experiment, with an m/z range of 100-1700. Monocharged precursors were excluded based on their position in the m/z-ion plane of mobility, and precursors that reached the target value of 20,000 a.u. were dynamically excluded for 0.4 min. The operational mode of TIMS-ToF, MS, and PASEF were controlled and synchronized using OtofControl 5.1 instrumental control software from Bruker Daltonik.
For data processing, the raw files were analyzed using MaxQuant software, version 2.4.0.0. Perseus analysis software, version 2.0.9.0 (Tyanova and Cox, 2018), was used to rigorously filter the proteomic results (90%), followed by a refinement using a script in the R programming language (
https://www.R-project.org/). Processed peak intensities were then further submitted to Metaboanalyst 5.0 [
12] to generate heatmaps after Pareto scaling, presenting only the 25 most significant proteins. Confirmed proteins retrieved from the Uniprot database for
Canis lupus familiaris (CANF) were displayed using their gene names, while predicted, inferred, uncertain or Ig-like proteins were indicated by their primary accession. Exclusive and shared proteins were visualized using Venn diagrams (
https://bioinformatics.psb.ugent.be/webtools/Venn/). Gene Ontology Enrichment Analysis [
13] was performed for the unique features found in each group, with outputs having a false discovery rate (FDR) < 0.05 considered enriched, and -log (p value) reported. To assess protein interaction or network formation, the protein list was submitted to the STRING database [
14]. Finally, unique proteins were cross-referenced in the Uniprot database, and squares were inserted in the figures to show only the confirmed, predicted, or inferred proteins while uncertain or Ig-like proteins were not shown.
Statistical Analysis
Statistical calculations were performed for SEVs using and GraphPad Prism software (v.10.2.3). All quantitative data were expressed either as the median with interquartile ranges or as means with standard errors. Quantitative significance was performed by One-way ANOVA for comparison between groups followed by Tukey’s post-hoc test. Comparison between two paired groups was determined by Student's t test. p <0.05 was considered statistically significant.
4. Discussion
Our findings demonstrate a diverse array of proteins present in each investigated condition. As depicted in
Figure 3, the heatmap reveals distinct patterns of protein expression between healthy dogs and those with benign and malignant tumors, suggesting significant alterations in the proteomic content of SEVs in different pathological states. Additionally, the Venn diagram highlights exclusive proteins found in each group, providing information on potential molecular markers associated with each condition. Gene ontology enrichment analysis reveals specific biological pathways and molecular functions that may be affected in malignant tumors, indicating potential therapeutic targets. Similarly,
Figure 4,
Figure 5 and
Figure 6 provide valuable information into changes in the proteomic profile of SEVs at different stages of the disease and treatment response. These results suggest that SEVs may play a significant role in the progression of canine mammary cancer and could serve as useful biomarkers for monitoring treatment efficacy and predicting patient prognosis.
EVs have emerged as pivotal mediators of intercellular communication and are increasingly recognized as potential source of biomarkers for liquid biopsy in various cancer types [
16]. They exhibit selectivity in tumor cell signaling and are continuously secreted into body fluids from the early stages of the disease [
16,
17]. Aguilera-Rojas and colleagues (2018) [
17] isolated SEVs from blood serum of dogs and reported a nanoparticle concentration (xE10/mL) ranging from 107.4 +/− 6.8 to 403.2 +/− 25.8 in non-cancer patients, and 225.6 +/− 10.4 to 500.4 +/− 76.4 in cancer patients, respectively. However, no significant difference was observed between non-cancerous and cancerous samples, despite notable variation within individual samples [
18]. Similarly, our results demonstrated that the SEVs concentration did not exhibit statistically significant difference in the canine cohorts, either before or after surgical excision of the primary tumor. In fact, the concentration before mastectomy was higher than after mastectomy, considering individual values. However, SEVs concentration in the remission, relapse, and control groups was even higher than that obtained for the mammary tumor groups. Like Aguilera-Rojas and colleagues (2018) [
17], we also noted considerable variation within individual samples, despite the lack of statistical significance. The constant influx of SEVs, the dynamic release-uptake dynamics of different cells, and the lack of precise characterization of SEVs' origin make it challenging to determine whether the quantity of SEVs in cancer patients differs significantly from that in healthy control patients.
On the contrary, scientific evidence suggests that the cargo of SEVs, comprising proteins, nucleic acids, and lipids, may be more relevant than merely their quantity or size in discerning between cancer and non-cancer patients. Their cargo could impact the signaling of the host's immune defenses, facilitating the spread of the neoplasm by priming distant environments for metastasis [
19,
20]. This aspect is particularly crucial for identifying biomarkers and gaining a better understanding of cancer pathophysiology. Consequently, the number of studies focusing on the proteomics of SEVs in patients diagnosed with mammary tumor has increased, with the potential for early diagnosis and longitudinal prognostic evaluation [
9].
In this study we observed that the protein characteristics of SEV’s cargo can be more discriminatory than their mere presence or absence. The relative abundance of proteins among the various groups did not exhibit significant differences, rendering it an insignificant factor within the study's context. However, the pertinent observation lies in the consistent detection of distinct proteins across all groups, despite certain groups sharing common proteins, presumably linked to normal cellular constituents of the species. This discernment has enabled the identification of potential biomarkers for diagnostic and prognostic purposes, and conceivably for guiding treatment modalities.
The proteins identified in blood samples from female dogs, particularly in the malignant mammary tumor group, encompass various roles in cancer biology. One of these very abundant proteins, inter-alpha-trypsin inhibitor heavy chain (ITIH2) contributes to matrix stability and may serve as a hyaluronan transporter or binding protein [
21,
22]. These proteins can be positively or negatively regulated, but there is already evidence that all members of the ITIH family play a significant role in cell malignant processes and tumor growth [
22]. Kopylov and colleagues (2020) [
22] studied the proteomics of plasma samples from colorectal cancer patients at different stages and healthy patients. They observed significant changes in the levels of inter-alpha-trypsin inhibitor heavy chains (such as ITIH2) due to their implication in tumor growth and the malignancy process [
22]. The literature describes that ITIH2 may present significant downregulation in solid tumors [
23], contrary to the results obtained in the study by Kopylov and colleagues (2020) [
22] and in our study, which observed a greater abundance of the ITIH2 protein in the malignant carcinoma group. According to Hamm (2008) [
20] ITIH2 expression was clearly detectable in normal breast epithelium, in hyperplastic gland epithelium, and in ductal carcinoma in-situ. However, in 44% of invasive breast carcinomas, ITIH2 expression was strongly reduced or completely lost, whereas 56% of invasive carcinomas maintained moderate to strong ITIH2 expression. Furthermore, strong expression of ITIH2 was highly significantly associated with the presence of estrogen receptor. Our results may reflect differences between canine and human carcinomas, and further studies are needed to prove or refute these findings. Also, van den Broek and colleagues (2010) [
23] observed significantly increased serum concentrations in the breast cancer group for ITIH4, significantly decreasing after surgery. The authors suggested a potential for ITIH4 in post-surgical breast cancer monitoring [
23]. Accordingly, van Winden and colleagues (2010) [
23] and Opstal-van Winden and colleagues (2011) [
24] discovered that isoform 1 of the heavy chain of inter-alpha trypsin inhibitor H4 (ITIH4) presented elevated values in pre-diagnosis breast cancer and were already altered up to three years before cancer detection. Corroborating the previous studies, the study by Yang and colleagues (2016) [
25] aimed to identify new serum peptide biomarkers for women with breast cancer (BC) through the analysis of the serum proteomic profile. Three peptide biomarkers were identified, including ITIH4, and the authors concluded that this protein has diagnostic and prognostic potential for breast cancer.
Another high abundant protein found mostly in the malign carcinoma group was phosphopyruvate hydratase or alpha enolase (ENO1), that plays a crucial role in cancer development by promoting cell proliferation, invasion and metastasis. Its overexpression has been linked to chemoresistance and poor prognosis [
26]. Song and colleagues (2015) [
27] analyzed the involvement of ENO1 in tumor progression and prognosis of human glioma. They concluded that overexpression of ENO1 was associated with glioma progression, while reducing ENO1 expression led to suppression of cell growth, migration and invasion progression by inactivating the PI3K/Akt pathway in glioma cells. In cancer metabolic reprogramming, ENO1 stimulates cancer cells to create energy largely through the breakdown of glucose in a non-oxidative way, rather than typical oxidative phosphorylation, which is known as the Warburg effect [
28,
29]. Indeed, mammary canine tumors undergo metabolic alterations [
30] that favors tumor growth. Chu and colleagues (2011) [
31] examined ENO1 expression by immunohistochemical staining and evaluated its importance in canine mammary carcinoma. ENO1 overexpression significantly correlated with shorter survival but was not associated with ER positivity in canine mammary carcinoma [
31]. The study of Zamani-Ahmadmahmudi and colleagues (2014) [
32] used serological proteomic analysis to detect autoantigens that provoke a humoral response in dogs with mammary tumors. They found four autoantigens, including ENO1, with significantly greater immunoreactivity in tumor samples than in normal samples, which were identified as biomarker candidates [
32].
Besides ITIH2, ITIH4 and ENO1, we also found eukaryotic translation elongation factor 2 (eEF2) as abundant protein. It’s a member of the GTP-binding translation elongation factor family that is essential for protein synthesis. In our study, this protein was abundantly found in the SEVs proteomic of the malignant carcinoma group. Eukaryotic translation elongation factor 2 (eEF2) and its kinase eEF2K regulates many cellular processes, such as protein synthesis, cellular differentiation and malignant transformation [
33]. Scientific evidence shows that eEF2K regulates the cell cycle, autophagy, apoptosis, angiogenesis, invasion and metastasis in various types of cancer. Expression of eEF2K promotes cancer cell survival, and the level of this protein is increased in many cancer cells to adapt them to microenvironmental conditions, including hypoxia, nutrient depletion, and acidosis [
34]. The study by Wang and colleagues (2019) [
35] investigated the role of autophagy and its regulator, eukaryotic elongation factor 2 kinase (eEF2K), in determining the biological nature of triple-negative breast cancers (TNBCs). eEF2K and autophagy play key roles in maintaining aggressive tumor behavior and chemoresistance in TNBC, and eEF2K silencing may be a new strategy for the treatment of TNBC.
Actin B (ACTB) was also found as abundant and exclusive protein in this study. Beta-actin (ACTB) is a cytoskeletal structural protein widely distributed in all eukaryotic cells and plays critical roles in cell migration and cell division. ACTB is generally upregulated in most tumor cells and tissues, being associated with the invasiveness and metastasis of cancers [
36]. Fang et al. (2019) [
37] evaluated the role of serum actin-4 as a biomarker for the diagnosis of BC, as well as the association between its levels and clinicopathological characteristics. They reported that serum ACTN4 level was upregulated in BC patients compared to healthy controls [
37]. Furthermore, high ACTN4 expression was significantly associated with clinical stage, tumor grade, and lymph node status [
37,
38]. Numerous clinical studies showed that changes in ACTN4 gene expression are correlated with aggressiveness, invasion, and metastasis in certain tumors [
39]. Wang and colleagues (2017) [
40] sought to determine the role and regulation of ACTN4 expression in human breast cancer metastasis under ellagic acid (EA)-based therapy. EA inhibited breast cancer growth and metastasis by directly targeting ACTN4 in vitro and in vivo. ACTN4 knockdown resulted in blocking malignant cell proliferation, colony formation and improving metastasis potency [
40]. Increased ACTN4 expression was directly associated with advanced cancer stage, increased incidence of metastases and short survival period global [
40].
Chung and colleagues (2014) [
41] used mass spectrometry to identify proteins differentially expressed in sera from women with breast cancer and healthy volunteers. From 51 protein peaks that were significantly up- or down-regulated, they obtained five protein peaks that showed positive association with large tumor size and lymph node involvement and could distinguish women with breast cancer with high sensitivity and specificity. Interestingly, transthyretin and beta2-glycoprotein were between those proteins. In agreement, our study also identified an abundance and exclusive peak for both proteins. However, previously, the study by Nasim et al. (2012) [
42] identified several acute phase proteins as potential biomarkers for breast carcinoma but reported that transthyretin was significantly downregulated in the sera of BC patients [
42]. Recently, Sharma et al. (2023) [
43] studied the action of organochlorine pesticides (OCP) and their relationship with the spread of breast cancer in Indian women, through the proteomic analysis of plasma from breast cancer patients. They found 17 dysregulated proteins, but transthyretin (TTR) was three times higher than in healthy controls, and docking and molecular dynamics studies revealed a competitive affinity between the pesticide endosulfan II and the thyroxine-binding site of TTR that could result in endocrine dysregulation and lead to breast cancer. The authors pointed to a putative role of TTR in OCP-mediated CM [
43].
Apolipoprotein H (APOH) or beta2-glycoprotein I (β2-GPI) is a plasma glycoprotein that has been implicated in a variety of physiological functions. Apolipoproteins (APOs) bind to lipids to form lipoproteins. By functioning as lipid carriers, apolipoproteins act as ligands for cell membrane receptors, enzyme cofactors and structural components of lipoproteins [
44,
45]. Several studies demonstrate the interaction of APOs with classical tumorigenesis pathways. Furthermore, dysregulation of APOs may indicate the occurrence and progression of cancer, thus serving as potential biomarkers for cancer patients [
46]. Besides Chung and colleagues [
41], Lee and colleagues (2023) [
47] also investigated the role of β2-GPI in tumor cells from mammary cancer patients and its correlation with tumor prognosis. β2-GPI expression was predominantly observed in the cells from breast cancer patients and significantly correlated with tumor stage and lymph node metastasis of breast cancer. High expression of β2-GPI was significantly correlated with better overall survival (OS) and disease-free survival (DFS). Their results corroborated Lin and colleagues previously results [
48].
Cytoskeletal rearrangement occurs in several cellular processes and involves a broad spectrum of proteins such as gelsolin superfamily proteins [
49]. Gelsolin (GSN) is a multifunctional actin-binding protein that is greatly decreased in many transformed cell lines in tumor tissues, including breast cancers. The study by Mielnicki et al. (1999) [
50] demonstrated epigenetic modification leading to downregulation of gelsolin expression in human breast cancer. Baig et al. (2013) [
49] investigated the germline mutations and expressional profile of gelsolin in human breast cancer tissues. Different types of mutations were observed in gelsolin coding regions and the transcript level was significantly lower in breast tumor tissues compared to control samples as well as in metastatic patients compared to disease-free patients at final follow-up. Downregulation of gelsolin was also described by Winston and colleagues (2001) [
51] in invasive mammary carcinomas. Stock et al. (2015) [
52] analyzed gelsolin mRNA levels and concluded that high gelsolin levels are associated with a better prognosis in ER(+) HER2(-) breast cancer and a reduction in tumor cell migration. However, gelsolin was also found as exclusive protein in this study, instead of downregulated. This finding was intriguing but considering the differences between canine and human species and the scarcity of studies in dogs, we hypothesize that the abundance of gelsolin protein in samples from dogs with malignant mammary tumors may be related to a better prognosis in these individuals, and not related to its function as a diagnostic biomarker. Therefore, it is likely that gelsolin may play a dual role in cancer genesis. For example, Rao and colleagues (2002) [
53] reported decreased gelsolin expression in the early stages of malignant transformation in urothelial carcinomas. However, there was an increase in gelsolin expression in the transition from non-invasive to invasive urothelial carcinomas, indicating a biphasic expression pattern. Therefore, the authors concluded that increased gelsolin expression could indicate the conversion of a superficial tumor into an invasive tumor. Corroborating our finding, Van den Abbeele and colleagues (2007) [
54] demonstrated that downregulation of gelsolin in several cancer cell types significantly reduces the invasive and motile properties of the cells, as well as cell aggregation, pointing to a role for gelsolin as tumor activator. Also, Chen et al. (2015) [
53] investigated the role of gelsolin (GSN) in epithelial-mesenchymal transition (EMT) in breast cancer cells in response to transforming growth factor beta 1 (TGF-β1). The results showed that TGF-β1 induced the expression of GSN and EMT and increased cell migration and invasion [
55]. GSN overexpression affected cell proliferation, cell cycle and migration, and modulated vimentin expression [
55]. The authors concluded that TGF-β1 promoted the demethylation of the GSN gene promoter in these cells, inducing epigenetic modifications that contributed to the progression of EMT in breast cancer cells [
55]. In the study by Zhang et al. (2020) [
56] the role of gelsolin (GSN) in hepatocellular carcinoma (HCC) was investigated. GSN was found to be overexpressed in HCC tissues and correlated with advanced tumor grade and poor prognosis, while GSN knockdown inhibited tumor cell migration and invasion and GSN overexpression had the opposite effect. The authors suggested that GSN promotes HCC progression by regulating epithelial-mesenchymal transition (EMT).
Interestingly, Vitamin D-binding protein (VDBP) or group-specific component (GC) was identified as an abundant protein in the remission group. We didn’t find literature about the relationship between cancer and VDBP, but Tagliabue and colleagues (2015) [
57] performed a meta-analysis of vitamin D binding protein and cancer risk. They included 28 independent studies relating to the following tumors: basal cell carcinoma, bladder, breast, colon-rectum, endometrium, liver, esophagus, stomach, melanoma, pancreas, prostate and kidney. The authors concluded they could not assert that DBP would be a marker of cancer remission, but the results suggest a trend towards a decreased risk of cancer in individuals with higher levels of DBP, which could indicate a potential protective role of DBP in cancer.