A Bayesian framework to study tumor subclone-specific expression by combining bulk DNA and single-cell RNA sequencing data.

Genetic and gene expression heterogeneity is an essential hallmark of many tumors, allowing the cancer to evolve and to develop resistance to treatment. Currently, the most commonly used data types for studying such heterogeneity are bulk tumor/normal whole-genome or whole-exome sequencing (WGS, WES); and single-cell RNA sequencing (scRNA-seq), respectively. However, tools are currently lacking to link genomic tumor subclonality with transcriptomic heterogeneity by integrating genomic and single-cell transcriptomic data collected from the same tumor. To address this gap, we developed scBayes, a Bayesian probabilistic framework that uses tumor subclonal structure inferred from bulk DNA sequencing data to determine the subclonal identity of cells from single-cell gene expression (scRNA-seq) measurements. Grouping together cells representing the same genetically defined tumor subclones allows comparison of gene expression across different subclones, or investigation of gene expression changes within the same subclone across time (i.e., progression, treatment response, or relapse) or space (i.e., at multiple metastatic sites and organs). We used simulated data sets, in silico synthetic data sets, as well as biological data sets generated from cancer samples to extensively characterize and validate the performance of our method, as well as to show improvements over existing methods. We show the validity and utility of our approach by applying it to published data sets and recapitulating the findings, as well as arriving at novel insights into cancer subclonal expression behavior in our own data sets. We further show that our method is applicable to a wide range of single-cell sequencing technologies including single-cell DNA sequencing as well as Smart-seq and 10x Genomics scRNA-seq protocols.

Circulating immune cell phenotype dynamics reflect the strength of tumor-immune cell interactions in patients during immunotherapy.

The extent to which immune cell phenotypes in the peripheral blood reflect within-tumor immune activity prior to and early in cancer therapy is unclear. To address this question, we studied the population dynamics of tumor and immune cells, and immune phenotypic changes, using clinical tumor and immune cell measurements and single-cell genomic analyses. These samples were serially obtained from a cohort of advanced gastrointestinal cancer patients enrolled in a trial with chemotherapy and immunotherapy. Using an ecological population model, fitted to clinical tumor burden and immune cell abundance data from each patient, we find evidence of a strong tumor-circulating immune cell interaction in responder patients but not in those patients that progress on treatment. Upon initiation of therapy, immune cell abundance increased rapidly in responsive patients, and once the peak level is reached tumor burden decreases, similar to models of predator-prey interactions; these dynamic patterns were absent in nonresponder patients. To interrogate phenotype dynamics of circulating immune cells, we performed single-cell RNA sequencing at serial time points during treatment. These data show that peripheral immune cell phenotypes were linked to the increased strength of patients' tumor-immune cell interaction, including increased cytotoxic differentiation and strong activation of interferon signaling in peripheral T cells in responder patients. Joint modeling of clinical and genomic data highlights the interactions between tumor and immune cell populations and reveals how variation in patient responsiveness can be explained by differences in peripheral immune cell signaling and differentiation soon after the initiation of immunotherapy.

Topsentinol L Trisulfate, a Marine Natural Product That Targets Basal-like and Claudin-Low Breast Cancers.

Patients diagnosed with basal-like breast cancer suffer from poor prognosis and limited treatment options. There is an urgent need to identify new targets that can benefit patients with basal-like and claudin-low (BL-CL) breast cancers. We screened fractions from our Marine Invertebrate Compound Library (MICL) to identify compounds that specifically target BL-CL breast cancers. We identified a previously unreported trisulfated sterol, i.e., topsentinol L trisulfate (TLT), which exhibited increased efficacy against BL-CL breast cancers relative to luminal/HER2+ breast cancer. Biochemical investigation of the effects of TLT on BL-CL cell lines revealed its ability to inhibit activation of AMP-activated protein kinase (AMPK) and checkpoint kinase 1 (CHK1) and to promote activation of p38. The importance of targeting AMPK and CHK1 in BL-CL cell lines was validated by treating a panel of breast cancer cell lines with known small molecule inhibitors of AMPK (dorsomorphin) and CHK1 (Ly2603618) and recording the increased effectiveness against BL-CL breast cancers as compared with luminal/HER2+ breast cancer. Finally, we generated a drug response gene-expression signature and projected it against a human tumor panel of 12 different cancer types to identify other cancer types sensitive to the compound. The TLT sensitivity gene-expression signature identified breast and bladder cancer as the most sensitive to TLT, while glioblastoma multiforme was the least sensitive.

Exploiting collateral sensitivity controls growth of mixed culture of sensitive and resistant cells and decreases selection for resistant cells in a cell line model.

BACKGROUND: CDK4/6 inhibitors such as ribociclib are becoming widely used targeted therapies in hormone-receptor-positive (HR+) human epidermal growth factor receptor 2-negative (HER2-) breast cancer. However, cancers can advance due to drug resistance, a problem in which tumor heterogeneity and evolution are key features. METHODS: Ribociclib-resistant HR+/HER2- CAMA-1 breast cancer cells were generated through long-term ribociclib treatment. Characterization of sensitive and resistant cells were performed using RNA sequencing and whole exome sequencing. Lentiviral labeling with different fluorescent proteins enabled us to track the proliferation of sensitive and resistant cells under different treatments in a heterogeneous, 3D spheroid coculture system using imaging microscopy and flow cytometry. RESULTS: Transcriptional profiling of sensitive and resistant cells revealed the downregulation of the G2/M checkpoint in the resistant cells. Exploiting this acquired vulnerability; resistant cells exhibited collateral sensitivity for the Wee-1 inhibitor, adavosertib (AZD1775). The combination of ribociclib and adavosertib achieved additional antiproliferative effect exclusively in the cocultures compared to monocultures, while decreasing the selection for resistant cells. CONCLUSIONS: Our results suggest that optimal antiproliferative effects in heterogeneous cancers can be achieved via an integrative therapeutic approach targeting sensitive and resistant cancer cell populations within a tumor, respectively.

Time- and dose-dependent gene expression analysis of macrophage response as a function of porosity of silica nanoparticles.

There is a limited amount of information available on gene expression regulation of macrophages in response to changing the time of exposure, concentration, and physicochemical properties of nanomaterials. In this study, RAW264.7 macrophages were treated with spherical nonporous and mesoporous silica nanoparticles of similar size at different incubation times and concentrations. RNA-sequencing was used to study transcriptional profiles. Bioinformatics analyses, functional annotation clustering, and network analyses were employed to understand signaling pathways of cellular response as a function of porosity, incubation time, and concentration. Porosity introduced drastic changes to the genomic response of macrophages at equitoxic concentrations and incubation times. Direct relations between increases in time and concentration with an increased number of differentially expressed genes were observed.

Genotoxicity of amorphous silica nanoparticles: Status and prospects.

Amorphous silica nanoparticles (SNPs) are widely used in biomedical applications and consumer products. Little is known, however, about their genotoxicity and potential to induce gene expression regulation. Despite recent efforts to study the underlying mechanisms of genotoxicity of SNPs, inconsistent results create a challenge. A variety of factors determine particle-cell interactions and underlying mechanisms. Further, high-throughput studies are required to carefully assess the impact of silica nanoparticle physicochemical properties on induction of genotoxic response in different cell lines and animal models. In this article, we review the strategies available for evaluation of genotoxicity of nanoparticles (NPs), survey current status of silica nanoparticle gene alteration and genotoxicity, discuss particle-mediated inflammation as a contributing factor to genotoxicity, identify existing gaps and suggest future directions for this research.

Differential gene expression patterns in vein regions susceptible versus resistant to neointimal hyperplasia.

Arteriovenous hemodialysis graft (AVG) stenosis results in thrombosis and AVG failure, but prevention of stenosis has been unsuccessful due in large part to our limited understanding of the molecular processes involved in neointimal hyperplasia (NH) formation. AVG stenosis develops chiefly as a consequence of highly localized NH formation in the vein-graft anastomosis region. Surprisingly, the vein region just downstream of the vein-graft anastomosis (herein termed proximal vein region) is relatively resistant to NH. We hypothesized that the gene expression profiles of the NH-prone and NH-resistant regions will be different from each other after graft placement, and analysis of their genomic profiles may yield potential therapeutic targets to prevent AVG stenosis. To test this, we evaluated the vein-graft anastomosis (NH-prone) and proximal vein (NH-resistant) regions in a porcine model of AVG stenosis with a porcine microarray. Gene expression changes in these two distinct vein regions, relative to the gene expression in unoperated control veins, were examined at early (5 days) and later (14 days) time points following graft placement. Global genomic changes were much greater in the NH-prone region than in the NH-resistant region at both time points. In the NH-prone region, genes related to regulation of cell proliferation and osteo-/chondrogenic vascular remodeling were most enriched among the significantly upregulated genes, and genes related to smooth muscle phenotype were significantly downregulated. These results provide insights into the spatial and temporal genomic modulation underlying NH formation in AVG and suggest potential therapeutic strategies to prevent and/or limit AVG stenosis.

Global gene expression analysis of macrophage response induced by nonporous and porous silica nanoparticles.

Little is known about the global gene expression profile of macrophages in response to changes in size and porosity of silica nanoparticles (SNPs). Spherical nonporous SNPs of two different diameters, and mesoporous spherical SNPs with comparable size were characterized. Reactive oxygen species, mitochondrial membrane potential, lysosome degradation capacity, and lysosome pH were measured to evaluate the influence of nonporous and mesoporous SNPs on mitochondrial and lysosomal function. RNA-sequencing was utilized to generate transcriptional profiles of RAW264.7 macrophages exposed to non-toxic SNP doses. DESeq2, limma, and BinReg2 software were used to analyze the data based on both unsupervised and supervised strategies to identify genes with greatest differences among NP treatments. Utilizing GATHER and DAVID software, possible induced pathways were studied. We found that mesoporous silica nanoparticles are capable of altering gene expression in macrophages at doses that do not elicit acute cytotoxicity, while gene transcription was minimally affected by nonporous SNPs.

Integrative analyses reveal signaling pathways underlying familial breast cancer susceptibility.

The signaling events that drive familial breast cancer (FBC) risk remain poorly understood. While the majority of genomic studies have focused on genetic risk variants, known risk variants account for at most 30% of FBC cases. Considering that multiple genes may influence FBC risk, we hypothesized that a pathway-based strategy examining different data types from multiple tissues could elucidate the biological basis for FBC. In this study, we performed integrated analyses of gene expression and exome-sequencing data from peripheral blood mononuclear cells and showed that cell adhesion pathways are significantly and consistently dysregulated in women who develop FBC. The dysregulation of cell adhesion pathways in high-risk women was also identified by pathway-based profiling applied to normal breast tissue data from two independent cohorts. The results of our genomic analyses were validated in normal primary mammary epithelial cells from high-risk and control women, using cell-based functional assays, drug-response assays, fluorescence microscopy, and Western blotting assays. Both genomic and cell-based experiments indicate that cell-cell and cell-extracellular matrix adhesion processes seem to be disrupted in non-malignant cells of women at high risk for FBC and suggest a potential role for these processes in FBC development.

Single Cell Analysis of Peripheral TB-Associated Granulomatous Lymphadenitis.

We successfully employed a single cell RNA sequencing (scRNA-seq) approach to describe the cells and the communication networks characterizing granulomatous lymph nodes of TB patients. When mapping cells from individual patient samples, clustered based on their transcriptome similarities, we uniformly identify several cell types that known to characterize human and non-human primate granulomas. Whether high or low Mtb burden, we find the T cell cluster to be one of the most abundant. Many cells expressing T cell markers are clearly quantifiable within this CD3 expressing cluster. Other cell clusters that are uniformly detected, but that vary dramatically in abundance amongst the individual patient samples, are the B cell, plasma cell and macrophage/dendrocyte and NK cell clusters. When we combine all our scRNA-seq data from our current 23 patients (in order to add power to cell cluster identification in patient samples with fewer cells), we distinguish T, macrophage, dendrocyte and plasma cell subclusters, each with distinct signaling activities. The sizes of these subclusters also varies dramatically amongst the individual patients. In comparing FNA composition we noted trends in which T cell populations and macrophage/dendrocyte populations were negatively correlated with NK cell populations. In addition, we also discovered that the scRNA-seq pipeline, designed for quantification of human cell mRNA, also detects Mtb RNA transcripts and associates them with their host cell's transcriptome, thus identifying individual infected cells. We hypothesize that the number of detected bacterial transcript reads provides a measure of Mtb burden, as does the number of Mtb-infected cells. The number of infected cells also varies dramatically in abundance amongst the patient samples. CellChat analysis identified predominating signaling pathways amongst the cells comprising the various granulomas, including many interactions between stromal or endothelial cells and the other component cells, such as Collagen, FN1 and Laminin,. In addition, other more selective communications pathways, including MIF, MHC-1, MHC-2, APP, CD 22, CD45, and others, are identified as originating or being received by individual immune cell components.

Transcriptional responses of human aortic endothelial cells to nanoconstructs used in biomedical applications.

Understanding the potential toxicities of manufactured nanoconstructs used for drug delivery and biomedical applications may help improve their safety. We sought to determine if surface-modified silica nanoparticles and poly(amido amine) dendrimers elicit genotoxic responses on vascular endothelial cells. The nanoconstructs utilized in this study had a distinct geometry (spheres vs worms) and surface charge, which were used to evaluate the contributions of these parameters to any potential adverse effects of these materials. Time-dependent cytotoxicity was found for surfaced-functionalized but geometrically distinct silica materials, while amine-terminated dendrimers displayed time-independent cytotoxicity and carboxylated dendrimers were nontoxic in our assays. Transcriptomic evaluation of human aortic endothelial cell (HAEC) responses indicated time-dependent gene induction following silica exposure, consisting of cell cycle gene repression and pro-inflammatory gene induction. However, the dendrimers did not induce genomic toxicity, despite displaying general cytotoxicity.

A chimeric p53 evades mutant p53 transdominant inhibition in cancer cells.

Because of the dominant negative effect of mutant p53, there has been limited success with wild-type (wt) p53 cancer gene therapy. Therefore, an alternative oligomerization domain for p53 was investigated to enhance the utility of p53 for gene therapy. The tetramerization domain of p53 was substituted with the coiled-coil (CC) domain from Bcr (breakpoint cluster region). Our p53 variant (p53-CC) maintains proper nuclear localization in breast cancer cells detected via fluorescence microscopy and shows a similar expression profile of p53 target genes as wt-p53. Additionally, similar tumor suppressor activities of p53-CC and wt-p53 were detected by terminal deoxynucleotidyl transferase dUTP nick end labeling (TUNEL), annexin-V, 7-aminoactinomycin D (7-AAD), and colony-forming assays. Furthermore, p53-CC was found to cause apoptosis in four different cancer cell lines, regardless of endogenous p53 status. Interestingly, the transcriptional activity of p53-CC was higher than wt-p53 in 3 different reporter gene assays. We hypothesized that the higher transcriptional activity of p53-CC over wt-p53 was due to the sequestration of wt-p53 by endogenous mutant p53 found in cancer cells. Co-immunoprecipitation revealed that wt-p53 does indeed interact with endogenous mutant p53 via its tetramerization domain, while p53-CC escapes this interaction. Therefore, we investigated the impact of the presence of a transdominant mutant p53 on tumor suppressor activities of wt-p53 and p53-CC. Overexpression of a potent mutant p53 along with wt-p53 or p53-CC revealed that, unlike wt-p53, p53-CC retains the same level of tumor suppressor activity. Finally, viral transduction of wt-p53 and p53-CC into a breast cancer cell line that harbors a tumor derived transdominant mutant p53 validated that p53-CC indeed evades sequestration and consequent transdominant inhibition by endogenous mutant p53.

Inhibition of TRPA1, Endoplasmic Reticulum Stress, Human Airway Epithelial Cell Damage, and Ectopic MUC5AC Expression by Vasaka (Adhatoda vasica; Malabar Nut) Tea.

This study tested whether a medicinal plant, Vasaka, typically consumed as a tea to treat respiratory malaise, could protect airway epithelial cells (AECs) from wood smoke particle-induced damage and prevent pathological mucus production. Wood/biomass smoke is a pneumotoxic air pollutant. Mucus normally protects the airways, but excessive production can obstruct airflow and cause respiratory distress. Vasaka tea pre- and co-treatment dose-dependently inhibited mucin 5AC (MUC5AC) mRNA induction by AECs treated with wood smoke particles. This correlated with transient receptor potential ankyrin-1 (TRPA1) inhibition, an attenuation of endoplasmic reticulum (ER) stress, and AEC damage/death. Induction of mRNA for anterior gradient 2, an ER chaperone/disulfide isomerase required for MUC5AC production, and TRP vanilloid-3, a gene that suppresses ER stress and wood smoke particle-induced cell death, was also attenuated. Variable inhibition of TRPA1, ER stress, and MUC5AC mRNA induction was observed using selected chemicals identified in Vasaka tea including vasicine, vasicinone, apigenin, vitexin, isovitexin, isoorientin, 9-oxoODE, and 9,10-EpOME. Apigenin and 9,10-EpOME were the most cytoprotective and mucosuppressive. Cytochrome P450 1A1 (CYP1A1) mRNA was also induced by Vasaka tea and wood smoke particles. Inhibition of CYP1A1 enhanced ER stress and MUC5AC mRNA expression, suggesting a possible role in producing protective oxylipins in stressed cells. The results provide mechanistic insights and support for the purported benefits of Vasaka tea in treating lung inflammatory conditions, raising the possibility of further development as a preventative and/or restorative therapy.

Major differences among chemopreventive organoselenocompounds in the sustained elevation of cytoprotective genes.

Cytoprotective enzyme elevation through the nuclear erythroid 2 p45-related factor 2 (Nrf2)-Kelch-like ECH-associated protein 1/antioxidant response element pathway has been promulgated for cancer prevention. This study compares the redox insult and sustained cytoprotective enzyme elevation by organoselenocompounds and sulforaphane (SF) in lung cells. SF elicited a rise in reactive oxygen species (ROS) and drop in glutathione (GSH) at 2 h; nuclear accumulation of Nrf2 at 4 h; and a GSH rebound and elevation in NAD(P)H quinone oxidoreductase (NQO1), thioredoxin reductase (TR1), and glutamate-cysteine ligase (GCL) at 24 h. Selenocystine (SECY) elicited a similar 24 h response, despite lesser earlier time-point changes. 2-Cyclohexylselenazolidine-4-carboxylic acid effects were similar to SECY's but with a larger Nrf2 change and the largest 24 h increase in GSH, GCL, TR1, and NQO1 of any compound investigated. Selenomethionine elicited a similar acute rise in ROS, but lesser depletion of GSH, no 4 h increase in nuclear Nrf2, only minor 24 h elevations in TR1 and NQO1, and a GCL elevation insufficient to elevate GSH.

Thioredoxin Reductase 1 Modulates Pigmentation and Photobiology of Murine Melanocytes in vivo.

Pigment-producing melanocytes overcome frequent oxidative stress in their physiological role of protecting the skin against the deleterious effects of solar UV irradiation. This is accomplished by the activity of several endogenous antioxidant systems, including the thioredoxin antioxidant system, in which thioredoxin reductase 1 (TR1) plays an important part. To determine whether TR1 contributes to the redox regulation of melanocyte homeostasis, we have generated a selective melanocytic Txnrd1-knockout mouse model (Txnrd1mel‒/‒), which exhibits a depigmentation phenotype consisting of variable amelanotic ventral spotting and reduced pigmentation on the extremities (tail tip, ears, and paws). The antioxidant role of TR1 was further probed in the presence of acute neonatal UVB irradiation, which stimulates melanocyte activation and introduces a spike in oxidative stress in the skin microenvironment. Interestingly, we observed a significant reduction in overall melanocyte count and proliferation in the absence of TR1. Furthermore, melanocytes exhibited an elevated level of UV-induced DNA damage in the form of 8-oxo-2'-deoxyguanosine after acute UVB treatment. We also saw an engagement of compensatory antioxidant mechanisms through increased nuclear localization of transcription factor NRF2. Altogether, these data indicate that melanocytic TR1 positively regulates melanocyte homeostasis and pigmentation during development and protects against UVB-induced DNA damage and oxidative stress.

MITF Is Regulated by Redox Signals Controlled by the Selenoprotein Thioredoxin Reductase 1.

TR1 and other selenoproteins have paradoxical effects in melanocytes and melanomas. Increasing selenoprotein activity with supplemental selenium in a mouse model of UV-induced melanoma prevents oxidative damage to melanocytes and delays melanoma tumor formation. However, TR1 itself is positively associated with progression in human melanomas and facilitates metastasis in melanoma xenografts. Here, we report that melanocytes expressing a microRNA directed against TR1 (TR1low) grow more slowly than control cell lines and contain significantly less melanin. This phenotype is associated with lower tyrosinase (TYR) activity and reduced transcription of tyrosinase-like protein-1 (TYRP1). Melanoma cells in which the TR1 gene (TXNRD1) was disrupted using Crispr/Cas9 showed more dramatic effects including the complete loss of the melanocyte-specific isoform of MITF; other MITF isoforms were unaffected. We provide evidence that TR1 depletion results in oxidation of MITF itself. This newly discovered mechanism for redox modification of MITF has profound implications for controlling both pigmentation and tumorigenesis in cells of the melanocyte lineage.

A `one-two punch' therapy strategy to target chemoresistance in estrogen receptor positive breast cancer.

Cancer cell phenotypes evolve during a tumor's treatment. In some cases, tumor cells acquire cancer stem cell-like (CSL) traits such as resistance to chemotherapy and diminished differentiation; therefore, targeting these cells may be therapeutically beneficial. In this study we show that in progressive estrogen receptor positive (ER+) metastatic breast cancer tumors, resistant subclones that emerge following chemotherapy have increased CSL abundance. Further, in vitro organoid growth of ER+ patient cancer cells also shows that chemotherapy treatment leads to increased abundance of ALDH+/CD44+ CSL cells. Chemotherapy induced CSL abundance is blocked by treatment with a pan-HDAC inhibitor, belinostat. Belinostat treatment diminished both mammosphere formation and size following chemotherapy, indicating a decrease in progenitor CSL traits. HDAC inhibitors specific to class IIa (HDAC4, HDAC5) and IIb (HDAC6) were shown to primarily reverse the chemo-resistant CSL state. Single-cell RNA sequencing analysis with patient samples showed that HDAC targets and MYC signaling were promoted by chemotherapy and inhibited upon HDAC inhibitor treatment. In summary, HDAC inhibition can block chemotherapy-induced drug resistant phenotypes with 'one-two punch' strategy in refractory breast cancer cells.

ZnO particulate matter requires cell contact for toxicity in human colon cancer cells.

There is ongoing concern regarding the toxicity of nanoparticles with sizes less than 100 nm as compared to larger particles of the same nominal substance. Two commercial ZnO types, one sold as a 8-10 nm powder and the other described as -325 mesh (<44 mum) powder, were evaluated in human colon-derived RKO cells. The powders had a volume-to-surface area ratio equivalent to 40 and 330 nm spheres, respectively. Both materials formed micrometer-sized agglomerates in cell culture media. The nanosized ZnO was more cytotoxic than the micrometer-sized ZnO with LC(50) values of 15 +/- 1 and 29 +/- 4 mug/cm(2), respectively. Transfer of Zn from the solid phase to the cell culture media in the presence of RKO cells was time- and concentration-dependent. However, direct particle-cell contact was required for RKO cell cytotoxicity, and the toxicity of particles was independent of the amount of soluble Zn in the cell culture media. The mechanism of cell death includes the disruption of mitochondrial function. Robust markers of apoptosis, Annexin V staining, loss of mitochondrial potential, and increased generation of superoxide were observed when cells were treated with ZnO particulate matter but not when treated with comparable concentration of a soluble Zn salt. Both ZnO samples induced similar mechanisms of toxicity, but there was a statistically significant increase in potency per unit mass with the smaller particles.

Thioredoxin reductase 1 ablation sensitizes colon cancer cells to methylseleninate-mediated cytotoxicity.

The relationship between selenium and cancer is complex because individuals with low serum selenium levels benefit from selenium supplementation, but those with high serum selenium levels are at increased risk for other diseases. This suggests that the use of selenocompounds might be limited to particular circumstances, such as adjuvant therapy. A contributor to this dichotomy may be the activity of certain selenium containing enzymes like the cytosolic thioredoxin reductase (TR1). We evaluated the cellular response to select selenocompounds that have anticancer activity when TR1 was attenuated by siRNA in RKO colon cancer cells. Methylseleninic acid (MSA), which is a substrate for TR1, enhanced cytotoxicity to colon cancer cells when TR1 was attenuated. MSA induced stress in the endoplasmic reticulum, as measured by GRP78 protein levels. However, this pathway did not appear to account for the change in cytotoxicity when TR1 was attenuated. Instead, knockdown of the cytosolic TR plus incubation with MSA increased autophagy, as measured by LC3B cleavage, and apoptosis, as measured by Annexin V and mitochondrial dysfunction. Therefore, the use of selenocompounds with anticancer activity, like MSA, might be utilized most effectively with agents that targets TR1 in chemotherapeutic applications.

Differential gene expression in primary human skin keratinocytes and fibroblasts in response to ionizing radiation.

Although skin is usually exposed during human exposures to ionizing radiation, there have been no thorough examinations of the transcriptional response of skin fibroblasts and keratinocytes to radiation. The transcriptional response of quiescent primary fibroblasts and keratinocytes exposed to from 10 cGy to 5 Gy and collected 4 h after treatment was examined. RNA was isolated and examined by microarray analysis for changes in the levels of gene expression. Exposure to ionizing radiation altered the expression of 279 genes across both cell types. Changes in RNA expression could be arranged into three main categories: (1) changes in keratinocytes but not in fibroblasts, (2) changes in fibroblasts but not in keratinocytes, and (3) changes in both. All of these changes were primarily of p53 target genes. Similar radiation-induced changes were induced in immortalized fibroblasts or keratinocytes. In separate experiments, protein was collected and analyzed by Western blotting for expression of proteins observed in microarray experiments to be overexpressed at the mRNA level. Both Q-PCR and Western blot analysis experiments validated these transcription changes. Our results are consistent with changes in the expression of p53 target genes as indicating the magnitude of cell responses to ionizing radiation.