A serum glycomics approach to breast cancer biomarkers

Because the glycosylation of proteins is known to change in tumor cells during the development of breast cancer, a glycomics approach is used here to find relevant biomarkers of breast cancer. These glycosylation changes are known to correlate with increasing tumor burden and poor prognosis. Current antibody-based immunochemical tests for cancer biomarkers of ovarian (CA125), breast (CA27.29 or CA15-3), pancreatic, gastric, colonic, and carcinoma (CA19-9) target highly glycosylated mucin proteins. However, these tests lack the specificity and sensitivity for use in early detection. This glycomics approach to find glycan biomarkers of breast cancer involves chemically cleaving oligosaccharides (glycans) from glycosylated proteins that are shed or secreted by breast cancer tumor cell lines. The resulting free glycan species are analyzed by MALDI-FT-ICR MS. Further structural analysis of the glycans can be performed in FTMS through the use of tandem mass spectrometry with infrared multiphoton dissociation. Glycan profiles were generated for each cell line and compared. These methods were then used to analyze sera obtained from a mouse model of breast cancer and a small number of serum samples obtained from human patients diagnosed with breast cancer or patients with no known history of breast cancer. In addition to the glycosylation changes detected in mice as mouse mammary tumors developed, glycosylation profiles were found to be sufficiently different to distinguish patients with cancer from those without. Although the small number of patient samples analyzed so far is inadequate to make any legitimate claims at this time, these promising but very preliminary results suggest that glycan profiles may contain distinct glycan biomarkers that may correspond to glycan "signatures of cancer."     

Identification of potential immunotherapy biomarkers for breast cancer by bioinformatics analysis

Breast cancer is a serious malignancy with a high incidence worldwide and a tendency to relapse. We used integrated bioinformatics analysis to identify potential biomarkers in breast carcinoma in the present study. Microarray data, 127breast tumor samples and 23 non-tumor samples, received from the Gene Expression Omnibus (GEO) dataset; 121 differentially expressed genes (DEGs) were selected. Functional analysis using DAVID revealed that these DEGs were highly gathered in endodermal cell differentiation and proteinaceous extracellular matrix. Five bioactive compounds (prostaglandin J2, tanespimycin, semustine, 5182598, and flunarizine) were identified using Connectivity Map. We used Cytoscape software and STRING dataset to structure a protein–protein interaction (PPI) network. The expression of CD24, MMP1, SDC1, and SPP1 was much higher in breast carcinoma tissue than in Para cancerous tissues analyzed by Gene Expression Profiling Interactive Analysis (GEPIA) and ONCOMINE. Overexpression ofCD24, MMP1, SDC1, and SPP1 indicated the poor prognosis in breast carcinoma patients analyzed by Kaplan–Meier (KM) Plotter. Immunohistochemistry microarray was used to further confirm that protein expression of CD24, MMP1, SDC1, and SPP1 was much higher in tumor sections than in Para cancerous tissues. Hub genes expression at the protein level was correlated tothe breast cancer subtype and grade. Furthermore, immunity analysis showed that CD24, MMP1, SDC1, and SPP1 were potentially associated with five immune cell types infiltration (CD8+ T cells, CD4+ T cells, neutrophils, macrophages,and dendritic cells) by TIMER. Thus, this study indicates potential biomarkers that could have applications in the development of immune therapy for breast cancer. However, further studies are required for verifying these results in vivo and vitro.     

Concurrent Gene Signatures for Han Chinese Breast Cancers

The interplay between copy number variation (CNV) and differential gene expression may be able to shed light on molecular process underlying breast cancer and lead to the discovery of cancer-related genes. In the current study, genes concurrently identified in array comparative genomic hybridization (CGH) and gene expression microarrays were used to derive gene signatures for Han Chinese breast cancers.We performed 23 array CGHs and 81 gene expression microarrays in breast cancer samples from Taiwanese women. Genes with coherent patterns of both CNV and differential gene expression were identified from the 21 samples assayed using both platforms. We used these genes to derive signatures associated with clinical ER and HER2 status and disease-free survival.Distributions of signature genes were strongly associated with chromosomal location: chromosome 16 for ER and 17 for HER2. A breast cancer risk predictive model was built based on the first supervised principal component from 16 genes (RCAN3, MCOLN2, DENND2D, RWDD3, ZMYM6, CAPZA1, GPR18, WARS2, TRIM45, SCRN1, CSNK1E, HBXIP, CSDE1, MRPL20, IKZF1, and COL20A1), and distinct survival patterns were observed between the high- and low-risk groups from the combined dataset of 408 microarrays. The risk score was significantly higher in breast cancer patients with recurrence, metastasis, or mortality than in relapse-free individuals (0.241 versus 0, P<0.001). The concurrent gene risk predictive model remained discriminative across distinct clinical ER and HER2 statuses in subgroup analysis. Prognostic comparisons with published gene expression signatures showed a better discerning ability of concurrent genes, many of which were rarely identifiable if expression data were pre-selected by phenotype correlations or variability of individual genes.We conclude that parallel analysis of CGH and microarray data, in conjunction with known gene expression patterns, can be used to identify biomarkers with prognostic values in breast cancer.     

Identification of novel high-frequency DNA methylation changes in breast cancer

Recent data have revealed that epigenetic alterations, including DNA methylation and chromatin structure changes, are among the earliest molecular abnormalities to occur during tumorigenesis. The inherent thermodynamic stability of cytosine methylation and the apparent high specificity of the alterations for disease may accelerate the development of powerful molecular diagnostics for cancer. We report a genome-wide analysis of DNA methylation alterations in breast cancer. The approach efficiently identified a large collection of novel differentially DNA methylated loci (approximately 200), a subset of which was independently validated across a panel of over 230 clinical samples. The differential cytosine methylation events were independent of patient age, tumor stage, estrogen receptor status or family history of breast cancer. The power of the global approach for discovery is underscored by the identification of a single differentially methylated locus, associated with the GHSR gene, capable of distinguishing infiltrating ductal breast carcinoma from normal and benign breast tissues with a sensitivity and specificity of 90% and 96%, respectively. Notably, the frequency of these molecular abnormalities in breast tumors substantially exceeds the frequency of any other single genetic or epigenetic change reported to date. The discovery of over 50 novel DNA methylation-based biomarkers of breast cancer may provide new routes for development of DNA methylation-based diagnostics and prognostics, as well as reveal epigenetically regulated mechanism involved in breast tumorigenesis.     

Explicit formulation of titration models for isothermal titration calorimetry

Isothermal titration calorimetry (ITC) produces a differential heat signal with respect to the total titrant concentration. This feature gives ITC excellent sensitivity for studying the thermodynamics of complex biomolecular interactions in solution. Currently, numerical methods for data fitting are based primarily on indirect approaches rooted in the usual practice of formulating biochemical models in terms of integrated variables. Here, a direct approach is presented wherein ITC models are formulated and solved as numerical initial value problems for data fitting and simulation purposes. To do so, the ITC signal is cast explicitly as a first-order ordinary differential equation (ODE) with total titrant concentration as independent variable and the concentration of a bound or free ligand species as dependent variable. This approach was applied to four ligand-receptor binding and homotropic dissociation models. Qualitative analysis of the explicit ODEs offers insights into the behavior of the models that would be inaccessible to indirect methods of analysis. Numerical ODEs are also highly compatible with regression analysis. Since solutions to numerical initial value problems are straightforward to implement on common computing platforms in the biochemical laboratory, this method is expected to facilitate the development of ITC models tailored to any experimental system of interest.     

Common plasma protein marker LCAT in aggressive human breast cancer and canine mammary tumor

Breast cancer is one of the most frequently diagnosed cancers. Although biomarkers are continuously being discovered, few specific markers, rather than classification markers, representing the aggressiveness and invasiveness of breast cancer are known. In this study, we used samples from canine mammary tumors in a comparative approach. We subjected 36 fractions of both canine normal and mammary tumor plasmas to high-performance quantitative proteomics analysis. Among the identified proteins, LCAT was selectively expressed in mixed tumor samples. With further MRM and Western blot validation, we discovered that the LCAT protein is an indicator of aggressive mammary tumors, an advanced stage of cancer, possibly highly metastatic. Interestingly, we also found that LCAT is overexpressed in high-grade and lymphnode-positive breast cancer in silico data. We also demonstrated that LCAT is highly expressed in the sera of advanced-stage human breast cancers within the same classification. In conclusion, we identified a possible common plasma protein biomarker, LCAT, that is highly expressed in aggressive human breast cancer and canine mammary tumor.     

Integrative “Omic” Analysis for Tamoxifen Sensitivity through Cell Based Models

It has long been observed that tamoxifen sensitivity varies among breast cancer patients. Further, ethnic differences of tamoxifen therapy between Caucasian and African American have also been reported. Since most studies have been focused on Caucasian people, we sought to comprehensively evaluate genetic variants related to tamoxifen therapy in African-derived samples. An integrative “omic” approach developed by our group was used to investigate relationships among endoxifen (an active metabolite of tamoxifen) sensitivity, SNP genotype, mRNA and microRNA expressions in 58 HapMap YRI lymphoblastoid cell lines. We identified 50 SNPs that associate with cellular sensitivity to endoxifen through their effects on 34 genes and 30 microRNA expression. Some of these findings are shared in both Caucasian and African samples, while others are unique in the African samples. Among gene/microRNA that were identified in both ethnic groups, the expression of TRAF1 is also correlated with tamoxifen sensitivity in a collection of 44 breast cancer cell lines. Further, knock-down TRAF1 and over-expression of hsa-let-7i confirmed the roles of hsa-let-7i and TRAF1 in increasing tamoxifen sensitivity in the ZR-75-1 breast cancer cell line. Our integrative omic analysis facilitated the discovery of pharmacogenomic biomarkers that potentially affect tamoxifen sensitivity.     

Differential expression analysis in ovarian cancer: A functional genomics and systems biology approach

BackgroundOvarian cancer is one of the rarest lethal oncologic diseases that have hardly any specific biomarkers. The availability of high-throughput genomic data and advancement in bioinformatics tools allow us to predict gene biomarkers and apply systems biology approaches to get better diagnosis, and prognosis of the disease with a tentative drug that may be repurposed.ObjectiveTo perform genome-wide association studies using microarray gene expression of ovarian cancer and identify gene biomarkers, construction and analyze networks, perform survival analysis, and drug interaction studies for better diagnosis, prognosis, and treatment of ovarian cancer.MethodThe gene expression profiles of both healthy and serous ovarian cancer epithelial samples were considered. We applied a series of bioinformatics methods and tools, including fold-change statistics for differential expression analysis, DisGeNET and NCBI-Gene databases for gene-disease association mapping, DAVID 6.8 for GO enrichment analysis, GeneMANIA for network construction, Cytoscape 3.8 with its plugins for network visualization, analysis, and module detection, the UALCAN for patient survival analysis, and PubChem, DrugBank and DGIdb for gene-drug interaction.ResultsWe identified 8 seed genes that were subjected for drug-gene interaction studies. Because of over-expression in all the four stages of ovarian cancer, we discern that genes HMGA1 and PSAT1 are potential therapeutic biomarkers for its diagnosis at an early stage (stage I). Our analysis suggests that there are 11 drugs common in the seed genes. However, hypermethylated seed genes HMGA1 and PSAT1 showcased a good interaction affinity with drugs cisplatin, cyclosporin, bisphenol A, progesterone, and sunitinib, and are crucial in the proliferation of ovarian cancer.ConclusionOur study reveals that HMGA1 and PSAT1 can be deployed for initial screening of ovarian cancer and drugs cisplatin, bisphenol A, cyclosporin, progesterone, and sunitinib are effective in curbing the epigenetic alteration.     

Discovery of pathway biomarkers from coupled proteomics and systems biology methods

Background

Breast cancer is worldwide the second most common type of cancer after lung cancer. Plasma proteome profiling may have a higher chance to identify protein changes between plasma samples such as normal and breast cancer tissues. Breast cancer cell lines have long been used by researches as model system for identifying protein biomarkers. A comparison of the set of proteins which change in plasma with previously published findings from proteomic analysis of human breast cancer cell lines may identify with a higher confidence a subset of candidate protein biomarker.

Results

In this study, we analyzed a liquid chromatography (LC) coupled tandem mass spectrometry (MS/MS) proteomics dataset from plasma samples of 40 healthy women and 40 women diagnosed with breast cancer. Using a two-sample t-statistics and permutation procedure, we identified 254 statistically significant, differentially expressed proteins, among which 208 are over-expressed and 46 are under-expressed in breast cancer plasma. We validated this result against previously published proteomic results of human breast cancer cell lines and signaling pathways to derive 25 candidate protein biomarkers in a panel. Using the pathway analysis, we observed that the 25 “activated” plasma proteins were present in several cancer pathways, including ‘Complement and coagulation cascades’, ‘Regulation of actin cytoskeleton’, and ‘Focal adhesion’, and match well with previously reported studies. Additional gene ontology analysis of the 25 proteins also showed that cellular metabolic process and response to external stimulus (especially proteolysis and acute inflammatory response) were enriched functional annotations of the proteins identified in the breast cancer plasma samples. By cross-validation using two additional proteomics studies, we obtained 86% and 83% similarities in pathway-protein matrix between the first study and the two testing studies, which is much better than the similarity we measured with proteins.

Conclusions

We presented a ‘systems biology’ method to identify, characterize, analyze and validate panel biomarkers in breast cancer proteomics data, which includes 1) t statistics and permutation process, 2) network, pathway and function annotation analysis, and 3) cross-validation of multiple studies. Our results showed that the systems biology approach is essential to the understanding molecular mechanisms of panel protein biomarkers.

     

Comprehensive analysis of gene expression profiles to identify differential prognostic factors of primary and metastatic breast cancer

Breast cancer accounts for nearly half of all cancer-related deaths in women worldwide. However, the molecular mechanisms that lead to tumour development and progression remain poorly understood and there is a need to identify candidate genes associated with primary and metastatic breast cancer progression and prognosis. In this study, candidate genes associated with prognosis of primary and metastatic breast cancer were explored through a novel bioinformatics approach. Primary and metastatic breast cancer tissues and adjacent normal breast tissues were evaluated to identify biomarkers characteristic of primary and metastatic breast cancer. The Cancer Genome Atlas-breast invasive carcinoma (TCGA-BRCA) dataset (ID: HS-01619) was downloaded using the mRNASeq platform. Genevestigator 8.3.2 was used to analyse TCGA-BRCA gene expression profiles between the sample groups and identify the differentially-expressed genes (DEGs) in each group. For each group, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were used to determine the function of DEGs. Networks of protein–protein interactions were constructed to identify the top hub genes with the highest degree of interaction. Additionally, the top hub genes were validated based on overall survival and immunohistochemistry using The Human Protein Atlas. Of the top 20 hub genes identified, four (KRT14, KIT, RAD51, and TTK) were considered as prognostic risk factors based on overall survival. KRT14 and KIT expression levels were upregulated while those of RAD51 and TTK were downregulated in patients with breast cancer. The four proposed candidate hub genes might aid in further understanding the molecular changes that distinguish primary breast tumours from metastatic tumours as well as help in developing novel therapeutics. Furthermore, they may serve as effective prognostic risk markers based on the strong correlation between their expression and patient overall survival.     

Immunoexpression Analysis and Prognostic Value of BLCAP in Breast Cancer

Bladder Cancer Associated Protein (BLCAP, formerly Bc10), was identified by our laboratory as being down-regulated in bladder cancer with progression. BLCAP is ubiquitously expressed in different tissues, and several studies have found differential expression of BLCAP in various cancer types, such as cervical and renal cancer, as well as human tongue carcinoma and osteosarcoma. Here we report the first study of the expression patterns of BLCAP in breast tissue. We analyzed by immunohistochemistry tissue sections of normal and malignant specimens collected from 123 clinical high-risk breast cancer patients within the Danish Center for Translational Breast Cancer Research (DCTB) prospective study dataset. The staining pattern, the distribution of the immunostaining, and its intensity were studied in detail. We observed weak immunoreactivity for BLCAP in mammary epithelial cells, almost exclusively localizing to the cytoplasm and found that levels of expression of BLCAP were generally higher in malignant cells as compared to normal cells. Quantitative IHC analysis of BLCAP expression in breast tissues confirmed this differential BLCAP expression in tumor cells, and we could establish, in a 62-patient sample matched cohort, that immunostaining intensity for BLCAP was increased in tumors relative to normal tissue, in more than 45% of the cases examined, indicating that BLCAP may be of value as a marker for breast cancer. We also analyzed BLCAP expression and prognostic value using a set of tissue microarrays comprising an independent cohort of 2,197 breast cancer patients for which we had follow-up clinical information.     

Characterization of MTAP Gene Expression in Breast Cancer Patients and Cell Lines

MTAP is a ubiquitously expressed gene important for adenine and methionine salvage. The gene is located at 9p21, a chromosome region often deleted in breast carcinomas, similar to CDKN2A, a recognized tumor suppressor gene. Several research groups have shown that MTAP acts as a tumor suppressor, and some therapeutic approaches were proposed based on a tumors´ MTAP status. We analyzed MTAP and CDKN2A gene (RT-qPCR) and protein (western-blotting) expression in seven breast cancer cell lines and evaluated their promoter methylation patterns to better characterize the contribution of these genes to breast cancer. Cytotoxicity assays with inhibitors of de novo adenine synthesis (5-FU, AZA and MTX) after MTAP gene knockdown showed an increased sensitivity, mainly to 5-FU. MTAP expression was also evaluated in two groups of samples from breast cancer patients, fresh tumors and paired normal breast tissue, and from formalin-fixed paraffin embedded (FFPE) core breast cancer samples diagnosed as Luminal-A tumors and triple negative breast tumors (TNBC). The difference of MTAP expression between fresh tumors and normal tissues was not statistically significant. However, MTAP expression was significantly higher in Luminal-A breast tumors than in TNBC, suggesting the lack of expression in more aggressive breast tumors and the possibility of using the new approaches based on MTAP status in TNBC.     

Novel comprehensive approach for accessible biomarker identification and absolute quantification from precious human tissues

The identification of specific biomarkers obtained directly from human pathological lesions remains a major challenge, because the amount of tissue available is often very limited. We have developed a novel, comprehensive, and efficient method permitting the identification and absolute quantification of potentially accessible proteins in such precious samples. This protein subclass comprises cell membrane associated and extracellular proteins, which are reachable by systemically deliverable substances and hence especially suitable for diagnosis and targeted therapy applications. To isolate such proteins, we exploited the ability of chemically modified biotin to label ex vivo accessible proteins and the fact that most of these proteins are glycosylated. This approach consists of three successive steps involving first the linkage of potentially accessible proteins to biotin molecules followed by their purification. The remaining proteins are then subjected to glycopeptide isolation. Finally, the analysis of the nonglycosylated peptides and their involvement in an in silico method increased the confident identification of glycoproteins. The value of the technique was demonstrated on human breast cancer tissue samples originating from 5 individuals. Altogether, the method delivered quantitative data on more than 400 potentially accessible proteins (per sample and replicate). In comparison to biotinylation or glycoprotein analysis alone, the sequential method significantly increased the number (≥30% and ≥50% respectively) of potentially therapeutically and diagnostically valuable proteins. The sequential method led to the identification of 93 differentially modulated proteins, among which several were not reported to be associated with the breast cancer. One of these novel potential biomarkers was CD276, a cell membrane-associated glycoprotein. The immunohistochemistry analysis showed that CD276 is significantly differentially expressed in a series of breast cancer lesions. Due to the fact that our technology is applicable to any type of tissue biopsy, it bears the ability to accelerate the discovery of new relevant biomarkers in a broad spectrum of pathologies.     

Low abundance protein enrichment for discovery of candidate plasma protein biomarkers for early detection of breast cancer

Molecular biomarkers of early stage breast cancer may improve the sensitivity and specificity of diagnosis. Plasma biomarkers have additional value in that they can be monitored with minimal invasiveness. Plasma biomarker discovery by genome-wide proteomic methods is impeded by the wide dynamic range of protein abundance and the heterogeneity of protein expression in healthy and disease populations which requires the analysis of a large number of samples. We addressed these issues through the development of a novel protocol that couples a combinatorial peptide ligand library protein enrichment strategy with isobaric label-based 2D LC-MS/MS for the identification of candidate biomarkers in high throughput. Plasma was collected from patients with stage I breast cancer or benign breast lesions. Low abundance proteins were enriched using a bead-based combinatorial library of hexapeptides. This resulted in the identification of 397 proteins, 22% of which are novel plasma proteins. Twenty-three differentially expressed plasma proteins were identified, demonstrating the effectiveness of the described protocol and defining a set of candidate biomarkers to be validated in independent samples. This work can be used as the basis for the design of properly powered investigations of plasma protein expression for biomarker discovery in larger cohorts of patients with complex disease.     

Integrated pipeline for mass spectrometry-based discovery and confirmation of biomarkers demonstrated in a mouse model of breast cancer

Despite their potential to impact diagnosis and treatment of cancer, few protein biomarkers are in clinical use. Biomarker discovery is plagued with difficulties ranging from technological (inability to globally interrogate proteomes) to biological (genetic and environmental differences among patients and their tumors). We urgently need paradigms for biomarker discovery. To minimize biological variation and facilitate testing of proteomic approaches, we employed a mouse model of breast cancer. Specifically, we performed LC-MS/MS of tumor and normal mammary tissue from a conditional HER2/Neu-driven mouse model of breast cancer, identifying 6758 peptides representing >700 proteins. We developed a novel statistical approach (SASPECT) for prioritizing proteins differentially represented in LC-MS/MS datasets and identified proteins over- or under-represented in tumors. Using a combination of antibody-based approaches and multiple reaction monitoring-mass spectrometry (MRM-MS), we confirmed the overproduction of multiple proteins at the tissue level, identified fibulin-2 as a plasma biomarker, and extensively characterized osteopontin as a plasma biomarker capable of early disease detection in the mouse. Our results show that a staged pipeline employing shotgun-based comparative proteomics for biomarker discovery and multiple reaction monitoring for confirmation of biomarker candidates is capable of finding novel tissue and plasma biomarkers in a mouse model of breast cancer. Furthermore, the approach can be extended to find biomarkers relevant to human disease.