A single-cell transcriptional gradient in human cutaneous memory T cells restricts Th17/Tc17 identity.

The homeostatic mechanisms that fail to restrain chronic tissue inflammation in diseases, such as psoriasis vulgaris, remain incompletely understood. We profiled transcriptomes and epitopes of single psoriatic and normal skin-resident T cells, revealing a gradated transcriptional program of coordinately regulated inflammation-suppressive genes. This program, which is sharply suppressed in lesional skin, strikingly restricts Th17/Tc17 cytokine and other inflammatory mediators on the single-cell level. CRISPR-based deactivation of two core components of this inflammation-suppressive program, ZFP36L2 and ZFP36, replicates the interleukin-17A (IL-17A), granulocyte macrophage-colony-stimulating factor (GM-CSF), and interferon gamma (IFNγ) elevation in psoriatic memory T cells deficient in these transcripts, functionally validating their influence. Combinatoric expression analysis indicates the suppression of specific inflammatory mediators by individual program members. Finally, we find that therapeutic IL-23 blockade reduces Th17/Tc17 cell frequency in lesional skin but fails to normalize this inflammatory-suppressive program, suggesting how treated lesions may be primed for recurrence after withdrawal of treatment.

Identification and validation of expressed HLA-binding breast cancer neoepitopes for potential use in individualized cancer therapy.

BACKGROUND: Therapeutic regimens designed to augment the immunological response of a patient with breast cancer (BC) to tumor tissue are critically informed by tumor mutational burden and the antigenicity of expressed neoepitopes. Herein we describe a neoepitope and cognate neoepitope-reactive T-cell identification and validation program that supports the development of next-generation immunotherapies. METHODS: Using GPS Cancer, NantOmics research, and The Cancer Genome Atlas databases, we developed a novel bioinformatic-based approach which assesses mutational load, neoepitope expression, human leukocyte antigen (HLA)-binding prediction, and in vitro confirmation of T-cell recognition to preferentially identify targetable neoepitopes. This program was validated by application to a BC cell line and confirmed using tumor biopsies from two patients with BC enrolled in the Tumor-Infiltrating Lymphocytes and Genomics (TILGen) study. RESULTS: The antigenicity and HLA-A2 restriction of the BC cell line predicted neoepitopes were determined by reactivity of T cells from HLA-A2-expressing healthy donors. For the TILGen subjects, tumor-infiltrating lymphocytes (TILs) recognized the predicted neoepitopes both as peptides and on retroviral expression in HLA-matched Epstein-Barr virus-lymphoblastoid cell line and BC cell line MCF-7 cells; PCR clonotyping revealed the presence of T cells in the periphery with T-cell receptors for the predicted neoepitopes. These high-avidity immune responses were polyclonal, mutation-specific and restricted to either HLA class I or II. Interestingly, we observed the persistence and expansion of polyclonal T-cell responses following neoadjuvant chemotherapy. CONCLUSIONS: We demonstrate our neoepitope prediction program allows for the successful identification of neoepitopes targeted by TILs in patients with BC, providing a means to identify tumor-specific immunogenic targets for individualized treatment, including vaccines or adoptively transferred cellular therapies.

Single-Cell Profiling Reveals Divergent, Globally Patterned Immune Responses in Murine Skin Inflammation.

Inflammatory response heterogeneity has impeded high-resolution dissection of diverse immune cell populations during activation. We characterize mouse cutaneous immune cells by single-cell RNA sequencing, after inducing inflammation using imiquimod and oxazolone dermatitis models. We identify 13 CD45+ subpopulations, which broadly represent most functionally characterized immune cell types. Oxazolone pervasively upregulates Jak2/Stat3 expression across T cells and antigen-presenting cells (APCs). Oxazolone also induces Il4/Il13 expression in newly infiltrating basophils, and Il4ra and Ccl24, most prominently in APCs. In contrast, imiquimod broadly upregulates Il17/Il22 and Ccl4/Ccl5. A comparative analysis of single-cell inflammatory transcriptional responses reveals that APC response to oxazolone is tightly restricted by cell identity, whereas imiquimod enforces shared programs on multiple APC populations in parallel. These global molecular patterns not only contrast immune responses on a systems level but also suggest that the mechanisms of new sources of inflammation can eventually be deduced by comparison to known signatures.

A risk-associated Active transcriptome phenotype expressed by histologically normal human breast tissue and linked to a pro-tumorigenic adipocyte population.

BACKGROUND: Previous studies have identified and validated a risk-associated Active transcriptome phenotype commonly expressed in the cancer-adjacent and histologically normal epithelium, stroma, and adipose containing peritumor microenvironment of clinically established invasive breast cancers, conferring a 2.5- to 3-fold later risk of dying from recurrent breast cancer. Expression of this Active transcriptome phenotype has not yet been evaluated in normal breast tissue samples unassociated with any benign or malignant lesions; however, it has been associated with increased peritumor adipocyte composition. METHODS: Detailed histologic and transcriptomic (RNAseq) analyses were performed on normal breast biopsy samples from 151 healthy, parous, non-obese (mean BMI = 29.60 ± 7.92) women, ages 27-66 who donated core breast biopsy samples to the Komen Tissue Bank, and whose average breast cancer risk estimate (Gail score) at the time of biopsy (1.27 ± 1.34) would not qualify them for endocrine prevention therapy. RESULTS: Full genome RNA sequencing (RNAseq) identified 52% (78/151) of these normal breast samples as expressing the Active breast phenotype. While Active signature genes were found to be most variably expressed in mammary adipocytes, donors with the Active phenotype had no difference in BMI but significantly higher Gail scores (1.46 vs. 1.18; p = 0.007). Active breast samples possessed 1.6-fold more (~ 80%) adipocyte nuclei, larger cross-sectional adipocyte areas (p < 0.01), and 0.5-fold fewer stromal and epithelial cell nuclei (p < 1e-6). Infrequent low-level expression of cancer gene hotspot mutations was detected but not enriched in the Active breast samples. Active samples were enriched in gene sets associated with adipogenesis and fat metabolism (FDR q ≤ 10%), higher signature scores for cAMP-dependent lipolysis known to drive breast cancer progression, white adipose tissue browning (Wilcoxon p < 0.01), and genes associated with adipocyte activation (leptin, adiponectin) and remodeling (CAV1, BNIP3), adipokine growth factors (IGF-1, FGF2), and pro-inflammatory fat signaling (IKBKG, CCL13). CONCLUSIONS: The risk-associated Active transcriptome phenotype first identified in cancer-adjacent breast tissues also occurs commonly in healthy women without breast disease who do not qualify for breast cancer chemoprevention, and independently of breast expressed cancer-associated mutations. The risk-associated Active phenotype appears driven by a pro-tumorigenic adipocyte microenvironment that can predate breast cancer development.

Transcriptomic silencing as a potential mechanism of treatment resistance.

Next-generation sequencing (NGS) has not revealed all the mechanisms underlying resistance to genomically matched drugs. Here, we performed in 1417 tumors whole-exome tumor (somatic)/normal (germline) NGS and whole-transcriptome sequencing, the latter focusing on a clinically oriented 50-gene panel in order to examine transcriptomic silencing of putative driver alterations. In this large-scale study, approximately 13% of the somatic single nucleotide variants (SNVs) were unexpectedly not expressed as RNA; 23% of patients had ≥1 nonexpressed SNV. SNV-bearing genes consistently transcribed were TP53, PIK3CA, and KRAS; those with lower transcription rates were ALK, CSF1R, ERBB4, FLT3, GNAS, HNF1A, KDR, PDGFRA, RET, and SMO. We also determined the frequency of tumor mutations being germline, rather than somatic, in these and an additional 462 tumors with tumor/normal exomes; 33.8% of germline SNVs within the gene panel were rare (not found after filtering through variant information domains) and at risk of being falsely reported as somatic. Both the frequency of silenced variant transcription and the risk of falsely identifying germline mutations as somatic/tumor related are important phenomena. Therefore, transcriptomics is a critical adjunct to genomics when interrogating patient tumors for actionable alterations, because, without expression of the target aberrations, there will likely be therapeutic resistance.

Overcoming hypoxia-induced functional suppression of NK cells.

BACKGROUND: Natural killer (NK) cells are immune cells capable of killing virally infected cells and tumor cells without the need for antigen stimulation. Tumors, however, can create a suppressive microenvironment that decreases NK function. A feature of many tumors is hypoxia (low oxygen perfusion), which has been previously shown to decrease NK function. A high affinity NK (haNK) cell has been engineered to express a high affinity CD16 receptor as well as internal interleukin (IL)-2 for increased antibody-dependent cellular cytotoxicity (ADCC) and activation, respectively. We sought to investigate the tolerance of NK cells versus haNK cells to hypoxia. METHODS: We exposed healthy donor (HD) NK and X-irradiated haNK cells to normoxia (20% oxygen) as well as hypoxia (0% oxygen) and investigated their ability to kill prostate, breast and lung tumor cell lines after 5 hours. We also used monoclonal antibodies cetuximab (anti-EGFR) or avelumab (antiprogrammed death-ligand 1) to investigate the effects of hypoxia on NK ADCC. Genomic and proteomic analyzes were done to determine the effect of hypoxia on the expression of factors important to NK cell function. RESULTS: While HD NK cell cytolytic abilities were markedly and significantly impaired under hypoxic conditions, haNK cells maintained killing capacity under hypoxic conditions. NK killing, serial killing and ADCC were maintained under hypoxia in haNK cells. IL-2 has been previously implicated in serial killing and perforin regeneration and thus the endogenous IL-2 produced by haNK cells is likely a driver of the maintained killing capacity of haNK cells under hypoxic conditions. Activation of signal transducer and activator of transcription 3 (STAT3) is not seen in haNKs under hypoxia but is significant in HD NK cells. Pharmaceutical activation of STAT3 in haNKs led to reduced killing, implicating active STAT3 in reduced NK cell function. CONCLUSIONS: In contrast to HD NK cells, haNK cells are resistant to acute hypoxia. The potent cytolytic function of haNK cells was maintained in an environment comparable to what would be encountered in a tumor. The data presented here provide an additional mechanism of action for haNK cells that are currently being evaluated in clinical trials for several tumor types.

A deep learning image-based intrinsic molecular subtype classifier of breast tumors reveals tumor heterogeneity that may affect survival.

BACKGROUND: Breast cancer intrinsic molecular subtype (IMS) as classified by the expression-based PAM50 assay is considered a strong prognostic feature, even when controlled for by standard clinicopathological features such as age, grade, and nodal status, yet the molecular testing required to elucidate these subtypes is not routinely performed. Furthermore, when such bulk assays as RNA sequencing are performed, intratumoral heterogeneity that may affect prognosis and therapeutic decision-making can be missed. METHODS: As a more facile and readily available method for determining IMS in breast cancer, we developed a deep learning approach for approximating PAM50 intrinsic subtyping using only whole-slide images of H&E-stained breast biopsy tissue sections. This algorithm was trained on images from 443 tumors that had previously undergone PAM50 subtyping to classify small patches of the images into four major molecular subtypes-Basal-like, HER2-enriched, Luminal A, and Luminal B-as well as Basal vs. non-Basal. The algorithm was subsequently used for subtype classification of a held-out set of 222 tumors. RESULTS: This deep learning image-based classifier correctly subtyped the majority of samples in the held-out set of tumors. However, in many cases, significant heterogeneity was observed in assigned subtypes across patches from within a single whole-slide image. We performed further analysis of heterogeneity, focusing on contrasting Luminal A and Basal-like subtypes because classifications from our deep learning algorithm-similar to PAM50-are associated with significant differences in survival between these two subtypes. Patients with tumors classified as heterogeneous were found to have survival intermediate between Luminal A and Basal patients, as well as more varied levels of hormone receptor expression patterns. CONCLUSIONS: Here, we present a method for minimizing manual work required to identify cancer-rich patches among all multiscale patches in H&E-stained WSIs that can be generalized to any indication. These results suggest that advanced deep machine learning methods that use only routinely collected whole-slide images can approximate RNA-seq-based molecular tests such as PAM50 and, importantly, may increase detection of heterogeneous tumors that may require more detailed subtype analysis.

Multiregion exome sequencing of ovarian immature teratomas reveals 2N near-diploid genomes, paucity of somatic mutations, and extensive allelic imbalances shared across mature, immature, and disseminated components.

Immature teratoma is a subtype of malignant germ cell tumor of the ovary that occurs most commonly in the first three decades of life, frequently with bilateral ovarian disease. Despite being the second most common malignant germ cell tumor of the ovary, little is known about its genetic underpinnings. Here we performed multiregion whole-exome sequencing to interrogate the genetic zygosity, clonal relationship, DNA copy number, and mutational status of 52 pathologically distinct tumor components from ten females with ovarian immature teratomas, with bilateral tumors present in five cases and peritoneal dissemination in seven cases. We found that ovarian immature teratomas are genetically characterized by 2N near-diploid genomes with extensive loss of heterozygosity and an absence of genes harboring recurrent somatic mutations or known oncogenic variants. All components within a single ovarian tumor (immature teratoma, mature teratoma with different histologic patterns of differentiation, and yolk sac tumor) were found to harbor an identical pattern of loss of heterozygosity across the genome, indicating a shared clonal origin. In contrast, the four analyzed bilateral teratomas showed distinct patterns of zygosity changes in the right versus left sided tumors, indicating independent clonal origins. All disseminated teratoma components within the peritoneum (including gliomatosis peritonei) shared a clonal pattern of loss of heterozygosity with either the right or left primary ovarian tumor. The observed genomic loss of heterozygosity patterns indicate that diverse meiotic errors contribute to the formation of ovarian immature teratomas, with 11 out of the 15 genetically distinct clones determined to result from nondisjunction errors during meiosis I or II. Overall, these findings suggest that copy-neutral loss of heterozygosity resulting from meiotic abnormalities may be sufficient to generate ovarian immature teratomas from germ cells.

Efficient Tumor Clearance and Diversified Immunity through Neoepitope Vaccines and Combinatorial Immunotherapy.

Progressive tumor growth is associated with deficits in the immunity generated against tumor antigens. Vaccines targeting tumor neoepitopes have the potential to address qualitative defects; however, additional mechanisms of immune failure may underlie tumor progression. In such cases, patients would benefit from additional immune-oncology agents targeting potential mechanisms of immune failure. This study explores the identification of neoepitopes in the MC38 colon carcinoma model by comparison of tumor to normal DNA and tumor RNA sequencing technology, as well as neoepitope delivery by both peptide- and adenovirus-based vaccination strategies. To improve antitumor efficacies, we combined the vaccine with a group of rationally selected immune-oncology agents. We utilized an IL15 superagonist to enhance the development of antigen-specific immunity initiated by the neoepitope vaccine, PD-L1 blockade to reduce tumor immunosuppression, and a tumor-targeted IL12 molecule to facilitate T-cell function within the tumor microenvironment. Analysis of tumor-infiltrating leukocytes demonstrated this multifaceted treatment regimen was required to promote the influx of CD8+ T cells and enhance the expression of transcripts relating to T-cell activation/effector function. Tumor-targeted IL12 resulted in a marked increase in clonality of T-cell repertoire infiltrating the tumor, which when sculpted with the addition of either a peptide or adenoviral neoepitope vaccine promoted efficient tumor clearance. In addition, the neoepitope vaccine induced the spread of immunity to neoepitopes expressed by the tumor but not contained within the vaccine. These results demonstrate the importance of combining neoepitope-targeting vaccines with a multifaceted treatment regimen to generate effective antitumor immunity.

Loss of ZNF750 in ocular and cutaneous sebaceous carcinoma.

BACKGROUND: Sebaceous carcinoma (SeC) is an uncommon malignancy arising from sebaceous glands of the conjunctiva and skin. Recurrent mutations in the ZNF750 were recently identified in ocular SeC. We assessed whether ZNF750 loss is a specific feature of ocular SeC or a general feature of sebaceous tumors. METHODS: Immunostaining for ZNF750 expression was performed in 54 benign and malignant sebocytic proliferations. Staining for ZNF750 was scored on a three-tier scale: positive (>75%), partially positive (5%-74%), and negative (<5%). RESULTS: ZNF750 expression was negative in 4/11 ocular SeC, and partially positive in 4/11 ocular SeC and 6/13 cutaneous SeC. No extraocular tumors were negative. No loss was found in sebaceous adenoma or sebaceous hyperplasia. In nine previously sequenced ocular SeCs, two lacked detectable somatic mutations in ZNF750, but showed complete loss of staining, indicating non-mutational inactivation of ZNF750. CONCLUSION: We show complete loss of the ZNF750 epidermal differentiation regulator in about half of ocular SeC, highlighting the most common genetic defect in this cancer type. Loss of ZNF750 expression is seen even in tumors without truncating mutations and reduced in many of the remaining ocular and cutaneous SeC. In contrast, no ZNF750 loss was detected in benign sebaceous proliferations.

A20 and ABIN1 Suppression of a Keratinocyte Inflammatory Program with a Shared Single-Cell Expression Signature in Diverse Human Rashes.

Genetic variation in the NF-κB inhibitors, ABIN1 and A20, increase risk for psoriasis. While critical for hematopoietic immune cell function, these genes are believed to additionally inhibit psoriasis by dampening inflammatory signaling in keratinocytes. We dissected ABIN1 and A20's regulatory role in human keratinocyte inflammation using an RNA sequencing-based comparative genomic approach. Here we show subsets of the IL-17 and tumor necrosis factor-α signaling pathways are robustly restricted by A20 overexpression. In contrast, ABIN1 overexpression inhibits these genes more modestly for IL-17, and weakly for tumor necrosis factor-α. Our genome-scale analysis also indicates that inflammatory program suppression appears to be the major transcriptional influence of A20/ABIN1 overexpression, without obvious influence on keratinocyte viability genes. Our findings thus enable dissection of the differing anti-inflammatory mechanisms of two distinct psoriasis modifiers, which may be directly exploited for therapeutic purposes. Importantly, we report that IL-17-induced targets of A20 show similar aberrant epidermal layer-specific transcriptional upregulation in keratinocytes from diseases as diverse as psoriasis, atopic dermatitis, and erythrokeratodermia variabilis, suggesting a contributory role for epidermal inflammation in a broad spectrum of rashes.

Transcriptional Programming of Normal and Inflamed Human Epidermis at Single-Cell Resolution.

Perturbations in the transcriptional programs specifying epidermal differentiation cause diverse skin pathologies ranging from impaired barrier function to inflammatory skin disease. However, the global scope and organization of this complex cellular program remain undefined. Here we report single-cell RNA sequencing profiles of 92,889 human epidermal cells from 9 normal and 3 inflamed skin samples. Transcriptomics-derived keratinocyte subpopulations reflect classic epidermal strata but also sharply compartmentalize epithelial functions such as cell-cell communication, inflammation, and WNT pathway modulation. In keratinocytes, ∼12% of assessed transcript expression varies in coordinate patterns, revealing undescribed gene expression programs governing epidermal homeostasis. We also identify molecular fingerprints of inflammatory skin states, including S100 activation in the interfollicular epidermis of normal scalp, enrichment of a CD1C+CD301A+ myeloid dendritic cell population in psoriatic epidermis, and IL1ßhiCCL3hiCD14+ monocyte-derived macrophages enriched in foreskin. This compendium of RNA profiles provides a critical step toward elucidating epidermal diseases of development, differentiation, and inflammation.

APOBEC mutation drives early-onset squamous cell carcinomas in recessive dystrophic epidermolysis bullosa.

Recessive dystrophic epidermolysis bullosa (RDEB) is a rare inherited skin and mucous membrane fragility disorder complicated by early-onset, highly malignant cutaneous squamous cell carcinomas (SCCs). The molecular etiology of RDEB SCC, which arises at sites of sustained tissue damage, is unknown. We performed detailed molecular analysis using whole-exome, whole-genome, and RNA sequencing of 27 RDEB SCC tumors, including multiple tumors from the same patient and multiple regions from five individual tumors. We report that driver mutations were shared with spontaneous, ultraviolet (UV) light-induced cutaneous SCC (UV SCC) and head and neck SCC (HNSCC) and did not explain the early presentation or aggressive nature of RDEB SCC. Instead, endogenous mutation processes associated with apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) deaminases dominated RDEB SCC. APOBEC mutation signatures were enhanced throughout RDEB SCC tumor evolution, relative to spontaneous UV SCC and HNSCC mutation profiles. Sixty-seven percent of RDEB SCC driver mutations was found to emerge as a result of APOBEC and other endogenous mutational processes previously associated with age, potentially explaining a >1000-fold increased incidence and the early onset of these SCCs. Human papillomavirus-negative basal and mesenchymal subtypes of HNSCC harbored enhanced APOBEC mutational signatures and transcriptomes similar to those of RDEB SCC, suggesting that APOBEC deaminases drive other subtypes of SCC. Collectively, these data establish specific mutagenic mechanisms associated with chronic tissue damage. Our findings reveal a cause for cancers arising at sites of persistent inflammation and identify potential therapeutic avenues to treat RDEB SCC.

Comprehensive genomic transcriptomic tumor-normal gene panel analysis for enhanced precision in patients with lung cancer.

A CMS approved test for lung cancer is based on tumor-only analysis of a targeted 35 gene panel, specifically excluding the use of the patient's normal germline tissue. However, this tumor-only approach increases the risk of mistakenly identifying germline single nucleotide polymorphisms (SNPs) as somatically-derived cancer driver mutations (false positives). 621 patients with 30 different cancer types, including lung cancer, were studied to compare the precision of tumor somatic variant calling in 35 genes using tumor-only DNA sequencing versus tumor-normal DNA plus RNA sequencing. When sequencing of lung cancer was performed using tumor genomes alone without normal germline controls, 94% of variants identified were SNPs and thus false positives. Filtering for common SNPs still resulted in as high as 48% false positive variant calling. With tumor-only sequencing, 29% of lung cancer patients had a false positive variant call in at least one of twelve genes with directly targetable drugs. RNA analysis showed 18% of true somatic variants were not expressed. Thus, sequencing and analysis of both normal germline and tumor genomes is necessary for accurate identification of molecular targets. Treatment decisions based on tumor-only analysis may result in the administration of ineffective therapies while also increasing the risk of negative drug-related side effects.

Cell of origin and mutation pattern define three clinically distinct classes of sebaceous carcinoma.

Sebaceous carcinomas (SeC) are cutaneous malignancies that, in rare cases, metastasize and prove fatal. Here we report whole-exome sequencing on 32 SeC, revealing distinct mutational classes that explain both cancer ontogeny and clinical course. A UV-damage signature predominates in 10/32 samples, while nine show microsatellite instability (MSI) profiles. UV-damage SeC exhibited poorly differentiated, infiltrative histopathology compared to MSI signature SeC (p = 0.003), features previously associated with dissemination. Moreover, UV-damage SeC transcriptomes and anatomic distribution closely resemble those of cutaneous squamous cell carcinomas (SCC), implicating sun-exposed keratinocytes as a cell of origin. Like SCC, this UV-damage subclass harbors a high somatic mutation burden with >50 mutations per Mb, predicting immunotherapeutic response. In contrast, ocular SeC acquires far fewer mutations without a dominant signature, but show frequent truncations in the ZNF750 epidermal differentiation regulator. Our data exemplify how different mutational processes convergently drive histopathologically related but clinically distinct cancers.

A gene signature for late distant metastasis in breast cancer identifies a potential mechanism of late recurrences.

INTRODUCTION: Breast cancer risk of recurrence is known to span 20 years, yet existing prognostic signatures are best at predicting early recurrences (≤ 5 years). There is a critical need to identify those patients at risk of late-relapse (>5 years), in order to select potential candidates for further treatment and to identify molecular targets for such treatment. METHODS: A total of 252 breast primary tumors were selected at the Netherlands Cancer Institute from a retrospective series of ER+, HER2- breast cancer patients with a follow-up of at least 10 years. Gene expression analysis was performed using Agilent 4x44K microarrays. Patients were classified in 3 groups: no relapse (M0); relapse before 5 years (M0-5) or after 5 years (M5-15). We assessed the correlation of clinico-pathological variables with late Distant Metastases (DM). We divided the patient series into a training set of untreated patients (n = 140) and a test set of treated patients (n = 112), to investigate whether a gene-signature or single genes could be identified for predicting late DM. Pathway level late DM correlates were identified using PARADIGM and DAVID. RESULTS: Of the clinico-pathologic variables tested, only lymph node status associated with late DM. A 241-gene signature developed on the NKI training set was able to classify M5-15 patients in the test set with a sensitivity of 77% and a specificity of 33% (AUC 0.654). This signature showed enrichment in genes involved in immune response and extracellular matrix. An alternative analysis of individual genes identified CH25H as an independent predictor of distant metastasis in our patient series. CONCLUSIONS: We identified a gene signature for late metastasis in breast cancer. Our data are consistent with a model in which suppressed anti-tumoral immunity enables dormant tumor cells to re-enter the cell cycle to form metastases in response to extrinsic events in the microenvironment.

Learning subgroup-specific regulatory interactions and regulator independence with PARADIGM.

UNLABELLED: High-dimensional '-omics' profiling provides a detailed molecular view of individual cancers; however, understanding the mechanisms by which tumors evade cellular defenses requires deep knowledge of the underlying cellular pathways within each cancer sample. We extended the PARADIGM algorithm (Vaske et al., 2010, Bioinformatics, 26, i237-i245), a pathway analysis method for combining multiple '-omics' data types, to learn the strength and direction of 9139 gene and protein interactions curated from the literature. Using genomic and mRNA expression data from 1936 samples in The Cancer Genome Atlas (TCGA) cohort, we learned interactions that provided support for and relative strength of 7138 (78%) of the curated links. Gene set enrichment found that genes involved in the strongest interactions were significantly enriched for transcriptional regulation, apoptosis, cell cycle regulation and response to tumor cells. Within the TCGA breast cancer cohort, we assessed different interaction strengths between breast cancer subtypes, and found interactions associated with the MYC pathway and the ER alpha network to be among the most differential between basal and luminal A subtypes. PARADIGM with the Naive Bayesian assumption produced gene activity predictions that, when clustered, found groups of patients with better separation in survival than both the original version of PARADIGM and a version without the assumption. We found that this Naive Bayes assumption was valid for the vast majority of co-regulators, indicating that most co-regulators act independently on their shared target. AVAILABILITY: http://paradigm.five3genomics.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The UCSC Interaction Browser: multidimensional data views in pathway context.

High-throughput data sets such as genome-wide protein-protein interactions, protein-DNA interactions and gene expression data have been published for several model systems, especially for human cancer samples. The University of California, Santa Cruz (UCSC) Interaction Browser (http://sysbio.soe.ucsc.edu/nets) is an online tool for biologists to view high-throughput data sets simultaneously for the analysis of functional relationships between biological entities. Users can access several public interaction networks and functional genomics data sets through the portal as well as upload their own networks and data sets for analysis. Users can navigate through correlative relationships for focused sets of genes belonging to biological pathways using a standard web browser. Using a new visual modality called the CircleMap, multiple 'omics' data sets can be viewed simultaneously within the context of curated, predicted, directed and undirected regulatory interactions. The Interaction Browser provides an integrative viewing of biological networks based on the consensus of many observations about genes and their products, which may provide new insights about normal and disease processes not obvious from any isolated data set.

Genomic aberrations in normal tissue adjacent to HER2-amplified breast cancers: field cancerization or contaminating tumor cells?

Field cancerization effects as well as isolated tumor cell foci extending well beyond the invasive tumor margin have been described previously to account for local recurrence rates following breast conserving surgery despite adequate surgical margins and breast radiotherapy. To look for evidence of possible tumor cell contamination or field cancerization by genetic effects, a pilot study (Study 1: 12 sample pairs) followed by a verification study (Study 2: 20 sample pairs) were performed on DNA extracted from HER2-positive breast tumors and matching normal adjacent mammary tissue samples excised 1-3 cm beyond the invasive tumor margin. High-resolution molecular inversion probe (MIP) arrays were used to compare genomic copy number variations, including increased HER2 gene copies, between the paired samples; as well, a detailed histologic and immunohistochemical (IHC) re-evaluation of all Study 2 samples was performed blinded to the genomic results to characterize the adjacent normal tissue composition bracketing the DNA-extracted samples. Overall, 14/32 (44 %) sample pairs from both studies produced genome-wide evidence of genetic aberrations including HER2 copy number gains within the adjacent normal tissue samples. The observed single-parental origin of monoallelic HER2 amplicon haplotypes shared by informative tumor-normal pairs, as well as commonly gained loci elsewhere on 17q, suggested the presence of contaminating tumor cells in the genomically aberrant normal samples. Histologic and IHC analyses identified occult 25-200 µm tumor cell clusters overexpressing HER2 scattered in more than half, but not all, of the genomically aberrant normal samples re-evaluated, but in none of the genomically normal samples. These genomic and microscopic findings support the conclusion that tumor cell contamination rather than genetic field cancerization represents the likeliest cause of local clinical recurrence rates following breast conserving surgery, and mandate caution in assuming the genomic normalcy of histologically benign appearing peritumor breast tissue.

Subtype and pathway specific responses to anticancer compounds in breast cancer.

Breast cancers are comprised of molecularly distinct subtypes that may respond differently to pathway-targeted therapies now under development. Collections of breast cancer cell lines mirror many of the molecular subtypes and pathways found in tumors, suggesting that treatment of cell lines with candidate therapeutic compounds can guide identification of associations between molecular subtypes, pathways, and drug response. In a test of 77 therapeutic compounds, nearly all drugs showed differential responses across these cell lines, and approximately one third showed subtype-, pathway-, and/or genomic aberration-specific responses. These observations suggest mechanisms of response and resistance and may inform efforts to develop molecular assays that predict clinical response.