Comprehensive analysis of microRNAs in breast cancer

  • Hong-Tai Chang1, 2Email author,

    Affiliated with

    • Sung-Chou Li3Email author,

      Affiliated with

      • Meng-Ru Ho4,

        Affiliated with

        • Hung-Wei Pan5,

          Affiliated with

          • Luo-Ping Ger5,

            Affiliated with

            • Ling-Yueh Hu6,

              Affiliated with

              • Shou-Yu Yu5,

                Affiliated with

                • Wen-Hsiung Li4, 7 and

                  Affiliated with

                  • Kuo-Wang Tsai5, 8Email author

                    Affiliated with

                    BMC Genomics201213(Suppl 7):S18

                    DOI: 10.1186/1471-2164-13-S7-S18

                    Published: 13 December 2012

                    Abstract

                    Background

                    MicroRNAs (miRNAs) are short noncoding RNAs (approximately 22 nucleotides in length) that play important roles in breast cancer progression by downregulating gene expression. The detailed mechanisms and biological functions of miRNA molecules in breast carcinogenesis have yet to be fully elucidated. This study used bioinformatics and experimental approaches to conduct detailed analysis of the dysregulated miRNAs, arm selection preferences, 3' end modifications, and position shifts in isoforms of miRNAs (isomiRs) in breast cancer.

                    Methods

                    Next-generation sequencing (NGS) data on breast cancer was obtained from the NCBI Sequence Read Archive (SRA). The miRNA expression profiles and isomiRs in normal breast and breast tumor tissues were determined by mapping the clean reads back to human miRNAs. Differences in miRNA expression and pre-miRNA 5p/3p arm usage between normal and breast tumor tissues were further investigated using stem-loop reverse transcription and real-time polymerase chain reaction.

                    Results

                    The analysis identified and confirmed the aberrant expression of 22 miRNAs in breast cancer. Results from pathway enrichment analysis further indicated that the aberrantly expressed miRNAs play important roles in breast carcinogenesis by regulating the mitogen-activated protein kinase (MAPK) signaling pathway. Data also indicated that the position shifts in isomiRs and 3' end modifications were consistent in breast tumor and adjacent normal tissues, and that 5p/3p arm usage of some miRNAs displayed significant preferences in breast cancer.

                    Conclusions

                    Expression pattern and arm selection of miRNAs are significantly varied in breast cancers through analyzing NGS data and experimental approach. These miRNA candidates have high potential to play critical roles in the progression of breast cancer and could potentially provide as targets for future therapy.

                    Background

                    Breast cancer is one of the major causes of cancer-related deaths worldwide and the most common cancer among women [1]. Metastatis to distant organs and lymph nodes represents a major problem, usually leading to high mortality. The investigation of breast cancer-associated genes for early detection or therapeutic targeting could potentially improve the survival rates of breast cancer patients.

                    MicroRNAs (miRNAs) are small RNA molecules with important regulatory functions in several physiological activities [2]. MicroRNAs are processed from primary transcripts (pri-miRNAs) in 2 maturation steps. First, the pri-miRNAs are processed by Drosha, forming the precursor miRNAs (pre-miRNAs), composed of a 5p arm, a 3p arm, and a terminal loop, approximately 70 nucleotides in length. Following the transport of pre-miRNAs to the cytoplasm by exportin 5, they are further processed by Dicer to release the terminal loop and the duplex (5p arm/3p arm), 22 nucleotides in length. The 5p arm/3p arm of the duplex is unwound at the end because of weaker hydrogen binding. The 5p or the 3p arm is selectively loaded into the RNA-induced silencing complex (RISC) and serves as mature miRNA [35]. Recent studies described a phenomenon in which RNA editing or nucleotide addition generated 3' end sequence variants of miRNAs [613]. Fernandez-Valverde et al. [6] reported that miR-282 and miR-312a are enriched for 3' adenosine additions during early embryonic development, which increases miRNA stability or enhances miRNA and mRNA interaction.

                    MicroRNAs exert their effects by repressing their target genes. They downregulate target gene expression by repressing translation or by degrading mRNAs. Previous studies reported that miRNAs play important roles in the oncogenesis pathway [1418]. The tumor-associated miRNAs were either tumor repressors, preferentially expressed in normal tissue, or onco-miRNAs, preferentially expressed in tumor tissue. These are aberrantly expressed in human breast cancer, including miR-9, miR-21, miR-31, miR-34a, miR-155, miR-200, miR-205, miR-206, and miR-335 [1922]. Although several studies have investigated the functions of miRNAs in breast tumors, these only included a small fraction of existing miRNAs [2326]. Using miRNA profiling approach, numerous breast cancer-associated miRNAs were identified [19, 2730]. Ryu et al. [27] identified 189 candidate novel microRNAs in human breast cancer cell lines by deep sequencing technology. Therefore, emerging NGS technologies can be used not only to identify novel miRNAs, but can also be applied in several miRNA-associated studies.

                    In the studies using NGS data for miRNA profiling, it is usually observed that miRNA sequence reads exist as isoforms, named isomiRs, with position and length shift compared with the reference miRNAs [31]. Recently, more and more studies worked on the isomiR issues, such as isomiR pattern preferences in specific libraries, target gene selection difference between different isomiRs and so on [6, 32]. Therefore, NGS data provides a good resource for miRNA expression profiling and isomiR related studies. In 2011, Farazi et al., [33] used NGS data to determine miRNA expression profiles in breast tissues with differing tumor malignancies. They focused on the relevance of specific miRNAs and the tumor malignancy type, without providing further experimental validation. The present study applied their NGS data to conduct analysis of miRNA-associated changes in breast cancer, including differential miRNA expression, position shifts in isomiRs, 3' end modifications, and arm selection preferences of pre-miRNAs.

                    Materials and methods

                    Collection and preprocessing of sequence reads

                    The small RNA transcriptome data of breast tumor (accession number: SRP006574) was downloaded from the NCBI Sequence Read Archive (SRA). These data included more than 200 samples and was classified into 2 libraries: normal and tumor (invasive ductal carcinoma). The initial sequence reads were subjected to 3' adaptor trimming to generate the clean reads, as described previously [10, 13]. For higher confidence, only the clean reads with read count ≥ 2 were included in further analysis.

                    Mapping clean reads to pre-miRNAs

                    MicroRNA expression profiles in different libraries were determined by mapping the clean reads back to human pre-miRNAs (miRBase 17). Several miRNA genes show high similarity (such as 68 mir-548 paralogous miRNAs in humans). This results in multiple ambiguous hits when mapping a read back to human miRNAs if variations are allowed. To eliminate the ambiguous mapping hits, no mismatch was allowed during the mapping procedure. Previous reports described observing nucleotide additions at the 3' ends of miRNAs [31, 3436] that could cause mismatches at the terminal of the mapping alignment. In order to follow the no-mismatch policy and keep the 3' end variation, like the method in Fernandez-Valverde's study [6], we trimmed out the terminal 3' end mismatch one by one until the perfect match reads were at least 18 nucleotides in length. By doing so, we can keep not only an at least 18-nt perfect alignment but also the 3' end variations.

                    Classifying non-miRNA reads into different data sets

                    The non-miRNA sequence reads were further classified into 9 classes by mapping to different data sets with bowtie [37], allowing a single nucleotide variation. The sequences of mRNAs and other ncRNAs were derived from NCBI RefSeq 47 [38]. The sequences of tRNAs was downloaded from the Genomic tRNA database [39]; sequences of rRNAs were downloaded from the SILVA database [40]. The sequences of snoRNAs, scaRNAs, and snRNAs were all downloaded from NONCODE [41]. The sequence reads not belonging to any of the described RNA classes were uploaded to RepeatMasker for identification of repeat elements and classified as unknown.

                    Samples and RNA extraction

                    Ten paired (tumor and adjacent normal) samples were collected from breast cancer patients receiving surgical operation at the Department of Surgery, Kaohsiung Veterans General Hospital. The total RNA of tissue was extracted using a TRIzol reagent (Invitrogen, USA), according to the instruction manual. Briefly, tissue samples were homogenized in 1 ml of TRIzol reagent and mixed with 0.2 ml chloroform to extract protein; RNA was precipitated using 0.5 ml isopropanol. The concentration, purity, and amount of total RNA were determined using a Nanodrop 1000 spectrophotometer (Nanodrop Technologies Inc., USA).

                    Stem-loop reverse transcription (RT) and real-time PCR

                    Reverse transcription primers were specifically designed for the examined miRNAs according to the methods reported by Chen et al., [42]. One microgram of total RNA was reverse transcribed in a stem-loop RT reaction with RT primers and SuperScript III Reverse Transcriptase according to the user's manual (Invitrogen, Carlsbad, CA, USA). The reaction was performed under the following incubation conditions: 30 min at 16°C, followed by 50 cycles of 20°C for 30 s, 42°C for 30 s, and 50°C for 1 s. The enzyme was subsequently inactivated by incubating at 85°C for 5 min. Real-time PCR reactions were performed using an miRNA-specific forward primer and a universal reverse primer with incubation at 94°C for 10 min, followed by 40 cycles of 94°C for 15s and 60°C for 32 s. Gene expression levels were detected using SYBR Green I assay (Applied Biosystems, Foster City, CA, USA), and miRNA expression levels were normalized to that of U6. The primer sequences for the examined miRNAs are listed in Additional File 1.

                    Pathway enrichment analysis

                    Human miRNA target gene data was downloaded from TargetScan 6.0. The target genes of differentially expressed miRNAs were extracted, then mapped onto KEGG pathways based on the enzyme commission (EC) numbers using the R package SubPathwayMiner v.3.1 [43]. The hypergeometric test was then performed to identify significantly enriched pathways and calculate the false positive discovery rate in FDR-corrected q-value.

                    Results and discussion

                    Analysis of miRNA sequence reads

                    After subjecting the downloaded small RNA reads to the 3' adaptor trimming procedure, the normal library contained approximately 2.7 million clean reads and the tumor library approximately 38.3 million clean reads (Table 1). Using the mapping criteria, the clean reads were mapped back to human pre-miRNAs (miRBase 17). Approximately 75% of the normal clean reads and 85% of the tumor clean reads belonged to miRNAs. The detection of a greater number of miRNA and pre-miRNAs in the tumor library than in the normal library reflected greater sequencing depth in the tumor library. The libraries contained an unequal number of initially used reads; therefore, comparisons in miRNA expression between the libraries were performed using the unit of transcript per million (TPM). There were no significant differences in the expression levels of most miRNAs between normal and tumor tissues (Additional File 2).
                    Table 1

                    Categories of sequence reads in the 2 breast libraries.

                    Library

                    Normal

                    Tumor

                    # clean reads

                    2 785 848

                    38 335 412

                    % miRNA reads

                    75.24%

                    84.79%

                    # detected pre-miRNAs

                    455

                    689

                    # detected miRNAs

                    631

                    906

                    # detected miRNAs at opposite arm

                    54

                    150

                    % miRNA reads with 3' end modification

                    17.65%

                    12.45%

                    According to miRBase 17 annotation, there are 1424 human pre-miRNAs, encoding 1733 mature miRNAs. These pre-miRNAs and mature miRNAs result in 1902 pairs of pre-miRNA/miRNA.

                    The clean reads not defined as miRNAs were further classified into an additional 9 classes by mapping to other transcript sequences. As shown in Table 2, in most categories, the clean read distributions in the 2 libraries were similar. However, compared with the normal library, the tumor library had a higher percentage of miRNA reads but lower percentages of other ncRNA and unknown reads. After further investigation, more than 85% of the unknown category (data not shown) in the 2 libraries could not be mapped back to the human genome (hg19). Our previous study showed that such unmappable reads could be derived from the exon-exon junctions of novel alternative splicing transcripts or from the transcripts of infection viruses [44]. These unmappable reads warrant further investigation.
                    Table 2

                    Categories of clean sequence reads in the 2 breast libraries

                    Category

                    miRNA

                    mRNA

                    tRNA

                    rRNA

                    snoRNA

                    scaRNA

                    snRNA

                    other ncRNA

                    repeat

                    unknown

                    Normal

                    75.24%

                    2.16%

                    0.47%

                    0.87%

                    0.84%

                    0.04%

                    0.12%

                    9.78%

                    0.0349%

                    10.45%

                    Tumor

                    84.79%

                    0.86%

                    0.37%

                    0.33%

                    0.35%

                    0.04%

                    0.14%

                    6.40%

                    0.0006%

                    6.70%

                    Differentially expressed miRNAs

                    A comparison of the miRNA expression levels in normal and tumor breast tissues provided a list of differentially expressed miRNAs. Similar to other high throughput technologies, NGS data is susceptible to background noise, producing biased results. Stem-loop RT-PCR was thus used to further validate the 23 miRNAs with the highest fold changes in expression levels between normal breast and breast tumor tissues. Stem-loop RT-PCR reportedly shows a high correlation with NGS technologies [10, 45]. In common with RT-PCR assays, an internal control is used for comparing the same genes in different samples. In NGS assays, no such internal control is used. The fold changes in the same miRNAs between normal and tumor breast tissues, as detected using PCR and NGS, were thus compared. Among the 23 evaluated miRNAs, the fold changes detected using PCR and NGS were highly correlated, with a Pearson correlation coefficient of 0.89 (Figure 1). This was a similar finding to those of previous studies.
                    http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-S7-S18/MediaObjects/12864_2012_4428_Fig1_HTML.jpg
                    Figure 1

                    Comparison of fold changes in expression determined using NGS analysis and PCR. The scatter plot shows that NGS and PCR results were highly correlated, with a Pearson correlation coefficient of 0.89

                    Among the 23 examined miRNAs, 15 are preferentially expressed in normal tissue, whereas 8 are preferentially expressed in tumor tissue (Table 3). The tissue preferences of miRNAs determined using NGS are consistent with those identified using PCR with only one exception: hsa-miR-423-3p. The present study, therefore, generated an authentic list of differentially expressed miRNAs in breast cancer using an experimental approach.
                    Table 3

                    Validation of differentially expressed miRNAs using PCR.

                    miRNA

                    NGS normal

                    NGS tumor

                    NGS fc

                    PCR fc

                    Consistent?

                    has-let-7b

                    32 594.1

                    16 107.1

                    0.49

                    0.46

                    Y

                    hsa-let-7c

                    7 710.1

                    2 439.1

                    0.32

                    0.44

                    Y

                    hsa-miR-22

                    69 532.1

                    24 576.1

                    0.35

                    0.24

                    Y

                    hsa-miR-125b

                    44 624.1

                    6 613.1

                    0.15

                    0.37

                    Y

                    hsa-miR-143

                    38 419.1

                    17 359.1

                    0.45

                    0.14

                    Y

                    hsa-miR-144

                    3 548.1

                    351.1

                    0.10

                    0.30

                    Y

                    hsa-miR-145

                    28 682.1

                    6 435.1

                    0.22

                    0.26

                    Y

                    hsa-miR-193a-5p

                    1 836.1

                    298.1

                    0.16

                    0.27

                    Y

                    hsa-miR-193b

                    2 469.1

                    1 058.1

                    0.43

                    0.30

                    Y

                    hsa-miR-199b-5p

                    1 608.1

                    1 055.1

                    0.66

                    0.83

                    Y

                    hsa-miR-320a

                    16 734.1

                    3 060.1

                    0.18

                    0.75

                    Y

                    hsa-miR-378

                    12 610.1

                    1 542.1

                    0.12

                    0.24

                    Y

                    hsa-miR-497

                    8 005.1

                    2 010.1

                    0.25

                    0.29

                    Y

                    hsa-miR-99a

                    31 360.1

                    5 548.1

                    0.18

                    0.16

                    Y

                    hsa-miR-141

                    5007.1

                    14 541.1

                    2.90

                    5.88

                    Y

                    hsa-miR-375

                    240.1

                    2 372.1

                    9.88

                    37.95

                    Y

                    hsa-miR-425

                    592.1

                    2 529.1

                    4.27

                    1.93

                    Y

                    hsa-miR-203

                    197.1

                    1 045.1

                    5.30

                    9.80

                    Y

                    hsa-miR-342-3p

                    350.1

                    1 331.1

                    3.80

                    3.43

                    Y

                    hsa-miR-182

                    316.1

                    1 653.1

                    5.23

                    11.86

                    Y

                    hsa-miR-200b

                    1 362.1

                    4 162.1

                    3.06

                    12.61

                    Y

                    hsa-miR-183

                    118.1

                    1 067.1

                    9.04

                    6.12

                    Y

                    hsa-miR-423-3p

                    2 446.1

                    1 993.1

                    0.81

                    5.96

                    N

                    NGS normal and NGS tumor denote miRNA relative read counts (TPM) in normal and tumor tissues. The NGS fold change (NGS fc) values are the quotients of miRNA read counts in normal tissue divided by those in tumor tissue. The PCR fold change (PCR fc) values are the quotients of miRNA expression levels, relative to U6 in normal tissue, divided by those in tumor tissue, averaged from 4 independent experiments

                    Enrichment analysis of miRNA-involved pathway

                    Following the identification of the differentially expressed miRNAs, the subsequent stage was to identify their functions as defined by their target genes. Several computational methods can be used to identify the putative target genes of miRNAs [4649]. However, these computational methods typically depend on the hydrodynamic stability of the miRNA/3'UTR duplex, and usually produce several false positive results. The most recently developed method is dependent on computational identification and also on the enrichment analysis of a target gene pathway [50]. The present study applied the same strategy and did the analysis on the miRNAs with the following features: with high expression alteration and high expression level in at least one library. So, we selected the union of the target genes of hsa-miR-141 and hsa-miR-200b for tumor-preferring miRNAs, and the union of the target genes of hsa-miR-22, hsa-miR-125b, and hsa-miR-99a for normal-preferring miRNAs. The 2 unions of genes were individually subjected to pathway enrichment analysis. The pathway enrichment analysis result showed that the target genes of tumor-preferring miRNAs are significantly enriched in the mitogen-activated protein kinase (MAPK) pathway, with a p-value of 2.1E-6 (Additional file 3); while, the target genes of normal-preferring miRNAs are significantly enriched also in the mitogen-activated protein kinase (MAPK) pathway, with a p-value of 2.4E-6 (Additional file 4).

                    The MAPK pathways are highly conserved kinase modules involved in fundamental cellular processes such as growth, proliferation, migration, and apoptosis [51]. Studies have identified that the miR-200 family (miR-200a, miR-200b, miR-200c, miR-141 and miR-429) is overexpressed in breast cancer, promoting breast cancer metastasis and drug resistance [25, 52]. Studies have also reported that overexpression of miR-200 in epithelial cell lines leads to an inhibition of transforming growth factor-β (TGF-β) and induction of epithelial mesenchymal transition (EMT) [53, 54]. These results indicated that the miR-200 family plays dual roles in modulating breast cancer metastasis through the regulation of the complex MAPK signaling pathway. However, miR-99a plays opposite roles in the regulation of cancer progression. Oneyama et al. [55] reported that overexpression of miR-99a led to the suppression c-Src-transformed cell growth, by controlling the mTOR/FGFR3 pathway in various human cancers. In epithelial NMUMG cells, however, miR-99a promoted proliferation and migration by regulating TGF-β-induced breast EMT [56]. Imbalance in the MAPK signaling pathway can, therefore, lead to promotion or inhibition of cancer cell progression. This data indicates that aberrant miRNA expression can result in the dysregulation of breast cancer cell proliferation, apoptosis, cell cycle, and migration through regulation of the component genes of the complex MAPK signaling pathway.

                    Previous studies usually focus on the regulation relationship between one miRNA versus one pathway. Since many miRNAs are simultaneous up- or down-regulation in the same tissue, we are curious if they simultaneously act on the same pathway. So, we pick up the target gene unions of up- or down-regulated miRNAs, followed by pathway enrichment analysis on the target genes of the same union. Our result showed that the simultaneously up-regulated or down-regulated miRNAs execute their functions by acting on the same MAPK pathway. Actually, we also did the pathway enrichment analysis as usual, one miRNA versus one pathway. And, we got the same conclusion that the simultaneously up- or down-regulated miRNAs simultaneously act on the same pathway.

                    MicroRNA 3' end additional nontemplate nucleotides in breast cancer

                    MicroRNAs reportedly undergo RNA editing or nucleotide addition at the 3' end, causing a mismatch at the termini of the mapping alignments [31, 3436]. As shown in Figure 2, the isomiRs of miR-511 have adenine (A) or uracil (U) added at their 3' ends. Therefore, the 3' end additional nontemplate nucleotides widely occur during the miRNA maturation process. Further investigation of the 3' end nucleotide modification events using the alternative mapping procedure showed that more than 12% of miRNA reads undergo 3' end nucleotide modification in all libraries (Table 1). As shown in Additional File 5, A and U accounted for approximately 80% of modification events. In addition, AA, UU, AU, cytosine (C), and guanine (G) contributed at least 1% of 3' end nucleotide modification events. The 3' end nucleotide modification event, therefore, prefers the A or U nucleotide. The same nucleotide modification preferences were displayed by normal breast and breast tumor tissues.
                    http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-S7-S18/MediaObjects/12864_2012_4428_Fig2_HTML.jpg
                    Figure 2

                    Mapping results of hsa-miR-511-1. As shown in 5p:16-36, hsa-miR-511-1 encodes mature miRNA at its 5p arm only and miRNA spans from nucleotide 16 to 36 of the hairpin. The integer values on the left denote the read count of each isomiR. The comma-separated values denote the position shift in the isomiR relative to the miRBase annotated positions (16 to 36). The nucleotides in lowercase denote the sequence fragments originating from the 3' modification event (Additional file 5)

                    Position shifts in isomiRs in breast cancer

                    It has been widely observed that miRNA exists as isoforms, or isomiRs, generated by a position shift during the maturation process [31]. The present study detected all isomiRs by mapping the clean reads back to pre-miRNAs. For example, hsa-miR-511 has 5 isomiRs, each with differing length. The expression of an miRNA can be derived by summarizing the read counts of its isomiRs (Figure 2). Morin et al. [31] showed that the miRBase reference miRNA (the isomiR with position shift 0,0 in Figure 2) is not always the most abundant isomiR. Among the detected miRNAs in libraries, approximately 55% are the most abundant isomiRs; most of the remaining miRNAs are the second or third most abundant isomiRs.

                    Although highly expressed miRNAs tend to have more isomiR types, the isomiR distribution between different libraries can differ [10, 13], indicating their diverse regulatory roles. The present study's data showed that position shifts occurred more frequently at the 3' end of miRNAs. Detailed analysis of the position shifts in miRNAs at both ends revealed that the 2 libraries displayed similar patterns in position shift (Figure 3). At the 5' end, position shift "0" dominated, contributing 95% of all miRNA reads; position shift "-1" and "1" accounted for only 1% and 4% of miRNA reads. At the 3' end, position shift "0" contributed 56% of all miRNA reads, remaining dominant. Other position shifts "-2," "-1," "1," and "2" accounted for 5%, 20%, 15%, and 1% of miRNA reads, respectively. In summary, most isomiRs displayed higher frequency position shifts at their 3' end during the miRNA maturation procedure. Additional file 6 displays the overall mapping results of the miRNAs of all libraries.
                    http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-S7-S18/MediaObjects/12864_2012_4428_Fig3_HTML.jpg
                    Figure 3

                    Position shifts at the 5' end and 3' end of miRNA reads. The position shifts at the miRNA 5' end and 3' end can be measured by comparing with the location of the miRBase reference miRNAs. (a) The position shift at the 5' end is less diverse and highly dominated by "0"; (b) the position shift at the 3' end is more diverse

                    Arm selection preferences in breast cancer

                    According to miRBase annotation, several pre-miRNAs encode mature miRNAs at both arms; however, several encode mature miRNAs at only one arm. With increasing NGS sequencing depth, additional miRNAs can be detected at the opposite arms of the pre-miRNAs, originally annotated to encode mature miRNAs at only one arm [10, 13]. This phenomenon was observed by mapping the clean reads back to pre-miRNAs (Table 1). Figure 2 presents one case in which hsa-miR-511-1 was annotated to encode mature miRNA at its 5p arm according to miRBase 17. However, additional miRNA was detected at its 3p arm. The newly detected 3p miRNA displayed higher expression levels than the original 5p miRNA. Squadrito et al. [57] reported that miR-511-3p is the major active strand of miR-511 in tumor-associated macrophages. In NGS analysis, results showed a higher abundance of miR-511-3p than miR-511-5p in breast cancer cells. To further investigate regulation of 5p and 3p arm selection of pre-miRNAs during breast cancer progression, the pre-miRNAs annotated to encode mature miRNAs at both arms were selected and their ratios of 5p arm to 3p arm expression were further compared in normal breast and breast tumor tissues. Based on hydrogen bonding theory, the selection preference between the 5p arm and 3p arm of pre-miRNA is an intrinsic characteristic of pre-miRNA. Data indicated that most of the examined pre-miRNAs were compatible with this theory, and the selection preferences of the 5p and 3p arms were consistent in normal and tumor tissues. In 17 of the pre-miRNAs, 5p/3p selection preference was not always consistent during breast cancer progression (Table 4). Stem-loop real-time PCR was used to further examine the expression of the 5p arm and 3p arm miRNAs of miR-324 and miR-455 in 10 paired (tumor and adjacent normal) breast cancer tissue samples. In Figure 4, the NGS data showed that the expression ratios of the 5p arm to 3p arm miRNAs of miR-324 and miR-455 were individually increased and decreased in the tumor, respectively, compared with in normal tissue. The 10 paired tissue samples provided consistent results. The ratio of 5p to 3p expression in miR-324 was increased significantly in breast cancer cells. These results showed that the miRNAs of samples did not display 5p or 3p preference consistently, suggesting the existence of another regulatory mechanism. This issue warrants further investigation.
                    Table 4

                    Arm selection preference of 5p arm and 3p arm miRNAs in normal breast and breast tumor tissues.

                    pre-miRNA

                    Location

                    N5p

                    N3p

                    T5p

                    T3p

                    hsa-miR-214

                    mi:30-51; MA:71-92

                    88

                    226

                    133

                    171

                    hsa-miR-576

                    5p:16-37; 3p:55-76

                    7

                    44

                    19

                    17

                    hsa-miR-154

                    MA:15-36; mi:51-72

                    21

                    8

                    7

                    14

                    hsa-miR-193a

                    5p:21-42; 3p:55-76

                    1 836

                    1 127

                    488

                    2 384

                    hsa-miR-296

                    5p:14-34; 3p:48-69

                    24

                    34

                    8

                    4

                    hsa-miR-361

                    5p:6-27; 3p:45-67

                    105

                    215

                    1 043

                    310

                    hsa-miR-324

                    5p:16-38; 3p:53-72

                    233

                    341

                    317

                    91

                    hsa-miR-339

                    5p:15-37; 3p:50-72

                    215

                    217

                    298

                    194

                    hsa-miR-493

                    mi:16-37; MA:57-78

                    2

                    5

                    22

                    9

                    hsa-miR-455

                    5p:16-37; 3p:54-74

                    99

                    72

                    107

                    308

                    hsa-miR-664

                    mi:11-34; MA:49-71

                    58

                    8

                    21

                    87

                    hsa-miR-212

                    3pM:71-91

                    22

                    7

                    7

                    6

                    hsa-miR-142

                    5p:16-36; 3p:52-74

                    914

                    1 102

                    8 661

                    10 074

                    hsa-miR-362

                    5p:5-28; 3p:42-63

                    174

                    34

                    145

                    114

                    hsa-miR-376a-1

                    mi:7-28; MA:44-64

                    41

                    24

                    9

                    13

                    hsa-miR-382

                    5pM:11-32

                    217

                    23

                    51

                    30

                    hsa-miR-151

                    5p:11-31; 3p:47-67

                    946

                    1 249

                    2 721

                    3 109

                    N5p, N3p, T5p, and T3p denote miRNA relative read counts (TPM) of 5p arm and 3p arm in normal and tumor tissues

                    http://static-content.springer.com/image/art%3A10.1186%2F1471-2164-13-S7-S18/MediaObjects/12864_2012_4428_Fig4_HTML.jpg
                    Figure 4

                    The 5p/3p arm selection of miR-324 and miR-455 in 10 breast cancer patients. (A) Ratios of 5p and 3p arms of miR-324 and miR-455 as observed from NGS data. (B) and (C) Ratios of 5p and 3p arms of miR-324 and miR-455 in the breast cancer tissues of 10 patients as evaluated using stem-loop quantitative PCR (left panels). Ratios of 5p and 3p arm usage are indicated in the boxplots in the right panels. All samples were assessed in triplicate and analyzed using Student's t test (p < 0.05 was considered significant)

                    Conclusions

                    The present study performed a series of sequence analysis to evaluate miRNA-associated changes in breast cancer, including miRNA expression, arm selection, 3' end modifications, and position shifts in isomiRs. We identified 22 differentially expressed miRNAs in normal breast and breast tumor tissue that might be involved in breast cancer progression through regulation of MAPK signaling. MicroRNAs widely displayed 3' end modifications and position shifts in isomiRs in breast cancer. However, no significant differences emerged between normal breast and breast tumor tissue during carcinogenesis. Arm usage of some miRNAs displayed significant preferences in breast cancer, suggesting that hydrogen bonding theory does not sufficiently explain 5p or 3p arm selection during carcinogenesis. Further investigation of the possible effects of arm selection of miRNAs on breast carcinogenesis is needed. The present study's findings provide insights into breast cancer that might facilitate the development of future cancer therapy.

                    Declarations

                    Acknowledgements

                    This work was supported by grants from Kaohsiung Veterans General Hospital (VGHKS 101-010 and VGHKS101-118) and National Sciences Council (NSC100-2314-B-075B-008)

                    This article has been published as part of BMC Genomics Volume 13 Supplement 7, 2012: Eleventh International Conference on Bioinformatics (InCoB2012): Computational Biology. The full contents of the supplement are available online at http://​www.​biomedcentral.​com/​bmcgenomics/​supplements/​13/​S7.

                    Authors’ Affiliations

                    (1)
                    Department of Surgery, Kaohsiung Veterans General Hospital
                    (2)
                    Department of Emergency, Kaohsiung Veterans General Hospital
                    (3)
                    Genomics Research Center, Academia Sinica
                    (4)
                    Biodiversity Research Center, Academia Sinica
                    (5)
                    Department of Medical Education and Research, Kaohsiung Veterans General Hospital
                    (6)
                    Institute of Biomedical Sciences, Academia Sinica
                    (7)
                    Department of Ecology and Evolution, University of Chicago
                    (8)
                    Department of Biotechnology, Tajen University

                    References

                    1. Gravgaard KH, Lyng MB, Laenkholm AV, Sokilde R, Nielsen BS, Litman T, Ditzel HJ: The miRNA-200 family and miRNA-9 exhibit differential expression in primary versus corresponding metastatic tissue in breast cancer. Breast cancer research and treatment 2012.
                    2. Yekta S, Shih IH, Bartel DP: MicroRNA-directed cleavage of HOXB8 mRNA. Science (New York, NY) 2004,304(5670):594–596.View Article
                    3. Newman MA, Hammond SM: Emerging paradigms of regulated microRNA processing. Genes & development 2010,24(11):1086–1092.View Article
                    4. Trabucchi M, Briata P, Filipowicz W, Rosenfeld MG, Ramos A, Gherzi R: How to control miRNA maturation? RNA biology 2009,6(5):536–540.PubMedView Article
                    5. Slezak-Prochazka I, Durmus S, Kroesen BJ, van den Berg A: MicroRNAs, macrocontrol: regulation of miRNA processing. RNA (New York, NY) 2010,16(6):1087–1095.View Article
                    6. Fernandez-Valverde SL, Taft RJ, Mattick JS: Dynamic isomiR regulation in Drosophila development. RNA (New York, NY) 2010,16(10):1881–1888.View Article
                    7. Guo L, Lu Z: Global expression analysis of miRNA gene cluster and family based on isomiRs from deep sequencing data. Comput Biol Chem 2010,34(3):165–171.PubMedView Article
                    8. Cloonan N, Wani S, Xu Q, Gu J, Lea K, Heater S, Barbacioru C, Steptoe AL, Martin HC, Nourbakhsh E, et al.: MicroRNAs and their isomiRs function cooperatively to target common biological pathways. Genome biology 2011,12(12):R126.PubMedView Article
                    9. Guo L, Yang Q, Lu J, Li H, Ge Q, Gu W, Bai Y, Lu Z: A comprehensive survey of miRNA repertoire and 3' addition events in the placentas of patients with pre-eclampsia from high-throughput sequencing. PloS one 2011,6(6):e21072.PubMedView Article
                    10. Li SC, Liao YL, Chan WC, Ho MR, Tsai KW, Hu LY, Lai CH, Hsu CN, Lin WC: Interrogation of rabbit miRNAs and their isomiRs. Genomics 2011,98(6):453–459.PubMedView Article
                    11. Guo L, Li H, Liang T, Lu J, Yang Q, Ge Q, Lu Z: Consistent isomiR expression patterns and 3' addition events in miRNA gene clusters and families implicate functional and evolutionary relationships. Molecular biology reports 2012,39(6):6699–6706.PubMedView Article
                    12. Humphreys DT, Hynes CJ, Patel HR, Wei GH, Cannon L, Fatkin D, Suter CM, Clancy JL, Preiss T: Complexity of murine cardiomyocyte miRNA biogenesis, sequence variant expression and function. PloS one 2012,7(2):e30933.PubMedView Article
                    13. Li SC, Liao YL, Ho MR, Tsai KW, Lai CH, Lin WC: miRNA arm selection and isomiR distribution in gastric cancer. BMC Genomics 2012,13(Suppl 1):S13.PubMedView Article
                    14. Chan SH, Wu CW, Li AF, Chi CW, Lin WC: miR-21 microRNA expression in human gastric carcinomas and its clinical association. Anticancer Res 2008,28(2A):907–911.PubMed
                    15. Chen Z, Zeng H, Guo Y, Liu P, Pan H, Deng A, Hu J: miRNA-145 inhibits non-small cell lung cancer cell proliferation by targeting c-Myc. J Exp Clin Cancer Res 2010, 29:151.PubMedView Article
                    16. Hausler SF, Keller A, Chandran PA, Ziegler K, Zipp K, Heuer S, Krockenberger M, Engel JB, Honig A, Scheffler M, et al.: Whole blood-derived miRNA profiles as potential new tools for ovarian cancer screening. Br J Cancer 2010,103(5):693–700.PubMedView Article
                    17. Theodore SC, Rhim JS, Turner T, Yates C: MiRNA 26a expression in a novel panel of African American prostate cancer cell lines. Ethn Dis 2010,20(1 Suppl 1):S1. 96–100PubMed
                    18. Yu S, Lu Z, Liu C, Meng Y, Ma Y, Zhao W, Liu J, Yu J, Chen J: miRNA-96 suppresses KRAS and functions as a tumor suppressor gene in pancreatic cancer. Cancer research 2010,70(14):6015–6025.PubMedView Article
                    19. Volinia S, Galasso M, Sana ME, Wise TF, Palatini J, Huebner K, Croce CM: Breast cancer signatures for invasiveness and prognosis defined by deep sequencing of microRNA. Proceedings of the National Academy of Sciences of the United States of America 2012,109(8):3024–3029.PubMedView Article
                    20. Yan LX, Huang XF, Shao Q, Huang MY, Deng L, Wu QL, Zeng YX, Shao JY: MicroRNA miR-21 overexpression in human breast cancer is associated with advanced clinical stage, lymph node metastasis and patient poor prognosis. RNA (New York, NY 2008,14(11):2348–2360.View Article
                    21. Iorio MV, Ferracin M, Liu CG, Veronese A, Spizzo R, Sabbioni S, Magri E, Pedriali M, Fabbri M, Campiglio M, et al.: MicroRNA gene expression deregulation in human breast cancer. Cancer research 2005,65(16):7065–7070.PubMedView Article
                    22. Kong W, Yang H, He L, Zhao JJ, Coppola D, Dalton WS, Cheng JQ: MicroRNA-155 is regulated by the transforming growth factor beta/Smad pathway and contributes to epithelial cell plasticity by targeting RhoA. Molecular and cellular biology 2008,28(22):6773–6784.PubMedView Article
                    23. Adams BD, Furneaux H, White BA: The micro-ribonucleic acid (miRNA) miR-206 targets the human estrogen receptor-alpha (ERalpha) and represses ERalpha messenger RNA and protein expression in breast cancer cell lines. Mol Endocrinol 2007,21(5):1132–1147.PubMedView Article
                    24. Cui W, Zhang Y, Hu N, Shan C, Zhang S, Zhang W, Zhang X, Ye L: miRNA-520b and miR-520e sensitize breast cancer cells to complement attack via directly targeting 3'UTR of CD46. Cancer Biol Ther 2010,10(3):232–241.PubMedView Article
                    25. Kastl L, Brown I, Schofield AC: miRNA-34a is associated with docetaxel resistance in human breast cancer cells. Breast cancer research and treatment 2012.
                    26. Tan Y, Zhang B, Wu T, Skogerbo G, Zhu X, Guo X, He S, Chen R: Transcriptional inhibiton of Hoxd4 expression by miRNA-10a in human breast cancer cells. BMC Mol Biol 2009, 10:12.PubMedView Article
                    27. Ryu S, Joshi N, McDonnell K, Woo J, Choi H, Gao D, McCombie WR, Mittal V: Discovery of novel human breast cancer microRNAs from deep sequencing data by analysis of pri-microRNA secondary structures. PloS one 2011,6(2):e16403.PubMedView Article
                    28. Bockmeyer CL, Christgen M, Muller M, Fischer S, Ahrens P, Langer F, Kreipe H, Lehmann U: MicroRNA profiles of healthy basal and luminal mammary epithelial cells are distinct and reflected in different breast cancer subtypes. Breast cancer research and treatment 2011,130(3):735–745.PubMedView Article
                    29. Buffa FM, Camps C, Winchester L, Snell CE, Gee HE, Sheldon H, Taylor M, Harris AL, Ragoussis J: microRNA-associated progression pathways and potential therapeutic targets identified by integrated mRNA and microRNA expression profiling in breast cancer. Cancer research 2011,71(17):5635–5645.PubMedView Article
                    30. Hannafon BN, Sebastiani P, de las Morenas A, Lu J, Rosenberg CL: Expression of microRNA and their gene targets are dysregulated in preinvasive breast cancer. Breast Cancer Res 2011,13(2):R24.PubMedView Article
                    31. Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, et al.: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res 2008,18(4):610–621.PubMedView Article
                    32. Fehniger TA, Wylie T, Germino E, Leong JW, Magrini VJ, Koul S, Keppel CR, Schneider SE, Koboldt DC, Sullivan RP, et al.: Next-generation sequencing identifies the natural killer cell microRNA transcriptome. Genome Res 2011,20(11):1590–1604.View Article
                    33. Farazi TA, Horlings HM, Ten Hoeve JJ, Mihailovic A, Halfwerk H, Morozov P, Brown M, Hafner M, Reyal F, van Kouwenhove M, et al.: MicroRNA sequence and expression analysis in breast tumors by deep sequencing. Cancer research 2011,71(13):4443–4453.PubMedView Article
                    34. Ebhardt HA, Tsang HH, Dai DC, Liu Y, Bostan B, Fahlman RP: Meta-analysis of small RNA-sequencing errors reveals ubiquitous post-transcriptional RNA modifications. Nucleic Acids Res 2009,37(8):2461–2470.PubMedView Article
                    35. Landgraf P, Rusu M, Sheridan R, Sewer A, Iovino N, Aravin A, Pfeffer S, Rice A, Kamphorst AO, Landthaler M, et al.: A mammalian microRNA expression atlas based on small RNA library sequencing. Cell 2007,129(7):1401–1414.PubMedView Article
                    36. Reid JG, Nagaraja AK, Lynn FC, Drabek RB, Muzny DM, Shaw CA, Weiss MK, Naghavi AO, Khan M, Zhu H, et al.: Mouse let-7 miRNA populations exhibit RNA editing that is constrained in the 5'-seed/cleavage/anchor regions and stabilize predicted mmu-let-7a:mRNA duplexes. Genome Res 2008,18(10):1571–1581.PubMedView Article
                    37. Langmead B, Trapnell C, Pop M, Salzberg SL: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome biology 2009,10(3):R25.PubMedView Article
                    38. Pruitt KD, Tatusova T, Klimke W, Maglott DR: NCBI Reference Sequences: current status, policy and new initiatives. Nucleic Acids Res 2009,37(Database issue):D32–36.PubMedView Article
                    39. Chan PP, Lowe TM: GtRNAdb: a database of transfer RNA genes detected in genomic sequence. Nucleic Acids Res 2009,37(Database issue):D93–97.PubMedView Article
                    40. Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glockner FO: SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res 2007,35(21):7188–7196.PubMedView Article
                    41. Liu C, Bai B, Skogerbo G, Cai L, Deng W, Zhang Y, Bu D, Zhao Y, Chen R: NONCODE: an integrated knowledge database of non-coding RNAs. Nucleic Acids Res 2005,33(Database issue):D112–115.PubMedView Article
                    42. Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, et al.: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res 2005,33(20):e179.PubMedView Article
                    43. Li C, Li X, Miao Y, Wang Q, Jiang W, Xu C, Li J, Han J, Zhang F, Gong B, et al.: SubpathwayMiner: a software package for flexible identification of pathways. Nucleic Acids Res 2009,37(19):e131.PubMedView Article
                    44. Li SC, Chan WC, Lai CH, Tsai KW, Hsu CN, Jou YS, Chen HC, Chen CH, Lin WC: UMARS: Un-MAppable Reads Solution. BMC Bioinformatics 2011,12(Suppl 1):S9.PubMedView Article
                    45. Chen SJ, Chen GH, Chen YH, Liu CY, Chang KP, Chang YS, Chen HC: Characterization of Epstein-Barr virus miRNAome in nasopharyngeal carcinoma by deep sequencing. PloS one 2010.,5(9):
                    46. Krek A, Grun D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, MacMenamin P, da Piedade I, Gunsalus KC, Stoffel M, et al.: Combinatorial microRNA target predictions. Nat Genet 2005,37(5):495–500.PubMedView Article
                    47. Rajewsky N: microRNA target predictions in animals. Nat Genet 2006, (38 Suppl):S8–13.
                    48. Rehmsmeier M, Steffen P, Hochsmann M, Giegerich R: Fast and effective prediction of microRNA/target duplexes. RNA (New York, NY 2004,10(10):1507–1517.View Article
                    49. Thadani R, Tammi MT: MicroTar: predicting microRNA targets from RNA duplexes. BMC Bioinformatics 2006,7(Suppl 5):S20.PubMedView Article
                    50. Hou J, Lin L, Zhou W, Wang Z, Ding G, Dong Q, Qin L, Wu X, Zheng Y, Yang Y, et al.: Identification of miRNomes in human liver and hepatocellular carcinoma reveals miR-199a/b-3p as therapeutic target for hepatocellular carcinoma. Cancer Cell 2011,19(2):232–243.PubMedView Article
                    51. Dhillon AS, Hagan S, Rath O, Kolch W: MAP kinase signalling pathways in cancer. Oncogene 2007,26(22):3279–3290.PubMedView Article
                    52. Korpal M, Ell BJ, Buffa FM, Ibrahim T, Blanco MA, Celia-Terrassa T, Mercatali L, Khan Z, Goodarzi H, Hua Y, et al.: Direct targeting of Sec23a by miR-200s influences cancer cell secretome and promotes metastatic colonization. Nature medicine 2011,17(9):1101–1108.PubMedView Article
                    53. Eades G, Yao Y, Yang M, Zhang Y, Chumsri S, Zhou Q: miR-200a regulates SIRT1 expression and epithelial to mesenchymal transition (EMT)-like transformation in mammary epithelial cells. The Journal of biological chemistry 2011,286(29):25992–26002.PubMedView Article
                    54. Gregory PA, Bracken CP, Smith E, Bert AG, Wright JA, Roslan S, Morris M, Wyatt L, Farshid G, Lim YY, et al.: An autocrine TGF-beta/ZEB/miR-200 signaling network regulates establishment and maintenance of epithelial-mesenchymal transition. Molecular biology of the cell 2011,22(10):1686–1698.PubMedView Article
                    55. Oneyama C, Ikeda J, Okuzaki D, Suzuki K, Kanou T, Shintani Y, Morii E, Okumura M, Aozasa K, Okada M: MicroRNA-mediated downregulation of mTOR/FGFR3 controls tumor growth induced by Src-related oncogenic pathways. Oncogene 2011,30(32):3489–3501.PubMedView Article
                    56. Turcatel G, Rubin N, El-Hashash A, Warburton D: MIR-99a and MIR-99b modulate TGF-beta induced epithelial to mesenchymal plasticity in normal murine mammary gland cells. PloS one 2012,7(1):e31032.PubMedView Article
                    57. Squadrito ML, Pucci F, Magri L, Moi D, Gilfillan GD, Ranghetti A, Casazza A, Mazzone M, Lyle R, Naldini L, et al.: miR-511–3p modulates genetic programs of tumor-associated macrophages. Cell reports 2012,1(2):141–154.PubMedView Article

                    Copyright

                    © Chang et al.; licensee BioMed Central Ltd. 2012

                    This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.