Temporally resolved proteomics identifies nidogen-2 as a cotarget in pancreatic cancer that modulates fibrosis and therapy response.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by increasing fibrosis, which can enhance tumor progression and spread. Here, we undertook an unbiased temporal assessment of the matrisome of the highly metastatic KPC (Pdx1-Cre, LSL-KrasG12D/+, LSL-Trp53R172H/+) and poorly metastatic KPflC (Pdx1-Cre, LSL-KrasG12D/+, Trp53fl/+) genetically engineered mouse models of pancreatic cancer using mass spectrometry proteomics. Our assessment at early-, mid-, and late-stage disease reveals an increased abundance of nidogen-2 (NID2) in the KPC model compared to KPflC, with further validation showing that NID2 is primarily expressed by cancer-associated fibroblasts (CAFs). Using biomechanical assessments, second harmonic generation imaging, and birefringence analysis, we show that NID2 reduction by CRISPR interference (CRISPRi) in CAFs reduces stiffness and matrix remodeling in three-dimensional models, leading to impaired cancer cell invasion. Intravital imaging revealed improved vascular patency in live NID2-depleted tumors, with enhanced response to gemcitabine/Abraxane. In orthotopic models, NID2 CRISPRi tumors had less liver metastasis and increased survival, highlighting NID2 as a potential PDAC cotarget.

Circulating Tumour DNA Biomarkers Associated with Outcomes in Metastatic Prostate Cancer Treated with Lutetium-177-PSMA-617.

Background: Lutetium-177-prostate-specific membrane antigen- 617 (Lu-PSMA) is an effective therapy for metastatic castration-resistant prostate cancer (mCRPC). However, treatment responses are heterogeneous despite stringent positron emission tomography (PET)-based imaging selection criteria. Molecularly based biomarkers have potential to refine patient selection and optimise outcomes. Objective: To identify circulating tumour DNA (ctDNA) features associated with treatment outcomes for men treated with Lu-PSMA. Design setting and participants: ctDNA from men treated with Lu-PSMA in combination with idronoxil for progressive mCRPC were analysed using an 85-gene customised sequencing assay. ctDNA fractions, molecular profiles, and the presence of alterations in aggressive-variant prostate cancer (AVPC) genes were analysed at baseline, cycle 3 and at disease progression. Intervention: Men received Lu-PSMA with idronoxil every 6 wk for up to six cycles. Outcome measurements and statistical analysis: Baseline and exit PSMA and fluorodeoxyglucose PET/computed tomography (CT) imaging was conducted at baseline and study exit. Single-photon emission CT (SPECT) scans were performed 24 h after Lu-PSMA. Blood samples were collected at baseline,cycle 3 and at disease progression. Cox proportional-hazards models were used to assess associations and derive hazard ratios (HRs) and confidence intervals (CIs) for associations between molecular factors, imaging features, and clinical outcomes. Results and limitations: Sixty samples from 32 men were sequenced (32 at baseline, 24 at cycle 3, four from patients with disease progression); two samples (baseline, on-treatment) from one individual were excluded from analysis owing to poor quality of the baseline sequencing data. Alterations in AVPC genes were associated with shorter prostate-specific antigen (PSA) progression-free survival (PFS) and overall survival (OS) in univariate (HR 3.4, 95% CI 1.5-7.7; p = 0.0036; and HR 3.3, 95% CI 1.4-7.7; p = 0.0063, respectively) and multivariate analyses (HR 4.8, 95% CI 1.8-13; p = 0.0014; and HR 4.1, 95% CI 1.6-11; p = 0.004). Conclusions: ctDNA alterations in AVPC genes were associated with shorter PSA PFS and OS among men treated with Lu-PSMA and intermittent idronoxil. These candidate molecular biomarkers warrant further study to determine whether they have predictive value and potential to guide synergistic combination strategies to enhance outcomes for men treated with Lu-PSMA for mCRPC. Patient summary: Certain DNA/gene changes detected in the blood of men with advanced prostate cancer were associated with shorter benefit from lutetium PSMA, a targeted radioactive therapy. This information may be useful in determining which men may benefit most from this treatment, but additional research is needed.

Neutrophil Conversion to a Tumor-Killing Phenotype Underpins Effective Microbial Therapy.

The inflammatory microenvironment of solid tumors creates a protumorigenic milieu that resembles chronic inflammation akin to a subverted wound healing response. Here, we investigated the effect of converting the tumor microenvironment from a chronically inflamed state to one of acute microbial inflammation by injecting microbial bioparticles directly into tumors. Intratumoral microbial bioparticle injection led to rapid and dramatic changes in the tumor immune composition, the most striking of which was a substantial increase in the presence of activated neutrophils. In situ photoconversion and intravital microscopy indicated that tumor neutrophils transiently switched from sessile producers of VEGF to highly motile neutrophils that clustered to make neutrophil-rich domains in the tumor. The neutrophil clusters remodeled tumor tissue and repressed tumor growth. Single-cell transcriptional analysis of microbe-stimulated neutrophils showed a profound shift in gene expression towards heightened activation and antimicrobial effector function. Microbe-activated neutrophils also upregulated chemokines known to regulate neutrophil and CD8+ T-cell recruitment. Microbial therapy also boosted CD8+ T-cell function and enhanced the therapeutic benefit of checkpoint inhibitor therapy in tumor-bearing mice and provided protection in a model of tumor recurrence. These data indicate that one of the major effector mechanisms of microbial therapy is the conversion of tumor neutrophils from a wound healing to an acutely activated cytotoxic phenotype, highlighting a rationale for broader deployment of microbial therapy in the treatment of solid cancers. SIGNIFICANCE: Intratumoral injection of microbial bioparticles stimulates neutrophil antitumor functions, suggesting pathways for optimizing efficacy of microbial therapies and paving the way for their broader utilization in the clinic.

Quantitative ctDNA Detection in Hepatoblastoma: Implications for Precision Medicine.

Hepatoblastoma is characterized by driver mutations in CTNNB1, making it an attractive biomarker for a liquid biopsy approach utilizing circulating tumor DNA (ctDNA). This prospective observational study sought to ascertain the feasibility of ctDNA detection in patients with hepatoblastoma and explore its associations with established clinical indicators and biomarkers, including serum Alpha-fetoprotein (AFP). We obtained 38 plasma samples and 17 tumor samples from 20 patients with hepatoblastoma. These samples were collected at various stages: 10 at initial diagnosis, 17 during neoadjuvant chemotherapy, 6 post-operatively, and 5 at disease recurrence. Utilizing a bespoke sequencing assay we developed called QUENCH, we identified single nucleotide variants and deletions in CTNNB1 ctDNA. Our study demonstrated the capability to quantitate ctDNA down to a variant allele frequency of 0.3%, achieving a sensitivity of 90% for patients at initial diagnosis, and a specificity of 100% at the patient level. Notably, ctDNA positivity correlated with tumor burden, and ctDNA levels exhibited associations with macroscopic residual disease and treatment response. Our findings provide evidence for the utility of quantitative ctDNA detection in hepatoblastoma management. Given the distinct detection targets, ctDNA and AFP-based stratification and monitoring approaches could synergize to enhance clinical decision-making. Further research is needed to elucidate the interplay between ctDNA and AFP and determine the optimal clinical applications for both methods in risk stratification and residual disease detection.

Integrated genomic analysis reveals mutated ELF3 as a potential gallbladder cancer vaccine candidate.

Gallbladder cancer (GBC) is an aggressive gastrointestinal malignancy with no approved targeted therapy. Here, we analyze exomes (n = 160), transcriptomes (n = 115), and low pass whole genomes (n = 146) from 167 gallbladder cancers (GBCs) from patients in Korea, India and Chile. In addition, we also sequence samples from 39 GBC high-risk patients and detect evidence of early cancer-related genomic lesions. Among the several significantly mutated genes not previously linked to GBC are ETS domain genes ELF3 and EHF, CTNNB1, APC, NSD1, KAT8, STK11 and NFE2L2. A majority of ELF3 alterations are frame-shift mutations that result in several cancer-specific neoantigens that activate T-cells indicating that they are cancer vaccine candidates. In addition, we identify recurrent alterations in KEAP1/NFE2L2 and WNT pathway in GBC. Taken together, these define multiple targetable therapeutic interventions opportunities for GBC treatment and management.

Inducers of the endothelial cell barrier identified through chemogenomic screening in genome-edited hPSC-endothelial cells.

The blood-retina barrier and blood-brain barrier (BRB/BBB) are selective and semipermeable and are critical for supporting and protecting central nervous system (CNS)-resident cells. Endothelial cells (ECs) within the BRB/BBB are tightly coupled, express high levels of Claudin-5 (CLDN5), a junctional protein that stabilizes ECs, and are important for proper neuronal function. To identify novel CLDN5 regulators (and ultimately EC stabilizers), we generated a CLDN5-P2A-GFP stable cell line from human pluripotent stem cells (hPSCs), directed their differentiation to ECs (CLDN5-GFP hPSC-ECs), and performed flow cytometry-based chemogenomic library screening to measure GFP expression as a surrogate reporter of barrier integrity. Using this approach, we identified 62 unique compounds that activated CLDN5-GFP. Among them were TGF-ß pathway inhibitors, including RepSox. When applied to hPSC-ECs, primary brain ECs, and retinal ECs, RepSox strongly elevated barrier resistance (transendothelial electrical resistance), reduced paracellular permeability (fluorescein isothiocyanate-dextran), and prevented vascular endothelial growth factor A (VEGFA)-induced barrier breakdown in vitro. RepSox also altered vascular patterning in the mouse retina during development when delivered exogenously. To determine the mechanism of action of RepSox, we performed kinome-, transcriptome-, and proteome-profiling and discovered that RepSox inhibited TGF-ß, VEGFA, and inflammatory gene networks. In addition, RepSox not only activated vascular-stabilizing and barrier-establishing Notch and Wnt pathways, but also induced expression of important tight junctions and transporters. Taken together, our data suggest that inhibiting multiple pathways by selected individual small molecules, such as RepSox, may be an effective strategy for the development of better BRB/BBB models and novel EC barrier-inducing therapeutics.

Peripheral T cell expansion predicts tumour infiltration and clinical response.

Despite the resounding clinical success in cancer treatment of antibodies that block the interaction of PD1 with its ligand PDL11, the mechanisms involved remain unknown. A major limitation to understanding the origin and fate of T cells in tumour immunity is the lack of quantitative information on the distribution of individual clonotypes of T cells in patients with cancer. Here, by performing deep single-cell sequencing of RNA and T cell receptors in patients with different types of cancer, we survey the profiles of various populations of T cells and T cell receptors in tumours, normal adjacent tissue, and peripheral blood. We find clear evidence of clonotypic expansion of effector-like T cells not only within the tumour but also in normal adjacent tissue. Patients with gene signatures of such clonotypic expansion respond best to anti-PDL1 therapy. Notably, expanded clonotypes found in the tumour and normal adjacent tissue can also typically be detected in peripheral blood, which suggests a convenient approach to patient identification. Analyses of our data together with several external datasets suggest that intratumoural T cells, especially in responsive patients, are replenished with fresh, non-exhausted replacement cells from sites outside the tumour, suggesting continued activity of the cancer immunity cycle in these patients, the acceleration of which may be associated with clinical response.

S100a4 upregulation in Pik3caH1047R;Trp53R270H;MMTV-Cre-driven mammary tumors promotes metastasis.

BACKGROUND: PIK3CA mutations are frequent in human breast cancer. Pik3caH1047R mutant expression in mouse mammary gland promotes tumorigenesis. TP53 mutations co-occur with PIK3CA mutations in human breast cancers. We previously generated a conditionally activatable Pik3caH1047R;MMTV-Cre mouse model and found a few malignant sarcomatoid (spindle cell) carcinomas that had acquired spontaneous dominant-negative Trp53 mutations. METHODS: A Pik3caH1047R;Trp53R270H;MMTV-Cre double mutant mouse breast cancer model was generated. Tumors were characterized by histology, marker analysis, transcriptional profiling, single-cell RNA-seq, and bioinformatics. Cell lines were developed from mutant tumors and used to identify and confirm genes involved in metastasis. RESULTS: We found Pik3caH1047R and Trp53R270H cooperate in driving oncogenesis in mammary glands leading to a shorter latency than either alone. Double mutant mice develop multiple histologically distinct mammary tumors, including adenocarcinoma and sarcomatoid (spindle cell) carcinoma. We found some tumors to be invasive and a few metastasized to the lung and/or the lymph node. Single-cell RNA-seq analysis of the tumors identified epithelial, stromal, myeloid, and T cell groups. Expression analysis of the metastatic tumors identified S100a4 as a top candidate gene associated with metastasis. Metastatic tumors contained a much higher percentage of epithelial-mesenchymal transition (EMT)-signature positive and S100a4-expressing cells. CRISPR/CAS9-mediated knockout of S100a4 in a metastatic tumor-derived cell line disrupted its metastatic potential indicating a role for S100a4 in metastasis. CONCLUSIONS: Pik3caH1047R;Trp53R270H;MMTV-Cre mouse provides a preclinical model to mimic a subtype of human breast cancers that carry both PIK3CA and TP53 mutations. It also allows for understanding the cooperation between the two mutant genes in tumorigenesis. Our model also provides a system to study metastasis and develop therapeutic strategies for PIK3CA/TP53 double-positive cancers. S100a4 found involved in metastasis in this model can be a potential diagnostic and therapeutic target.

An Empirical Approach Leveraging Tumorgrafts to Dissect the Tumor Microenvironment in Renal Cell Carcinoma Identifies Missing Link to Prognostic Inflammatory Factors.

By leveraging tumorgraft (patient-derived xenograft) RNA-sequencing data, we developed an empirical approach, DisHet, to dissect the tumor microenvironment (eTME). We found that 65% of previously defined immune signature genes are not abundantly expressed in renal cell carcinoma (RCC) and identified 610 novel immune/stromal transcripts. Using eTME, genomics, pathology, and medical record data involving >1,000 patients, we established an inflamed pan-RCC subtype (IS) enriched for regulatory T cells, natural killer cells, TH1 cells, neutrophils, macrophages, B cells, and CD8+ T cells. IS is enriched for aggressive RCCs, including BAP1-deficient clear-cell and type 2 papillary tumors. The IS subtype correlated with systemic manifestations of inflammation such as thrombocytosis and anemia, which are enigmatic predictors of poor prognosis. Furthermore, IS was a strong predictor of poor survival. Our analyses suggest that tumor cells drive the stromal immune response. These data provide a missing link between tumor cells, the TME, and systemic factors.Significance: We undertook a novel empirical approach to dissect the renal cell carcinoma TME by leveraging tumorgrafts. The dissection and downstream analyses uncovered missing links between tumor cells, the TME, systemic manifestations of inflammation, and poor prognosis. Cancer Discov; 8(9); 1142-55. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 1047.

Comprehensive genomic analysis identifies pathogenic variants in maturity-onset diabetes of the young (MODY) patients in South India.

BACKGROUND: Maturity-onset diabetes of the young (MODY) is an early-onset, autosomal dominant form of non-insulin dependent diabetes. Genetic diagnosis of MODY can transform patient management. Earlier data on the genetic predisposition to MODY have come primarily from familial studies in populations of European origin. METHODS: In this study, we carried out a comprehensive genomic analysis of 289 individuals from India that included 152 clinically diagnosed MODY cases to identify variants in known MODY genes. Further, we have analyzed exome data to identify putative MODY relevant variants in genes previously not implicated in MODY. Functional validation of MODY relevant variants was also performed. RESULTS: We found MODY 3 (HNF1A; 7.2%) to be most frequently mutated followed by MODY 12 (ABCC8; 3.3%). They together account for ~ 11% of the cases. In addition to known MODY genes, we report the identification of variants in RFX6, WFS1, AKT2, NKX6-1 that may contribute to development of MODY. Functional assessment of the NKX6-1 variants showed that they are functionally impaired. CONCLUSIONS: Our findings showed HNF1A and ABCC8 to be the most frequently mutated MODY genes in south India. Further we provide evidence for additional MODY relevant genes, such as NKX6-1, and these require further validation.

Recurrent Loss of NFE2L2 Exon 2 Is a Mechanism for Nrf2 Pathway Activation in Human Cancers.

The Nrf2 pathway is frequently activated in human cancers through mutations in Nrf2 or its negative regulator KEAP1. Using a cell-line-derived gene signature for Nrf2 pathway activation, we found that some tumors show high Nrf2 activity in the absence of known mutations in the pathway. An analysis of splice variants in oncogenes revealed that such tumors express abnormal transcript variants from the NFE2L2 gene (encoding Nrf2) that lack exon 2, or exons 2 and 3, and encode Nrf2 protein isoforms missing the KEAP1 interaction domain. The Nrf2 alterations result in the loss of interaction with KEAP1, Nrf2 stabilization, induction of a Nrf2 transcriptional response, and Nrf2 pathway dependence. In all analyzed cases, transcript variants were the result of heterozygous genomic microdeletions. Thus, we identify an alternative mechanism for Nrf2 pathway activation in human tumors and elucidate its functional consequences.

Untangling the brain's neuroinflammatory and neurodegenerative transcriptional responses.

A common approach to understanding neurodegenerative disease is comparing gene expression in diseased versus healthy tissues. We illustrate that expression profiles derived from whole tissue RNA highly reflect the degenerating tissues' altered cellular composition, not necessarily transcriptional regulation. To accurately understand transcriptional changes that accompany neuropathology, we acutely purify neurons, astrocytes and microglia from single adult mouse brains and analyse their transcriptomes by RNA sequencing. Using peripheral endotoxemia to establish the method, we reveal highly specific transcriptional responses and altered RNA processing in each cell type, with Tnfr1 required for the astrocytic response. Extending the method to an Alzheimer's disease model, we confirm that transcriptomic changes observed in whole tissue are driven primarily by cell type composition, not transcriptional regulation, and identify hundreds of cell type-specific changes undetected in whole tissue RNA. Applying similar methods to additional models and patient tissues will transform our understanding of aberrant gene expression in neurological disease.

Comprehensive genomic analysis of malignant pleural mesothelioma identifies recurrent mutations, gene fusions and splicing alterations.

We analyzed transcriptomes (n = 211), whole exomes (n = 99) and targeted exomes (n = 103) from 216 malignant pleural mesothelioma (MPM) tumors. Using RNA-seq data, we identified four distinct molecular subtypes: sarcomatoid, epithelioid, biphasic-epithelioid (biphasic-E) and biphasic-sarcomatoid (biphasic-S). Through exome analysis, we found BAP1, NF2, TP53, SETD2, DDX3X, ULK2, RYR2, CFAP45, SETDB1 and DDX51 to be significantly mutated (q-score ≥ 0.8) in MPMs. We identified recurrent mutations in several genes, including SF3B1 (∼2%; 4/216) and TRAF7 (∼2%; 5/216). SF3B1-mutant samples showed a splicing profile distinct from that of wild-type tumors. TRAF7 alterations occurred primarily in the WD40 domain and were, except in one case, mutually exclusive with NF2 alterations. We found recurrent gene fusions and splice alterations to be frequent mechanisms for inactivation of NF2, BAP1 and SETD2. Through integrated analyses, we identified alterations in Hippo, mTOR, histone methylation, RNA helicase and p53 signaling pathways in MPMs.

Complex regulation of ADAR-mediated RNA-editing across tissues.

BACKGROUND: RNA-editing is a tightly regulated, and essential cellular process for a properly functioning brain. Dysfunction of A-to-I RNA editing can have catastrophic effects, particularly in the central nervous system. Thus, understanding how the process of RNA-editing is regulated has important implications for human health. However, at present, very little is known about the regulation of editing across tissues, and individuals. RESULTS: Here we present an analysis of RNA-editing patterns from 9 different tissues harvested from a single mouse. For comparison, we also analyzed data for 5 of these tissues harvested from 15 additional animals. We find that tissue specificity of editing largely reflects differential expression of substrate transcripts across tissues. We identified a surprising enrichment of editing in intronic regions of brain transcripts, that could account for previously reported higher levels of editing in brain. There exists a small but remarkable amount of editing which is tissue-specific, despite comparable expression levels of the edit site across multiple tissues. Expression levels of editing enzymes and their isoforms can explain some, but not all of this variation. CONCLUSIONS: Together, these data suggest a complex regulation of the RNA-editing process beyond transcript expression levels.

A peptide filtering relation quantifies MHC class I peptide optimization.

Major Histocompatibility Complex (MHC) class I molecules enable cytotoxic T lymphocytes to destroy virus-infected or cancerous cells, thereby preventing disease progression. MHC class I molecules provide a snapshot of the contents of a cell by binding to protein fragments arising from intracellular protein turnover and presenting these fragments at the cell surface. Competing fragments (peptides) are selected for cell-surface presentation on the basis of their ability to form a stable complex with MHC class I, by a process known as peptide optimization. A better understanding of the optimization process is important for our understanding of immunodominance, the predominance of some T lymphocyte specificities over others, which can determine the efficacy of an immune response, the danger of immune evasion, and the success of vaccination strategies. In this paper we present a dynamical systems model of peptide optimization by MHC class I. We incorporate the chaperone molecule tapasin, which has been shown to enhance peptide optimization to different extents for different MHC class I alleles. Using a combination of published and novel experimental data to parameterize the model, we arrive at a relation of peptide filtering, which quantifies peptide optimization as a function of peptide supply and peptide unbinding rates. From this relation, we find that tapasin enhances peptide unbinding to improve peptide optimization without significantly delaying the transit of MHC to the cell surface, and differences in peptide optimization across MHC class I alleles can be explained by allele-specific differences in peptide binding. Importantly, our filtering relation may be used to dynamically predict the cell surface abundance of any number of competing peptides by MHC class I alleles, providing a quantitative basis to investigate viral infection or disease at the cellular level. We exemplify this by simulating optimization of the distribution of peptides derived from Human Immunodeficiency Virus Gag-Pol polyprotein.

Characterisation of microRNA expression in post-natal mouse mammary gland development.

BACKGROUND: The differential expression pattern of microRNAs (miRNAs) during mammary gland development might provide insights into their role in regulating the homeostasis of the mammary epithelium. Our aim was to analyse these regulatory functions by deriving a comprehensive tissue-specific combined miRNA and mRNA expression profile of post-natal mouse mammary gland development.We measured the expression of 318 individual murine miRNAs by bead-based flow-cytometric profiling of whole mouse mammary glands throughout a 16-point developmental time course, including juvenile, puberty, mature virgin, gestation, lactation, and involution stages. In parallel whole-genome mRNA expression data were obtained. RESULTS: One third (n = 102) of all murine miRNAs analysed were detected during mammary gland development. MicroRNAs were represented in seven temporally co-expressed clusters, which were enriched for both miRNAs belonging to the same family and breast cancer-associated miRNAs. Global miRNA and mRNA expression was significantly reduced during lactation and the early stages of involution after weaning. For most detected miRNA families we did not observe systematic changes in the expression of predicted targets. For miRNA families whose targets did show changes, we observed inverse patterns of miRNA and target expression. The data sets are made publicly available and the combined expression profiles represent an important community resource for mammary gland biology research. CONCLUSION: MicroRNAs were expressed in likely co-regulated clusters during mammary gland development. Breast cancer-associated miRNAs were significantly enriched in these clusters. The mechanism and functional consequences of this miRNA co-regulation provide new avenues for research into mammary gland biology and generate candidates for functional validation.

AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma.

Invasive urothelial cell carcinoma (UCC) is characterized by increased chromosomal instability and follows an aggressive clinical course in contrast to non-invasive disease. To identify molecular processes that confer and maintain an aggressive malignant phenotype, we used a high-throughput genome-wide approach to interrogate a cohort of high and low clinical risk UCC tumors. Differential expression analyses highlighted cohesive dysregulation of critical genes involved in the G(2)/M checkpoint in aggressive UCC. Hierarchical clustering based on DNA Damage Response (DDR) genes separated tumors according to a pre-defined clinical risk phenotype. Using array-comparative genomic hybridization, we confirmed that the DDR was disrupted in tumors displaying high genomic instability. We identified DNA copy number gains at 20q13.2-q13.3 (AURKA locus) and determined that overexpression of AURKA accompanied dysregulation of DDR genes in high risk tumors. We postulated that DDR-deficient UCC tumors are advantaged by a selective pressure for AURKA associated override of M phase barriers and confirmed this in an independent tissue microarray series. This mechanism that enables cancer cells to maintain an aggressive phenotype forms a rationale for targeting AURKA as a therapeutic strategy in advanced stage UCC.

MicroRNA expression profiling of human breast cancer identifies new markers of tumor subtype.

BACKGROUND: MicroRNAs (miRNAs), a class of short non-coding RNAs found in many plants and animals, often act post-transcriptionally to inhibit gene expression. RESULTS: Here we report the analysis of miRNA expression in 93 primary human breast tumors, using a bead-based flow cytometric miRNA expression profiling method. Of 309 human miRNAs assayed, we identify 133 miRNAs expressed in human breast and breast tumors. We used mRNA expression profiling to classify the breast tumors as luminal A, luminal B, basal-like, HER2+ and normal-like. A number of miRNAs are differentially expressed between these molecular tumor subtypes and individual miRNAs are associated with clinicopathological factors. Furthermore, we find that miRNAs could classify basal versus luminal tumor subtypes in an independent data set. In some cases, changes in miRNA expression correlate with genomic loss or gain; in others, changes in miRNA expression are likely due to changes in primary transcription and or miRNA biogenesis. Finally, the expression of DICER1 and AGO2 is correlated with tumor subtype and may explain some of the changes in miRNA expression observed. CONCLUSION: This study represents the first integrated analysis of miRNA expression, mRNA expression and genomic changes in human breast cancer and may serve as a basis for functional studies of the role of miRNAs in the etiology of breast cancer. Furthermore, we demonstrate that bead-based flow cytometric miRNA expression profiling might be a suitable platform to classify breast cancer into prognostic molecular subtypes.