Dual targeting of MAPK and PI3K pathways unlocks redifferentiation of Braf-mutated thyroid cancer organoids

Thyroid cancer is the most common endocrine malignancy and several genetic events have been described to promote the development of thyroid carcinogenesis. Besides the effects of specific mutations on thyroid cancer development, the molecular mechanisms controlling tumorigenesis, tumor behavior, and drug resistance are still largely unknown. Cancer organoids have been proposed as a powerful tool to study aspects related to tumor development and progression and appear promising to test individual responses to therapies. Here, using mESC-derived thyroid organoids, we developed a BrafV637E-inducible model able to recapitulate the features of papillary thyroid cancer in vitro. Overexpression of the murine BrafV637E mutation, equivalent to BrafV600E in humans, rapidly triggers to MAPK activation, cell dedifferentiation, and disruption of follicular organization. BrafV637E-expressing organoids show a transcriptomic signature for p53, focal adhesion, ECM-receptor interactions, EMT, and inflammatory signaling pathways. Finally, PTC-like thyroid organoids were used for drug screening assays. The combination of MAPK and PI3K inhibitors reversed BrafV637E oncogene-promoted cell dedifferentiation while restoring thyroid follicle organization and function in vitro. Our results demonstrate that pluripotent stem cells-derived thyroid cancer organoids can mimic tumor development and features while providing an efficient tool for testing novel targeted therapies.

Transplantable human thyroid organoids generated from embryonic stem cells to rescue hypothyroidism

The thyroid gland captures iodide in order to synthesize hormones that act on almost all tissues and are essential for normal growth and metabolism. Low plasma levels of thyroid hormones lead to hypothyroidism, which is one of the most common disorder in humans and is not always satisfactorily treated by lifelong hormone replacement. Therefore, in addition to the lack of in vitro tractable models to study human thyroid development, differentiation and maturation, functional human thyroid organoids could pave the way to explore new therapeutic approaches. Here we report the generation of transplantable thyroid organoids derived from human embryonic stem cells capable of restoring plasma thyroid hormone in athyreotic mice as a proof of concept for future therapeutic development.

Thyroid cancer under the scope of emerging technologies

The vast majority of thyroid cancers originate from follicular cells. We outline outstanding issues at each step along the path of cancer patient care, from prevention to post-treatment follow-up and highlight how emerging technologies will help address them in the coming years. Three directions will dominate the coming technological landscape. Genomics will reveal tumoral evolutionary history and shed light on how these cancers arise from the normal epithelium and the genomics alteration driving their progression. Transcriptomics will gain cellular and spatial resolution providing a full account of intra-tumor heterogeneity and opening a window on the microenvironment supporting thyroid tumor growth. Artificial intelligence will set morphological analysis on an objective quantitative ground laying the foundations of a systematic thyroid tumor classification system. It will also integrate into unified representations the molecular and morphological perspectives on thyroid cancer.

Dynamic Cancer Cell Heterogeneity: Diagnostic and Therapeutic Implications

Though heterogeneity of cancers is recognized and has been much discussed in recent years, the concept often remains overlooked in different routine examinations. Indeed, in clinical or biological articles, reviews, and textbooks, cancers and cancer cells are generally presented as evolving distinct entities rather than as an independent heterogeneous cooperative cell population with its self-oriented biology. There are, therefore, conceptual gaps which can mislead the interpretations/diagnostic and therapeutic approaches. In this short review, we wish to summarize and discuss various aspects of this dynamic evolving heterogeneity and its biological, pathological, clinical, diagnostic, and therapeutic implications, using thyroid carcinoma as an illustrative example.

Fluorescent Multiplex Immunohistochemistry Coupled With Other State-Of-The-Art Techniques to Systematically Characterize the Tumor Immune Microenvironment

Our expanding knowledge of the interactions between tumor cells and their microenvironment has helped to revolutionize cancer treatments, including the more recent development of immunotherapies. Immune cells are an important component of the tumor microenvironment that influence progression and treatment responses, particularly to the new immunotherapies. Technological advances that help to decipher the complexity and diversity of the tumor immune microenvironment (TIME) are increasingly used in translational research and biomarker studies. Current techniques that facilitate TIME evaluation include flow cytometry, multiplex bead-based immunoassays, chromogenic immunohistochemistry (IHC), fluorescent multiplex IHC, immunofluorescence, and spatial transcriptomics. This article offers an overview of our representative data, discusses the application of each approach to studies of the TIME, including their advantages and challenges, and reviews the potential clinical applications. Flow cytometry and chromogenic and fluorescent multiplex IHC were used to immune profile a HER2+ breast cancer, illustrating some points. Spatial transcriptomic analysis of a luminal B breast tumor demonstrated that important additional insight can be gained from this new technique. Finally, the development of a multiplex panel to identify proliferating B cells, Tfh, and Tfr cells on the same tissue section demonstrates their co-localization in tertiary lymphoid structures.

Single-cell transcriptome analysis reveals thyrocyte diversity in the zebrafish thyroid gland

The thyroid gland regulates growth and metabolism via production of thyroid hormone in follicles composed of thyrocytes. So far, thyrocytes have been assumed to be a homogenous population. To uncover heterogeneity in the thyrocyte population and molecularly characterize the non-thyrocyte cells surrounding the follicle, we developed a single-cell transcriptome atlas of the region containing the zebrafish thyroid gland. The 6249-cell atlas includes profiles of thyrocytes, blood vessels, lymphatic vessels, immune cells, and fibroblasts. Further, the thyrocytes show expression heterogeneity, including bimodal expression of the transcription factor pax2a. To validate thyrocyte heterogeneity, we generated a CRISPR/Cas9-based pax2a knock-in line that monitors pax2a expression in the thyrocytes. A population of pax2a-low mature thyrocytes interspersed in individual follicles can be distinguished. We corroborate heterogeneity within the thyrocyte population using RNA sequencing of pax2a-high and pax2a-low thyrocytes, which demonstrates 20% differential expression in transcriptome between the two subpopulations. Our results identify and validate transcriptional differences within the presumed homogenous thyrocyte population.

Dissecting the Role of Thyrotropin in the DNA Damage Response in Human Thyrocytes after 131I, γ Radiation and H2O2

BACKGROUND

The early molecular events in human thyrocytes after 131I exposure have not yet been unravelled. Therefore, we investigated the role of TSH in the 131I-induced DNA damage response and gene expression in primary cultured human thyrocytes.

METHODS

Following exposure of thyrocytes, in the presence or absence of TSH, to 131I (β radiation), γ radiation (3 Gy), and hydrogen peroxide (H2O2), we assessed DNA damage, proliferation, and cell-cycle status. We conducted RNA sequencing to profile gene expression after each type of exposure and evaluated the influence of TSH on each transcriptomic response.

RESULTS

Overall, the thyrocyte responses following exposure to β or γ radiation and to H2O2 were similar. However, TSH increased 131I-induced DNA damage, an effect partially diminished after iodide uptake inhibition. Specifically, TSH increased the number of DNA double-strand breaks in nonexposed thyrocytes and thus predisposed them to greater damage following 131I exposure. This effect most likely occurred via Gα q cascade and a rise in intracellular reactive oxygen species (ROS) levels. β and γ radiation prolonged thyroid cell-cycle arrest to a similar extent without sign of apoptosis. The gene expression profiles of thyrocytes exposed to β/γ radiation or H2O2 were overlapping. Modulations in genes involved in inflammatory response, apoptosis, and proliferation were observed. TSH increased the number and intensity of modulation of differentially expressed genes after 131I exposure.

CONCLUSIONS

TSH specifically increased 131I-induced DNA damage probably via a rise in ROS levels and produced a more prominent transcriptomic response after exposure to 131I.

BRCA gene mutations do not shape the extent and organization of tumor infiltrating lymphocytes in triple negative breast cancer

This study investigated the prevalence of TIL subpopulations, TLS, PD-1 and PD-L1 in tumors from TNBC patients harboring wild-type or mutated BRCA1 or BRCA2 germline genes. This TNBC cohort included 85% TIL-positive (≥10%) tumors with 21% classified as TIL^hi (≥50%). Interestingly, the BRCA^mut group had a significantly higher incidence of TIL^pos tumors compared to the BRCA^wt group (P = 0.037). T cells were dominant in the infiltrate but no statistically significant differences were detected between BRCA^wt and BRCA^mut for CD3⁺, CD4⁺ and CD8⁺ T cells or CD20⁺ B cells. TLS were detected in 74% of tumors but again no significant differences between the BRCA groups. PD-1 expression was observed in 33% and PD-L1 in 53% (any cell, cut-off ≥1%) tumors for the entire TNBC cohort. PD-1 expression correlated with PD-L1 and both with TIL and TLS but was not associated with BRCA mutational status. Our analyses reveal that BRCA^wt and BRCA^mut TNBC are similar except for a significant increase of TIL^pos tumors in the BRCA^mut group. While BRCA gene mutations may not directly drive immune infiltration, the greater number of TIL^pos tumors could signal greater immunogenicity in this group.

Factors contributing to the resistance of the thyrocyte to hydrogen peroxide

We studied the mechanism that may explain the relative resistance of thyrocytes to H₂O₂ compared to other cell types. Ability to degrade H₂O₂, glutathione peroxidase (GPx) activity, heme oxygenase-1 (HO-1) expression, cell survival and capacity to repair DNA damage after H₂O₂ exposure or irradiation were measured in human thyrocytes in primary culture and compared to the values obtained in human T-cells and different cell lines. Compared to other cell types, thyrocytes presented a low mortality rate after H₂O₂ exposure, rapidly degraded extracellular H₂O₂ and presented a high basal seleno-dependent GPx activity. Only in thyrocytes, H₂O₂ up-regulated GPx activity and expression of HO-1 mRNA. These effects were not reproduced by irradiation. DNA damage caused by H₂O₂ was more slowly repaired than that caused by irradiation and not repaired at all in T-cells. Our study demonstrates that the thyrocyte has specific protective mechanisms against H₂O₂ and its mutagenic effects.

Abstract 3801: Functional relevance of A-to-I RNA editing on the immune response in breast cancer

INTRODUCTION: RNA editing is a post-transcriptional modification that changes double stranded RNA sequences (dsRNA). The most common way of editing in humans is the Adenosine-to-Inosine type, which is catalyzed by members of the ADAR enzymes family. A-to-I RNA editing has been shown to be essential for normal tissue development but also plays a crucial role in many cancers. We have notably shown that in breast cancer the editing frequency is increased in tumors as compared to normal tissue and correlates with ADAR expression. We also showed that type I IFN response and ADAR DNA copy number together explain 53% of ADAR expression in breast cancers. Moreover, the implication of ADAR-mediated RNA editing in preventing the sensing of endogenous dsRNA leading to inhibition of the immune response has been recently reported. ADAR-mediated RNA editing may play a crucial role in tumor progression by modulating immune response through the dsRNA sensing pathway. Here, we aim to assess the influence of ADAR-mediated RNA editing in the sensing of endogenous dsRNA and the immune response in breast cancer.

MATERIAL & METHODS: ADAR expression was modulated in MDAMB231 breast cancer cell line using shRNA lentiviral transduction and lentiviral plasmid transduction for ADAR inhibition and ADAR wild-type and editing-deficient mutant overexpression respectively. ADAR expression was assessed by RTqPCR/Western blotting. dsRNA level was evaluated by immunofluorescence using the dsRNA-specific J2 antibody (Scicons).

RESULTS: Poly I:C transfection, mimicking dsRNAs, as well as IFN Type I exposure induced the immune response in BC cell lines, in particular it increases the expression of RIG-I, MDA5 and IRF7 involved in the dsRNA sensing pathway. DNMTi treatment, known as an inducer of dsRNA level, also activated cytosolic dsRNA response in these cells in a dose dependent manner. In order to assess the role of ADAR-mediated RNA editing in the activation of the dsRNA sensing pathway, we exposed MDAMB231 cell lines in which ADAR expression was modulated to DNMTi. A wide range of doses of DNMTi, from 0.3 to 50μM, and time of exposure, from 3h to 72h, were tested. The cytosolic dsRNA pathway was induced in a dose dependent manner after DNMTi exposure. However, no effect of ADAR modulation on the activation of the pathway was observed. No difference in dsRNA levels was observed after activation of the immune pathway by DNMTi or IFN exposure, nor after ADAR modulation.

CONCLUSION & PERSPECTIVES: dsRNA pathway activation by DNMTi, poly I:C transfection or IFN Type I treatment is not dependent of ADAR in MDAMB231 cell lines. Moreover, in contrast to data in colorectal cell lines, our data suggest that DNMTi exposure does not increase cytosolic dsRNA level.

Citation Format: Floriane Dupont, Florian Clatot, David Venet, Véronique Kruys, Debora Fumagalli, Vincent Detours, Françoise Rothé, Christos Sotiriou. Functional relevance of A-to-I RNA editing on the immune response in breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3801.

Genomic hotspots but few recurrent fusion genes in breast cancer

The advent of next generation sequencing technologies has boosted the interest in exploring the role of fusion genes in the development and progression of solid tumors. In breast cancer, most of the detected gene fusions seem to be "passenger" events while the presence of recurrent and driver fusions is still under study. We performed RNA sequencing in 55 well-characterized breast cancer samples and 10 adjacent normal breast tissues, complemented by an analysis of SNP array data. We explored the presence of fusion genes and defined their association with breast cancer subtypes, clinical-pathologic characteristics and copy number aberrations. Overall, 370 fusions were detected across the majority of the samples. HER2+ samples had significantly more fusions than triple negative and luminal subtypes. The number of fusions was correlated with histological grade, Ki67 and tumor size. Clusters of fusion genes were observed across the genome and a significant correlation of fusions with copy number aberrations and more specifically amplifications was also revealed. Despite the large number of fusion events, only a few were recurrent, while recurrent individual genes forming fusions with different partners were also detected including the estrogen receptor 1 gene in the previously detected ESR1-CCDC170 fusion. Overall we detected novel gene fusion events while we confirmed previously reported fusions. Genomic hotspots of fusion genes, differences between subtypes and small number of recurrent fusions are the most relevant characteristics of these events in breast cancer. Further investigation is necessary to comprehend the biological significance of these fusions.

Distinctive Desmoplastic 3D Morphology Associated With BRAFV600E in Papillary Thyroid Cancers

CONTEXT

Although 60% of papillary thyroid carcinomas are BRAFV600E mutant (PTCV600E), the increased aggressiveness of these cancers is still debated.

OBJECTIVE

For PTCV600E we aimed to further characterize the extent of the stroma and its activation, the three-dimensional (3D) tumor-stroma interface, and the proliferation rates of tumor and stromal fibroblasts.

DESIGN

We analyzed exomes, transcriptomes, and images of 364 papillary thyroid carcinoma (PTCs) from The Cancer Genome Atlas (TCGA), including 211 PTCV600E; stained 22 independent PTCs for BRAFV600E and Ki67; sequenced the exomes and stained BRAFV600E in 5 primary tumor blocks and 4 nodal metastases from one patient with PTCV600E; and reconstructed the 3D volumes of one tumor and one metastatic block at histological resolution.

RESULTS

In TCGA, BRAFV600E was associated with higher expression of proliferation markers and lower expression of thyroid differentiation markers, independently of tumor purity. Moreover, PTCV600E, in line with their overall lower purity, also had higher expression of fibroblast- and T cell-associated genes and presented more fibrosis. Tumor cells that appeared disconnected on two-dimensional histological slices were revealed to be part of a unique tumor component in the 3D reconstructed microvolumes, and they formed a surprisingly complex connected space, infiltrating a proliferative stroma. Finally, in our PTC set, both stromal fibroblasts and tumor cells presented higher proliferation rates in PTCV600E.

CONCLUSIONS

Our results support the increased aggressiveness associated with BRAFV600E in PTC and shed light on the important role of the stroma in tumor expansion. The greater and more active fibrotic component predicts better efficiency of combined targeted treatments, as previously proposed for melanomaV600E.

miRNA expression and function in thyroid carcinomas: a comparative and critical analysis and a model for other cancers

As in many cancer types, miRNA expression profiles and functions have become an important field of research on non-medullary thyroid carcinomas, the most common endocrine cancers. This could lead to the establishment of new diagnostic tests and new cancer therapies. However, different studies showed important variations in their research strategies and results. In addition, the action of miRNAs is poorly considered as a whole because of the use of underlying dogmatic truncated concepts. These lead to discrepancies and limits rarely considered. Recently, this field has been enlarged by new miRNA functional and expression studies. Moreover, studies using next generation sequencing give a new view of general miRNA differential expression profiles of papillary thyroid carcinoma. We analyzed in detail this literature from both physiological and differential expression points of view. Based on explicit examples, we reviewed the progresses but also the discrepancies and limits trying to provide a critical approach of where this literature may lead. We also provide recommendations for future studies. The conclusions of this systematic analysis could be extended to other cancer types.

Immune Checkpoint Molecules on Tumor-Infiltrating Lymphocytes and Their Association with Tertiary Lymphoid Structures in Human Breast Cancer

There is an exponentially growing interest in targeting immune checkpoint molecules in breast cancer (BC), particularly in the triple-negative subtype where unmet treatment needs remain. This study was designed to analyze the expression, localization, and prognostic role of PD-1, PD-L1, PD-L2, CTLA-4, LAG3, and TIM3 in primary BC. Gene expression analysis using the METABRIC microarray dataset found that all six immune checkpoint molecules are highly expressed in basal-like and HER2-enriched compared to the other BC molecular subtypes. Flow cytometric analysis of fresh tissue homogenates from untreated primary tumors show that PD-1 is principally expressed on CD4⁺ or CD8⁺ T cells and CTLA-4 is expressed on CD4⁺ T cells. The global proportion of PD-L1⁺, PD-L2⁺, LAG3⁺, and TIM3⁺ tumor-infiltrating lymphocytes (TIL) was low and detectable in only a small number of tumors. Immunohistochemically staining fixed tissues from the same tumors was employed to score TIL and tertiary lymphoid structures (TLS). PD-L1⁺, PD-L2⁺, LAG3⁺, and TIM3⁺ cells were detected in some TLS in a pattern that resembles secondary lymphoid organs. This observation suggests that TLS are important sites of immune activation and regulation, particularly in tumors with extensive baseline immune infiltration. Significantly improved overall survival was correlated with PD-1 expression in the HER2-enriched and PD-L1 or CTLA-4 expression in basal-like BC. PD-1 and CTLA-4 proteins were most frequently detected on TIL, which supports the correlations observed between their gene expression and improved long-term outcome in basal-like and HER2-enriched BC. PD-L1 expression by tumor or immune cells is uncommon in BC. Overall, the data presented here distinguish PD-1 as a marker of T cell activity in both the T and B cell areas of BC associated TLS. We found that immune checkpoint molecule expression parallels the extent of TIL and TLS, although there is a noteworthy amount of heterogeneity between tumors even within the same molecular subtype. These data indicate that assessing the levels of immune checkpoint molecule expression in an individual patient has important implications for the success of therapeutically targeting them in BC.

Abstract P1-05-02: Genomic, transcriptomic and immune features of breast cancer according to the patient's body mass index at diagnosis

This abstract was withdrawn by the authors.

A research note regarding "Variation in cancer risk among tissues can be explained by the number of stem cell divisions"

Tomasetti and Vogelstein argued that 2/3 of human cancers are due to 'bad luck' and that "primary prevention measures [against cancer] are not likely to be very effective". We demonstrate that their calculations for hepatocellular carcinomas overlooked a major subset of these cancers proven to be preventable through vaccination. The problem, which is not limited to hepatocellular carcinoma, arises from the general reliance of their analysis on average incidences in the United States and the omission of incidences in specific risk groups.

Gene signature of the post-Chernobyl papillary thyroid cancer

Purpose

Following the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC.

Methods

We investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized.

Results

There were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test, p < 0.01), with 300 differently expressed probe sets (p < 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the BRAF mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (PPME1, HDAC11, SOCS7, CIC, THRA, ERBB2, PPP1R9A, HDGF, RAD51AP1, and CDK1) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples.

Conclusion

Significant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure.

Electronic supplementary material

The online version of this article (doi:10.1007/s00259-015-3303-3) contains supplementary material, which is available to authorized users.

New global analysis of the microRNA transcriptome of primary tumors and lymph node metastases of papillary thyroid cancer

Background

Papillary Thyroid Cancer (PTC) is the most prevalent type of endocrine cancer. Its incidence has rapidly increased in recent decades but little is known regarding its complete microRNA transcriptome (miRNome). In addition, there is a need for molecular biomarkers allowing improved PTC diagnosis.

Methods

We performed small RNA deep-sequencing of 3 PTC, their matching normal tissues and lymph node metastases (LNM). We designed a new bioinformatics framework to handle each aspect of the miRNome: whole expression profiles, isomiRs distribution, non-templated additions distributions, RNA-editing or mutation. Results were validated experimentally by qRT-PCR on normal samples, tumors and LNM from 14 independent patients and in silico using the dataset from The Cancer Genome Atlas (small RNA deepsequencing of 59 normal samples, 495 PTC, and 8 LNM).

Results

We performed small RNA deep-sequencing of 3 PTC, their matching normal tissues and lymph node metastases (LNM). We designed a new bioinformatics framework to handle each aspect of the miRNome: whole expression profiles, isomiRs distribution, non-templated additions distributions, RNA-editing or mutation. Results were validated experimentally by qRT-PCR on normal samples, tumors and LNM from 14 independent patients and in silico using the dataset from The Cancer Genome Atlas (small RNA deep-sequencing of 59 normal samples, 495 PTC, and 8 LNM). We confirmed already described up-regulations of microRNAs in PTC, such as miR-146b-5p or miR-222-3p, but we also identified down-regulated microRNAs, such as miR-7-5p or miR-30c-2-3p. We showed that these down-regulations are linked to the tumorigenesis process of thyrocytes. We selected the 14 most down-regulated microRNAs in PTC and we showed that they are potential biomarkers of PTC samples. Nevertheless, they can distinguish histological classical variants and follicular variants of PTC in the TCGA dataset. In addition, 12 of the 14 down-regulated microRNAs are significantly less expressed in aggressive PTC compared to non-aggressive PTC. We showed that the associated aggressive expression profile is mainly due to the presence of the BRAF V600E mutation. In general, primary tumors and LNM presented similar microRNA expression profiles but specific variations like the down-regulation of miR-7-2-3p and miR-30c-2-3p in LNM were observed. Investigations of the 5p-to-3p arm expression ratios, non-templated additions or isomiRs distributions revealed no major implication in PTC tumorigenesis process or LNM appearance.

Conclusions

Our results showed that down-regulated microRNAs can be used as new potential common biomarkers of PTC and to distinguish main subtypes of PTC. MicroRNA expressions can be linked to the development of LNM of PTC. The bioinformatics framework that we have developed can be used as a starting point for the global analysis of any microRNA deep-sequencing data in an unbiased way.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2082-3) contains supplementary material, which is available to authorized users.

Abstract S4-02: Principles governing A-to-I RNA editing in breast cancer transcriptome

Background: Messenger RNA (mRNA) is the target of a series of post-transcriptional modifications that can affect its structure and stability, one of the most relevant being RNA editing. The most common form of RNA editing in humans is of the A-to-I type and is catalyzed by the adenosine deaminase acting on RNA (ADAR) family of enzymes. Currently, little is known about how RNA editing operates in cancer. The main objectives of this study were to investigate and characterize the extent of A-to-I RNA editing in breast cancer (BC) and to define the principles governing the editing process in this as well as other cancers.

Material and Methods: The exome and transcriptome of 58 BC samples representing the four main known subtypes, namely TN, HER2+, luminal A and luminal B, and 10 matched normal samples were sequenced using the Illumina platform. For the same series, gene expression and copy number profiles were obtained using the Affymetrix platform. RNA-DNA single nucleotide differences (RDDs) were called using a pipeline in line with the most updated bioinformatics tools and validated in an independent cohort of 15 BC samples and breast cell lines.

Results: Overall, we detected 16,027 RDDs present in one or more samples, with all possible base changes represented. Among these, 560 RDDs were located in Alu regions and were all of the A-to-I type consistent with the notion that A-to-I editing occurs predominantly in forward-facing Alu forming double stranded RNA duplexes processed by ADAR.

We found that the same sites were edited in normal and tumor breast tissues. However, the editing frequency was significantly higher in tumors compared to matched normal breast tissues. Moreover, high editing frequencies were observed in samples in which more editing sites were detectable and/or in which ADAR expression was high.

We identified two key factors that independently determine ADAR expression and therefore A-to-I RNA editing in breast and the majority of other human cancers: 1) the type-I interferon response in tumors and 2) gains in ADAR copy number. The mean editing frequency was found to be significantly correlated with the expression of STAT1 and other type I interferon target genes, both in our patient series, in a large pool of BC datasets and a panel of normal and breast cancer cell lines treated with interferon (IFN) α, β, or γ for 1, 2 and 5 days in vitro. Moreover, the association between editing and STAT1 expression or ADAR amplification was validated in 19 additional cancer types obtained from the TCGA dataset.

Conclusions: Our work, which represents the largest survey so far on RNA editing in BC, shows that A-to-I editing is a pervasive and well-controlled phenomenon in cancer that can drive aberrant transcriptome expression in breast and potentially the vast majority of cancers. Moreover, it suggests that the immune response can profoundly impact the transcriptome sequence in tumor cells and thereby influence the internal mechanisms governing their behavior.

Citation Format: Debora Fumagalli, David Gacquer, Françoise Rothé, Anne Lefort, Frederick Libert, David N Brown, Naima Kheddoumi, Adam Shlien, Tomasz Konopka, Roberto Salgado, Denis Larsimont, Kornelia Polyak, Karen Willard-Gallo, Christine Desmedt, Martine Piccart, Marc Abramowicz, Peter J Campbell, Vincent Detours, Christos Sotiriou. Principles governing A-to-I RNA editing in breast cancer transcriptome [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr S4-02.

Intratumor heterogeneity and clonal evolution in an aggressive papillary thyroid cancer and matched metastases

The contribution of intratumor heterogeneity to thyroid metastatic cancers is still unknown. The clonal relationships between the primary thyroid tumors and lymph nodes (LN) or distant metastases are also poorly understood. The objective of this study was to determine the phylogenetic relationships between matched primary thyroid tumors and metastases. We searched for non-synonymous single-nucleotide variants (nsSNVs), gene fusions, alternative transcripts, and loss of heterozygosity (LOH) by paired-end massively parallel sequencing of cDNA (RNA-Seq) in a patient diagnosed with an aggressive papillary thyroid cancer (PTC). Seven tumor samples from a stage IVc PTC patient were analyzed by RNA-Seq: two areas from the primary tumor, four areas from two LN metastases, and one area from a pleural metastasis (PLM). A large panel of other thyroid tumors was used for Sanger sequencing screening. We identified seven new nsSNVs. Some of these were early events clonally present in both the primary PTC and the three matched metastases. Other nsSNVs were private to the primary tumor, the LN metastases and/or the PLM. Three new gene fusions were identified. A novel cancer-specific KAZN alternative transcript was detected in this aggressive PTC and in dozens of additional thyroid tumors. The PLM harbored an exclusive whole-chromosome 19 LOH. We have presented the first, to our knowledge, deep sequencing study comparing the mutational spectra in a PTC and both LN and distant metastases. This study has yielded novel findings concerning intra-tumor heterogeneity, clonal evolution and metastases dissemination in thyroid cancer.

No significant viral transcription detected in whole breast cancer transcriptomes

Background

Studies evaluating the presence of viral sequences in breast cancer (BC), including various strains of human papillomavirus and human herpes virus, have yielded conflicting results. Most were based on RT-PCR and in situ hybridization.

Methods

In this report we searched for expressed viral sequences in 58 BC transcriptomes using five distinct in silico methods. In addition, we complemented our RNA sequencing results with exome sequencing, PCR and immunohistochemistry (IHC) analyses. A control sample was used to test our in silico methods.

Results

All of the computational methods correctly detected viral sequences in the control sample. We identified a small number of viral sequences belonging to human herpesvirus 4 and 6 and Merkel cell polyomavirus. The extremely low expression levels—two orders of magnitude lower than in a typical hepatitis B virus infection in hepatocellular carcinoma—did not suggest active infections. The presence of viral elements was confirmed in sample-matched exome sequences, but could not be confirmed by PCR or IHC.

Conclusions

Our results show that no viral sequences are expressed in significant amounts in the BC investigated. The presence of non-transcribed viral DNA cannot be excluded.

Transfer of clinically relevant gene expression signatures in breast cancer: from Affymetrix microarray to Illumina RNA-Sequencing technology

Background

Microarrays have revolutionized breast cancer (BC) research by enabling studies of gene expression on a transcriptome-wide scale. Recently, RNA-Sequencing (RNA-Seq) has emerged as an alternative for precise readouts of the transcriptome. To date, no study has compared the ability of the two technologies to quantify clinically relevant individual genes and microarray-derived gene expression signatures (GES) in a set of BC samples encompassing the known molecular BC’s subtypes. To accomplish this, the RNA from 57 BCs representing the four main molecular subtypes (triple negative, HER2 positive, luminal A, luminal B), was profiled with Affymetrix HG-U133 Plus 2.0 chips and sequenced using the Illumina HiSeq 2000 platform. The correlations of three clinically relevant BC genes, six molecular subtype classifiers, and a selection of 21 GES were evaluated.

Results

16,097 genes common to the two platforms were retained for downstream analysis. Gene-wise comparison of microarray and RNA-Seq data revealed that 52% had a Spearman’s correlation coefficient greater than 0.7 with highly correlated genes displaying significantly higher expression levels. We found excellent correlation between microarray and RNA-Seq for the estrogen receptor (ER; r_s = 0.973; 95% CI: 0.971-0.975), progesterone receptor (PgR; r_s = 0.95; 0.947-0.954), and human epidermal growth factor receptor 2 (HER2; r_s = 0.918; 0.912-0.923), while a few discordances between ER and PgR quantified by immunohistochemistry and RNA-Seq/microarray were observed. All the subtype classifiers evaluated agreed well (Cohen’s kappa coefficients >0.8) and all the proliferation-based GES showed excellent Spearman correlations between microarray and RNA-Seq (all r_s >0.965). Immune-, stroma- and pathway-based GES showed a lower correlation relative to prognostic signatures (all r_s >0.6).

Conclusions

To our knowledge, this is the first study to report a systematic comparison of RNA-Seq to microarray for the evaluation of single genes and GES clinically relevant to BC. According to our results, the vast majority of single gene biomarkers and well-established GES can be reliably evaluated using the RNA-Seq technology.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1008) contains supplementary material, which is available to authorized users.

Systems biology of cancer: entropy, disorder, and selection-driven evolution to independence, invasion and "swarm intelligence"

Our knowledge of the biology of solid cancer has greatly progressed during the last few years, and many excellent reviews dealing with the various aspects of this biology have appeared. In the present review, we attempt to bring together these subjects in a general systems biology narrative. It starts from the roles of what we term entropy of signaling and noise in the initial oncogenic events, to the first major transition of tumorigenesis: the independence of the tumor cell and the switch in its physiology, i.e., from subservience to the organism to its own independent Darwinian evolution. The development after independence involves a constant dynamic reprogramming of the cells and the emergence of a sort of collective intelligence leading to invasion and metastasis and seldom to the ultimate acquisition of immortality through inter-individual infection. At each step, the probability of success is minimal to infinitesimal, but the number of cells possibly involved and the time scale account for the relatively high occurrence of tumorigenesis and metastasis in multicellular organisms.

Comparative analysis of the thyrocytes and T cells: responses to H2O2 and radiation reveals an H2O2-induced antioxidant transcriptional program in thyrocytes

CONTEXT

Radiation is an established cause of thyroid cancer, and growing evidence supports a role for hydrogen peroxide (H2O2) in spontaneous thyroid carcinogenesis. Little is known about the molecular programs activated by these agents in thyrocytes.

OBJECTIVE

The purpose of this study was to compare the responses of thyrocytes and T cells to H2O2 and radiation.

METHODS

We profiled the DNA damage and cell death induced by γ-radiation (0.1-5 Gy) and H2O2 (0.0025-0.3 mM) in primary human thyrocytes and T cells. We next prepared thyroid and T-cell primary cultures from 8 donors operated for noncancerous thyroid pathological conditions and profiled their genome-wide transcriptional response 4 hours after (1) exposure to 1-Gy radiation, (2) treatment with H2O2 and (3) no treatment. Two H2O2 concentrations were investigated, calibrated in each cell type to elicit levels of single- and double-strand breaks equivalent to 1-Gy γ-radiation.

RESULTS

Although thyrocytes and T cells had comparable radiation responses, 3- to 10-fold more H2O2 was needed to induce detectable DNA damage in thyrocytes. At H2O2 and radiation doses inducing double-strand breaks, cell death occurred after 24 hours in T cells but not in thyrocytes. The transcriptional responses of thyrocytes and T cells to radiation were similar, involving DNA repair and cell death genes. In addition to this transcriptional program, H2O2 also up-regulated antioxidant genes in thyrocytes, including glutathione peroxidases and heme oxygenase at the double-strand breaks-inducing concentration. In contrast, a transcriptional storm involving thousands of genes was raised in T cells. Finally, we showed that inhibiting glutathione peroxidases activity increased the DNA damaging effect of H2O2 in thyrocytes.

CONCLUSION

We propose that high H2O2 production in thyrocytes is matched with specific transcriptionally regulated antioxidant protection.

TriageTools: tools for partitioning and prioritizing analysis of high-throughput sequencing data

High-throughput sequencing is becoming a popular research tool but carries with it considerable costs in terms of computation time, data storage and bandwidth. Meanwhile, some research applications focusing on individual genes or pathways do not necessitate processing of a full sequencing dataset. Thus, it is desirable to partition a large dataset into smaller, manageable, but relevant pieces. We present a toolkit for partitioning raw sequencing data that includes a method for extracting reads that are likely to map onto pre-defined regions of interest. We show the method can be used to extract information about genes of interest from DNA or RNA sequencing samples in a fraction of the time and disk space required to process and store a full dataset. We report speedup factors between 2.6 and 96, depending on settings and samples used. The software is available at http://www.sourceforge.net/projects/triagetools/.

A measure of the signal-to-noise ratio of microarray samples and studies using gene correlations

Background

The quality of gene expression data can vary dramatically from platform to platform, study to study, and sample to sample. As reliable statistical analysis rests on reliable data, determining such quality is of the utmost importance. Quality measures to spot problematic samples exist, but they are platform-specific, and cannot be used to compare studies.

Results

As a proxy for quality, we propose a signal-to-noise ratio for microarray data, the “Signal-to-Noise Applied to Gene Expression Experiments”, or SNAGEE. SNAGEE is based on the consistency of gene-gene correlations. We applied SNAGEE to a compendium of 80 large datasets on 37 platforms, for a total of 24,380 samples, and assessed the signal-to-noise ratio of studies and samples. This allowed us to discover serious issues with three studies. We show that signal-to-noise ratios of both studies and samples are linked to the statistical significance of the biological results.

Conclusions

We showed that SNAGEE is an effective way to measure data quality for most types of gene expression studies, and that it often outperforms existing techniques. Furthermore, SNAGEE is platform-independent and does not require raw data files. The SNAGEE R package is available in BioConductor.

InSilico DB genomic datasets hub: an efficient starting point for analyzing genome-wide studies in GenePattern, Integrative Genomics Viewer, and R/Bioconductor

Genomics datasets are increasingly useful for gaining biomedical insights, with adoption in the clinic underway. However, multiple hurdles related to data management stand in the way of their efficient large-scale utilization. The solution proposed is a web-based data storage hub. Having clear focus, flexibility and adaptability, InSilico DB seamlessly connects genomics dataset repositories to state-of-the-art and free GUI and command-line data analysis tools. The InSilico DB platform is a powerful collaborative environment, with advanced capabilities for biocuration, dataset sharing, and dataset subsetting and combination. InSilico DB is available from https://insilicodb.org.

Piecewise polynomial representations of genomic tracks

Genomic data from micro-array and sequencing projects consist of associations of measured values to chromosomal coordinates. These associations can be thought of as functions in one dimension and can thus be stored, analyzed, and interpreted as piecewise-polynomial curves. We present a general framework for building piecewise polynomial representations of genome-scale signals and illustrate some of its applications via examples. We show that piecewise constant segmentation, a typical step in copy-number analyses, can be carried out within this framework for both array and (DNA) sequencing data offering advantages over existing methods in each case. Higher-order polynomial curves can be used, for example, to detect trends and/or discontinuities in transcription levels from RNA-seq data. We give a concrete application of piecewise linear functions to diagnose and quantify alignment quality at exon borders (splice sites). Our software (source and object code) for building piecewise polynomial models is available at http://sourceforge.net/projects/locsmoc/.

Most random gene expression signatures are significantly associated with breast cancer outcome

Bridging the gap between animal or in vitro models and human disease is essential in medical research. Researchers often suggest that a biological mechanism is relevant to human cancer from the statistical association of a gene expression marker (a signature) of this mechanism, that was discovered in an experimental system, with disease outcome in humans. We examined this argument for breast cancer. Surprisingly, we found that gene expression signatures—unrelated to cancer—of the effect of postprandial laughter, of mice social defeat and of skin fibroblast localization were all significantly associated with breast cancer outcome. We next compared 47 published breast cancer outcome signatures to signatures made of random genes. Twenty-eight of them (60%) were not significantly better outcome predictors than random signatures of identical size and 11 (23%) were worst predictors than the median random signature. More than 90% of random signatures >100 genes were significant outcome predictors. We next derived a metagene, called meta-PCNA, by selecting the 1% genes most positively correlated with proliferation marker PCNA in a compendium of normal tissues expression. Adjusting breast cancer expression data for meta-PCNA abrogated almost entirely the outcome association of published and random signatures. We also found that, in the absence of adjustment, the hazard ratio of outcome association of a signature strongly correlated with meta-PCNA (R² = 0.9). This relation also applied to single-gene expression markers. Moreover, >50% of the breast cancer transcriptome was correlated with meta-PCNA. A corollary was that purging cell cycle genes out of a signature failed to rule out the confounding effect of proliferation. Hence, it is questionable to suggest that a mechanism is relevant to human breast cancer from the finding that a gene expression marker for this mechanism predicts human breast cancer outcome, because most markers do. The methods we present help to overcome this problem.

Proving that research findings from in vitro or animal models are relevant to human diseases is a major bottleneck in medical science. Hundreds of researchers have suggested the human relevance of oncogenic mechanisms from the statistical association between gene expression markers of these mechanisms and disease outcome. Such evidence has become easier to obtain recently with the advent of microarray screens and of large public-domain genome-wide expression datasets with patient follow-up. We demonstrated that in breast cancer any set of 100 genes or more selected at random has a 90% chance to be significantly associated with outcome. Thus, investigators are bound to find an association however whimsical their marker is. For example, we could establish outcome associations for a signature of postprandial laughter and a signature of social defeat in mice. Association was not stronger than expected at random for 28 (60%) of 47 published breast cancer signatures. The odds of association are 5–17% with random single gene markers—a finding relevant to older breast cancer studies. We explained these results by showing that much of the breast cancer transcriptome is correlated with proliferation, which integrates most prognostic information in this disease.

Genes expressed in specific areas of the human fetal cerebral cortex display distinct patterns of evolution

The developmental mechanisms through which the cerebral cortex increased in size and complexity during primate evolution are essentially unknown. To uncover genetic networks active in the developing cerebral cortex, we combined three-dimensional reconstruction of human fetal brains at midgestation and whole genome expression profiling. This novel approach enabled transcriptional characterization of neurons from accurately defined cortical regions containing presumptive Broca and Wernicke language areas, as well as surrounding associative areas. We identified hundreds of genes displaying differential expression between the two regions, but no significant difference in gene expression between left and right hemispheres. Validation by qRTPCR and in situ hybridization confirmed the robustness of our approach and revealed novel patterns of area- and layer-specific expression throughout the developing cortex. Genes differentially expressed between cortical areas were significantly associated with fast-evolving non-coding sequences harboring human-specific substitutions that could lead to divergence in their repertoires of transcription factor binding sites. Strikingly, while some of these sequences were accelerated in the human lineage only, many others were accelerated in chimpanzee and/or mouse lineages, indicating that genes important for cortical development may be particularly prone to changes in transcriptional regulation across mammals. Genes differentially expressed between cortical regions were also enriched for transcriptional targets of FoxP2, a key gene for the acquisition of language abilities in humans. Our findings point to a subset of genes with a unique combination of cortical areal expression and evolutionary patterns, suggesting that they play important roles in the transcriptional network underlying human-specific neural traits.

Human cancer cell lines: Experimental models for cancer cells in situ? For cancer stem cells?

Established human cancer cell lines are routinely used as experimental models for human cancers. Their validity for such use is analyzed and discussed, with particular focus on thyroid tumors. Although cell lines retain some properties of the cells of origin, from the points of view of their genetics, epigenetics and gene expression, they show clear differences in these properties compared to in vivo tumors. This can be explained by a prior selection of initiating cells and a Darwinian evolution in vitro. The properties of the cell lines are compared to those of the postulated cancer stem cells and their use as models in this regard are discussed. Furthermore, other proper and possible uses of the cell lines are discussed.

Gene expression in RET/PTC3 and E7 transgenic mouse thyroids: RET/PTC3 but not E7 tumors are partial and transient models of human papillary thyroid cancers

We studied gene expression profiles in two mouse models of human thyroid carcinoma: the Tg-RET/PTC3 (RP3) and Tg-E7 mice. RP3 fusion gene is the most frequent mutation found in the first wave post-Chernobyl papillary thyroid cancers (PTCs). E7 is an oncoprotein derived from the human papillomavirus 16 responsible for most cervical carcinoma in women. Both transgenic mice develop thyroid hyperplasia followed by solid differentiated carcinoma in older animals. To understand the different steps leading to carcinoma, we analyzed thyroid gene expression in both strains at different ages by microarray technology. Important biological processes were differentially regulated in the two tumor types. In E7 thyroids, cell cycle was the most up-regulated process, an observation consistent with the huge size of these tumors. In RP3 thyroids, contrary to E7 tumors, several human PTC characteristics were observed: overexpression of many immune-related genes, regulation of human PTC markers, up-regulation of EGF-like growth factors and significant regulation of angiogenesis and extracellular matrix remodeling-related genes. However, similarities were incomplete; they did not concern the overall gene expression and were not conserved in old animals. Therefore, RP3 tumors are partial and transient models of human PTC. They constitute a good model, especially in young animals, to study the respective role of the biological processes shared with human PTC and will allow testing drugs targeting these validated variables.

Gene expression profiles of post-Chernobyl thyroid cancers

PURPOSE OF REVIEW

We discuss new evidence supporting the existence of a susceptibility to develop cancer following radiation exposure that is variable in the general population and could be measurable from gene expression.

RECENT FINDINGS

Microarray analysis of spontaneous and post-Chernobyl thyroid cancers has uncovered gene expression radiation signatures, one of which could be related to the putative cause of these tumors and to a DNA repair pathway. A gene expression signature distinguishes the lymphocytes drawn from parents of children with retinoblastoma and the lymphocytes of parents of healthy children. The first are more radiosensitive. A familial clustering pattern is observed in radiation-induced meningiomas.

SUMMARY

The existence of a susceptibility to develop radiation-induced cancer would explain why only a minority of the population most heavily exposed to radiation following the Chernobyl disaster developed a cancer. The possibility of measuring this susceptibility from gene expression has a number of implications for research, medicine and radioprotection.

[Transcription markers of survival in breast neoplasms measuring the signals from the proliferation]

We show that proliferation-related signals are omnipresent in the breast cancer transcriptome. As a result, many transcriptional signatures generated at random are valuable for the prognosis of disease-free survival: despite their biological rationale, 30-60% of published prognostic signatures are not significantly better. We propose a mathematical transformation, the super PCNA decovolution, which removes proliferation-related signals from tumours transcriptional profiles. Both random and published signatures loose nearly all their prognostic value after removal of these signals.

Human thyroid tumor cell lines derived from different tumor types present a common dedifferentiated phenotype

Cell lines are crucial to elucidate mechanisms of tumorigenesis and serve as tools for cancer treatment screenings. Therefore, careful validation of whether these models have conserved properties of in vivo tumors is highly important. Thyrocyte-derived tumors are very interesting for cancer biology studies because from one cell type, at least five histologically characterized different benign and malignant tumor types can arise. To investigate whether thyroid tumor-derived cell lines are representative in vitro models, characteristics of eight of those cell lines were investigated with microarrays, differentiation markers, and karyotyping. Our results indicate that these cell lines derived from differentiated and undifferentiated tumor types have evolved in vitro into similar phenotypes with gene expression profiles the closest to in vivo undifferentiated tumors. Accordingly, the absence of expression of most thyrocyte-specific genes, the nonresponsiveness to thyrotropin, as well as their large number of chromosomal abnormalities, suggest that these cell lines have acquired characteristics of fully dedifferentiated cells. They represent the outcome of an adaptation and evolution in vitro, which questions the reliability of these cell lines as models for differentiated tumors. However, they may represent useful models for undifferentiated cancers, and by their comparison with differentiated cells, can help to define the genes involved in the differentiation/dedifferentiation process. The use of any cell line as a model for a cancer therefore requires prior careful and thorough validation for the investigated property.

Roles of hydrogen peroxide in thyroid physiology and disease

CONTEXT

The long-lived thyroid cell generates, for the synthesis of thyroid hormones, important amounts of H2O2 that are toxic in other cell types. This review analyzes the protection mechanisms of the cell and the pathological consequences of disorders of this system.

EVIDENCE ACQUISITION

The literature on H2O2 generation and disposal, thyroid hormone synthesis, and their control in the human thyroid is analyzed.

EVIDENCE SYNTHESIS

In humans, H2O2 production by dual-oxidases and consequently thyroid hormone synthesis by thyroperoxidase are controlled by the phospholipase C-Ca2+-diacylglycerol arm of TSH receptor action. H2O2 in various cell types, and presumably in thyroid cells, is a signal, a mitogen, a mutagen, a carcinogen, and a killer. The various protection mechanisms of the thyroid cell against H2O2 are analyzed. They include the separation of the generating enzymes (dual-oxidases), their coupling to thyroperoxidase in a proposed complex, the thyroxisome, and H2O2 degradation systems.

CONCLUSIONS

It is proposed that various pathologies can be explained, at least in part, by overproduction and lack of degradation of H2O2 (tumorigenesis, myxedematous cretinism, and thyroiditis) and by failure of the H2O2 generation or its positive control system (congenital hypothyroidism).

Genome-wide gene expression profiling suggests distinct radiation susceptibilities in sporadic and post-Chernobyl papillary thyroid cancers

Papillary thyroid cancers (PTCs) incidence dramatically increased in the vicinity of Chernobyl. The cancer-initiating role of radiation elsewhere is debated. Therefore, we searched for a signature distinguishing radio-induced from sporadic cancers. Using microarrays, we compared the expression profiles of PTCs from the Chernobyl Tissue Bank (CTB, n=12) and from French patients with no history of exposure to ionising radiations (n=14). We also compared the transcriptional responses of human lymphocytes to the presumed aetiological agents initiating these tumours, γ-radiation and H₂O₂. On a global scale, the transcriptomes of CTB and French tumours are indistinguishable, and the transcriptional responses to γ-radiation and H₂O₂ are similar. On a finer scale, a 118 genes signature discriminated the γ-radiation and H₂O₂ responses. This signature could be used to classify the tumours as CTB or French with an error of 15–27%. Similar results were obtained with an independent signature of 13 genes involved in homologous recombination. Although sporadic and radio-induced PTCs represent the same disease, they are distinguishable with molecular signatures reflecting specific responses to γ-radiation and H₂O₂. These signatures in PTCs could reflect the susceptibility profiles of the patients, suggesting the feasibility of a radiation susceptibility test.

Gene expression and the biological phenotype of papillary thyroid carcinomas

The purpose of this paper is to correlate the molecular phenotype of papillary thyroid carcinoma (PTC) to their biological pathology. We hybridized 26 PTC on microarrays and showed that nearly 44% of the transcriptome was regulated in these tumors. We then combined our data set with two published PTC microarray studies to produce a platform- and study-independent list of PTC-associated genes. We further confirmed the mRNA regulation of 15 genes from this list by quantitative reverse transcription-PCR. Analysis of this list with statistical tools led to several conclusions: (1) there is a change in cell population with an increased expression of genes involved in the immune response, reflecting lymphocyte infiltration in the tumor compared to the normal tissue. (2) The c-jun N-terminal kinase pathway is activated by overexpression of its components. (3) The activation of ERKK1/2 by genetic alterations is supplemented by activation of the epidermal growth factor but not of the insulin-like growth factor signaling pathway. (4) There is a downregulation of immediate early genes. (5) We observed an overexpression of many proteases in accordance with tumor remodeling, and suggested a probable role of S100 proteins and annexin A2 in this process. (6) Numerous overexpressed genes favor the hypothesis of a collective migration mode of tumor cells.

Gene expression in thyroid autonomous adenomas provides insight into their physiopathology

The purpose of this study was to use the microarray technology to define expression profiles characteristic of thyroid autonomous adenomas and relate these findings to physiological mechanisms. Experiments were performed on a series of separated adenomas and their normal counterparts on Micromax cDNA microarrays covering 2400 genes (analysis I), and on a pool of adenomatous tissues and their corresponding normal counterparts using microarrays of 18,000 spots (analysis II). Results for genes present on the two arrays corroborated and several gene regulations previously determined by Northern blotting or microarrays in similar lesions were confirmed. Five overexpressed and 24 underexpressed genes were also confirmed by real-time RT-PCR in some of the samples used for microarray analysis, and in additional tumor specimens. Our results show: (1) a change in the cell populations of the tumor, with a marked decrease in lymphocytes and blood cells and an increase in endothelial cells. The latter increase would correspond to the establishment of a close relation between thyrocytes and endothelial cells and is related to increased N-cadherin expression. It explains the increased blood flow in the tumor; (2) a homogeneity of tumor samples correlating with their common physiopathological mechanism: the constitutive activation of the thyrotropin (TSH)/cAMP cascade; (3) a low proportion of regulated genes consistent with the concept of a minimal deviation tumor; (4) a higher expression of genes coding for specific functional proteins, consistent with the functional hyperactivity of the tumors; (5) an increase of phosphodiesterase gene expression which explains the relatively low cyclic AMP levels measured in these tumors; (6) an overexpression of antiapoptotic genes and underexpression of proapoptotic genes compatible with their low apoptosis rate; (7) an overexpression of N-cadherin and downregulation of caveolins, which casts doubt about the use of these expressions as markers for malignancy.

Gene expression in human thyrocytes and autonomous adenomas reveals suppression of negative feedbacks in tumorigenesis

The cAMP signaling pathway regulates growth of many cell types, including somatotrophs, thyrocytes, melanocytes, ovarian follicular granulosa cells, adrenocortical cells, and keratinocytes. Mutations of partners from the cAMP signaling cascade are involved in tumor formation. Thyroid-stimulating hormone (TSH) receptor and Gsalpha activating mutations have been detected in thyroid autonomous adenomas, Gsalpha mutations in growth hormone-secreting pituitary adenomas, and PKAR1A mutations in Carney complex, a multiple neoplasia syndrome. To gain more insight into the role of cAMP signaling in tumor formation, human primary cultures of thyrocytes were treated for different times (1.5, 3, 16, 24, and 48 h) with TSH to characterize modulations in gene expression using cDNA microarrays. This kinetic study showed a clear difference in expression, early (1.5 and 3 h) and late (16-48 h) after the onset of TSH stimulation. This result suggests a progressive sequential process leading to a change of cell program. The gene expression profile of the long-term stimulated cultures resembled the autonomous adenomas, but not papillary carcinomas. The molecular phenotype of the adenomas thus confirms the role of long-term stimulation of the TSH-cAMP cascade in the pathology. TSH induced a striking up-regulation of different negative feedback modulators of the cAMP cascade, presumably insuring the one-shot effect of the stimulus. Some were down- or nonregulated in adenomas, suggesting a loss of negative feedback control in the tumors. These results suggest that in tumorigenesis, activation of proliferation pathways may be complemented by suppression of multiple corresponding negative feedbacks, i.e., specific tumor suppressors.

Absence of a specific radiation signature in post-Chernobyl thyroid cancers

Thyroid cancers have been the main medical consequence of the Chernobyl accident. On the basis of their pathological features and of the fact that a large proportion of them demonstrate RET-PTC translocations, these cancers are considered as similar to classical sporadic papillary carcinomas, although molecular alterations differ between both tumours. We analysed gene expression in post-Chernobyl cancers, sporadic papillary carcinomas and compared to autonomous adenomas used as controls. Unsupervised clustering of these data did not distinguish between the cancers, but separates both cancers from adenomas. No gene signature separating sporadic from post-Chernobyl PTC (chPTC) could be found using supervised and unsupervised classification methods although such a signature is demonstrated for cancers and adenomas. Furthermore, we demonstrate that pooled RNA from sporadic and chPTC are as strongly correlated as two independent sporadic PTC pools, one from Europe, one from the US involving patients not exposed to Chernobyl radiations. This result relies on cDNA and Affymetrix microarrays. Thus, platform-specific artifacts are controlled for. Our findings suggest the absence of a radiation fingerprint in the chPTC and support the concept that post-Chernobyl cancer data, for which the cancer-causing event and its date are known, are a unique source of information to study naturally occurring papillary carcinomas.

Integration and cross-validation of high-throughput gene expression data: comparing heterogeneous data sets

Data analysis--not data production--is becoming the bottleneck in gene expression research. Data integration is necessary to cope with an ever increasing amount of data, to cross-validate noisy data sets, and to gain broad interdisciplinary views of large biological data sets. New Internet resources may help researchers to combine data sets across different gene expression platforms. However, noise and disparities in experimental protocols strongly limit data integration. A detailed review of four selected studies reveals how some of these limitations may be circumvented and illustrates what can be achieved through data integration.

Clustering patterns of cytotoxic T-lymphocyte epitopes in human immunodeficiency virus type 1 (HIV-1) proteins reveal imprints of immune evasion on HIV-1 global variation

The human cytotoxic T-lymphocyte (CTL) response to human immunodeficiency virus type 1 (HIV-1) has been intensely studied, and hundreds of CTL epitopes have been experimentally defined, published, and compiled in the HIV Molecular Immunology Database. Maps of CTL epitopes on HIV-1 protein sequences reveal that defined epitopes tend to cluster. Here we integrate the global sequence and immunology databases to systematically explore the relationship between HIV-1 amino acid sequences and CTL epitope distributions. CTL responses to five HIV-1 proteins, Gag p17, Gag p24, reverse transcriptase (RT), Env, and Nef, have been particularly well characterized in the literature to date. Through comparing CTL epitope distributions in these five proteins to global protein sequence alignments, we identified distinct characteristics of HIV amino acid sequences that correlate with CTL epitope localization. First, experimentally defined HIV CTL epitopes are concentrated in relatively conserved regions. Second, the highly variable regions that lack epitopes bear cumulative evidence of past immune escape that may make them relatively refractive to CTLs: a paucity of predicted proteasome processing sites and an enrichment for amino acids that do not serve as C-terminal anchor residues. Finally, CTL epitopes are more highly concentrated in alpha-helical regions of proteins. Based on amino acid sequence characteristics, in a blinded fashion, we predicted regions in HIV regulatory and accessory proteins that would be likely to contain CTL epitopes; these predictions were then validated by comparison to new sets of experimentally defined epitopes in HIV-1 Rev, Tat, Vif, and Vpr.

Prediction of proteasome cleavage motifs by neural networks

We present a predictive method that can simulate an essential step in the antigen presentation in higher vertebrates, namely the step involving the proteasomal degradation of polypeptides into fragments which have the potential to bind to MHC Class I molecules. Proteasomal cleavage prediction algorithms published so far were trained on data from in vitro digestion experiments with constitutive proteasomes. As a result, they did not take into account the characteristics of the structurally modified proteasomes--often called immunoproteasomes--found in cells stimulated by gamma-interferon under physiological conditions. Our algorithm has been trained not only on in vitro data, but also on MHC Class I ligand data, which reflect a combination of immunoproteasome and constitutive proteasome specificity. This feature, together with the use of neural networks, a non-linear classification technique, make the prediction of MHC Class I ligand boundaries more accurate: 65% of the cleavage sites and 85% of the non-cleavage sites are correctly determined. Moreover, we show that the neural networks trained on the constitutive proteasome data learns a specificity that differs from that of the networks trained on MHC Class I ligands, i.e. the specificity of the immunoproteasome is different than the constitutive proteasome. The tools developed in this study in combination with a predictor of MHC and TAP binding capacity should give a more complete prediction of the generation and presentation of peptides on MHC Class I molecules. Here we demonstrate that such an approach produces an accurate prediction of the CTL the epitopes in HIV Nef. The method is available at www.cbs.dtu.dk/services/NetChop/.

Evolutionary and immunological implications of contemporary HIV-1 variation

Evolutionary modelling studies indicate less than a century has passed since the most recent common ancestor of the HIV-1 pandemic strains and, in that time frame, an extraordinarily diverse viral population has developed. HIV-1 employs a multitude of schemes to generate variants: accumulation of base substitutions, insertions and deletions, addition and loss of glycosylation sites in the envelope protein, and recombination. A comparison between HIV and influenza virus illustrates the extraordinary scale of HIV variation, and underscores the importance of exploring innovative HIV vaccine strategies. Deeper understanding of the implications of variation for both antibody and T-cell responses may help in the effort to rationally design vaccines that stimulate broad cross-reactivity. The impact of HIV-1 variation on host immune response is reviewed in this context.

The paradox of alloreactivity and self MHC restriction: quantitative analysis and statistics

Although 1-24% of T cells are alloreactive, i.e., respond to MHC molecules encoded by a foreign haplotype, it is generally believed that T cells cannot recognize foreign peptides binding foreign MHC molecules. We show using a quantitative model that, if T cell selection and activation are affinity-driven, then an alloreactivity of 1-24% is incompatible with the textbook notion that self MHC restriction is absolute. If an average of 1% of clones are alloreactive, then according to our model, at most 20-fold more clones should, on average, be activated by antigens presented on self MHC than by antigens presented on foreign MHC. This ratio is at best 5 if alloreactivity is 5%. These results describe average properties of the murine immune system, but not the outcome of individual experiments. Using supercomputer technology, we simulated 100,000 MHC restriction experiments. Although the average restriction ratio was 7.1, restriction was absolute in 10% of the simulated experiments, greater than 100, although not absolute, in 29%, and below 6 in 24%. This extreme variability agrees with experimental estimates. Our analysis suggests that alloreactivity and average self MHC restriction both cannot be high, but that a low average restriction level is compatible with high levels in a significant number of experiments.

Deriving quantitative constraints on T cell selection from data on the mature T cell repertoire

The T cell repertoire is shaped in the thymus through positive and negative selection. Thus, data about the mature repertoire may be used to infer information on how TCR generation and selection operate. Assuming that T cell selection is affinity driven, we derive the quantitative constraints that the parameters driving these processes must fulfill to account for the experimentally observed levels of alloreactivity, self MHC restriction and the frequency of cells recognizing a given foreign Ag. We find that affinity-driven selection is compatible with experimental estimates of these latter quantities only if 1) TCRs see more peptide residues than MHC polymorphic residues, 2) the majority of positively selected clones are deleted by negative selection, 3) between 1 and 3.6 clonal divisions occur on average in the thymus after completion of TCR rearrangement, and 4) selection is driven by 103-105 self peptides.

A quantitative theory of affinity-driven T cell repertoire selection

Binding of the T cell antigen receptor (TCR) to peptides presented on molecules encoded by major histocompatibility complex (MHC) genes is the key event driving T cell development and activation. Selection of the T cell repertoire in the thymus involves two steps. First, positive selection promotes the survival of cells binding thymic self-MHC-peptide complexes with sufficient affinity. The resulting repertoire is self-MHC restricted: it recognizes foreign peptides presented on self, but not foreign MHC. Second, negative selection deletes cells which may be potentially harmful because their receptors interact with self-MHC-peptide complexes with too high an affinity. The mature repertoire is also highly alloreactive: a large fraction of T cells respond to tissues harboring foreign MHC. We derive mathematical expressions giving the frequency of alloreactivity, the level of self-MHC restriction, and the fraction of the repertoire activated by a foreign peptide, as a function of the parameters driving the generation and selection of the repertoire: self-MHC and self-peptide diversity, the stringencies of positive and negative selection, and the number of peptide and MHC polymorphic residues that contribute to T cell receptor binding. Although the model is based on a simplified digit string representation of receptors, all the parameters but one relate directly to experimentally determined quantities. The only parameter without a biological counterpart has no effect on the model's behavior besides a trivial and easily preventable discretization effect. We further analyse the role of the MHC and peptide contribution to TCR binding, and find that their relative, rather than absolute value, is important in shaping the mature repertoire. This result makes it possible to adopt different physical interpretations for the digit string formalism. We also find that the alloreactivity level can be inferred directly from data on the stringency of selection, and that, in agreement with recent experiments, it is not affected by thymic selection.