Single-nucleotide polymorphisms (SNP) associated with polygenetic disorders, such as breast cancer (BC), can create, destroy, or modify microRNA (miRNA) binding sites; however, the extent to which SNPs interfere with miRNA gene regulation and affect cancer susceptibility remains largely unknown. We hypothesize that disruption of miRNA target binding by SNPs is a widespread mechanism relevant to cancer susceptibility. To test this, we analyzed SNPs known to be associated with BC risk, in silico and in vitro, for their ability to modify miRNA binding sites and miRNA gene regulation and referred to these as target SNPs. We identified rs1982073-TGFB1 and rs1799782-XRCC1 as target SNPs, whose alleles could modulate gene expression by differential interaction with miR-187 and miR-138, respectively. Genome-wide bioinformatics analysis predicted ∼64% of transcribed SNPs as target SNPs that can modify (increase/decrease) the binding energy of putative miRNA::mRNA duplexes by >90%. To assess whether target SNPs are implicated in BC susceptibility, we conducted a case-control population study and observed that germline occurrence of rs799917-BRCA1 and rs334348-TGFR1 significantly varies among populations with different risks of developing BC. Luciferase activity of target SNPs, allelic variants, and protein levels in cancer cell lines with different genotypes showed differential regulation of target genes following overexpression of the two interacting miRNAs (miR-638 and miR-628-5p). Therefore, we propose that transcribed target SNPs alter miRNA gene regulation and, consequently, protein expression, contributing to the likelihood of cancer susceptibility, by a novel mechanism of subtle gene regulation. Cancer Res; 70(7); 2789–98
- Single-nucleotide polymorphism
- tumor susceptibility
MicroRNAs (miRNA) are a family of endogenous, short noncoding RNAs that modulate posttranscriptional gene regulation. They exert their regulatory role on protein-coding gene (PCG) expression by binding to either full or partial complementary sequences primarily in the 3′ untranslated region (UTR), but also inside the coding sequence (CDS) and the 5′ UTR, of the corresponding mRNAs. Eventually, this affects mRNA stability and translation (1–3). The role of miRNAs in human cancer pathogenesis has been well established by the identification of genetic alterations in miRNA loci, miRNA expression signatures that define different neoplastic phenotypes, and numerous oncogenes and tumor suppressor genes as miRNA targets (4).
The discovery of germline variants that influence cancer susceptibility is currently the subject of intense research (5). The most common variations are single-nucleotide polymorphisms (SNP), which occur approximately once every 100 to 300 bp in the human genome and constitute an immense source of information for unraveling the heterogeneity of human cancer pathogenesis, clinical course, and response to treatment. Several research groups have studied cancer risk association with SNPs located in candidate genes, and it was estimated that around 50,000 to 250,000 SNPs, mostly located in or around the 25,000 PCGs, have a putative biological effect (6). SNPs can affect protein function by changing the amino acid sequences (nonsynonymous SNP) or by perturbing their regulation [e.g., affecting promoter activity (7), splicing process (8), and DNA and pre-mRNA conformation]. When SNPs occur in 3′ UTRs, they may interfere with mRNA stability and translation by altering polyadenylation, protein::mRNA, and miRNA::mRNA regulatory interactions. Recent studies have found 3′ UTR SNPs that affect PCG expression via miRNA gene regulation in different diseases (9–14). These experimental evidences led us to hypothesize that the disruption of miRNA target binding by SNPs located in the 5′ UTRs, CDS, or 3′ UTRs is a widespread mechanism leading to cancer susceptibility and initiation. In support of our hypothesis, it was recently shown by Chin and colleagues (15) that a kRAS variant within the let-7 target site increased the risk for non–small cell lung carcinoma among moderate smokers. Moreover, because the 3′ UTRs of PCGs were insufficiently screened for mutations/polymorphisms, it is possible that the extent of such abnormalities might be much greater than initially predicted.
Breast cancer (BC) is one of the most common female malignancies, with more than 1 million new cases diagnosed every year (16). It is widely accepted that the majority of BC risk might be specified by the combined product of multiple low-penetrance alleles (see ref. 17 for review); however, the multiplicity of variants identified thus far alone cannot account for ∼80% of familial BC cases that are unrelated to high-penetrance BC susceptibility genes. The role of polymorphic variants located in BC-relevant genes in tumor susceptibility has been extensively addressed (Supplementary Table S1; refs. 18, 19). However, none of the reports investigating the pathogenetic implications of genomic variation on BC susceptibility has evaluated the role of SNPs in miRNA::mRNA gene regulation and its effect in BC development.
Here, we analyzed SNPs associated with BC susceptibility for their ability to affect miRNA binding sites and miRNA::mRNA gene regulation. We identified SNP-dependent miRNA interactions that might explain the pathogenetic relevance of known BC-associated SNPs. We performed a genome-wide scan for SNPs inside PCG transcribed regions (3′ UTR, CDS, and 5′ UTR) that could potentially alter miRNA binding (referred to as target SNPs). Finally, we illustrated the functional consequences of target SNPs on miRNA regulation of protein expression and their possible clinical implications due to differences in frequency distribution among familial BCs, sporadic BCs, and controls.
Materials and Methods
We obtained mature human miRNA sequences (total of 885) from National Center for Biotechnology Information (NCBI) build 36 database (version 13.0). We downloaded SNP data from HapMap (Release 21a). We retrieved 5′ UTR, 3′ UTR, and CDS mRNA genomic locations from the University of California at Santa Cruz (20) known gene table (human genome assembly, NCBI build 36, hg18).
miRNA::SNP interaction analyses
For each SNP, we retrieved two sequences, centered on each allele with 25 nucleotides flanking both sides. We used the miRNA target prediction program miRanda (2) with two different cutoffs (score ≥80, MFE ≤−16 kcal/mol and score ≥50, MFE ≤−5 kcal/mol) to calculate minimum free energy (MFE) for all the possible miRNA::SNP centered sequences. For all the predicted interactions, we computed the allele-dependent MFE changes. Based on the distribution of MFE changes, we identified the 20 and 80 percentile threshold values for miRNA::target SNP identification (Fig. 1). Based on these calculations, we termed SNP alleles as either active or nonactive alleles, the active allele being the variant that induces a decrease in MFE and a stronger miRNA binding with the target.
Patient sample collection
Blood samples from patients with familial and sporadic BC and from control subjects were collected at the Istituto Nazionale Tumori (Milan, Italy) with their informed consent. Total genomic DNA was purified from peripheral blood leukocytes using the QIAamp DNA Mini kit (Qiagen).
Selection of SNPs associated with BC
We performed a literature search by use of the MedLine/PubMed database (21) covering the time between January 2006 and December 2008 to retrieve articles reporting association studies between SNPs and BC risk. We selected SNPs that were confirmed to be associated with BC risk in one or more articles and/or from large-scale genome-wide association studies. Supplementary Table S1 shows the initial list of SNPs use for subsequent miRanda analysis.
DNA sequencing and SNP cloning
PCR amplification of the transcribed SNP-containing region was performed on genomic DNA (Platinum Taq High Fidelity, Invitrogen). Primers are available on request. Sequences were performed using the BigDye Terminator Reaction Chemistry v3.1 on Applied Biosystems 3730 DNA analyzers (Applied Biosystems). All sequence analyses and alignments were performed with the SeqmanPro program, Lasergene version 7.1 (DNASTAR). PCR-amplified SNP-containing regions carrying either the active or the nonactive alleles were XbaI cloned into the 3′ UTR of the pGL3-control vector (Promega).
Cell cultures, transfection, and immunoblotting
All cell lines were obtained from the American Type Culture Collection. Cells were transfected using Lipofectamine 2000 (Invitrogen). Total protein extracts were prepared in 0.5% NP40 lysis buffer. For immunoblotting, protein extracts were separated in 4% to 20% SDS-PAGE (Criterion Precast Gel, Bio-Rad) and transferred to nitrocellulose membranes (Bio-Rad). Antibodies for TGFB, TGFBR1, and XRCC1 were from Cell Signaling; antibodies for BRCA1 were from Calbiochem (clone MS110); and antibodies for vinculin were from Santa Cruz Biotechnology (clone N-19). Bands were quantified using GelDoc XR software (Bio-Rad).
MCF7 cells (200 × 105 per 24-well plates) were cotransfected with 0.4 μg of pGL3 (luciferase), 0.08 μg of pRLTK (Renilla; Promega), and 50 nmol/L of precursor miRNA molecules or scrambled negative control (Ambion). Thirty hours after transfection, cells were lysed in 100 μL of passive lysis buffer according to the dual luciferase reporter assay protocol (Promega), and luciferase activity was measured with Veritas luminometer (Turner BioSytems).
Fisher's exact test was used to determine the association of SNPs and patient disease status. Cochran-Armitage trend test was used to investigate whether there was a trend in binomial proportions of cancer patients (familiar + sporadic, familiar, or sporadic, respectively) across levels of the different genotypes for each SNP. The univariate and multivariate logistic regression models were fit to evaluate the association of SNP genotype with patient state status.
Additional methods, including patient information, statistical analysis, RNA extraction, retrotranscription, and real-time PCR, are available in Supplementary Materials and Methods.
Analysis of BC susceptibility alleles for their ability to affect miRNA binding
To test our hypothesis that SNPs can affect BC susceptibility by altering miRNA::mRNA binding, we took advantage of genetic variants that are already known to be associated with BC (see Materials and Methods and Supplementary Table S1). For each of these SNPs, we retrieved two sequences 51 nucleotides long, one centered on the common allele and the other centered on the variant allele. We scanned the retrieved sequences with miRanda (2) in search of miRNA::mRNA target sites and calculated the MFE changes induced by the allelic variants (Supplementary Table S2). We observed that the allelic variants can either increase or decrease the MFE of the corresponding RNA duplexes, leading to either a stronger or a weaker miRNA::mRNA binding, respectively. This mechanism can also lead to either creation of a new binding site or destruction of an existing target site. We focused only on those miRNAs (included in Supplementary Table S2) that were expressed in BC samples and cell lines (Supplementary Fig. S1; data not shown) and selected three candidate target SNPs for functional validation, for a total of 16 miRNA::SNP interactions (Supplementary Table S3). Two of the 16 interactions chosen were negative controls, in which the interacting miRNA binding was not modified by the SNP allelic variants. To verify whether the selected BC-associated SNPs can actually affect miRNA binding, we performed in vitro luciferase reporter assays. For each SNP, we produced two pGL3-SNP constructs (Fig. 2A) carrying either the active or the nonactive allele and cotransfected them in parallel with the predicted interacting miRNAs or the scrambled negative control. rs1982073 inside TGFB1 and rs1799782 inside XRCC1 (Table 1) showed statistically significant effect on miR-187 and miR-138 activity, respectively (Table 2; Fig. 2B). The rs1982073 [C] active allele of TGFB1 was predicted by miRanda analysis to create a new target site for miR-187, which was absent with the more common [T] variant when stringent settings were used. Luciferase assay showed statistically significant suppressive effect of miR-187 on the construct carrying the [C] active allele, which is completely absent with the [T] nonactive allele (40% decrease in luciferase activity, P = 0.006; Fig. 2B, left). This agrees with the major line of evidence assigning an inhibitory role of miRNAs on gene expression. Conversely, according to miRanda prediction, the rs1799782-XRCC1 [T] variant was the active allele increasing the binding energy of miR-138 compared with the [C] variant. When miR-138 was assayed with pGL3-rs1799782-XRCC1 constructs, we observed a significant increase in the luciferase activity with the [T] allele compared with the [C] allele (34% increase, P = 0.0003; Fig. 2B, right). This clearly indicates a stabilizing role, rather than an inhibitory role, for miR-138 stronger binding with the rs1799782-XRCC1 [T] variant.
Subsequently, we tested the effect of these two target SNPs on endogenous TGFB1 and XRCC1 protein levels in cancer cell lines carrying different SNP genotypes by transfecting with the interacting miRNAs. miR-187, consistent with the luciferase results, preferentially downregulated protein levels of the TGFB1 gene carrying the [CC] rs1982073 active alleles, whereas miR-187 had an intermediate and opposite effect on the [TC] and [TT] rs1982073 genotypes, respectively (Fig. 2C, left). Similarly, miR-138 showed a stabilizing effect when overexpressed in a cell line that was heterozygous for rs1799782-XRCC1. Whereas [CC] carriers displayed decreased XRCC1 protein levels, a cell line carrying the [CT] genotype displayed increased XRCC1 levels, in agreement with the luciferase results (Fig. 2C, right). Among the 25 different cancer cell lines we screened, we did not find any rs1799782-XRCC1 [TT] homozygous carrier, supporting the notion that [TT] carriers are protected from developing cancer (22). Taken together, these results suggest that for ∼15% of the investigated transcribed SNPs (2 of 14 target SNP interactions predicted to be allele dependent and tested in our analysis; Fig. 2; Supplementary Table S3; data not shown), their known association with BC can be biologically due to interference with miRNA binding and therefore miRNA gene regulation. The fact that only one seventh of the target SNP interactions could be biologically confirmed is not surprising, given the high rate of false positives of all miRNA target prediction programs. Although miRanda program is known to produce more false-positive predictions than other existing programs that give preference to seed complementarity or interspecies target conservation, it has a higher sensitivity. We, therefore, choose miRanda to identify noncanonical miRNA targets (i.e., with low seed complementarity) such as rs1982073-TGFB1 and rs1799782-XRCC1 that would have been otherwise overlooked by using prediction programs focusing only on 3′ UTR and considering extensive 5′ complementarity (data not shown). Overall, from our in silico predictions, we validated target SNP interactions that were predicted to have a wide range of MFE change (ranging from 8% to 42%) with no direct correlation between MFE and probability of a true miRNA::target SNP interaction. Therefore, for subsequent analyses, we filtered the predicted target interactions with a MFE change of at least 8%, at which, in our in vitro experiments, a significant biological effect was detected. Moreover, at this cutoff, through a genome-wide target SNP analyses (see below), we noticed a clear peak in the MFE change distribution, further supporting a possible biological significance for MFE changes >8% (Supplementary Fig. S2).
Identification and characterization of target SNPs that disrupt miRNA::mRNA interaction
Of 3,839,363 SNPs comprised in the human HapMap database (Release 21a), we identified 90,985 SNPs that are located in mRNA regions (i.e., transcribed SNPs) and grouped them according to their genomic locations (5′ UTRs, CDS, and 3′ UTRs). We identified candidate target SNPs that can potentially create, destroy, or modify miRNA binding sites due to allele-specific MFE changes through an integrated bioinformatics approach (Fig. 1; Materials and Methods) that includes the predictions of miRanda (2). In Supplementary Tables S4 to S6, we present the top 100 miRNA::SNP allele-specific MFE changes for SNPs located in 3′ UTR, CDS, and 5′ UTR, respectively (additional target SNP interactions are available on request). We observed highly similar bimodal distribution of MFE change for 3′ UTR, 5′ UTR, and CDS (Supplementary Fig. S2) and considered all SNPs that induce at least an 8% MFE change as target SNPs. Using this strategy (Fig. 1) and stringent miRanda settings (see Materials and Methods), we found that 41%, 6%, and 53% of target SNPs were located inside 5′ UTR, CDS, and 3′ UTR, respectively. In the MFE change distribution of all transcribed SNPs, the lower 2.5 and the higher 97.5 percentile corresponded to minimum 94% MFE change. By using 94% cutoff, we found that ∼64% of transcribed SNPs are predicted to disrupt (increase/decrease) the MFE of putative miRNA::mRNA duplexes. In summary, through this genome-wide bioinformatics approach, we identified a set of putative target SNPs that modify the MFE of miRNA::mRNA binding and eventually influence the interaction and the regulatory function of miRNAs with their targets.
Additionally, we performed a similar analysis on the 17 SNPs that Saunders and colleagues (23) identified inside experimentally verified miRNA target sites with the corresponding interacting miRNAs (Supplementary Table S7). We observed that 78% of these SNPs were competent to change the miRNA MFE by at least 8%. Therefore, variant alleles inside experimentally verified miRNA target sites can change the status of miRNA::mRNA interactions and affect gene expression, consistent with our genome-wide bioinformatics predictions.
Target SNP distribution in BC and control populations
We sequenced a panel of BC patients (166 sporadic and 169 familial BRCA1- and BRCA2-negative probands) and controls (186; see Supplementary Materials and Methods for patient information) for the germline presence of selected target SNPs located inside cancer-relevant genes (Table 3). Supplementary Table S8 summarizes the distribution of target SNP genotypes and their association with BC risk in Caucasian cases (sporadic and familial BC) and controls. The logistic regression analysis indicated that rs799917-BRCA1 [T] carriers were more likely associated with BC [all cases versus controls: univariate analysis odds ratio (OR) for [CT] carriers of 1.57 (P = 0.03), and still significant in multivariate analysis with OR of 1.86 (P = 0.046), and univariate analysis OR for [TT] carriers of 1.95 (P = 0.03)] and particularly with sporadic BC [sporadic cases versus controls: univariate analysis OR for CT carriers of 1.98 (P = 0.007) and OR for [TT] carriers of 2.81 (P = 0.003); Table 4]. Furthermore, the analysis suggested that rs334348-TGFBR1 [AG] carriers were more likely to have BC (all cases versus controls: univariate analysis OR, 1.69; P = 0.048) and particularly familiar BC (familiar versus control: univariate analysis OR, 2.2; P = 0.005; familiar versus sporadic: univariate analysis OR, 1.99; P = 0.01), and this association was also significant in the multivariate analysis (familiar versus control: multivariate analysis OR, 2.67; P = 0.002; familiar versus sporadic: multivariate analysis OR, 2.17; P = 0.02; Table 4). Despite the small size of the case-control study population, this evidence suggests that specific target SNPs have differential distribution among familial BCs, sporadic BCs, and controls.
Validation of target SNP predictions
To validate the computational predictions and the biological relevance of target SNPs, we first carried out in vitro luciferase reporter assays. The identified target SNPs that are located in BRCA1, TGFBR1, and XIAP (Table 3) affected significantly (P < 0.05) the pGL3-SNP luciferase activity by the predicted interacting miRNAs (Fig. 3A; Supplementary Fig. S3). In all three cases tested, miRNAs displayed a higher repressive effect when interacting with the active allele. Next, we tested the effect of the two target SNPs that we found to be associated with BC, rs799917-BRCA1 and rs334348-TGFBR1, on endogenous BRCA1 and TGFBR1 protein levels. We transfected the interacting miRNAs (miR-638 and miR-628-5p) in cancer cell lines carrying different target SNP genotypes. Following miR-638 overexpression, seven of nine cell lines (in which BRCA1 was clearly detected to perform correct band quantification) showed a reduction in BRCA1 protein levels (Fig. 3B, left; data not shown). By subgrouping cell lines responsive to miR-638–induced BRCA1 reduction according to rs799917 genotype, we observed that the presence of the [CC] genotype was responsible for a stronger reduction of BRCA1 protein levels (61% reduction versus 79%, P = 0.047; Fig. 3C, left). Additionally, we excluded the possibility that miR-638–dependent BRCA1 reduction was due to a conserved target site inside BRCA1 3′ UTR predicted by TargetScan via luciferase and BRCA1 overexpression experiments (Supplementary Fig. S4), further confirming the functional significance of the target SNP interaction site. Following similar approach, we evaluated the effect of rs334348 on miR-628-5p regulation of TGFBR1. By overexpressing miR-628-5p in different cell lines, we observed that miR-628-5p behaves as a true repressor of TGFBR1 (Fig. 3B, right; data not shown) in a cell-specific manner and that its repressor activity on TGFBR1 protein levels is dependent on the rs334348 variant inside the TGFBR1 3′ UTR miRNA target sequence (80% versus 50% reduction with the [AA] and the [CC] genotypes, respectively, P = 0.006; Fig. 3C, right). The observed significant effects on protein modulation by miRNA in the presence of the active or nonactive alleles (Figs. 2 and 3; Supplementary Fig. S3) support our hypothesis that genetic variants inside miRNA targets can actually disrupt miRNA gene regulation and affect protein expression, eventually leading to tumor susceptibility.
In the present study, we identified a new pathogenetic mechanism to explain the association of certain SNPs with BC susceptibility by using the latest knowledge on posttranscriptional gene regulation by miRNAs. Evidence of the effect of SNPs on miRNA binding derived from our integrative analysis delineates a novel role for SNPs in gene expression modulation. Supporting this theory, a recent work suggested that differences in SNP allele frequency among ethnic groups account for differences in gene expression (24). Two of the 11 cis SNPs identified by Spielman and colleagues (24) were actually transcribed SNPs located in 3′ UTRs, and by bioinformatics analysis, we found that both of the SNPs were qualified to induce MFE absolute changes >8% with multiple putative interacting miRNAs (Supplementary Table S9). Therefore, we speculate that cis SNPs modulate phenotypic gene expression diversities, at least in part, through alteration of miRNA target binding capability, ultimately leading to differences in the susceptibility to complex genetic diseases, such as BC. Notably, among the target SNPs identified as differentially expressed in BC, we found rs334348 located in the 3′ UTR of TGFBR1. The association of this SNP with germline allele-specific expression of TGFBR1 was recently found, and it was shown to confer an increased risk of colorectal cancer (25). This phenotype could be explained, for example, by an altered miRNA interaction and therefore strongly supports our hypothesis that SNPs affecting miRNA function contribute to allele-specific protein expression and consequently play a role in tumor susceptibility. A recent study by Kim and Bartel (26) based on transcriptome analyses of heterozygous mouse strains showed that polymorphic miRNA target sites determine differences in gene expression, which further supports our hypothesis.
The focus of this study is to address the effect of transcribed SNPs on miRNA::mRNA interactions and thereby the regulation of gene expression by these SNPs with possible implications in tumorigenesis. We found that numerous target SNPs can effectively interfere with miRNA target recognition, and the distribution of these SNPs significantly varies among BC and control populations, suggesting a potential role in BC susceptibility. It is interesting to note that the BC-associated rs1982073-TGFB1 variant had been previously linked with low expression levels of TGFB1 during the initial phases of tumorigenesis, whereas higher expression levels were found in more advanced metastatic BC stages (27, 28). This switch in TGFB1 expression could be due to concomitant changes of miR-187 expression levels and could be favored by carrying the [T] variant. Additionally, the described model, together with TGFB1 protumorigenic or antitumorigenic activities in different phases of tumorigenesis, can explain the conflicting results observed for rs1982073-TGFB1 association with BC. Another interesting observation is that the stronger binding occurring between miRNA and active allele does unequivocally mean repression; in fact, as shown for rs1799782-XRCC1, miR-187 stabilizes XRCC1 expression when carrying the active allele. This finding is in line with other evidence showing that miRNAs can positively modulate the protein expression (29, 30). Furthermore, by searching the complete PCG transcripts (5′ UTR, CDS, and 3′ UTR) for miRNA target SNPs, we have identified miR-638 as regulator of BRCA1 via binding to its target site inside CDS (to our knowledge, the first miRNA regulating BRCA1 published thus far). Although a conserved interaction site was predicted by TargetScan in the 3′ UTR of BRCA1, we found that it is not a functional site in the conditions we tested. Recently, Tay and colleagues (3) also proved the existence and abundance of miRNA targets located inside the CDS of PCGs, confirming the possibility that miR-638 regulates BRCA1 in a similar way. Although large case-control studies have not reached a unifying conclusion on the relevance to BC susceptibility of the miRNA target SNP rs7799917 inside BRCA1 CDS (31, 32), its association with tumor susceptibility has been recently highlighted in a SNP analysis of DNA repair genes in glioblastomas multiforme samples (33). In the present study, we found that the rs7799917 [T] allele is associated with a weaker miR-638–dependent BRCA1 reduction and is linked with BC (Table 4). Further investigations are under way to explain this apparent inconsistency. We speculate that this is due to concomitant variations during tumor development of miR-638 expression levels. Moreover, we observed that low-frequency CDS mutations (both synonymous and nonsynonymous; data not shown) might exert on miRNA gene regulation effects similar to target SNPs (34). This provides an innovative insight on the importance of synonymous mutations, which in the future should be considered and investigated with renewed attention.
Supported by the data presented here, we highlight the possibility that both genetic variants and a miRNA expression patterns can account for gene expression differences among distinct populations. Therefore, the newly recognized class of target SNPs may contribute to cancer susceptibility not per se (or as tags for specific haplotypes; refs. 31, 32) but in concert with miRNA expression patterns, which are becoming more easily available in this new era of miRNA profiling.
In summary, considering the multifactorial model of BC susceptibility, we propose that transcribed allelic variants can alter miRNA PCG regulation and thus increase miRNA abnormalities. These contribute to tumor susceptibility and tumor general risk by a mechanism of subtle gene regulation still not fully appreciated. As the lists of miRNAs and target SNPs keep growing, the need to confirm the effect of target SNPs on gene expression becomes more pressing and will eventually help us achieve a deeper insight into many variables of the molecular pathways relevant to the development of BC and other human malignancies.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
We thank Gustavo Baldassarre for the kind gift of BRCA1-expressing vector and Ellen Ko for English editing.
Grant Support: G.A. Calin is supported as a Fellow at The University of Texas M.D. Anderson Research Trust, as a Fellow of The University of Texas System Regents Research Scholar, and by the Ladjevardian Regents Research Scholar Fund. This study was supported in part by Breast Cancer Specialized Program of Research Excellence Developmental Award to G.A. Calin. R.V. Davuluri holds a Philadelphia Healthcare Trust Endowed Chair Position; research in his laboratory is supported by the Philadelphia Healthcare Trust, National Human Genome Research Institute/NIH, American Cancer Society, and grant # P30CA010815.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: Supplementary data for this article are available at Cancer Research Online (http://cancerres.aacrjournals.org/).
- Received September 24, 2009.
- Revision received December 30, 2009.
- Accepted January 19, 2010.
- ©2010 American Association for Cancer Research.