Identification of metastasis-associated breast cancer genes using a high-resolution whole genome profiling approach

A Focus on the Synergy of Radiomics and RNA Sequencing in Breast Cancer

Radiological imaging is currently employed as the most effective technique for screening, diagnosis, and follow up of patients with breast cancer (BC), the most common type of tumor in women worldwide. However, the introduction of the omics sciences such as metabolomics, proteomics, and molecular genomics, have optimized the therapeutic path for patients and implementing novel information parallel to the mutational asset targetable by specific clinical treatments. Parallel to the "omics" clusters, radiological imaging has been gradually employed to generate a specific omics cluster termed "radiomics". Radiomics is a novel advanced approach to imaging, extracting quantitative, and ideally, reproducible data from radiological images using sophisticated mathematical analysis, including disease-specific patterns, that could not be detected by the human eye. Along with radiomics, radiogenomics, defined as the integration of "radiology" and "genomics", is an emerging field exploring the relationship between specific features extracted from radiological images and genetic or molecular traits of a particular disease to construct adequate predictive models. Accordingly, radiological characteristics of the tissue are supposed to mimic a defined genotype and phenotype and to better explore the heterogeneity and the dynamic evolution of the tumor over the time. Despite such improvements, we are still far from achieving approved and standardized protocols in clinical practice. Nevertheless, what can we learn by this emerging multidisciplinary clinical approach? This minireview provides a focused overview on the significance of radiomics integrated by RNA sequencing in BC. We will also discuss advances and future challenges of such radiomics-based approach.

Genomics applied to the treatment of breast cancer

Breast cancer remains a major health issue in the world with 1.7 million new cases in 2012 worldwide. It is the second cause of death from cancer in western countries. Genomics have started to modify the treatment of breast cancer, and the developments should become more and more significant, especially in the present era of treatment personalization and with the implementation of new technologies. With molecular signatures, genomics enabled a de-escalation of chemotherapy and personalized treatments of localized forms of estrogen-dependent breast cancers. Genomics can also make a real contribution to constitutional genetics, so as to identify mutations in a panel of candidate genes. In this review, we will discuss the contributions of genomics applied to the treatment of breast cancer, whether already validated contributions or possible future applications linked to research data.

Transcriptome profiling revealed multiple genes and ECM-receptor interaction pathways that may be associated with breast cancer

Background

Exploration of the genes with abnormal expression during the development of breast cancer is essential to provide a deeper understanding of the mechanisms involved. Transcriptome sequencing and bioinformatics analysis of invasive ductal carcinoma and paracancerous tissues from the same patient were performed to identify the key genes and signaling pathways related to breast cancer development.

Methods

Samples of breast tumor tissue and paracancerous breast tissue were obtained from 6 patients. Sequencing used the Illumina HiSeq platform. All. Only perfectly matched clean reads were mapped to the reference genome database, further analyzed and annotated based on the reference genome information. Differentially expressed genes (DEGs) were identified using the DESeq R package (1.10.1) and DEGSeq R package (1.12.0). Using KOBAS software to execute the KEGG bioinformatics analyses, enriched signaling pathways of DEGs involved in the occurrence of breast cancer were determined. Subsequently, quantitative real time PCR was used to verify the accuracy of the expression profile of key DEGs from the RNA-seq result and to explore the expression patterns of novel cancer-related genes on 8 different clinical individuals.

Results

The transcriptomic sequencing results showed 937 DEGs, including 487 upregulated and 450 downregulated genes in the breast cancer specimens. Further quantitative gene expression analysis was performed and captured 252 DEGs (201 downregulated and 51 upregulated) that showed the same differential expression pattern in all libraries. Finally, 6 upregulated DEGs (CST2, DRP2, CLEC5A, SCD, KIAA1211, DTL) and 6 downregulated DEGs (STAC2, BTNL9, CA4, CD300LG, GPIHBP1 and PIGR), were confirmed in a quantitative real time PCR comparison of breast cancer and paracancerous breast tissues from 8 clinical specimens. KEGG analysis revealed various pathway changes, including 20 upregulated and 21 downregulated gene enrichment pathways. The extracellular matrix–receptor (ECM-receptor) interaction pathway was the most enriched pathway: all genes in this pathway were DEGs, including the THBS family, collagen and fibronectin. These DEGs and the ECM-receptor interaction pathway may perform important roles in breast cancer.

Conclusion

Several potential breast cancer-related genes and pathways were captured, including 7 novel upregulated genes and 76 novel downregulated genes that were not found in other studies. These genes are related to cell proliferation, movement and adhesion. They may be important for research into breast cancer mechanisms, particularly CST2 and CA4. A key signaling pathway, the ECM-receptor interaction signal pathway, was also identified as possibly involved in the development of breast cancer.

Electronic supplementary material

The online version of this article (10.1186/s11658-019-0162-0) contains supplementary material, which is available to authorized users.

Silencing Tspan1 inhibits migration and invasion, and induces the apoptosis of human pancreatic cancer cells

Pancreatic cancer (PCC) is one of the most dangerous types of tumor as it is very difficult to treat and its 5-year survival rate is <6%. To date, there have been no effective therapeutic strategies to treat PCC, thus, novel effective therapeutic methods are required. Tetraspanin 1 (Tspan1) is a novel member of the tetraspanins superfamily and is highly expressed in a variety of types of cancer, including gastric, hepatocellular and colonic carcinomas. However, the detailed functional role of Tspan1 in pancreatic cancer cells is still unclear and further investigation is required to uncover its therapeutic potential for the treatment of different tumor types. The purpose of the present study was to investigate the expression of Tspan1 in human PCC tissues and cells, and explore the effect of Tspan1 silencing on invasion, migration, cell survival and apoptosis in human PCC to clarify its function. Expression levels of Tspan1 were analyzed in human pancreatic cancer tissues and the cell lines Capan-2 and SW1990 using immunohistochemistry staining, reverse transcription-quantitative polymerase chain reaction and western blotting. The effects of downregulation of Tspan1 expression on cell survival, apoptosis, invasion and migration were investigated viaTspan1-small interfering (si)RNA transfection into human PCC cell lines. The results indicated that Tspan1 expression was increased in human PCC tissues compared with the adjacent normal pancreatic tissues. Tspan1 was highly expressed in the human PCC cell lines Capan-2 and SW1990 when compared with the normal pancreatic cell line HPC-Y5. In addition, transfection with siRNA-targeting Tspan1 significantly reduced cell migration and invasion, and increased the cell apoptosis of Capan-2 and SW1990. The present findings highlighted the important role of Tspan1 in human PCC cell migration, invasion and apoptosis. Thus, Tspan1 RNA interference may serve as a potential therapeutic strategy to treat human PCC.

Decreased TSPAN1 promotes prostate cancer progression and is a marker for early biochemical recurrence after radical prostatectomy

Patients with prostate cancer (PCa) have a variable prognosis. It is challenging to recognize the progressive disease. In this study, we focused on TSPAN1, a new member of the tetraspanin family. Its expression was decreased in progressive PCa and was an independent prognosis factor of biochemical recurrence after radical prostatectomy. In vitro, knockdown and overexpression of TSPAN1 in PCa cell lines showed that TSPAN1 could inhibit cell proliferation and migration. TSPAN1 was positive related to PTEN in both clinical specimen and mouse models. The combination of these two markers could increase their prognosis value especially in low risk patients. In vitro TSPAN1 knockdown resulted in increased Akt phosphorylation and caused evident cell cycle transition from G1 to S phase. Our data suggests that TSPAN1 is a valuable marker to recognize more progressive PCa.

Array comparative genomic hybridization of 18 pancreatic ductal adenocarcinomas and their autologous metastases

BACKGROUND

Mortality rates of pancreatic cancer remain high, which is mainly due to advanced disease and metastasis. We hypothesized that genomic copy number alterations are enriched in metastatic cells compared to autologous primary tumors, which may inform on cancer-related pathways possibly serving as potential targets for specific therapies. We investigated 18 pancreatic ductal adenocarcinomas, including 39 lymph node and 5 distant metastases after surgical resection. Analysis was performed with array-based comparative genomic hybridization (aCGH).

RESULTS

Metastases acquire a higher frequency of copy number alterations with the highest in distant metastasis (median = 42, lymph node metastases: median = 23, primary tumors: median = 17). In lymph node metastases, gains were prevalent on chromosome bands 8q11.23-q24.3, 12q14.1, 17p12.1, 21q22.12, and losses on 3p21.31, 4p14, 8p23.3-p11.21,17p12-11.2. Genes on amplified regions are involved in cancer-related pathways such as WNT-signaling, also involved in metastasis.

CONCLUSIONS

Pancreatic cancers show a high degree of intratumor heterogeneity, which could lead to resistance of chemotherapy and worse outcome. ACGH analysis reveals regions preferentially gained or lost in synchronous metastases encoding for genes involved in cancer-related pathways, which could lead to novel therapeutic opportunities.

The cancer-associated cell migration protein TSPAN1 is under control of androgens and its upregulation increases prostate cancer cell migration

Cell migration drives cell invasion and metastatic progression in prostate cancer and is a major cause of mortality and morbidity. However the mechanisms driving cell migration in prostate cancer patients are not fully understood. We previously identified the cancer-associated cell migration protein Tetraspanin 1 (TSPAN1) as a clinically relevant androgen regulated target in prostate cancer. Here we find that TSPAN1 is acutely induced by androgens, and is significantly upregulated in prostate cancer relative to both normal prostate tissue and benign prostate hyperplasia (BPH). We also show for the first time, that TSPAN1 expression in prostate cancer cells controls the expression of key proteins involved in cell migration. Stable upregulation of TSPAN1 in both DU145 and PC3 cells significantly increased cell migration and induced the expression of the mesenchymal markers SLUG and ARF6. Our data suggest TSPAN1 is an androgen-driven contributor to cell survival and motility in prostate cancer.

Expression of Tspan-1 gene in patients with advanced gastric cancer

The present study investigated the correlations of the Tspan-1 gene expression with the clinical characteristics and survival prognoses of patients with advanced gastric cancer. A total of 150 patients with advanced gastric cancer were enrolled in the present study, of whom 84 were at stage II and 66 were at stage III according to the tumor node metastasis (TNM) staging; the immunohistochemical staining method and the semi-quantitative PCR method were used to detect the positive expression rates and mRNA relative expression levels of Tspan-1, vascular endothelial growth factor (VEGF), E-cadherin and N-cadherin. The positive expression rates of Tspan-1, VEGF, E-cadherin and N-cadherin were 58.0% (87 patients), 50.0% (75 patients), 28.0% (42 patients) and 53.3% (80 patients), respectively. The positive expressions and mRNA levels of Tspan-1, VEGF, E-cadherin and N-cadherin were not correlated with sex or age (P>0.05), but associated with the cancer state (stage II or stage III) and maximum tumor diameter (P<0.05). With the increase of stage and tumor diameter, the positive rates and mRNA levels of Tspan-1, VEGF and N-cadherin were increased, while those of E-cadherin were decreased. Among patients with stage II/III advanced gastric cancer, those with positive expression of Tspan-1, VEGF and N-cadherin had lower median survival time and survival rates than patients with negative expressions, while patients with positive expression of E-cadherin had higher median survival time and survival rate than those with negative expression (P<0.05). The high expression of Tspan-1 gene is associated with the TNM staging of advanced gastric cancer and the tumor diameter, influences the survival prognosis, and may involve the processes of angiogenesis and epithelial-mesenchymal transition.

Increased risk of brain metastases in women with breast cancer and p16 expression in metastatic lymph-nodes

PURPOSE

Metastatic breast cancer is a leading cause of mortality in women, partly on account of brain metastases. However, the mechanisms by which cancer cells cross the blood-brain barrier remain undeciphered. Most molecular studies predicting metastatic risk have been performed on primary breast cancer samples. Here we studied metastatic lymph-nodes from patients with breast cancers to identify markers associated with the occurrence of brain metastases.

RESULTS

Transcriptomic analyses identified CDKN2A/p16 as a gene potentially associated with brain metastases.

MATERIALS AND METHODS

Fifty-two patients with HER2-overexpressing or triple-negative breast carcinoma with lymph nodes and distant metastases were included in this study. Transcriptomic analyses were performed on laser-microdissected tumor cells from 28 metastatic lymph-nodes. Supervised analyses compared the transcriptomic profiles of women who developed brain metastases and those who did not. As a validation series, we studied metastatic lymph-nodes from 24 other patients.Immunohistochemistry investigations showed that p16 mean scores were significantly higher in patients with brain metastases than in patients without (7.4 vs. 1.7 respectively, p < 0.01). This result was confirmed on the validation series. Multivariate analyses showed that the p16 score was the only variable positively associated with the risk of brain metastases (p = 0.01).With the same threshold of 5 for p16 scores using a Cox model, overall survival was shorter in women with a p16 score over 5 in both series.

CONCLUSIONS

The risk of brain metastases in women with HER2-overexpressing or triple-negative breast cancer could be better assessed by studying p16 protein expression on surgically removed axillary lymph-nodes.

Molecular Concordance Between Primary Breast Cancer and Matched Metastases

Clinical management of breast cancer is increasingly personalized and based on molecular profiling. Often, primary tumors are used as proxies for systemic disease at the time of recurrence. However, recent studies have revealed substantial discordances between primary tumors and metastases, both with respect to traditional clinical treatment targets and on the genomic and transcriptomic level. With the increasing use of molecularly targeted therapy, discordance of actionable molecular targets between primary tumors and recurrences can result in nonoptimal treatment or unnecessary side effects. The purpose of this review is to illuminate the extent of cancer genome evolution through disease progression and the degree of molecular concordance between primary breast cancers and matched metastases. We present an overview of the most prominent studies investigating the expression of endocrine receptors, transcriptomics, and genome aberrations in primary tumors and metastases. In conclusion, biopsy of metastatic lesions at recurrence of breast cancer is encouraged to provide optimal treatment of the disease. Furthermore, molecular profiling of metastatic tissue provides invaluable mechanistic insight into the biology underlying metastatic progression and has the potential to identify novel, potentially druggable, drivers of progression.

Integrating genetics and epigenetics in breast cancer: biological insights, experimental, computational methods and therapeutic potential

BACKGROUND

Development of human cancer can proceed through the accumulation of different genetic changes affecting the structure and function of the genome. Combined analyses of molecular data at multiple levels, such as DNA copy-number alteration, mRNA and miRNA expression, can clarify biological functions and pathways deregulated in cancer. The integrative methods that are used to investigate these data involve different fields, including biology, bioinformatics, and statistics.

RESULTS

These methodologies are presented in this review, and their implementation in breast cancer is discussed with a focus on integration strategies. We report current applications, recent studies and interesting results leading to the identification of candidate biomarkers for diagnosis, prognosis, and therapy in breast cancer by using both individual and combined analyses.

CONCLUSION

This review presents a state of art of the role of different technologies in breast cancer based on the integration of genetics and epigenetics, and shares some issues related to the new opportunities and challenges offered by the application of such integrative approaches.

Tetraspanin 1 is involved in survival, proliferation and carcinogenesis of pancreatic cancer

Pancreatic cancer (PCC) is one of the most difficult cancers to treat and the 10th leading cause of cancer-related death in worldwide. Studies have demonstrated that the tetraspanin 1 (Tspan1) is overexpressed in various cancers and may be a potential therapeutic strategy for the treatment of different cancers. However, the possible role of Tspan1 in PCC is still unknown. In the present study, our data revealed that the increased Tspan1 in PCC tissues was associated with the clinicopathological features and survival rate of PCC patient. We also investigated the effects of Tspan1 gene knockdown on the biological behavior of human PCC. The expression of Tspan1 (detected by immunohistochemistry, qRT-PCR and western blot analysis) derived from human PCC tissues and cell lines (AsPC-1 and PANC-1), were significantly elevated compared with those of the control (P<0.05). Transfection with siRNA-targeting Tspan1 significantly decreased proliferation, increased the apoptosis and reduced migration and invasion of AsPC-1 and PANC-1 cells. The present study demonstrated that Tspan1 plays an important role in PCC carcinogenic progression, including migration and invasion. The siRNA targeting of Tspan1 may be a potential therapeutic strategy for the treatment of PCC.

Transcriptional Profiling of Breast Cancer Metastases Identifies Liver Metastasis-Selective Genes Associated with Adverse Outcome in Luminal A Primary Breast Cancer

PURPOSE

The complete molecular basis of the organ-specificity of metastasis is elusive. This study aimed to provide an independent characterization of the transcriptional landscape of breast cancer metastases with the specific objective to identify liver metastasis-selective genes of prognostic importance following primary tumor diagnosis.

EXPERIMENTAL DESIGN

A cohort of 304 women with advanced breast cancer was studied. Associations between the site of recurrence and clinicopathologic features were investigated. Fine-needle aspirates of metastases (n = 91) were subjected to whole-genome transcriptional profiling. Liver metastasis-selective genes were identified by significance analysis of microarray (SAM) analyses and independently validated in external datasets. Finally, the prognostic relevance of the liver metastasis-selective genes in primary breast cancer was tested.

RESULTS

Liver relapse was associated with estrogen receptor (ER) expression (P = 0.002), luminal B subtype (P = 0.01), and was prognostic for an inferior postrelapse survival (P = 0.01). The major variation in the transcriptional landscape of metastases was also associated with ER expression and molecular subtype. However, liver metastases displayed unique transcriptional fingerprints, characterized by downregulation of extracellular matrix (i.e., stromal) genes. Importantly, we identified a 17-gene liver metastasis-selective signature, which was significantly and independently prognostic for shorter relapse-free (P < 0.001) and overall (P = 0.001) survival in ER-positive tumors. Remarkably, this signature remained independently prognostic for shorter relapse-free survival (P = 0.001) among luminal A tumors.

CONCLUSIONS

Extracellular matrix (stromal) genes can be used to partition breast cancer by site of relapse and may be used to further refine prognostication in ER positive primary breast cancer.

Molecular signatures of lymph node status by intrinsic subtype: gene expression analysis of primary breast tumors from patients with and without metastatic lymph nodes

Background

Identification of a gene expression signature in primary breast tumors that could classify patients by lymph node status would allow patients to avoid the morbidities of surgical disruption of the lymph nodes. Attempts to identify such a signature have, to date, been unsuccessful. Because breast tumor subtypes have unique molecular characteristics and different sites of metastasis, molecular signatures for lymph node involvement may vary by subtype.

Methods

Gene expression data was generated from HG U133A 2.0 arrays for 135 node positive and 210 node negative primary breast tumors. Intrinsic subtype was assigned using the BreastPRS. Differential gene expression analysis was performed using one-way ANOVA using lymph node status as the variable with a False-discovery rate <0.05, to define significance.

Results

Luminal A tumors were most common (51%) followed by basal-like (27%), HER2-enriched (14%) luminal B (7%) and normal-like (1%). Basal-like and luminal A tumors were less likely to have metastatic lymph nodes (35% and 37%, respectively) compared to luminal B or HER2-enriched (52% and 51%, respectively). No differentially expressed genes associated with lymph node status were detected when all tumors were considered together or within each subtype.

Conclusions

Gene expression patterns from the primary tumor are not able to stratify patients by lymph node status. Although the primary breast tumor may influence tumor cell dissemination, once metastatic cells enter the lymphatics, it is likely that characteristics of the lymph node microenvironment, such as establishment of a pre-metastatic niche and release of pro-survival factors, determine which cells are able to colonize. The inability to utilize molecular profiles from the primary tumor to determine lymph node status suggest that other avenues of investigation, such as how systemic factors including diminished immune response or genetic susceptibility contribute to metastasis, may be critical in the development of tools for non-surgical assessment of lymph node status with a corresponding reduction in downstream sequelae associated with disruption of the lymphatics.

SCCRO3 (DCUN1D3) antagonizes the neddylation and oncogenic activity of SCCRO (DCUN1D1)

The activity of cullin-RING type ubiquitination E3 ligases is regulated by neddylation, a process analogous to ubiquitination that culminates in covalent attachment of the ubiquitin-like protein Nedd8 to cullins. As a component of the E3 for neddylation, SCCRO/DCUN1D1 plays a key regulatory role in neddylation and, consequently, cullin-RING ligase activity. The essential contribution of SCCRO to neddylation is to promote nuclear translocation of the cullin-ROC1 complex. The presence of a myristoyl sequence in SCCRO3, one of four SCCRO paralogues present in humans that localizes to the membrane, raises questions about its function in neddylation. We found that although SCCRO3 binds to CAND1, cullins, and ROC1, it does not efficiently bind to Ubc12, promote cullin neddylation, or conform to the reaction processivity paradigms, suggesting that SCCRO3 does not have E3 activity. Expression of SCCRO3 inhibits SCCRO-promoted neddylation by sequestering cullins to the membrane, thereby blocking its nuclear translocation. Moreover, SCCRO3 inhibits SCCRO transforming activity. The inhibitory effects of SCCRO3 on SCCRO-promoted neddylation and transformation require both an intact myristoyl sequence and PONY domain, confirming that membrane localization and binding to cullins are required for in vivo functions. Taken together, our findings suggest that SCCRO3 functions as a tumor suppressor by antagonizing the neddylation activity of SCCRO.

Lack of genomic heterogeneity at high-resolution aCGH between primary breast cancers and their paired lymph node metastases

Lymph-node metastasis (LNM) predict high recurrence rates in breast cancer patients. Systemic treatment aims to eliminate (micro)metastatic cells. However decisions regarding systemic treatment depend largely on clinical and molecular characteristics of primary tumours. It remains, however, unclear to what extent metastases resemble the cognate primary breast tumours, especially on a genomic level, and as such will be eradicated by the systemic therapy chosen. In this study we used high-resolution aCGH to investigate DNA copy number differences between primary breast cancers and their paired LNMs. To date, no recurrent LNM-specific genomic aberrations have been identified using array comparative genomic hybridization (aCGH) analysis. In our study we employ a high-resolution platform and we stratify on different breast cancer subtypes, both aspects that might have underpowered previously performed studies.To test the possibility that genomic instability in triple-negative breast cancers (TNBCs) might cause increased random and potentially also recurrent copy number aberrations (CNAs) in their LNMs, we studied 10 primary TNBC–LNM pairs and 10 ER-positive (ER+) pairs and verified our findings adding additionally 5 TNBC-LNM and 22 ER+-LNM pairs. We found that all LNMs clustered nearest to their matched tumour except for two cases, of which one was due to the presence of two distinct histological components in one tumour. We found no significantly altered CNAs between tumour and their LNMs in the entire group or in the subgroups. Within the TNBC subgroup, no absolute increase in CNAs was found in the LNMs compared to their primary tumours, suggesting that increased genomic instability does not lead to more CNAs in LNMs. Our findings suggest a high clonal relationship between primary breast tumours and its LNMs, at least prior to treatment, and support the use of primary tumour characteristics to guide adjuvant systemic chemotherapy in breast cancer patients.

Up-regulation of vitamin B1 homeostasis genes in breast cancer

An increased carbon flux and exploitation of metabolic pathways for the rapid generation of biosynthetic precursors is a common phenotype observed in breast cancer. To support this metabolic phenotype, cancer cells adaptively regulate the expression of glycolytic enzymes and nutrient transporters. However, activity of several enzymes involved in glucose metabolism requires an adequate supply of cofactors. In particular, vitamin B1 (thiamine) is utilized as an essential cofactor for metabolic enzymes that intersect at critical junctions within the glycolytic network. Intracellular availability of thiamine is facilitated by the activity of thiamine transporters and thiamine pyrophosphokinase-1 (TPK-1). Therefore, the objective of this study was to establish if the cellular determinants regulating thiamine homeostasis differ between breast cancer and normal breast epithelia. Employing cDNA arrays of breast cancer and normal breast epithelial tissues, SLC19A2, SLC25A19 and TPK-1 were found to be significantly up-regulated. Similarly, up-regulation was also observed in breast cancer cell lines compared to human mammary epithelial cells. Thiamine transport assays and quantitation of intracellular thiamine and thiamine pyrophosphate established a significantly greater extent of thiamine transport and free thiamine levels in breast cancer cell lines compared to human mammary epithelial cells. Overall, these findings demonstrate an adaptive response by breast cancer cells to increase cellular availability of thiamine.

Mapping of three genetic determinants of susceptibility to estrogen-induced mammary cancer within the Emca8 locus on rat chromosome 5

The ACI rat model of 17β-estradiol (E2)-induced mammary cancer has gained wide use in the study of breast cancer etiology, prevention, and genetics. Emca8, a QTL that determines susceptibility to E2-induced mammary cancer, was previously mapped to rat chromosome 5 (RNO5) in an intercross between resistant Brown Norway (BN) and susceptible ACI rats. In this study, a panel of congenic rat strains, each of which carries BN alleles across a defined segment of RNO5 on the ACI genetic background, was generated and used to map more precisely the Emca8 determinants of mammary cancer susceptibility. Three distinct genetic determinants were localized within Emca8, and two of these were mapped to intervals of less than 15 megabases. Emca8.1 harbors Cdkn2a, Cdkn2b, and other genes and is orthologous to the 9p21 breast cancer locus identified in genome-wide and candidate gene association studies. Emca8.2 harbors Cdkn2c and other genes and is orthologous to the 1p32 locus in humans that is frequently deleted in breast cancers. Both Emca8.1 and Emca8.2 harbor copy number variants that are orthologous to copy number variant regions in humans. Gene expression profiles were defined for mammary tissues from E2-treated ACI and ACI.BN-Emca8 rats to define the impact of Emca8 on gene expression and identify differentially expressed genes residing within Emca8.1 and Emca8.2. This study further illustrates the relevance of the ACI rat model of E2-induced mammary cancer for identifying novel genetic determinants of breast cancer susceptibility and defining the mechanisms through which estrogens contribute to breast cancer development. Cancer Prev Res; 6(1); 59-69. ©2012 AACR.

Differential copy number aberrations in novel candidate genes associated with progression from in situ to invasive ductal carcinoma of the breast

Only a minority of intraductal carcinomas of the breast give rise to stromally invasive disease. We microdissected 206 paraffin blocks representing 116 different cases of low-grade ductal carcinoma in situ (DCIS). Fifty-five were pure DCIS (PD) cases without progression to invasive carcinoma. Sixty-one cases had a small invasive component. DNA was extracted from microdissected sections and hybridized to high-density bacterial artificial chromosome arrays. Array comparative genomic hybridization analysis of 118 hybridized DNA samples yielded data on 69 samples that were suitable for further statistical analysis. This cohort included 20 pure DCIS cases, 25 mixed DCIS (MD), and 24 mixed invasive carcinoma samples. PD cases had a higher frequency of DNA copy number changes than MD cases, and the latter had similar DNA profiles compared to paired invasive carcinomas. Copy number changes on 13 chromosomal arms occurred at different rates in PD versus MD lesions. Eight of 19 candidate genes residing at those loci were confirmed to have differential copy number changes by quantitative PCR. NCOR2/SMRT and NR4A1 (both on 12q), DYNLRB2 (16q), CELSR1, UPK3A, and ST13 (all on 22q) were more frequently amplified in PD. Moreover, NCOR2, NR4A1, and DYNLRB2 showed more frequent copy number losses in MD. GRAP2 (22q) was more often amplified in MD, whereas TAF1C (16q) was more commonly deleted in PD. A multigene model comprising these candidate genes discriminated between PD and MD lesions with high accuracy. These findings suggest that the propensity to invade the stroma may be encoded in the genome of intraductal carcinomas.