Design, construction, and analysis of cell line arrays and tissue microarrays for gene expression analysis

Cell line array (CMA) and tissue microarray (TMA) technologies are high-throughput methods for analysing both the abundance and distribution of gene expression in a panel of cell lines or multiple tissue specimens in an efficient and cost-effective manner. The process is based on Kononen's method of extracting a cylindrical core of paraffin-embedded donor tissue and inserting it into a recipient paraffin block. Donor tissue from surgically resected paraffin-embedded tissue blocks, frozen needle biopsies or cell line pellets can all be arrayed in the recipient block. The representative area of interest is identified and circled on a haematoxylin and eosin (H&E)-stained section of the donor block. Using a predesigned map showing a precise spacing pattern, a high density array of up to 1,000 cores of cell pellets and/or donor tissue can be embedded into the recipient block using a tissue arrayer from Beecher Instruments. Depending on the depth of the cell line/tissue removed from the donor block 100-300 consecutive sections can be cut from each CMA/TMA block. Sections can be stained for in situ detection of protein, DNA or RNA targets using immunohistochemistry (IHC), fluorescent in situ hybridisation (FISH) or mRNA in situ hybridisation (RNA-ISH), respectively. This chapter provides detailed methods for CMA/TMA design, construction and analysis with in-depth notes on all technical aspects including tips to deal with common pitfalls the user may encounter. © Springer Science+Business Media, LLC 2011.

Oxidative stress induced lung cancer and COPD : opportunities for epigenetic therapy

Reactive oxygen species (ROS) form as a natural by-product of the normal metabolism of oxygen and play important roles within the cell. Under normal circumstances the cell is able to maintain an adequate homeostasis between the formation of ROS and its removal through particular enzymatic pathways or via antioxidants. If however, this balance is disturbed a situation called oxidative stress occurs. Critically, oxidative stress plays important roles in the pathogenesis of many diseases, including cancer. Epigenetics is a process where gene expression is regulated by heritable mechanisms that do not cause any direct changes to the DNA sequence itself, and disruption of epigenetic mechanisms has important implications in disease. Evidence is emerging that histone deacetylases (HDACs) play decisive roles in regulating important cellular oxidative stress pathways including those involved with sensing oxidative stress and those involved with regulating the cellular response to oxidative stress. In particular aberrant regulation of these pathways by HDACs may play critical roles in cancer progression. In this review we discuss the current evidence linking epigenetics and oxidative stress and cancer, using chronic obstructive pulmonary disease and non-small cell lung cancer to illustrate the importance of epigenetics on these pathways within these disease settings. © 2009 Foundation for Cellular and Molecular Medicine/Blackwell Publishing Ltd.

Expression and prognostic significance of a panel of tissue hypoxia markers in head-and-neck squamous cell carcinomas

Purpose: To investigate the expression pattern of hypoxia-induced proteins identified as being involved in malignant progression of head-and-neck squamous cell carcinoma (HNSCC) and to determine their relationship to tumor pO 2 and prognosis. Methods and Materials: We performed immunohistochemical staining of hypoxia-induced proteins (carbonic anhydrase IX [CA IX], BNIP3L, connective tissue growth factor, osteopontin, ephrin A1, hypoxia inducible gene-2, dihydrofolate reductase, galectin-1, IκB kinase β, and lysyl oxidase) on tumor tissue arrays of 101 HNSCC patients with pretreatment pO 2 measurements. Analysis of variance and Fisher's exact tests were used to evaluate the relationship between marker expression, tumor pO 2, and CA IX staining. Cox proportional hazard model and log-rank tests were used to determine the relationship between markers and prognosis. Results: Osteopontin expression correlated with tumor pO 2 (Eppendorf measurements) (p = 0.04). However, there was a strong correlation between lysyl oxidase, ephrin A1, and galectin-1 and CA IX staining. These markers also predicted for cancer-specific survival and overall survival on univariate analysis. A hypoxia score of 0-5 was assigned to each patient, on the basis of the presence of strong staining for these markers, whereby a higher score signifies increased marker expression. On multivariate analysis, increasing hypoxia score was an independent prognostic factor for cancer-specific survival (p = 0.015) and was borderline significant for overall survival (p = 0.057) when adjusted for other independent predictors of outcomes (hemoglobin and age). Conclusions: We identified a panel of hypoxia-related tissue markers that correlates with treatment outcomes in HNSCC. Validation of these markers will be needed to determine their utility in identifying patients for hypoxia-targeted therapy. © 2007 Elsevier Inc. All rights reserved.

Gene expression analysis of diagnostic biopsies predicts pathological response to neoadjuvant chemoradiotherapy of esophageal cancer

OBJECTIVE: This study explored gene expression differences in predicting response to chemoradiotherapy in esophageal cancer. PURPOSE:: A major pathological response to neoadjuvant chemoradiation is observed in about 40% of esophageal cancer patients and is associated with favorable outcomes. However, patients with tumors of similar histology, differentiation, and stage can have vastly different responses to the same neoadjuvant therapy. This dichotomy may be due to differences in the molecular genetic environment of the tumor cells. BACKGROUND DATA: Diagnostic biopsies were obtained from a training cohort of esophageal cancer patients (13), and extracted RNA was hybridized to genome expression microarrays. The resulting gene expression data was verified by qRT-PCR. In a larger, independent validation cohort (27), we examined differential gene expression by qRT-PCR. The ability of differentially-regulated genes to predict response to therapy was assessed in a multivariate leave-one-out cross-validation model. RESULTS: Although 411 genes were differentially expressed between normal and tumor tissue, only 103 genes were altered between responder and non-responder tumor; and 67 genes differentially expressed >2-fold. These included genes previously reported in esophageal cancer and a number of novel genes. In the validation cohort, 8 of 12 selected genes were significantly different between the response groups. In the predictive model, 5 of 8 genes could predict response to therapy with 95% accuracy in a subset (74%) of patients. CONCLUSIONS: This study has identified a gene microarray pattern and a set of genes associated with response to neoadjuvant chemoradiation in esophageal cancer. The potential of these genes as biomarkers of response to treatment warrants further investigation. Copyright © 2009 by Lippincott Williams & Wilkins.

Combination therapy with gefitinib and rofecoxib in patients with platinum-pretreated relapsed non-small-cell lung cancer

Purpose: In non-small-cell lung cancer (NSCLC), the epidermal growth factor receptor (EGFR) and cyclooxygenase-2 (COX-2) play major roles in tumorigenesis. This phase I/II study evaluated combined therapy with the EGFR tyrosine kinase inhibitor (TKI) gefitinib and the COX-2 inhibitor rofecoxib in platinum-pretreated, relapsed, metastatic NSCLC (n = 45). Patients and Methods: Gefitinib 250 mg/d was combined with rofecoxib (dose escalated from 12.5 to 25 to 50 mg/d through three cohorts, each n = 6). Because the rofecoxib maximum-tolerated dose was not reached, the 50 mg/d cohort was expanded for efficacy evaluation (n = 33). Results: Among the 42 assessable patients, there was one complete response (CR) and two partial responses (PRs) and 12 patients with stable disease (SD); disease control rate was 35.7% (95% CI, 21.6% to 52.0%). Median time to tumor progression was 55 days (95% CI, 47 to 70 days), and median survival was 144 days (95% CI, 103 to 190 days). In a pilot study, matrix-assisted laser desorption/ionization (MALDI) proteomics analysis of baseline serum samples could distinguish patients with an objective response from those with SD or progressive disease (PD), and those with disease control (CR, PR, and SD) from those with PD. The regimen was generally well tolerated, with predictable toxicities including skin rash and diarrhea. Conclusion: Gefitinib combined with rofecoxib provided disease control equivalent to that expected with single-agent gefitinib and was generally well tolerated. Baseline serum proteomics may help identify those patients most likely to benefit from EGFR TKIs. © 2007 by American Society of Clinical Oncology.

Randomized phase II study of cetuximab plus cisplatin/vinorelbine compared with cisplatin/vinorelbine alone as first-line therapy in EGFR-expressing advanced non-small-cell lung cancer

Background: The Lung Cancer Cetuximab Study is an open-label, randomized phase II pilot study of cisplatin and vinorelbine combined with the epidermal growth factor receptor (EGFR)-targeted monoclonal antibody cetuximab versus cisplatin and vinorelbine alone, in patients with advanced EGFR-expressing, non-small-cell lung cancer (NSCLC). End points of the study are activity, safety and pharmacokinetics. Patients and methods: Following randomization, for a maximum of eight cycles, patients received three-weekly cycles of cisplatin (80 mg/m2, day 1) and vinorelbine (25 mg/m2 on days 1 and 8) alone or following cetuximab treatment (initial dose 400 mg/m, followed by 250 mg/m2 weekly thereafter). Results: Eighty-six patients were randomly allocated to the study (43 per arm). Confirmed response rates were 28% in the cisplatin/vinorelbine arm (A) and 35% in the cetuximab plus cisplatin/vinorelbine arm (B). Median progression-free survival (PFS) was 4.6 months in arm A and 5.0 months in arm B, with PFS rates at 12 months of 0% and 15%, respectively. Median survival was 7.3 months in arm A and 8.3 months in arm B. The 24-month survival rates were 0% and 16%, respectively. The cetuximab combination was well tolerated. Conclusion: In the first-line treatment of advanced NSCLC, the combination of cetuximab plus cisplatin/vinorelbine demonstrated an acceptable safety profile and the potential to improve activity over cisplatin/vinorelbine alone. © 2007 European Society for Medical Oncology.

Transparency tools in gene patenting for informing policy and practice

In the recent decision Association for Molecular Pathology v. Myriad Genetics1, the US Supreme Court held that naturally occurring sequences from human genomic DNA are not patentable subject matter. Only certain complementary DNAs (cDNA), modified sequences and methods to use sequences are potentially patentable. It is likely that this distinction will hold for all DNA sequences, whether animal, plant or microbial2. However, it is not clear whether this means that other naturally occurring informational molecules, such as polypeptides (proteins) or polysaccharides, will also be excluded from patents. The decision underscores a pressing need for precise analysis of patents that disclose and reference genetic sequences, especially in the claims. Similarly, data sets, standards compliance and analytical tools must be improved—in particular, data sets and analytical tools must be made openly accessible—in order to provide a basis for effective decision making and policy setting to support biological innovation. Here, we present a web-based platform that allows such data aggregation, analysis and visualization in an open, shareable facility. To demonstrate the potential for the extension of this platform to global patent jurisdictions, we discuss the results of a global survey of patent offices that shows that much progress is still needed in making these data freely available for aggregation in the first place.

Phase II trial of Oxaliplatin and 5-Fluorouracil/Leucovorin combination in epithelial ovarian carcinoma relapsing within 2 years of platinum-based therapy

Objective To evaluate the efficacy and toxicity of Oxaliplatin and 5-Fluorouracil (5-FU)/Leucovorin (LV) combination in ovarian cancer relapsing within 2 years of prior platinum-based chemotherapy in a phase II trial. Methods Eligible patients had at least one prior platinum-based chemotherapy regimen, elevated CA-125 ≥ 60 IU/l, radiological evidence of disease progression and adequate hepatic, renal and bone marrow function. Patients with raised CA-125 levels alone as marker of disease relapse were not eligible. Oxaliplatin (85 mg/m 2) was given on day 1, and 5-Fluorouracil (370 mg/m 2) and Leucovorin (30 mg) was given on days 1 and 8 of a 14-day cycle. Results Twenty-seven patients were enrolled. The median age was 57 years (range 42-74 years). The median platinum-free interval (PFI) was 5 months (range 0-17 months) with only 30% of patients being platinum sensitive (PFI > 6 months). Six patients (22%) had two prior regimens of chemotherapy. A total of 191 cycles were administered (median 7; range 2-12). All patients were evaluable for toxicity. The following grade 3/4 toxicities were noted: anemia 4%; neutropenia 15%; thrombocytopenia 11%; neurotoxicity 8%; lethargy 4%; diarrhea 4%; hypokalemia 11%; hypomagnesemia 11%. Among 27 enrolled patients, 20 patients were evaluable for response by WHO criteria and 25 patients were evaluable by Rustin's CA-125 criteria. The overall response rate (RR) by WHO criteria was 30% (95% CI: 15- 52) [three complete responses (CRs) and three partial responses (PRs)]. The CA-125 response rate was 56% (95% CI: 37-73). Significantly, a 25% (95% CI: 9-53) radiological and a 50% (95% CI: 28-72) CA-125 response rate were noted in platinum resistant patients (PFI < 6 months). The median response duration was 4 months (range 3-12) and the median overall survival was 10 months. Conclusion Oxaliplatin and 5-Fluorouracil/ Leucovorin combination has a good safety profile and is active in platinum-pretreated advanced epithelial ovarian cancer. © 2004 Elsevier Inc. All rights reserved.

Targeting amino acid transport in metastatic castration-resistant prostate cancer : effects on cell cycle, cell growth, and tumor development

Background L-type amino acid transporters (LATs) uptake neutral amino acids including L-leucine into cells, stimulating mammalian target of rapamycin complex 1 signaling and protein synthesis. LAT1 and LAT3 are overexpressed at different stages of prostate cancer, and they are responsible for increasing nutrients and stimulating cell growth. Methods We examined LAT3 protein expression in human prostate cancer tissue microarrays. LAT function was inhibited using a leucine analog (BCH) in androgen-dependent and -independent environments, with gene expression analyzed by microarray. A PC-3 xenograft mouse model was used to study the effects of inhibiting LAT1 and LAT3 expression. Results were analyzed with the Mann-Whitney U or Fisher exact tests. All statistical tests were two-sided. Results LAT3 protein was expressed at all stages of prostate cancer, with a statistically significant decrease in expression after 4–7 months of neoadjuvant hormone therapy (4–7 month mean = 1.571; 95% confidence interval = 1.155 to 1.987 vs 0 month = 2.098; 95% confidence interval = 1.962 to 2.235; P = .0187). Inhibition of LAT function led to activating transcription factor 4–mediated upregulation of amino acid transporters including ASCT1, ASCT2, and 4F2hc, all of which were also regulated via the androgen receptor. LAT inhibition suppressed M-phase cell cycle genes regulated by E2F family transcription factors including critical castration-resistant prostate cancer regulatory genes UBE2C, CDC20, and CDK1. In silico analysis of BCH-downregulated genes showed that 90.9% are statistically significantly upregulated in metastatic castration-resistant prostate cancer. Finally, LAT1 or LAT3 knockdown in xenografts inhibited tumor growth, cell cycle progression, and spontaneous metastasis in vivo. Conclusion Inhibition of LAT transporters may provide a novel therapeutic target in metastatic castration-resistant prostate cancer, via suppression of mammalian target of rapamycin complex 1 activity and M-phase cell cycle genes.

Association between xenobiotic gene polymorphisms and non-Hodgkin's lymphoma risk

The last four decades have seen a significant increase in the incidence of non-Hodgkin's lymphoma (NHL) as a possible result of increasing environmental carcinogen exposure, particularly pesticides and solvents. Based on the increasing evidence for an association between carcinogen exposure-related cancer risk and xenobiotic gene polymorphisms, we have undertaken a case-control study of xenobiotic gene polymorphisms in individuals with a diagnosis of NHL. Polymorphisms of six xenobiotic genes (CYP1A1, GSTT1, GSTM1, PON1, NAT1, NAT2) were characterized in 169 individuals with NHL and 205 normal controls using polymerase chain reaction-based methods. Polymorphic frequencies were compared using Fisher's exact tests, and odds ratios for NHL risk were calculated. Among the NHL group, the incidence of GSTT1 null and PON1 BB genotypes were significantly increased compared with controls, 34% vs 14%, and 24% vs 11% respectively. Adjusted odds ratios calculated from multivariate analyses demonstrated that GSTT1 null conferred a fourfold increase in NHL risk (OR = 4.27; 95% CI, 2.40-7.61, P < 0.001) and PON1 BB a 2.9-fold increase (OR = 2.92; 95% CI, 1.49-5.72, P = 0.002). Furthermore, GSTT1 null combined with PON1 BB or GSTM1 null conferred an additional risk of NHL. This is the first time that a PON1 gene polymorphism has been shown to be associated with cancer risk. We conclude that the two polymorphisms, GSTT1 null and PON1 BB, are common genetic traits that pose low individual risk but may be important determinants of overall population NHL risk, particularly among groups exposed to NHL-related carcinogens.

Different in vitro activation methods for latent transforming growth factors (TGF)–β : considerable exogenous factors to promote higher mesenchymal-origin cell proliferation in a bioprocessing platform

Regenerative medicine includes two efficient techniques, namely tissue-engineering and cell-based therapy in order to repair tissue damage efficiently. Most importantly, huge numbers of autologous cells are required to deal these practices. Nevertheless, primary cells, from autologous tissue, grow very slowly while culturing in vitro; moreover, they lose their natural characteristics over prolonged culturing period. Transforming growth factors-beta (TGF-β) is a ubiquitous protein found biologically in its latent form, which prevents it from eliciting a response until conversion to its active form. In active form, TGF-β acts as a proliferative agent in many cell lines of mesenchymal origin in vitro. This article reviews on some of the important activation methods-physiochemical, enzyme-mediated, non-specific protein interaction mediated, and drug-induced- of TGF-β, which may be established as exogenous factors to be used in culturing medium to obtain extensive proliferation of primary cells.

Investigation of brain derived neurotrophic factor (BDNF) gene variants in migraine

A number of observations have suggested that brain derived neurotrophic factor (BDNF) plays a role in migraine pathophysiology. This study investigates whether variants in the BDNF gene are associated with migraine in an Australian case-control population. Background. Brain derived neurotrophic factor (BDNF) has an important role in neural growth, development and survival in the central nervous system and is an important modulator of central and peripheral pain responses. Variants in BDNF, in particular the functional Val66Met polymorphism (rs6265), have been found to be associated with a number of psychiatric disorders, cognitive function and obesity. As BDNF has been found to be differentially expressed in a number of aspects related to migraine, we tested for association between single nucleotide polymorphisms (SNPs) in BDNF and migraine. Methods. Five SNPs in the BDNF locus (rs1519480, rs6265, rs712507, rs2049046 and rs12273363) were genotyped initially in a cohort of 277 migraine cases, including 172 diagnosed with migraine with aura (MA) and 105 with migraine without aura (MO), and 277 age- and sex-matched controls. Three of these SNPs (rs6265, rs2049046 and rs12273363) were subsequently genotyped in a second cohort of 580 migraineurs, including 473 diagnosed with MA and 105 with O, and 580 matched controls. Results. – BDNF SNPs rs1519480, rs6265, rs712507 and rs12273363 were not significantly associated with migraine. However, rs2049046 showed a significant association with migraine, and in particular, MA in the first cohort. In the second cohort, although an increase in the rs2049046 T-allele frequency was observed in migraine cases, and in both MA and MO subgroups, it was not significantly different from controls. Analysis of data combined from both cohorts for rs2049046 showed significant differences in the genotypic and allelic distributions for this marker in both migraine and the MA sub-group. Conclusion. This study confirmed previous studies that the functional BDNF SNP rs6265 (Val66Met) is not associated with migraine. However, we found that rs2049046, which resides at the 5’ end of 3 one the BDNF transcripts, may be associated with migraine, suggesting that further investigations of this SNP may be warranted.

The relationship between BCOM1 gene variants and macular pigment optical density in persons with and without age-related macular degeneration

Background: Recent evidence indicates that gene variants related to carotenoid metabolism play a role in the uptake of macular pigments lutein (L) and zeaxanthine (Z). Moreover, these pigments are proposed to reduce the risk for advanced age-related macular degeneration (AMD). This study provides the initial examination of the relationship between the gene variants related to carotenoid metabolism, macular pigment optical density (MPOD) and their combined expression in healthy humans and patients with AMD. Participants and Methods: Forty-four participants were enrolled from a general population and a private practice including 20 healthy participants and 24 patients with advanced (neovascular) AMD. Participants were genotyped for the three single nucleotide polymorphisms (SNPs) upstream from BCMO1, rs11645428, rs6420424 and rs6564851 that have been shown to either up or down regulate beta-carotene conversion efficiency in the plasma. MPOD was determined by heterochromatic flicker photometry. Results: Healthy participants with the rs11645428 GG genotype, rs6420424 AA genotype and rs6564851 GG genotype all had on average significantly lower MPOD compared to those with the other genotypes (p < 0.01 for all three comparisons). When combining BCMO1 genotypes reported to have “high” (rs11645428 AA/rs6420424 GG/rs6564851 TT) and “low” (rs11645428 GG/rs6420424 AA/rs6564851 GG) beta-carotene conversion efficiency, we demonstrate clear differences in MPOD values (p<0.01). In patients with AMD there were no significant differences in MPOD for any of the three BCMO1 gene variants. Conclusion: In healthy participants MPOD levels can be related to high and low beta-carotene conversion BCMO1 genotypes. Such relationships were not found in patients with advanced neovascular AMD, indicative of additional processes influencing carotenoid uptake, possibly related to other AMD susceptibility genes. Our findings indicate that specific BCMO1 SNPs should be determined when assessing the effects of carotenoid supplementation on macular pigment and that their expression may be influenced by retinal disease.

Direct repression of MYB by ZEB1 suppresses proliferation and epithelial gene expression during epithelial-to-mesenchymal transition of breast cancer cells

Introduction Epithelial-to-mesenchymal transition (EMT) promotes cell migration and is important in metastasis. Cellular proliferation is often downregulated during EMT, and the reverse transition (MET) in metastases appears to be required for restoration of proliferation in secondary tumors. We studied the interplay between EMT and proliferation control by MYB in breast cancer cells. Methods MYB, ZEB1, and CDH1 expression levels were manipulated by lentiviral small-hairpin RNA (shRNA)-mediated knockdown/overexpression, and verified with Western blotting, immunocytochemistry, and qRT-PCR. Proliferation was assessed with bromodeoxyuridine pulse labeling and flow cytometry, and sulforhodamine B assays. EMT was induced with epidermal growth factor for 9 days or by exposure to hypoxia (1% oxygen) for up to 5 days, and assessed with qRT-PCR, cell morphology, and colony morphology. Protein expression in human breast cancers was assessed with immunohistochemistry. ZEB1-MYB promoter binding and repression were determined with Chromatin Immunoprecipitation Assay and a luciferase reporter assay, respectively. Student paired t tests, Mann–Whitney, and repeated measures two-way ANOVA tests determined statistical significance (P < 0.05). Results Parental PMC42-ET cells displayed higher expression of ZEB1 and lower expression of MYB than did the PMC42-LA epithelial variant. Knockdown of ZEB1 in PMC42-ET and MDA-MB-231 cells caused increased expression of MYB and a transition to a more epithelial phenotype, which in PMC42-ET cells was coupled with increased proliferation. Indeed, we observed an inverse relation between MYB and ZEB1 expression in two in vitro EMT cell models, in matched human breast tumors and lymph node metastases, and in human breast cancer cell lines. Knockdown of MYB in PMC42-LA cells (MYBsh-LA) led to morphologic changes and protein expression consistent with an EMT. ZEB1 expression was raised in MYBsh-LA cells and significantly repressed in MYB-overexpressing MDA-MB-231 cells, which also showed reduced random migration and a shift from mesenchymal to epithelial colony morphology in two dimensional monolayer cultures. Finally, we detected binding of ZEB1 to MYB promoter in PMC42-ET cells, and ZEB1 overexpression repressed MYB promoter activity. Conclusions This work identifies ZEB1 as a transcriptional repressor of MYB and suggests a reciprocal MYB-ZEB1 repressive relation, providing a mechanism through which proliferation and the epithelial phenotype may be coordinately modulated in breast cancer cells.

Assessment of gene expression of intracellular calcium channels, pumps and exchangers with epidermal growth factor-induced epithelial-mesenchymal transition in a breast cancer cell line

Background Epithelial-mesenchymal transition (EMT) is a process implicated in cancer metastasis that involves the conversion of epithelial cells to a more mesenchymal and invasive cell phenotype. In breast cancer cells EMT is associated with altered store-operated calcium influx and changes in calcium signalling mediated by activation of cell surface purinergic receptors. In this study, we investigated whether MDA-MB-468 breast cancer cells induced to undergo EMT exhibit changes in mRNA levels of calcium channels, pumps and exchangers located on intracellular calcium storing organelles, including the Golgi, mitochondria and endoplasmic reticulum (ER). Methods Epidermal growth factor (EGF) was used to induce EMT in MDA-MB-468 breast cancer cells. Serum-deprived cells were treated with EGF (50 ng/mL) for 12 h and gene expression was assessed using quantitative RT-PCR. Results and conclusions These data reveal no significant alterations in mRNA levels of the Golgi calcium pump secretory pathway calcium ATPases (SPCA1 and SPCA2), or the mitochondrial calcium uniporter (MCU) or Na+/Ca2+ exchanger (NCLX). However, EGF-induced EMT was associated with significant alterations in mRNA levels of specific ER calcium channels and pumps, including (sarco)-endoplasmic reticulum calcium ATPases (SERCAs), and inositol 1,4,5-trisphosphate receptor (IP3R) and ryanodine receptor (RYR) calcium channel isoforms. The most prominent change in gene expression between the epithelial and mesenchymal-like states was RYR2, which was enriched 45-fold in EGF-treated MDA-MB-468 cells. These findings indicate that EGF-induced EMT in breast cancer cells may be associated with major alterations in ER calcium homeostasis.