1
|
Grabski IN, Vito RD, Trippa L, Parmigiani G. Bayesian combinatorial MultiStudy factor analysis. Ann Appl Stat 2023; 17:2212-2235. [PMID: 37786772 PMCID: PMC10543692 DOI: 10.1214/22-aoas1715] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2023]
Abstract
Mutations in the BRCA1 and BRCA2 genes are known to be highly associated with breast cancer. Identifying both shared and unique transcript expression patterns in blood samples from these groups can shed insight into if and how the disease mechanisms differ among individuals by mutation status, but this is challenging in the high-dimensional setting. A recent method, Bayesian Multi-Study Factor Analysis (BMSFA), identifies latent factors common to all studies (or equivalently, groups) and latent factors specific to individual studies. However, BMSFA does not allow for factors shared by more than one but less than all studies. This is critical in our context, as we may expect some but not all signals to be shared by BRCA1-and BRCA2-mutation carriers but not necessarily other high-risk groups. We extend BMSFA by introducing a new method, Tetris, for Bayesian combinatorial multi-study factor analysis, which identifies latent factors that any combination of studies or groups can share. We model the subsets of studies that share latent factors with an Indian Buffet Process, and offer a way to summarize uncertainty in the sharing patterns using credible balls. We test our method with an extensive range of simulations, and showcase its utility not only in dimension reduction but also in covariance estimation. When applied to transcript expression data from high-risk families grouped by mutation status, Tetris reveals the features and pathways characterizing each group and the sharing patterns among them. Finally, we further extend Tetris to discover groupings of samples when group labels are not provided, which can elucidate additional structure in these data.
Collapse
Affiliation(s)
- Isabella N Grabski
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
| | - Roberta De Vito
- Department of Biostatistics and Data Science Initiative, Brown University, Providence, RI
| | - Lorenzo Trippa
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA
| | - Giovanni Parmigiani
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA
| |
Collapse
|
2
|
Wiggins GAR, Black MA, Dunbier A, Morley-Bunker AE, Pearson JF, Walker LC. Increased gene expression variability in BRCA1-associated and basal-like breast tumours. Breast Cancer Res Treat 2021; 189:363-375. [PMID: 34287743 PMCID: PMC8357684 DOI: 10.1007/s10549-021-06328-y] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2021] [Accepted: 07/07/2021] [Indexed: 11/21/2022]
Abstract
Purpose Inherited variants in the cancer susceptibility genes, BRCA1 and BRCA2 account for up to 5% of breast cancers. Multiple gene expression studies have analysed gene expression patterns that maybe associated with BRCA12 pathogenic variant status; however, results from these studies lack consensus. These studies have focused on the differences in population means to identified genes associated with BRCA1/2-carriers with little consideration for gene expression variability, which is also under genetic control and is a feature of cellular function. Methods We measured differential gene expression variability in three of the largest familial breast cancer datasets and a 2116 breast cancer meta-cohort. Additionally, we used RNA in situ hybridisation to confirm expression variability of EN1 in an independent cohort of more than 500 breast tumours. Results BRCA1-associated breast tumours exhibited a 22.8% (95% CI 22.3–23.2) increase in transcriptome-wide gene expression variability compared to BRCAx tumours. Additionally, 40 genes were associated with BRCA1-related breast cancers that had ChIP-seq data suggestive of enriched EZH2 binding. Of these, two genes (EN1 and IGF2BP3) were significantly variable in both BRCA1-associated and basal-like breast tumours. RNA in situ analysis of EN1 supported a significant (p = 6.3 × 10−04) increase in expression variability in BRCA1-associated breast tumours. Conclusion Our novel results describe a state of increased gene expression variability in BRCA1-related and basal-like breast tumours. Furthermore, genes with increased variability may be driven by changes in DNA occupancy of epigenetic effectors. The variation in gene expression is replicable and led to the identification of novel associations between genes and disease phenotypes. Supplementary Information The online version contains supplementary material available at 10.1007/s10549-021-06328-y.
Collapse
Affiliation(s)
- George A R Wiggins
- Department of Pathology and Biomedical Science, University of Otago, Christchurch, New Zealand
| | - Michael A Black
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - Anita Dunbier
- Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - Arthur E Morley-Bunker
- Department of Pathology and Biomedical Science, University of Otago, Christchurch, New Zealand
| | | | - John F Pearson
- Department of Pathology and Biomedical Science, University of Otago, Christchurch, New Zealand.,Biostatistics and Computational Biology Unit, University of Otago, Christchurch, New Zealand
| | - Logan C Walker
- Department of Pathology and Biomedical Science, University of Otago, Christchurch, New Zealand.
| |
Collapse
|
3
|
Kothari C, Clemenceau A, Ouellette G, Ennour-Idrissi K, Michaud A, C.-Gaudreault R, Diorio C, Durocher F. TBC1D9: An Important Modulator of Tumorigenesis in Breast Cancer. Cancers (Basel) 2021; 13:3557. [PMID: 34298771 PMCID: PMC8304074 DOI: 10.3390/cancers13143557] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 07/05/2021] [Accepted: 07/12/2021] [Indexed: 01/02/2023] Open
Abstract
Triple-negative breast cancer (TNBC) is a major concern among the different subtypes of breast cancer (BC) due to the lack of effective treatment. In a previous study by our group aimed at understanding the difference between TNBC and non-TNBC tumors, we identified the gene TBC1 domain family member 9 (TBC1D9), the expression of which was lower in TNBC as compared to non-TNBC tumors. In the present study, analysis of TBC1D9 expression in TNBC (n = 58) and non-TNBC (n = 25) patient tumor samples validated that TBC1D9 expression can differentiate TNBC (low) from non-TNBC (high) samples and that expression of TBC1D9 was inversely correlated with grade and proliferative index. Moreover, we found that downregulation of the TBC1D9 gene decreases the proliferation marginally in non-TNBC and was associated with increased migratory and tumorigenic potential in both TNBC and luminal BC cell lines. This increase was mediated by the upregulation of ARL8A, ARL8B, PLK1, HIF1α, STAT3, and SPP1 expression in TBC1D9 knockdown cells. Our results suggest that TBC1D9 expression might limit tumor aggressiveness and that it has a differential expression in TNBC vs. non-TNBC tumors.
Collapse
Affiliation(s)
- Charu Kothari
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada; (C.K.); (A.C.); (G.O.); (R.C.-G.)
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
| | - Alisson Clemenceau
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada; (C.K.); (A.C.); (G.O.); (R.C.-G.)
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
| | - Geneviève Ouellette
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada; (C.K.); (A.C.); (G.O.); (R.C.-G.)
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
| | - Kaoutar Ennour-Idrissi
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
- Département de Biologie Moléculaire, de Biochimie Médicale et de Pathologie, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada
- Département de Médecine Sociale et Préventive, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada
| | - Annick Michaud
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
| | - René C.-Gaudreault
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada; (C.K.); (A.C.); (G.O.); (R.C.-G.)
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
- Laboratoire de Chimie Médicinale, l’Hôpital Saint-François d’Assise, Université Laval, Québec City, QC G1L 3L5, Canada
| | - Caroline Diorio
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
- Département de Médecine Sociale et Préventive, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada
- Centre des Maladies du Sein, Hôpital du Saint-Sacrement, Québec City, QC G1S 4L8, Canada
| | - Francine Durocher
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec City, QC G1T 1C2, Canada; (C.K.); (A.C.); (G.O.); (R.C.-G.)
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec City, QC G1V 4G2, Canada; (K.E.-I.); (A.M.); (C.D.)
| |
Collapse
|
4
|
Is Carboxypeptidase B1 a Prognostic Marker for Ductal Carcinoma In Situ? Cancers (Basel) 2021; 13:cancers13071726. [PMID: 33917306 PMCID: PMC8038727 DOI: 10.3390/cancers13071726] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2021] [Revised: 03/30/2021] [Accepted: 04/02/2021] [Indexed: 12/17/2022] Open
Abstract
Simple Summary Ductal carcinoma in situ (DCIS) is an early-stage breast cancer (BC), in which tumor cells are growing in a localized duct of the mammary gland. DCIS is considered a precursor disease for invasive BC and, therefore, treated as soon as it is identified. However, low-grade DCIS can be confused with atypical ductal hyperplasia, which is not a malignant lesion, leading to unnecessary surgery in around 70% of women with suspected DCIS. On the other hand, if left untreated, a DCIS has the potential to progress to IDC. In this retrospective study, we identified a gene signature, carboxypeptidase B1 (CPB1), the expression of which could help differentiate DCIS from an ADH lesion and DCIS that may progress to an invasive BC. Abstract Ductal carcinoma in situ (DCIS) is considered a non-obligatory precursor for invasive ductal carcinoma (IDC). Around 70% of women with atypical ductal hyperplasia (ADH) undergo unnecessary surgery due to the difficulty in differentiating ADH from low-grade DCIS. If untreated, 14–60% of DCIS progress to IDC, highlighting the importance of identifying a DCIS gene signature. Human transcriptome data of breast tissue samples representing each step of BC progression were analyzed and high expression of carboxypeptidase B1 (CPB1) expression strongly correlated with DCIS. This was confirmed by quantitative PCR in breast tissue samples and cell lines model. High CPB1 expression correlated with better survival outcome, and mRNA level was highest in DCIS than DCIS adjacent to IDC and IDC. Moreover, loss of CPB1 in a DCIS cell line led to invasive properties associated with activation of HIF1α, FN1, STAT3 and SPP1 and downregulation of SFRP1 and OS9. The expression of CPB1 could predict 90.1% of DCIS in a cohort consisting of DCIS and IDC. We identified CPB1, a biomarker that helps differentiate DCIS from ADH or IDC and in predicting if a DCIS is likely to progress to IDC, thereby helping clinicians in their decisions.
Collapse
|
5
|
Kothari C, Diorio C, Durocher F. Gene signatures of breast cancer development and the potential for novel targeted treatments. Pharmacogenomics 2021; 21:157-161. [PMID: 31967517 DOI: 10.2217/pgs-2019-0158] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- Charu Kothari
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec, QC G1V 0A6, Canada.,Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec, QC G1V 0A6, Canada
| | - Caroline Diorio
- Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec, QC G1V 0A6, Canada.,Département de Médecine Sociale et Préventive, Faculté de Médecine, Université Laval, Québec, QC G1V 0A6, Canada
| | - Francine Durocher
- Département de Médecine Moléculaire, Faculté de Médecine, Université Laval, Québec, QC G1V 0A6, Canada.,Centre de Recherche sur le Cancer, Centre de Recherche du CHU de Québec-Université Laval, Québec, QC G1V 0A6, Canada
| |
Collapse
|
6
|
Genome-Wide Gene Expression Analyses of BRCA1- and BRCA2-Associated Breast and Ovarian Tumours. Cancers (Basel) 2020; 12:cancers12103015. [PMID: 33081408 PMCID: PMC7603076 DOI: 10.3390/cancers12103015] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Revised: 09/28/2020] [Accepted: 10/14/2020] [Indexed: 12/13/2022] Open
Abstract
Germline pathogenic variants in BRCA1 and BRCA2 increase cumulative lifetime risk up to 75% for breast cancer and 76% for ovarian cancer. Genetic testing for BRCA1 and BRCA2 pathogenic variants has become an important part of clinical practice for cancer risk assessment and for reducing individual risk of developing cancer. Genetic testing can produce three outcomes: positive (a pathogenic variant), uninformative (no pathogenic variant) and uncertain significance (a variant of unknown clinical significance). More than one third of BRCA1 and BRCA2 variants identified have been classified as variants of uncertain significance, presenting a challenge for clinicians. To address this important clinical challenge, a number of studies have been undertaken to establish a gene expression phenotype for pathogenic BRCA1 and BRCA2 variant carriers in several diseased and normal tissues. However, the consistency of gene expression phenotypes described in studies has been poor. To determine if gene expression analysis has been a successful approach for variant classification, we describe the design and comparability of 23 published gene expression studies that have profiled cells from BRCA1 and BRCA2 pathogenic variant carriers. We show the impact of advancements in expression-based technologies, the importance of developing larger study cohorts and the necessity to better understand variables affecting gene expression profiles across different tissue types.
Collapse
|
7
|
Kothari C, Ouellette G, Labrie Y, Jacob S, Diorio C, Durocher F. Identification of a gene signature for different stages of breast cancer development that could be used for early diagnosis and specific therapy. Oncotarget 2018; 9:37407-37420. [PMID: 30647841 PMCID: PMC6324778 DOI: 10.18632/oncotarget.26448] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2018] [Accepted: 11/26/2018] [Indexed: 12/16/2022] Open
Abstract
Breast cancer (BC) is a heterogeneous disease where the survival rate of patients decreases with progression of the disease. BC usually has a linear progression, classified into normal/benign, atypical ductal hyperplasia (ADH), ductal carcinoma in situ (DCIS), and invasive ductal carcinoma (IDC). This study aimed to identify gene signature for each of these subgroups. We performed human transcriptome array analysis on 5 patient samples from each Normal, ADH, IDC and DCIS and 2 replicates of MCF10A cell line representative of each subgroup. We identified SFRP1 and snoRNAs (especially SNORD115 and SNORD114) as the initial regulators of cancer progression, accompanied by significant changes in extracellular matrix organization. Tumor progression to the IDC stage showed upregulation of tumor promoting genes responsible for increased invasion, inflammation, survival in stress environment and metastasis. The gene signatures identified in this study could represent potential biomarkers for each subgroup of breast cancer progression, which could assist in early diagnosis of breast cancer progression as well as treatment interventions. Moreover, these gene signatures could serve in discovery of specific targeted therapies for each subgroup.
Collapse
Affiliation(s)
- Charu Kothari
- Département de médecine moléculaire, Faculté de médecine, Université Laval, Québec, Canada
- Centre de recherche sur le cancer, Centre de recherche du CHU de Québec-Université Laval, Québec, Canada
| | - Geneviève Ouellette
- Département de médecine moléculaire, Faculté de médecine, Université Laval, Québec, Canada
- Centre de recherche sur le cancer, Centre de recherche du CHU de Québec-Université Laval, Québec, Canada
| | - Yvan Labrie
- Département de médecine moléculaire, Faculté de médecine, Université Laval, Québec, Canada
- Centre de recherche sur le cancer, Centre de recherche du CHU de Québec-Université Laval, Québec, Canada
| | - Simon Jacob
- Centre de recherche sur le cancer, Centre de recherche du CHU de Québec-Université Laval, Québec, Canada
- Laboratoire de pathologie, Hôpital du Saint-Sacrement, CHU de Québec-Université Laval, Québec, Canada
| | - Caroline Diorio
- Centre de recherche sur le cancer, Centre de recherche du CHU de Québec-Université Laval, Québec, Canada
- Département de médecine sociale et préventive, Faculté de médecine, Université Laval, Québec, Canada
| | - Francine Durocher
- Département de médecine moléculaire, Faculté de médecine, Université Laval, Québec, Canada
- Centre de recherche sur le cancer, Centre de recherche du CHU de Québec-Université Laval, Québec, Canada
| |
Collapse
|
8
|
Abstract
Background Inflammation is a major player in breast cancer (BC) progression. Allograft-inflammatory factor-1 (AIF1) is a crucial mediator in the inflammatory response. AIF1 reportedly plays a role in BC, but the mechanism remains to be elucidated. We identified two AIF1 isoforms, AIF1v1 and AIF1v3, which were differentially expressed between affected and unaffected sisters from families with high risk of BC with no deleterious BRCA1/BRCA2 mutations (BRCAX). We investigated potential functions of AIFv1/v3 in BC of varying severity and breast adipose tissue by evaluating their expression, and association with metabolic and clinical parameters of BC patients. Methods AIF1v1/v3 expression was determined in BC tissues and cell lines using quantitative real-time PCR. Potential roles and mechanisms were examined in the microenvironment (fibroblasts, adipose tissue, monocytes and macrophages), inflammatory response (cell reaction in BC subgroups), and metabolism [treatment with docosahexaenoic acid (DHA)]. Association of AIF1 transcript expression with clinical factors was determined by Spearman’s rank correlation. Bioinformatics analyses were performed to characterize transcripts. Results AIF1v1/v3 were mostly expressed in the less severe BC samples, and their expression appeared to originate from the tumor microenvironment. AIF1 isoforms had different expression rates and sources in breast adipose tissue; lymphocytes mostly expressed AIF1v1 while activated macrophages mainly expressed AIF1v3. Bioinformatics analysis revealed major structural differences suggesting distinct functions in BC progression. Lymphocytes were the most infiltrating cells in breast tumors and their number correlated with AIF1v1 adipose expression. Furthermore, DHA supplementation significantly lowered the expression of AIF1 isoforms in BRCAX cell lines. Finally, the expression of AIF1 isoforms in BC and breast adipose tissue correlated with clinical parameters of BC patients. Conclusions Results strongly suggest that AIF1v1 as much as AIF1v3 play a major role in the crosstalk between BC and infiltrating immune cells mediating tumor progression, implying their high potential as target molecules for BC diagnostic, prognostication and treatment. Electronic supplementary material The online version of this article (10.1186/s12935-018-0663-3) contains supplementary material, which is available to authorized users.
Collapse
|