Abstract
Background/Aim: Until now, little emphasis has been placed on the protein expression profile of male breast cancer (MBC) tumors, due to the rarity of the disease. The present study aimed to identify a proteomic pattern that is characteristic for malignant male breast tissue epithelium. Materials and Methods: The protein content of four male breast tumors and corresponding adjacent healthy (control) tissues was analyzed by high-throughput nano-liquid chromatography-MS/MS technology. Results: A total of 2,352 proteins were identified, that correspond to 1,249 single gene products, with diverse biological roles. Of those, a panel of 119 differentially expressed tissue proteins was identified in MBC samples compared to controls; 90 were found to be over-expressed in MBC tissues, while 29 were down-regulated. Concurrently, 844 proteins were detected only in MBC tumors and 197 were expressed exclusively in control mammary samples. Conclusion: Differential proteomic expression was found in MBC tissue, leading to improved understanding of MBC pathology and highlighting the need for personalized management of male patients.
Male breast cancer (MBC) is an uncommon disease, accounting for fewer than 1% of all diagnosed breast carcinoma cases (1). The estimates of the American Cancer Society, that predict merely 2,650 new cases of invasive MBC and 530 deaths from the disease in the United States for 2021, indicate its limited impact on male morbidity and mortality (2). Despite the fact that this entity is rare, the incidence of MBC is on the rise over the past two decades (3), attracting increasing interest in terms of deciphering its biology, etiopathogenesis, clinical presentation and unique management. Still, because of the low number of diagnoses, MBC remains understudied, prospective research-based progress is hindered, and decision-making for patients and physicians is based on observational retrospective studies (4, 5).
Owning to the scarcity of research in the field, clinical data on the management of MBC are usually extrapolated from trials on its female counterpart; this approach is becoming less and less acceptable due to the inherent discrepancies between the two entities. Characteristically, the largest retrospective collection and pathological review of MBC tumor samples so far, showed a lack of association between histological grading and outcome, which may be linked to the different distribution of disease subtypes in men compared with women (6). For example, men are much more likely to express the estrogen (ER) or androgen (AR) receptors than women, and less likely to over-express HER2 (6, 7). Lobular tumors are also much less frequent in men despite being relatively common in women, whereas most MBC cases are subtyped as either Luminal A–like or Luminal-B-like/HER-2-negative (6, 8). The genomic landscape of the disease further highlights the differences between sexes, with a recent publication reporting less frequent 16q losses, PIK3CA mutations, and TP53 mutations than those seen in ER-positive/HER2-negative female breast cancers, but more frequent somatic mutations in genes involved in DNA repair pathways (9). These findings do not only highlight the biologic heterogeneity between male and female patients, but also indicate the need for a better tailored approach for optimal male breast cancer classification and prognostication.
Current efforts are focusing on the identification of unique subgroups of breast malignancies in men, which differ from the intrinsic subtypes women harbor, mainly in the context of gene expression studies (10). Characteristically, gene expression profiling approaches have geared towards the stratification of MBC into two distinct subgroups that occur only among males and that differ from the well-established intrinsic subtypes of breast cancer in women; one group displays high expression of immune response genes and genes associated with ER signaling, whereas the other subgroup shows up-regulation of genes involved in activation of invasion/metastasis and induction of angiogenesis, such as HOXB7 (11). However, genes lack any catalytic functions and instead exercise their effects through proteins, the biomolecules that direct cellular processes and disease progression (12). Given the weak link between genomic alterations and protein functionality, due to biological processes such as RNA alternative splicing and post-translational modifications, using targeted proteomics in clinical settings allows us to gain a deeper molecular knowledge of breast cancer and shed light on biomarker research (13). Recent advances in proteomic technologies, such as High Performance Liquid Chromatography-tandem Mass Spectrometry (HPLC-MS/MS), enable us to combine high separation capacity and strong qualitative ability of proteins in biological samples (14).
In the current study, we employed HPLC-MS/MS to examine protein alterations in MBC tissue samples. HPLC-MS/MS is proposed as a powerful analytical tool for deciphering protein expression alterations as a function of breast cancer progression and, according to our knowledge, has not yet been utilized for the proteomic profiling of male breast cancer tumors. Our aim was to identify a panel of proteins that are differentially regulated in male breast tumors, by comparing protein profiles between tumor and adjacent non-tumor tissues. Unravelling the altered protein expression in malignant tissue could help researchers better understand the pathogenesis of male breast cancer and contribute towards the identification of more accurate markers that have the potential to be clinically meaningful.
Materials and Methods
Patients and tissue samples collection. Four male patients with sporadic breast cancer (MBC) attending the First Propaedeutic Surgical Department at Hippokration General Hospital in Athens, Greece, were enrolled in this study. Participation was voluntary and informed consent was provided by each participant. MBC was diagnosed based on a combination of standard clinical, radiological, and histological criteria (15). Breast tissue samples were collected from each patient at the time of surgery and stored at –80°C until use. The control group consisted of an equal number of adjacent healthy tissue samples, which were extracted from each patient during surgical treatment and preserved at the same storage conditions as the malignant tumors. The study was approved by the Ethical Committee of the participating hospital. The main demographic and clinical characteristics of the MBC patients are presented in Table I.
Clinical characteristics of the four male breast cancer patients.
Tissue sample preparation and peptide generation. The fresh frozen tissue samples were treated with lysis buffer, as previously described (16). Following lysis, tissue disruption was carried out, first by mechanical pressure and then by tip-sonication for three cycles of 18s duration each, under an amplification of 38%. The subsequent step involved centrifugation for 20 min at 13,000 × g, in order to discard the pellet which contained organelle debris and other insoluble materials. The Bradford protein assay was applied for total protein quantification of the supernatant. Then, 200 μg of protein from each tissue sample was treated with 7 M urea buffer and 80 mM triethylamonium bicarbonate. Next, protein alkylation was performed applying 55 mM iodoacetamide solution for 30 min at room temperature in dark conditions. Finally, a trypsin solution (Roche Diagnostics, Hoffman-La-Roche, Basel, Switzerland) at a final concentration of 500 ng/μl was added to all samples, and after thorough mixing, the digestion took place overnight, at 36°C, in a humidified environment. The following day, 5% (vol/vol) acetic acid was added to end trypsinization, and the peptide-containing mixtures were lyophilized in a vacuum concentrator for approximately 1 h, until they were completely dried. The resulting powder was dissolved in 0.1% formic acid in distilled water and LC-MS/MS analysis followed.
One-dimensional nano liquid chromatography (1D-nanoLC) separation and MS/MS analysis. MS analysis was performed as previously described (17). The peptides that derived from the previous steps were analyzed using an LTQ Orbitrap Elite coupled with a Dionex Ultimate-3000 nanoHPLC system (Thermo Scientific, Rockford, IL, USA). Specifically, one-dimensional liquid chromatography (1D-LC) separation took place after loading our samples subsequently onto two C-18 pre-columns (Thermo Scientific, PepMap® RSLC, 100 Å, 3-μm-bead-packed 15-cm column and 2-μm-bead-packed 50-cm column). The mobile phases A and B were 0.1% formic acid and 99% acetonitrile in water, respectively, at a flow rate of 0.3 μl/min. The gradient elution profile for all runs was 360 min for 2%-30% of buffer B (99.9% acetonitrile, 0.1% formic acid) at a constant 300 nl/min flow rate. Notably, to achieve maximum peptide recovery, we ended up testing several linear and non-linear gradient profiles, aiming at the optimal gradient, according to the material under study.
In the following step, an Orbitrap Elite mass spectrometer fitted with a nano-spray source was employed to collect mass spectra; data were acquired in accordance to a data-dependent acquisition strategy. Full-scan data were collected with a resolution of 60,000 and a maximum injection time of 120 ms. Tandem mass spectrometry for the 20 most intense ions per survey scan was performed, and peptide fragmentation was achieved at a high-energy collision dissociation fragmentation mode, set at a normalized collision energy of 36 NSE%. The MS/MS spectra of the resulting fragments were acquired with a resolving power of 15,000 and a maximum integration time of 120 ms. Measurements were performed using m/z 445.120025 as lock mass. Notably, to prevent repetitive selection of the same peptide, dynamic exclusion was employed within 45 s.
Data analysis. The raw data files that arose were processed with the Proteome Discoverer software (Thermo Scientific), using the SequestHT (version 28.0, Thermo Scientific) search engine applying the Homo sapiens UniProtKB/Swiss-Prot database. MS/MS searches were performed according to specific parameters; these included a 10 ppm parent peptide mass tolerance, a 0.05 fragment ion tolerance, up to two missed cleavage points for trypsin, cysteine methylation as a fixed modification, and oxidation of methionine as a variable modification. Peptide spectral matches (PSM) were considered valid at a 1% false discovery rate (FDR) (q-value <0.01). Further peptide filtering was done based on Xcorr versus peptide charge values (percolator maximum Delta Cn was 0.05). Of note, the minimum length of acceptable identified peptides was set as six amino acids.
Protein annotation and network analysis. Information on the identified proteins’ annotation were retrieved via the Proteome Discoverer software. Gene Ontology (GO) terms regarding biological process and pathway was used for the construction of protein classifications, while they were retrieved through PANTHER (Protein Analysis Through Evolutionary Relationships) version 16, which is an integrated knowledgebase of evolutionary and functional relationships between protein-coding genes (18). The potential interactions across previously detected proteins were studied with the STRING webserver.
Results
Focusing on deciphering the proteome of male breast cancer tumors at a high-resolution scale, four male breast cancer tumor and four adjacent normal tissue samples from the same patients were analyzed using untargeted tandem mass spectrometry. A total of 2,352 proteins were identified in the analyzed tissue samples, that correspond to 1,249 single gene products with diverse biological roles (Supplementary Tables I-IV). Of those, a panel of 119 differentially expressed tissue proteins was identified in MBC samples compared to controls; 90 were found to be over-expressed in MBC tissues, while 29 were down-regulated in male breast tumors. Concurrently, 844 proteins were detected exclusively in the MBC group and 197 were only expressed in the control samples.
The pathways which engage the differentially expressed proteins identified in our samples were studied using the STRING database (Figure 1, Figure 2, Figure 3 and Figure 4). These figures depict a non-directed graph with the functional interactions in the dataset presented as predicted by STRING. Each node represents a protein, and an edge between protein X and Y represents predicted functional interaction among these proteins. This interplay can be viewed as one protein being regulated by another, a genetic or physical interaction, the participation in the same biological process, or the common regulation of both proteins by a third molecule, among others. Thicker network lines demonstrate stronger protein relation as well as neighboring positions.
The figure depicts the interactions of the 90 proteins identified as being over-expressed in male breast cancer (MBC) samples. Each node represents a protein, Edges represent functional interaction between the proteins. Thicker network lines demonstrate stronger protein relation as well as neighboring positions. Protein abbreviations used are presented in Supplementary Table I.
STRING network of the 29 proteins that were down-regulated in male breast tumors. Network nodes represent proteins and edges represent protein-protein associations.
A total of 844 proteins were detected only in the male breast cancer (MBC) group. Due to the large number of identified molecules, the STRING network of these proteins was set at a minimum required interaction score of highest confidence (0.900), with no more than five interactors for each protein shown, while the disconnected nodes in the network were hidden.
The figure depicts the interactions of the 197 proteins identified as being exclusively expressed in our control samples. Network nodes represent proteins. Thicker network lines demonstrate stronger protein relation as well as neighboring positions.
For this interpretation, a multi-level comprehensive analysis based on GO categorization was performed using the PANTHER classification system, which sorts the proteins into respective groups based on their biological process, and their involvement in different signaling pathways. Concerning the 90 over-expressed proteins in MBC tumors, as presented in Figure 5, it is of great interest that the four major groups involved molecules related to cellular process (52.9%), metabolic process proteins (28.7%), biological regulation (17.2%), and members of the response to stimulus process (14.9%), factors that are known to be aberrant in cancer. Regarding the analysis of the 29 down-regulated proteins, as illustrated in Figure 5, the majority of them are involved in cellular process (37.9%), biological regulation (31.0%), metabolic process (31.0%), and response to stimulus (27.6%). Additionally, as shown in Figure 6, most of the 844 proteins identified as being exclusively expressed in MBC tumors participate in cellular processes (55.5%), metabolic processes (34.7%), and biological regulation (22.2%). Finally, the correlation of the 197 proteins identified exclusively in the control samples to biological processes, according to the GO categorization, revealed that half of the proteins identified were involved in cellular process (53.6%), followed by a significant number of proteins taking part in the metabolic process (29.3%) and the biological regulation (28.7%) (Figure 6).
Biological process of differentially expressed proteins, either over-expressed in male breast cancer (MBC) tissues (blue color) or down-regulated in MBC tissues (orange color).
Biological process of proteins exclusively expressed in male breast cancer (MBC) (blue color) or CONTROL (orange color).
Regarding the PANTHER pathway analysis, as illustrated in Figure 7 and Figure 8, several predominant pathways were identified in each group of proteins. Specifically, 10.3% of the over-expressed and 17.2% of the down-regulated proteins in MBC corresponded to proteins found in the integrin signaling pathway, demonstrating the importance of the crosstalk between cancer cells and their microenvironment that triggers a variety of critical signaling cues (Figure 7). With regards to the proteins found exclusively in MBC samples, 2.6% were involved in the Huntington disease pathway, 2.2% in the Wnt signaling pathway, 2.2% in the inflammation mediated by chemokine and cytokine signaling pathway, and 1.9% in the integrin signaling pathway (Figure 8). Finally, a significant number of proteins found exclusively in the control samples were implicated in the integrin signaling pathway (3.3%), in cytoskeletal regulation by Rho GTPase (2.8%), in inflammation mediated by chemokine and cytokine signaling pathway (2.8%), as well as in the heterotrimeric G-protein signaling pathway-Gi alpha and Gs alpha mediated pathway (2.8%) (Figure 8).
Pathway analysis of differentially expressed proteins in male breast cancer (MBC).
Pathway analysis of proteins exclusively expressed in male breast cancer (MBC) tissues (blue color) or CONTROL tissue samples (orange color). This figure depicts the signal transduction pathways in which more than three of the identified proteins are involved in.
Based on the interactions of differentially expressed proteins in cancerous and non-cancerous tissues, distinct protein groups with different biological processes are observed. The potential interactions across upregulated proteins in MBC, revealed several proteins in the center of the network that demonstrated stronger relations (Figure 1). These proteins are mainly involved in glycolysis (ENO1, TPI1, PKM), heat stress response (ANXA2) and chaperonin-mediated protein folding (HSPA5, HSPA8, HSPD1, VCP, PDIA3, PDIA6), protein biosynthesis (EEF1A1), cytoskeleton organization (PLEC), DNA synthesis (PHB), cell division (TUBB, TUBA1C) and ATP synthesis (ATP5F1A, ATP5F1B). Concerning the interactions among proteins down-regulated in MBC samples, the densest area presented in Figure 2, involved molecules associated with the integrin- and growth factor-signaling pathway (ILK), cell adhesion (VCL, MYH10), myelination (MBP), cytoskeleton organization (SPTAN1, SPTBN1, FLNA), and epithelial cell differentiation (TAGLN). Additionally, the network of the interactions among the proteins exclusively identified in MBC (Figure 3) revealed distinct subgroups involved in protein biosynthesis and transport, proteasome assembly, mRNA processing, cell adhesion, endocytosis, cell division, and mitochondrial electron transport chain. Finally, the functional relationships analysis of the proteins expressed only in control samples, resulted in a network with distinct areas encompassing proteins involved in cell adhesion (LAMC1, LAMA4, LAMA5, WVF, NID1, ICAM1, MCAM, CD99), in glycolysis/glycogen metabolism/glyconeogenesis (PYGL, PYGB, ALDOC, HK2, BPGM, GPD1), and in translation/protein biosynthesis/proteasome assembly (RPL19, RPL29, EIF3H, ADRM1, PSMA6) (Figure 4).
Discussion
In the current study, we attempted to characterize the tissue proteomic profile of patients with MBC using high-resolution mass spectrometric technology. The analysis identified a total of 2,352 proteins, belonging to 1,249 single gene products, that are involved in a multitude of biological processes. Several proteins (119 molecules) were differentially expressed in patients with MBC compared to controls. Of those, 90 were found to be upregulated in all tissue samples of the MBC group, while 29 were down-regulated in male breast tumors. Additionally, 844 proteins were detected only in samples from MBC tumors and 197 were exclusively expressed in the control samples. The comparison revealed outstanding differences between MBC samples and adjacent “healthy” tissue from the same patient, giving us insight into the perturbations of pathways within tumor cells and establishing an important foundation for future male breast cancer-related discoveries using differential proteomics profiling.
MS technologies have enabled increased proteome coverage rather than using the traditional approach of translating gene expression data. Especially in the context of clinical proteomics, cancerous tissue proteomic analysis provides the most accurate reflection of the tumor’s physiological state (19-22). However, most protein MS-based breast cancer screening strategies refer to female patients and although possible biomarkers have been identified, few have been validated so as to achieve clinical application (23-26). Regarding protein-expression alterations in male breast tumors, the only published study using MALDI-TOF mass spectrometric technology revealed over-expression of structural proteins, glycolytic enzymes, stress related proteins, enzymes involved in the synthesis of AMP, heat-shock proteins, and others (27). However, further advancements in the HPLC-MS/MS platform, which exhibits increased sensitivity and reliability in identifying proteins in complex mixtures, have not yet been utilized for the proteomic profiling of cancerous male breast regions that are obtained during surgery.
Breast cancer has traditionally been considered a female-specific disease, presenting approximately one hundred times more often in women compared to men (28). To date, no proteomic approach has identified the molecular mechanisms involved in gender-specific tumorigenesis – much is still to be learned. The hypothesis is that by looking at gender-related proteins, physiological differences between men and women could be revealed. Characteristically, none of the differentially expressed proteins that we identified are encoded by genes on the Y chromosome. On the other hand, a small number of over-expressed proteins in MBC are products of genes on the single X chromosome in males (Table II). Of these, the cytoskeletal protein MOES has been suggested as a potential marker for poor prognosis in females with ER-positive breast cancer treated with tamoxifen (29). HCD2 is a lipid metabolism- mitochondrial-related protein that has been found to be selectively upregulated in mammosphere-forming cells from breast tumor cell lines that represent the population of cancer stem cells (CSCs) that might lead to the failure of routine chemotherapy and the recurrence of the disease (30). NONO is a DNA-binding protein found in a number of breast tumor cell lines, and has been proposed as a significant independent prognostic factor for patients since its knockdown substantially represses the proliferation of breast cancer cells, acting through the pre-mRNA splicing of cell proliferation-related genes (31). Interestingly, NONO exhibits gender expression differences since higher expression levels in females have been reported (32). On the other hand, FLNA is an X-linked protein that we found to be down-regulated in MBC tissues. Although most of the studies point to a positive correlation between expression of FLNA and breast cancer development and progression (33), some studies are in line with our results, showing that FLNA expression is negatively correlated with the metastatic potential of breast cancer lines (34, 35). This is probably because FLNA is able to regulate focal adhesion disassembly and suppress breast cancer cell migration and invasion (35). These gender differences could have important clinical implications for treatment decision making and might eventually enable us to make better treatment choices, individualized for each patient.
Differentially expressed proteins in male breast cancer (MBC) encoded by the X chromosome genes.
Concerning the limitations of our study that deserve mention, it is important to advert that we analyzed a small number of tumor samples due to the rarity of the disease, deeming our results merely indicative and in need of further confirmation in a larger cohort of patients. Despite the small number of samples analyzed and thus our inability to draw safe conclusions on the detailed proteomic pattern of MBC tumors, we obtained high-throughput reproducible results that literature is currently lacking.
In conclusion, in this work we detected the differentially expressed protein in MBC tumors, using high resolution MS. The data generated highlight the important contribution of MS-based techniques towards deciphering the molecular basis of male breast malignancies. While some strides have been made to address the need for male-specific information, the challenges still exist and need to be tackled from diverse perspectives and on a larger scale. The in-depth proteomic profiling of this rare population might ensure that MBC patients receive the most appropriate treatment in respect to the gender specific differences in cancer biology and management.
Footnotes
Authors’ Contributions
GTT conceptualized the project and the methodology. EZ, with the support of AMK, analyzed the data, generated the tables and figures, and prepared the original draft under the supervision of FZ, MG and GTT. SCP and AKA performed the experimental procedures. All authors provided critical feedback, contributed to the manuscript, and approved the final version in accordance with criteria established by the International Committee of Medical Journal Editors (ICMJE).
↵Supplementary Material
The Supplementary Tables I to IV can be accessed online at: https://drive.google.com/drive/folders/1EErNPuxByG0SeQPV6g6g_LG-ghfzJj2Y?usp=sharing
Conflicts of Interest
The Authors declare no conflicts of interest regarding this study.
- Received November 12, 2021.
- Revision received December 29, 2021.
- Accepted December 31, 2021.
- Copyright © 2022, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved