Human-Specific NOTCH2NL Genes Expand Cortical Neurogenesis through Delta/Notch Regulation.

The cerebral cortex underwent rapid expansion and increased complexity during recent hominid evolution. Gene duplications constitute a major evolutionary force, but their impact on human brain development remains unclear. Using tailored RNA sequencing (RNA-seq), we profiled the spatial and temporal expression of hominid-specific duplicated (HS) genes in the human fetal cortex and identified a repertoire of 35 HS genes displaying robust and dynamic patterns during cortical neurogenesis. Among them NOTCH2NL, human-specific paralogs of the NOTCH2 receptor, stood out for their ability to promote cortical progenitor maintenance. NOTCH2NL promote the clonal expansion of human cortical progenitors, ultimately leading to higher neuronal output. At the molecular level, NOTCH2NL function by activating the Notch pathway through inhibition of cis Delta/Notch interactions. Our study uncovers a large repertoire of recently evolved genes active during human corticogenesis and reveals how human-specific NOTCH paralogs may have contributed to the expansion of the human cortex.

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.

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.

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.

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.

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.

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

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; rs = 0.973; 95% CI: 0.971-0.975), progesterone receptor (PgR; rs = 0.95; 0.947-0.954), and human epidermal growth factor receptor 2 (HER2; rs = 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 rs >0.965). Immune-, stroma- and pathway-based GES showed a lower correlation relative to prognostic signatures (all rs >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.

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.

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(2) = 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.

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.

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.

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.

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.

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.

Factors contributing to the resistance of the thyrocyte to hydrogen peroxide.

We studied the mechanism that may explain the relative resistance of thyrocytes to H2O2 compared to other cell types. Ability to degrade H2O2, glutathione peroxidase (GPx) activity, heme oxygenase-1 (HO-1) expression, cell survival and capacity to repair DNA damage after H2O2 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 H2O2 exposure, rapidly degraded extracellular H2O2 and presented a high basal seleno-dependent GPx activity. Only in thyrocytes, H2O2 up-regulated GPx activity and expression of HO-1 mRNA. These effects were not reproduced by irradiation. DNA damage caused by H2O2 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 H2O2 and its mutagenic effects.

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 TILhi (≥50%). Interestingly, the BRCAmut group had a significantly higher incidence of TILpos tumors compared to the BRCAwt group (P = 0.037). T cells were dominant in the infiltrate but no statistically significant differences were detected between BRCAwt and BRCAmut 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 BRCAwt and BRCAmut TNBC are similar except for a significant increase of TILpos tumors in the BRCAmut group. While BRCA gene mutations may not directly drive immune infiltration, the greater number of TILpos tumors could signal greater immunogenicity in this group.

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.