Joel S. Parker

Repeated observation of breast tumor subtypes in independent gene expression data sets

Characteristic patterns of gene expression measured by DNA microarrays have been used to classify tumors into clinically relevant subgroups. In this study, we have refined the previously defined subtypes of breast tumors that could be distinguished by their distinct patterns of gene expression. A total of 115 malignant breast tumors were analyzed by hierarchical clustering based on patterns of expression of 534 “intrinsic” genes and shown to subdivide into one basal-like, one ERBB2-overexpressing, two luminal-like, and one normal breast tissue-like subgroup. The genes used for classification were selected based on their similar expression levels between pairs of consecutive samples taken from the same tumor separated by 15 weeks of neoadjuvant treatment. Similar cluster analyses of two published, independent data sets representing different patient cohorts from different laboratories, uncovered some of the same breast cancer subtypes. In the one data set that included information on time to development of distant metastasis, subtypes were associated with significant differences in this clinical feature. By including a group of tumors from BRCA1 carriers in the analysis, we found that this genotype predisposes to the basal tumor subtype. Our results strongly support the idea that many of these breast tumor subtypes represent biologically distinct disease entities.

Ki67 Index, HER2 Status, and Prognosis of Patients With Luminal B Breast Cancer

Background Gene expression profiling of breast cancer has identified two biologically distinct estrogen receptor (ER)-positive subtypes of breast cancer: luminal A and luminal B. Luminal B tumors have higher proliferation and poorer prognosis than luminal A tumors. In this study, we developed a clinically practical immunohistochemistry assay to distinguish luminal B from luminal A tumors and investigated its ability to separate tumors according to breast cancer recurrence-free and disease-specific survival. Methods Tumors from a cohort of 357 patients with invasive breast carcinomas were subtyped by gene expression profile. Hormone receptor status, HER2 status, and the Ki67 index (percentage of Ki67-positive cancer nuclei) were determined immunohistochemically. Receiver operating characteristic curves were used to determine the Ki67 cut point to distinguish luminal B from luminal A tumors. The prognostic value of the immunohistochemical assignment for breast cancer recurrence-free and disease-specific survival was investigated with an independent tissue microarray series of 4046 breast cancers by use of Kaplan–Meier curves and multivariable Cox regression. Results Gene expression profiling classified 101 (28%) of the 357 tumors as luminal A and 69 (19%) as luminal B. The best Ki67 index cut point to distinguish luminal B from luminal A tumors was 13.25%. In an independent cohort of 4046 patients with breast cancer, 2847 had hormone receptor–positive tumors. When HER2 immunohistochemistry and the Ki67 index were used to subtype these 2847 tumors, we classified 1530 (59%, 95% confidence interval [CI] = 57% to 61%) as luminal A, 846 (33%, 95% CI = 31% to 34%) as luminal B, and 222 (9%, 95% CI = 7% to 10%) as luminal–HER2 positive. Luminal B and luminal–HER2-positive breast cancers were statistically significantly associated with poor breast cancer recurrence-free and disease-specific survival in all adjuvant systemic treatment categories. Of particular relevance are women who received tamoxifen as their sole adjuvant systemic therapy, among whom the 10-year breast cancer–specific survival was 79% (95% CI = 76% to 83%) for luminal A, 64% (95% CI = 59% to 70%) for luminal B, and 57% (95% CI = 47% to 69%) for luminal–HER2 subtypes. Conclusion Expression of ER, progesterone receptor, and HER2 proteins and the Ki67 index appear to distinguish luminal A from luminal B breast cancer subtypes.

Phenotypic and molecular characterization of the claudin-low intrinsic subtype of breast cancer

IntroductionIn breast cancer, gene expression analyses have defined five tumor subtypes (luminal A, luminal B, HER2-enriched, basal-like and claudin-low), each of which has unique biologic and prognostic features. Here, we comprehensively characterize the recently identified claudin-low tumor subtype.MethodsThe clinical, pathological and biological features of claudin-low tumors were compared to the other tumor subtypes using an updated human tumor database and multiple independent data sets. These main features of claudin-low tumors were also evaluated in a panel of breast cancer cell lines and genetically engineered mouse models.ResultsClaudin-low tumors are characterized by the low to absent expression of luminal differentiation markers, high enrichment for epithelial-to-mesenchymal transition markers, immune response genes and cancer stem cell-like features. Clinically, the majority of claudin-low tumors are poor prognosis estrogen receptor (ER)-negative, progesterone receptor (PR)-negative, and epidermal growth factor receptor 2 (HER2)-negative (triple negative) invasive ductal carcinomas with a high frequency of metaplastic and medullary differentiation. They also have a response rate to standard preoperative chemotherapy that is intermediate between that of basal-like and luminal tumors. Interestingly, we show that a group of highly utilized breast cancer cell lines, and several genetically engineered mouse models, express the claudin-low phenotype. Finally, we confirm that a prognostically relevant differentiation hierarchy exists across all breast cancers in which the claudin-low subtype most closely resembles the mammary epithelial stem cell.ConclusionsThese results should help to improve our understanding of the biologic heterogeneity of breast cancer and provide tools for the further evaluation of the unique biology of claudin-low tumors and cell lines.

A custom microarray platform for analysis of microRNA gene expression.

MicroRNAs are short, noncoding RNA transcripts that post-transcriptionally regulate gene expression. Several hundred microRNA genes have been identified in Caenorhabditis elegans, Drosophila, plants and mammals. MicroRNAs have been linked to developmental processes in C. elegans, plants and humans and to cell growth and apoptosis in Drosophila. A major impediment in the study of microRNA function is the lack of quantitative expression profiling methods. To close this technological gap, we have designed dual-channel microarrays that monitor expression levels of 124 mammalian microRNAs. Using these tools, we observed distinct patterns of expression among adult mouse tissues and embryonic stem cells. Expression profiles of staged embryos demonstrate temporal regulation of a large class of microRNAs, including members of the let-7 family. This microarray technology enables comprehensive investigation of microRNA expression, and furthers our understanding of this class of recently discovered noncoding RNAs.

RNA polymerase is poised for activation across the genome

Regulation of gene expression is integral to the development and survival of all organisms. Transcription begins with the assembly of a pre-initiation complex at the gene promoter, followed by initiation of RNA synthesis and the transition to productive elongation. In many cases, recruitment of RNA polymerase II (Pol II) to a promoter is necessary and sufficient for activation of genes. However, there are a few notable exceptions to this paradigm, including heat shock genes and several proto-oncogenes, whose expression is attenuated by regulated stalling of polymerase elongation within the promoter-proximal region. To determine the importance of polymerase stalling for transcription regulation, we carried out a genome-wide search for Drosophila melanogaster genes with Pol II stalled within the promoter-proximal region. Our data show that stalling is widespread, occurring at hundreds of genes that respond to stimuli and developmental signals. This finding indicates a role for regulation of polymerase elongation in the transcriptional responses to dynamic environmental and developmental cues.

Molecular classification of head and neck squamous cell carcinomas using patterns of gene expression.

The prognostication of head and neck squamous cell carcinoma (HNSCC) is largely based upon the tumor size and location and the presence of lymph node metastases. Here we show that gene expression patterns from 60 HNSCC samples assayed on cDNA microarrays allowed categorization of these tumors into four distinct subtypes. These subtypes showed statistically significant differences in recurrence-free survival and included a subtype with a possible EGFR-pathway signature, a mesenchymal-enriched subtype, a normal epithelium-like subtype, and a subtype with high levels of antioxidant enzymes. Supervised analyses to predict lymph node metastasis status were approximately 80% accurate when tumor subsite and pathological node status were considered simultaneously. This work represents an important step toward the identification of clinically significant biomarkers for HNSCC.

A Comparison of PAM50 Intrinsic Subtyping with Immunohistochemistry and Clinical Prognostic Factors in Tamoxifen-Treated Estrogen Receptor-Positive Breast Cancer

Purpose: To compare clinical, immunohistochemical (IHC), and gene expression models of prognosis applicable to formalin-fixed, paraffin-embedded blocks in a large series of estrogen receptor (ER)–positive breast cancers from patients uniformly treated with adjuvant tamoxifen. Experimental Design: Quantitative real-time reverse transcription-PCR (qRT-PCR) assays for 50 genes identifying intrinsic breast cancer subtypes were completed on 786 specimens linked to clinical (median follow-up, 11.7 years) and IHC [ER, progesterone receptor (PR), HER2, and Ki67] data. Performance of predefined intrinsic subtype and risk-of-relapse scores was assessed using multivariable Cox models and Kaplan-Meier analysis. Harrells C-index was used to compare fixed models trained in independent data sets, including proliferation signatures. Results: Despite clinical ER positivity, 10% of cases were assigned to nonluminal subtypes. qRT-PCR signatures for proliferation genes gave more prognostic information than clinical assays for hormone receptors or Ki67. In Cox models incorporating standard prognostic variables, hazard ratios for breast cancer disease-specific survival over the first 5 years of follow-up, relative to the most common luminal A subtype, are 1.99 [95% confidence interval (CI), 1.09-3.64] for luminal B, 3.65 (95% CI, 1.64-8.16) for HER2-enriched subtype, and 17.71 (95% CI, 1.71-183.33) for the basal-like subtype. For node-negative disease, PAM50 qRT-PCR–based risk assignment weighted for tumor size and proliferation identifies a group with >95% 10-year survival without chemotherapy. In node-positive disease, PAM50-based prognostic models were also superior. Conclusion: The PAM50 gene expression test for intrinsic biological subtype can be applied to large series of formalin-fixed, paraffin-embedded breast cancers, and gives more prognostic information than clinical factors and IHC using standard cut points. Clin Cancer Res; 16(21); 5222–32. ©2010 AACR.

microRNA expression in the prefrontal cortex of individuals with schizophrenia and schizoaffective disorder

BackgroundmicroRNAs (miRNAs) are small, noncoding RNA molecules that are now thought to regulate the expression of many mRNAs. They have been implicated in the etiology of a variety of complex diseases, including Tourettes syndrome, Fragile × syndrome, and several types of cancer.ResultsWe hypothesized that schizophrenia might be associated with altered miRNA profiles. To investigate this possibility we compared the expression of 264 human miRNAs from postmortem prefrontal cortex tissue of individuals with schizophrenia (n = 13) or schizoaffective disorder (n = 2) to tissue of 21 psychiatrically unaffected individuals using a custom miRNA microarray. Allowing a 5% false discovery rate, we found that 16 miRNAs were differentially expressed in prefrontal cortex of patient subjects, with 15 expressed at lower levels (fold change 0.63 to 0.89) and 1 at a higher level (fold change 1.77) than in the psychiatrically unaffected comparison subjects. The expression levels of 12 selected miRNAs were also determined by quantitative RT-PCR in our lab. For the eight miRNAs distinguished by being expressed at lower microarray levels in schizophrenia samples versus comparison samples, seven were also expressed at lower levels with quantitative RT-PCR.ConclusionThis study is the first to find altered miRNA profiles in postmortem prefrontal cortex from schizophrenia patients.

Ezh2 Orchestrates Gene Expression for the Stepwise Differentiation of Tissue-Specific Stem Cells

Although in vitro studies of embryonic stem cells have identified polycomb repressor complexes (PRCs) as key regulators of differentiation, it remains unclear as to how PRC-mediated mechanisms control fates of multipotent progenitors in developing tissues. Here, we show that an essential PRC component, Ezh2, is expressed in epidermal progenitors but diminishes concomitant with embryonic differentiation and with postnatal decline in proliferative activity. We show that Ezh2 controls proliferative potential of basal progenitors by repressing the Ink4A-Ink4B locus and tempers the developmental rate of differentiation by preventing premature recruitment of AP1 transcriptional activator to the structural genes that are required for epidermal differentiation. Together, our studies reveal that PRCs control epigenetic modifications temporally and spatially in tissue-restricted stem cells. They maintain their proliferative potential and globally repressing undesirable differentiation programs while selectively establishing a specific terminal differentiation program in a stepwise fashion.

A biochemical approach to identifying microRNA targets

Identifying the downstream targets of microRNAs (miRNAs) is essential to understanding cellular regulatory networks. We devised a direct biochemical method for miRNA target discovery that combined RNA-induced silencing complex (RISC) purification with microarray analysis of bound mRNAs. Because targets of miR-124a have been analyzed, we chose it as our model. We honed our approach both by examining the determinants of stable binding between RISC and synthetic target RNAs in vitro and by determining the dependency of both repression and RISC coimmunoprecipitation on miR-124a seed sites in two of its well characterized targets in vivo. Examining the complete spectrum of miR-124 targets in 293 cells yielded both a set that were down-regulated at the mRNA level, as previously observed, and a set whose mRNA levels were unaffected by miR-124a. Reporter assays validated both classes, extending the spectrum of mRNA targets that can be experimentally linked to the miRNA pathway.