Erratum to: Modeling precision treatment of breast cancer

During the type-setting of the final version of the article [1] some of the additional files were swapped. The correct files are republished in this Erratum.

ATM suppresses SATB1-induced malignant progression in breast epithelial cells

SATB1 drives metastasis when expressed in breast tumor cells by radically reprogramming gene expression. Here, we show that SATB1 also has an oncogenic activity to transform certain non-malignant breast epithelial cell lines. We studied the non-malignant MCF10A cell line, which is used widely in the literature. We obtained aliquots from two different sources (here we refer to them as MCF10A-1 and MCF10A-2), but found them to be surprisingly dissimilar in their responses to oncogenic activity of SATB1. Ectopic expression of SATB1 in MCF10A-1 induced tumor-like morphology in three-dimensional cultures, led to tumor formation in immunocompromised mice, and when injected into tail veins, led to lung metastasis. The number of metastases correlated positively with the level of SATB1 expression. In contrast, SATB1 expression in MCF10A-2 did not lead to any of these outcomes. Yet DNA copy-number analysis revealed that MCF10A-1 is indistinguishable genetically from MCF10A-2. However, gene expression profiling analysis revealed that these cell lines have significantly divergent signatures for the expression of genes involved in oncogenesis, including cell cycle regulation and signal transduction. Above all, the early DNA damage-response kinase, ATM, was greatly reduced in MCF10A-1 cells compared to MCF10A-2 cells. We found the reason for reduction to be phenotypic drift due to long-term cultivation of MCF10A. ATM knockdown in MCF10A-2 and two other non-malignant breast epithelial cell lines, 184A1 and 184B4, enabled SATB1 to induce malignant phenotypes similar to that observed for MCF10A-1. These data indicate a novel role for ATM as a suppressor of SATB1-induced malignancy in breast epithelial cells, but also raise a cautionary note that phenotypic drift could lead to dramatically different functional outcomes.

Cross-platform pathway-based analysis identifies markers of response to the PARP inhibitor olaparib

Poly(ADP-ribose) polymerase (PARP) is an enzyme involved in DNA repair. PARP inhibitors can act as chemosensitizers, or operate on the principle of synthetic lethality when used as single agent. Clinical trials have shown drugs in this class to be promising for BRCA mutation carriers. We postulated that inability to demonstrate response in non-BRCA carriers in which BRCA is inactivated by other mechanisms or with deficiency in homologous recombination for DNA repair is due to lack of molecular markers that define a responding subpopulation. We identified candidate markers for this purpose for olaparib (AstraZeneca) by measuring inhibitory effects of nine concentrations of olaparib in 22 breast cancer cell lines and identifying features in transcriptional and genome copy number profiles that were significantly correlated with response. We emphasized in this discovery process genes involved in DNA repair. We found that the cell lines that were sensitive to olaparib had a significant lower copy number of BRCA1 compared to the resistant cell lines (p value 0.012). In addition, we discovered seven genes from DNA repair pathways whose transcriptional levels were associated with response. These included five genes (BRCA1, MRE11A, NBS1, TDG, and XPA) whose transcript levels were associated with resistance and two genes (CHEK2 and MK2) whose transcript levels were associated with sensitivity. We developed an algorithm to predict response using the seven-gene transcription levels and applied it to 1,846 invasive breast cancer samples from 8 U133A/plus 2 (Affymetrix) data sets and found that 8–21 % of patients would be predicted to be responsive to olaparib. A similar response frequency was predicted in 536 samples analyzed on an Agilent platform. Importantly, tumors predicted to respond were enriched in basal subtype tumors. Our studies support clinical evaluation of the utility of our seven-gene signature as a predictor of response to olaparib.

Temporal dissection of tumorigenesis in primary cancers

Timely intervention for cancer requires knowledge of its earliest genetic aberrations. Sequencing of tumors and their metastases reveals numerous abnormalities occurring late in progression. A means to temporally order aberrations in a single cancer, rather than inferring them from serially acquired samples, would define changes preceding even clinically evident disease. We integrate DNA sequence and copy number information to reconstruct the order of abnormalities as individual tumors evolve for 2 separate cancer types. We detect vast, unreported expansion of simple mutations sharply demarcated by recombinative loss of the second copy of TP53 in cutaneous squamous cell carcinomas (cSCC) and serous ovarian adenocarcinomas, in the former surpassing 50 mutations per megabase. In cSCCs, we also report diverse secondary mutations in known and novel oncogenic pathways, illustrating how such expanded mutagenesis directly promotes malignant progression. These results reframe paradigms in which TP53 mutation is required later, to bypass senescence induced by driver oncogenes.

Immutable functional attributes of histologic grade revealed by context-independent gene expression in primary breast cancer cells

Inherent cancer phenotypes that are independent of fluctuating cross-talk with the surrounding tissue matrix are highly desirable candidates for targeting tumor cells. Our novel study design uses epithelial cell lines derived from low versus high histologic grade primary breast cancer to effectively diminish the breadth of transient variability generated within the tumor microenvironment of the host, revealing a "paracrine-independent expression of grade-associated" (PEGA) gene signature. PEGA members extended beyond "proliferation-driven" signatures commonly associated with aggressive, high-grade breast cancer. The calcium-binding protein S100P was prominent among PEGA genes overexpressed in high-grade tumors. A three-member fingerprint of S100P-correlated genes, consisting of GPRC5A, FXYD3, and PYCARD, conferred poor outcome in multiple breast cancer data sets, irrespective of estrogen receptor status but dependent on tumor size (P < 0.01). S100P silencing markedly diminished coregulated gene transcripts and reversed aggressive tumor behavior. Exposure to pathway-implicated agents, including the calmodulin inhibitor N-(6-aminohexyl)-5-chloro-1-naphthalenesulfonamide, phenothiazine, and chlorpromazine, resulted in rapid apoptotic cell death in high-grade tumor cells resistant to the chemotherapeutic drug cisplatin. This is the first comprehensive study describing molecular phenotypes intimately associated with histologic grade whose expression remains relatively fixed despite an unavoidably changing environment to which tumor cells are invariably exposed.

Molecular mapping of Rxp conditioning reaction to bacterial pustule in soybean

The Rxp locus in soybean [Glycine max (L.) Merr.] that conditions reaction to bacterial pustule was mapped by simple sequence repeat (SSR) marker analysis. A population of 116 F4-derived lines from a cross between the resistant parent Young and the susceptible parent PI 416937 was used for mapping. The Rxp locus was mapped 3.9 cM from Satt372 and 12.4 cM from Satt014 on linkage group D2. Linkage associations were confirmed by identifying a close association between the SSR genotype at each locus identified as flanking Rxp and the bacterial pustule reaction of individual lines derived from a population different from the one used for mapping. A molecular pedigree analysis showed that bacterial pustule-resistant cultivars inherited the resistance gene rxp from the ancestral cultivar CNS based on their consistent genotypic pattern at flanking marker loci. Based on the results of the study, marker-assisted selection for rxp would be very effective.