Abstract 1320: The landscape of somatic substitutions in the repetitive genome across cancer types

Somatic mutation calling from bulk DNA sequencing is a complex problem susceptible to elevated false positive rates. High mapping quality is considered an important feature of reliable variant calls. At the resolution of short reads, ~10% of the genome displays high sequence similarity with at least one other genomic region and is assigned low mapping quality by alignment algorithms. These low-mapping-quality regions represent recurrent blind spots for mutation callers, which discard many of the variants they harbor, overlooking true biological variation.

Here, we developed a pipeline to call substitutions in the low-mapping-quality genome. We used a published thesaurus approach to annotate the variant positions with their high-similarity links. We trained a classifier to emulate high-quality consensus calls made in unique regions using 20 features unrelated to mapping quality, reaching ~95% accuracy in those regions. In an independent sample more than 90% of the thesaurus calls were validated through linked-read sequencing. We then applied the classifier to all candidate substitutions of 2,658 cancer whole genomes from the PCAWG/ICGC consortium including variants falling in low-mapability regions.

We retrieve hidden thesaurus variants genome-wide in ~6% of the genome, including genic, coding, and promoter regions. Thesaurus calls are directly proportional in numbers to somatic calls falling in the low-mapping-quality genome and share a similar trinucleotide context spectrum. Rescuing these mutations reveal hidden signal in known cancer genes, including PIK3CA, and excess of mutations genome-wide in promoter, untranslated, and coding regions of many other genes. We also find potential excess of non-synonymous mutations, including in genes from the TRIM and POTE families, having been previously implicated in multiple cancer types.

Altogether, we developed a pipeline to call somatic substitutions in the low-mapping-quality genome and uncovered hidden somatic changes along the genomes of human cancers. In the future, this pipeline could be extended to indels and structural variants, and applied to the study of de novo germline variants.

Citation Format: Maxime Tarabichi, Jonas Demeulemeester, Annelien Verfaillie, Peter Van Loo, Tomasz Konopka. The landscape of somatic substitutions in the repetitive genome across cancer types [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1320.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

The evolutionary history of 2,658 cancers

Cancer develops through a process of somatic evolution1,2. Sequencing data from a single biopsy represent a snapshot of this process that can reveal the timing of specific genomic aberrations and the changing influence of mutational processes3. Here, by whole-genome sequencing analysis of 2,658 cancers as part of the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA)4, we reconstruct the life history and evolution of mutational processes and driver mutation sequences of 38 types of cancer. Early oncogenesis is characterized by mutations in a constrained set of driver genes, and specific copy number gains, such as trisomy 7 in glioblastoma and isochromosome 17q in medulloblastoma. The mutational spectrum changes significantly throughout tumour evolution in 40% of samples. A nearly fourfold diversification of driver genes and increased genomic instability are features of later stages. Copy number alterations often occur in mitotic crises, and lead to simultaneous gains of chromosomal segments. Timing analyses suggest that driver mutations often precede diagnosis by many years, if not decades. Together, these results determine the evolutionary trajectories of cancer, and highlight opportunities for early cancer detection.

A community effort to create standards for evaluating tumor subclonal reconstruction

Tumor DNA sequencing data can be interpreted by computational methods that analyze genomic heterogeneity to infer evolutionary dynamics. A growing number of studies have used these approaches to link cancer evolution with clinical progression and response to therapy. Although the inference of tumor phylogenies is rapidly becoming standard practice in cancer genome analyses, standards for evaluating them are lacking. To address this need, we systematically assess methods for reconstructing tumor subclonality. First, we elucidate the main algorithmic problems in subclonal reconstruction and develop quantitative metrics for evaluating them. Then we simulate realistic tumor genomes that harbor all known clonal and subclonal mutation types and processes. Finally, we benchmark 580 tumor reconstructions, varying tumor read depth, tumor type and somatic variant detection. Our analysis provides a baseline for the establishment of gold-standard methods to analyze tumor heterogeneity.

Investigation of somatic CNVs in brains of synucleinopathy cases using targeted SNCA analysis and single cell sequencing

Synucleinopathies are mostly sporadic neurodegenerative disorders of partly unexplained aetiology, and include Parkinson's disease (PD) and multiple system atrophy (MSA). We have further investigated our recent finding of somatic SNCA (α-synuclein) copy number variants (CNVs, specifically gains) in synucleinopathies, using Fluorescent in-situ Hybridisation for SNCA, and single-cell whole genome sequencing for the first time in a synucleinopathy. In the cingulate cortex, mosaicism levels for SNCA gains were higher in MSA and PD than controls in neurons (> 2% in both diseases), and for MSA also in non-neurons. In MSA substantia nigra (SN), we noted SNCA gains in > 3% of dopaminergic (DA) neurons (identified by neuromelanin) and neuromelanin-negative cells, including olig2-positive oligodendroglia. Cells with CNVs were more likely to have α-synuclein inclusions, in a pattern corresponding to cell categories mostly relevant to the disease: DA neurons in Lewy-body cases, and other cells in the striatonigral degeneration-dominant MSA variant (MSA-SND). Higher mosaicism levels in SN neuromelanin-negative cells may correlate with younger onset in typical MSA-SND, and in cingulate neurons with younger death in PD. Larger sample sizes will, however, be required to confirm these putative findings. We obtained genome-wide somatic CNV profiles from 169 cells from the substantia nigra of two MSA cases, and pons and putamen of one. These showed somatic CNVs in ~ 30% of cells, with clonality and origins in segmental duplications for some. CNVs had distinct profiles based on cell type, with neurons having a mix of gains and losses, and other cells having almost exclusively gains, although control data sets will be required to determine possible disease relevance. We propose that somatic SNCA CNVs may contribute to the aetiology and pathogenesis of synucleinopathies, and that genome-wide somatic CNVs in MSA brain merit further study.

2019 John Insall Award: Fructosamine is a better glycaemic marker compared with glycated haemoglobin (HbA1C) in predicting adverse outcomes following total knee arthroplasty: a prospective multicentre study

AIMS

The best marker for assessing glycaemic control prior to total knee arthroplasty (TKA) remains unknown. The purpose of this study was to assess the utility of fructosamine compared with glycated haemoglobin (HbA1c) in predicting early complications following TKA, and to determine the threshold above which the risk of complications increased markedly.

PATIENTS AND METHODS

This prospective multi-institutional study evaluated primary TKA patients from four academic institutions. Patients (both diabetics and non-diabetics) were assessed using fructosamine and HbA1c levels within 30 days of surgery. Complications were assessed for 12 weeks from surgery and included prosthetic joint infection (PJI), wound complication, re-admission, re-operation, and death. The Youden's index was used to determine the cut-off for fructosamine and HbA1c associated with complications. Two additional cut-offs for HbA1c were examined: 7% and 7.5% and compared with fructosamine as a predictor for complications.

RESULTS

Overall, 1119 patients (441 men, 678 women) were included in the study. Fructosamine level of 293 µmol/l was identified as the optimal cut-off associated with complications. Patients with high fructosamine (> 293 µmol/l) were 11.2 times more likely to develop PJI compared with patients with low fructosamine (p = 0.001). Re-admission and re-operation rates were 4.2 and 4.5 times higher in patients with fructosamine above the threshold (p = 0.005 and p = 0.019, respectively). One patient (1.7%) from the elevated fructosamine group died compared with one patient (0.1%) in the normal fructosamine group (p = 0.10). These complications remained statistically significant in multiple regression analysis. Unlike fructosamine, all three cut-offs for HbA1c failed to show a significant association with complications.

CONCLUSION

Fructosamine is a valid and an excellent predictor of complications following TKA. It better reflects the glycaemic control, has greater predictive power for adverse events, and responds quicker to treatment compared with HbA1c. These findings support the screening of all patients undergoing TKA using fructosamine and in those with a level above 293 µmol/l, the risk of surgery should be carefully weighed against its benefit. Cite this article: Bone Joint J 2019;101-B(7 Supple C):3-9.

Undifferentiated Sarcomas Develop through Distinct Evolutionary Pathways

Undifferentiated sarcomas (USARCs) of adults are diverse, rare, and aggressive soft tissue cancers. Recent sequencing efforts have confirmed that USARCs exhibit one of the highest burdens of structural aberrations across human cancer. Here, we sought to unravel the molecular basis of the structural complexity in USARCs by integrating DNA sequencing, ploidy analysis, gene expression, and methylation profiling. We identified whole genome duplication as a prevalent and pernicious force in USARC tumorigenesis. Using mathematical deconvolution strategies to unravel the complex copy-number profiles and mutational timing models we infer distinct evolutionary pathways of these rare cancers. In addition, 15% of tumors exhibited raised mutational burdens that correlated with gene expression signatures of immune infiltration, and good prognosis.

Abstract 218: The evolutionary history of 2,658 cancers

Cancer develops through a continuous process of somatic evolution. Whole genome sequencing provides a snapshot of the tumor genome at the point of sampling, however, the data can contain information that permits the reconstruction of a tumor's evolutionary past.

Here, we apply such life history analyses on an unprecedented scale, to a set of 2,658 tumors spanning 39 cancer types. We estimated the timing of large chromosomal gains during tumor evolution, by comparing the rates of doubled to non-doubled point mutations within gained regions. Although we find that such events typically occur in the second half of clonal evolution, we also observe distinctive and early chromosomal gains in some cancer types, such as gains of chromosomes 7, 19 and 20 in glioblastoma, and isochromosome 17q in medulloblastoma. By integrating these results with the qualitative timing of individual driver mutations, we obtained an overall ranking, from early to late, of frequent somatic events per cancer type, which both identified novel patterns of tumor evolution, and incorporated additional detail into known models, such as the progression of APC-KRAS-TP53 in colorectal cancer proposed by Vogelstein and Fearon.

To estimate how mutational processes acting on the tumor genome change over time, we classified mutations in each sample according to three broad time periods (early clonal, late clonal, and subclonal), and quantified the activity of mutational signatures in each period. Most mutational processes appear to remain remarkably constant, however, certain signatures show clear and consistent changes during clonal evolution. Particularly, mutational signatures associated with exposure to carcinogens, such as smoking and UV light, tend to decrease over time. In contrast, signatures associated with defective endogenous processes, such as APOBEC mutagenesis and defective double strand break repair, show an increase between early and late phases of tumor evolution.

Making use of clock-like mutational signatures, we converted mutational time estimates for large events, such as whole genome duplication (WGD), and the emergence of the most recent common ancestor (MRCA), into real time estimates, which allowed us to combine our analyses into overall timelines of cancer evolution, per tumor type. For example, the typical timeline of ovarian adenocarcinoma development shows that early tumor evolution is characterized by mutations in TP53, and widespread genome instability, with WGD events taking place on average 8 years prior to diagnosis. In later stages of evolution, signatures of defective repair processes increase, and the MRCA emerges on average 1 year before diagnosis.

Taken together, these data reveal the common and divergent evolutionary trajectories available to a cancer, which might be crucial in understanding specific tumor biology, and in providing new opportunities for early detection and cancer prevention.

Citation Format: Clemency Jolly, Moritz Gerstung, Ignaty Leshchiner, Stefan C. Dentro, Santiago Gonzalez, Thomas J. Mitchell, Yulia Rubanova, Pavana Anur, Daniel Rosebrock, Kaixian Yu, Maxime Tarabichi, Amit Deshwar, Jeff Wintersinger, Kortine Kleinheinz, Ignacio Vásquez-García, Kerstin Haase, Subhajit Sengupta, Geoff Macintyre, Salem Malikic, Nilgun Donmez, Dimitri G. Livitz, Mark Cmero, Jonas Demeulemeester, Steve Schumacher, Yu Fan, Xiaotong Yao, Juhee Lee, Matthias Schlesner, Paul C. Boutros, David D. Bowtell, Hongtu Zhu, Gad Getz, Marcin Imielinski, Rameen Beroukhim, S Cenk Sahinalp, Yuan Ji, Martin Peifer, Florian Markowetz, Ville Mustonen, Ke Juan, Wenyi Wang, Quaid D. Morris, Paul T. Spellman, David C. Wedge, Peter Van Loo, PCAWG Evolution and Heterogeneity Working Group. The evolutionary history of 2,658 cancers [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 218.

Recurrent rearrangements of FOS and FOSB define osteoblastoma

The transcription factor FOS has long been implicated in the pathogenesis of bone tumours, following the discovery that the viral homologue, v-fos, caused osteosarcoma in laboratory mice. However, mutations of FOS have not been found in human bone-forming tumours. Here, we report recurrent rearrangement of FOS and its paralogue, FOSB, in the most common benign tumours of bone, osteoblastoma and osteoid osteoma. Combining whole-genome DNA and RNA sequences, we find rearrangement of FOS in five tumours and of FOSB in one tumour. Extending our findings into a cohort of 55 cases, using FISH and immunohistochemistry, provide evidence of ubiquitous mutation of FOS or FOSB in osteoblastoma and osteoid osteoma. Overall, our findings reveal a human bone tumour defined by mutations of FOS and FOSB.

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.

Can next generation sequencing play a role in detecting pathogens in synovial fluid?

Aims

The diagnosis of periprosthetic joint infection can be difficult due to the high rate of culture-negative infections. The aim of this study was to assess the use of next-generation sequencing for detecting organisms in synovial fluid.

Materials and Methods

In this prospective, single-blinded study, 86 anonymized samples of synovial fluid were obtained from patients undergoing aspiration of the hip or knee as part of the investigation of a periprosthetic infection. A panel of synovial fluid tests, including levels of C-reactive protein, human neutrophil elastase, total neutrophil count, alpha-defensin, and culture were performed prior to next-generation sequencing.

Results

Of these 86 samples, 30 were alpha-defensin-positive and culture-positive (Group I), 24 were alpha-defensin-positive and culture-negative (Group II) and 32 were alpha-defensin-negative and culture-negative (Group III). Next-generation sequencing was concordant with 25 results for Group I. In four of these, it detected antibiotic resistant bacteria whereas culture did not. In another four samples with relatively low levels of inflammatory biomarkers, culture was positive but next-generation sequencing was negative. A total of ten samples had a positive next-generation sequencing result and a negative culture. In five of these, alpha-defensin was positive and the levels of inflammatory markers were high. In the other five, alpha-defensin was negative and the levels of inflammatory markers were low. While next-generation sequencing detected several organisms in each sample, in most samples with a higher probability of infection, there was a predominant organism present, while in those presumed not to be infected, many organisms were identified with no predominant organism.

Conclusion

Pathogens causing periprosthetic infection in both culture-positive and culture-negative samples of synovial fluid could be identified by next-generation sequencing. Cite this article: Bone Joint J 2018;100-B:127–33.

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.

The human bitter taste receptor T2R38 is broadly tuned for bacterial compounds

T2R38 has been shown to be a specific bacterial detector implicated in innate immune defense mechanism of human upper airway. Several clinical studies have demonstrated that this receptor is associated with the development of chronic rhinosinusitis (CRS). T2R38 was previously reported to bind to homoserine lactones (HSL), quorum sensing molecules specific of Pseudomonas Aeruginosa and other gram negative species. Nevertheless, these bacteria are not the major pathogens found in CRS. Here we report on the identification of bacterial metabolites acting as new agonists of T2R38 based on a single cell calcium imaging study. Two quorum sensing molecules (Agr D1 thiolactone from Staphylococcus Aureus and CSP-1 from Streptococcus Pneumoniae) and a list of 32 bacterial metabolites from pathogens frequently implicated in CRS were tested. First, we observed that HSL failed to activate T2R38 in our experimental system, but that the dimethylsulfoxide (DMSO), used as a solvent for these lactones may, by itself, account for the agonistic effect previously described. Secondly, we showed that both Agr D1 thiolactone and CSP-1 are inactive but that at least 7 bacterial metabolites (acetone, 2-butanone, 2-pentanone, 2-methylpropanal, dimethyl disulfide, methylmercaptan, γ-butyrolactone) are able to specifically trigger this receptor. T2R38 is thus much more broadly tuned for bacterial compounds than previously thought.

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.

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.

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.

Profiling of olfactory receptor gene expression in whole human olfactory mucosa

Olfactory perception is mediated by a large array of olfactory receptor genes. The human genome contains 851 olfactory receptor gene loci. More than 50% of the loci are annotated as nonfunctional due to frame-disrupting mutations. Furthermore haplotypic missense alleles can be nonfunctional resulting from substitution of key amino acids governing protein folding or interactions with signal transduction components. Beyond their role in odor recognition, functional olfactory receptors are also required for a proper targeting of olfactory neuron axons to their corresponding glomeruli in the olfactory bulb. Therefore, we anticipate that profiling of olfactory receptor gene expression in whole human olfactory mucosa and analysis in the human population of their expression should provide an opportunity to select the frequently expressed and potentially functional olfactory receptors in view of a systematic deorphanization. To address this issue, we designed a TaqMan Low Density Array (Applied Biosystems), containing probes for 356 predicted human olfactory receptor loci to investigate their expression in whole human olfactory mucosa tissues from 26 individuals (13 women, 13 men; aged from 39 to 81 years, with an average of 67±11 years for women and 63±12 years for men). Total RNA isolation, DNase treatment, RNA integrity evaluation and reverse transcription were performed for these 26 samples. Then 384 targeted genes (including endogenous control genes and reference genes specifically expressed in olfactory epithelium for normalization purpose) were analyzed using the same real-time reverse transcription PCR platform. On average, the expression of 273 human olfactory receptor genes was observed in the 26 selected whole human olfactory mucosa analyzed, of which 90 were expressed in all 26 individuals. Most of the olfactory receptors deorphanized to date on the basis of sensitivity to known odorant molecules, which are described in the literature, were found in the expressed olfactory receptors gene set.

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.

5-aza-2'-deoxycytidine has minor effects on differentiation in human thyroid cancer cell lines, but modulates genes that are involved in adaptation in vitro

BACKGROUND

In thyroid cancer, the lack of response to specific treatment, for example, radioactive iodine, can be caused by a loss of differentiation characteristics of tumor cells. It is hypothesized that this loss is due to epigenetic modifications. Therefore, drugs releasing epigenetic repression have been proposed to reverse this silencing.

METHODS

We investigated which genes were reinduced in dedifferentiated human thyroid cancer cell lines when treated with the demethylating agent 5-aza-2'-deoxycytidine (5-AzadC) and the histone deacetylase inhibitors trichostatin A (TSA) and suberoylanilide hydroxamic acid, by using reverse transcriptase-polymerase chain reaction and microarrays. These results were compared to the expression patterns in in vitro human differentiated thyrocytes and in in vivo dedifferentiated thyroid cancers. In addition, the effects of 5-AzadC on DNA quantities and cell viability were investigated.

RESULTS

Among the canonical thyroid differentiation markers, most were not, or only to a minor extent, re-expressed by 5-AzadC, whether or not combined with TSA or forskolin, an inducer of differentiation in normal thyrocytes. Furthermore, 5-AzadC-modulated overall mRNA expression profiles showed only few commonly regulated genes compared to differentiated cultured primary thyrocytes. In addition, most of the commonly strongly 5-AzadC-induced genes in cell lines were either not regulated or upregulated in anaplastic thyroid carcinomas. Further analysis of which genes were induced by 5-AzadC showed that they were involved in pathways such as apoptosis, antigen presentation, defense response, and cell migration. A number of these genes had similar expression responses in 5-AzadC-treated nonthyroid cell lines.

CONCLUSIONS

Our results suggest that 5-AzadC is not a strong inducer of differentiation in thyroid cancer cell lines. Under the studied conditions and with the model used, 5-AzadC treatment does not appear to be a potential redifferentiation treatment for dedifferentiated thyroid cancer. However, this may reflect primarily the inadequacy of the model rather than that of the treatment. Moreover, the observation that 5-AzadC negatively affected cell viability in cell lines could still suggest a therapeutic opportunity. Some of the genes that were modulated by 5-AzadC were also induced in nonthyroid cancer cell lines, which might be explained by an epigenetic modification resulting in the adaptation of the cell lines to their culture conditions.

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/.