Prediction of clusters of miRNA binding sites in mRNA candidate genes of breast cancer subtypes

The development of breast cancer (BC) subtypes is controlled by distinct sets of candidate genes, and the expression of these genes is regulated by the binding of their mRNAs with miRNAs. Predicting miRNA associations and target genes is thus essential when studying breast cancer. The MirTarget program identifies the initiation of miRNA binding to mRNA, the localization of miRNA binding sites in mRNA regions, and the free energy from the binding of all miRNA nucleotides with mRNA. Candidate gene mRNAs have clusters (miRNA binding sites with overlapping nucleotide sequences). mRNAs of EPOR, MAZ and NISCH candidate genes of the HER2 subtype have clusters, and there are four clusters in mRNAs of MAZ, BRCA2 and CDK6 genes. Candidate genes of the triple-negative subtype are targets for multiple miRNAs. There are 11 sites in CBL mRNA, five sites in MMP2 mRNA, and RAB5A mRNA contains two clusters in each of the three sites. In SFN mRNA, there are two clusters in three sites, and one cluster in 21 sites. Candidate genes of luminal A and B subtypes are targets for miRNAs: there are 21 sites in FOXA1 mRNA and 15 sites in HMGA2 mRNA. There are clusters of five sites in mRNAs of ITGB1 and SOX4 genes. Clusters of eight sites and 10 sites are identified in mRNAs of SMAD3 and TGFB1 genes, respectively. Organizing miRNA binding sites into clusters reduces the proportion of nucleotide binding sites in mRNAs. This overlapping of miRNA binding sites creates a competition among miRNAs for a binding site. From 6,272 miRNAs studied, only 29 miRNAs from miRBase and 88 novel miRNAs had binding sites in clusters of target gene mRNA in breast cancer. We propose using associations of miRNAs and their target genes as markers in breast cancer subtype diagnosis.