In Silico Identification of a BRCA1:miR-29:DNMT3 Axis Involved in the Control of Hormone Receptors in BRCA1-Associated Breast Cancers

Abstract

Germline inactivating mutations in the BRCA1 gene lead to an increased lifetime risk of ovarian and breast cancer (BC). Most BRCA1-associated BC are triple-negative tumors (TNBC), aggressive forms of BC characterized by a lack of expression of estrogen and progesterone hormone receptors (HR) and HER2. How BRCA1 inactivation may favor the development of such a specific BC phenotype remains to be elucidated. To address this question, we focused on the role of miRNAs and their networks in mediating BRCA1 functions. miRNA, mRNA, and methylation data were retrieved from the BRCA cohort of the TCGA project. The cohort was divided into a discovery set (Hi-TCGA) and a validation set (GA-TCGA) based on the platform used for miRNA analyses. The METABRIC, GSE81002, and GSE59248 studies were used as additional validation data sets. BCs were differentiated into BRCA1-like and non-BRCA1-like based on an established signature of BRCA1 pathway inactivation. Differential expression of miRNAs, gene enrichment analysis, functional annotation, and methylation correlation analyses were performed. The miRNAs downregulated in BRCA1-associated BC were identified by comparing the miRNome of BRCA1-like with non-BRCA1-like tumors from the Hi-TCGA discovery cohort. miRNAs:gene-target anticorrelation analyses were then performed. The target genes of miRNAs downregulated in the Hi-TCGA series were enriched in the BRCA1-like tumors from the GA-TCGA and METABRIC validation data sets. Functional annotation of these genes revealed an over-representation of several biological processes ascribable to BRCA1 activity. The enrichment of genes related to DNA methylation was particularly intriguing, as this is an aspect of BRCA1 functions that has been poorly explored. We then focused on the miR-29:DNA methyltransferase network and showed that the miR-29 family, which was downregulated in BRCA1-like tumors, was associated with poor prognosis in these BCs and inversely correlated with the expression of the DNA methyltransferases DNMT3A and DNMT3B. This, in turn, correlated with the methylation extent of the promoter of HR genes. These results suggest that BRCA1 may control the expression of HR via a miR-29:DNMT3:HR axis and that disruption of this network may contribute to the receptor negative phenotype of tumors with dysfunctional BRCA1.