Thyroid hormone receptor α1: a novel regulator of thyroid cancer cell differentiation.

Thyroid hormone receptor α1 (TRα1) mediates the genomic actions of thyroid hormone (T3). The biology of TRα1 in growth and development has been well studied, but the functional role of TRα1 in cancers remains to be elucidated. Analysis of the human thyroid cancer database of The Cancer Genome Atlas (TCGA) showed that THRA gene expression is lost in highly dedifferentiated anaplastic thyroid cancer (ATC). We, therefore, explored the effects of TRα1 on the progression of ATC. We stably expressed TRα1 in two human ATC cell lines, THJ-11T (11T-TRα1 #2, #7, and #8) and THJ-16T (16T-TRα1 #3, #4, and #8) cells. We found that the expressed TRα1 inhibited ATC cell proliferation and induced apoptosis. TCGA data showed that THRA gene expression was best correlated with the paired box gene 8 (PAX8). Consistently, we found that the PAX8 expression was barely detectable in parental 11T and 16T cells. However, PAX8 gene expression was elevated in 11T- and 16T-TRα1-expressing cells at the mRNA and protein levels. Using various molecular analyses, we found that TRα1 directly regulated the expression of the PAX8 gene. Single-cell transcriptomic analyses (scRNA-seq) demonstrated that TRα1 functions as a transcription factor through multiple signaling pathways to suppress tumor growth. Importantly, scRNA-seq analysis showed that TRα1-induced PAX8, via its transcription program, shifts the cell landscape of ATC toward a differentiated state. The present studies suggest that TRα1 is a newly identified regulator of thyroid differentiation and could be considered as a potential therapeutic target to improve the outcome of ATC patients.

Morphological and Functional Colonic Defects Caused by a Mutated Thyroid Hormone Receptor α.

Background: Mutations of thyroid hormone receptor α (TRα1) result in resistance to thyroid hormone (RTHα), exhibiting symptoms of retarded growth, delayed bone maturation, anemia, and severe constipation. Using a mouse model of RTHα (Thra1PV/+ mouse), we aimed at understanding the molecular basis underlying the severe constipation observed in patients. Methods: The Thra1PV/+ mouse expresses a strong dominant negative mutant, PV, which has lost T3 binding and transcription activity. Thra1PV/+ mouse faithfully reproduces growth abnormalities and anemia as shown in RTHα patients and therefore is a valid model to examine causes of severe constipation in patients. We used histopathological analysis, confocal fluorescence imaging, transmission electron microscopy (TEM), and gene expression profiles to comprehensively analyze the colonic abnormalities of Thra1PV/+ mouse. Results: We found a significant increase in colonic transit time and decrease stool water content in Thra1PV/+ mouse, mimicking constipation as found in patients. Histopathological analysis showed expanded lamina propria filled with interstitium fluid between crypt columns, enlarged muscularis mucosa, and increased content of collagen in expanded submucosa. The TEM analysis revealed shorter muscle fibers with wider gap junctions between muscle cells, fewer caveolae, and hypoplastic interstitial cells of Cajal (ICC) in the rectal smooth muscles of Thra1PV/+ mice. These abnormal histological manifestations suggested defective intercellular transfer of small molecules, electrolytes, and signals for communication among muscles cells, validated by Lucifer Yellow transferring assays. Expression of key smooth muscle contractility regulators, such as calmodulin, myosin light-chain kinase, and phosphorylated myosin light chain, was markedly lower, and c-KIT signaling in ICC was attenuated, resulting in decreased contractility of the rectal smooth muscles of Thra1PV/+ mice. Collectively, these abnormal histopathological alterations and diminished contractility regulators led to the constipation exhibited in patients. Conclusions: This is the first demonstration that TRα1 mutants could act to cause abnormal rectum smooth muscle organization, defects in intercellular exchange of small molecules, and decreased expression of contractility regulators to weaken the contractility of rectal smooth muscles. These findings provide new insights into the molecular basis underlying constipation found in RTHα patients.

Regulation of cancer stem cell activity by thyroid hormone receptor ß.

Increasing numbers of cancer stem cell markers have been recently identified. It is not known, however, whether a member of the nuclear receptor superfamily, thyroid hormone receptor ß (TRß), can function to regulate cancer stem cell (CSC) activity. Using anaplastic thyroid cancer cells (ATC) as a model, we highlight the role of TRß in CSC activity. ATC is one of the most aggressive solid cancers in humans and is resistant to currently available therapeutics. Recent studies provide evidence that CSC activity underlies aggressiveness and therapeutic resistance of ATC. Here we show that TRß inhibits CSC activity by suppressing tumor-sphere formation of human ATC cells and their tumor-initiating capacity. TRß suppresses the expression of CSC regulators, including ALDH, KLF2, SOX2, b-catenin, and ABCG2, in ATC cell-induced xenograft tumors. Single-cell transcriptomic analysis shows that TRß reduces CSC population in ATC-induced xenograft tumors. Analysis of The Cancer Genome Atlas (TCGA) database demonstrates that the inhibition of CSC capacity by TRß contributes to favorable clinical outcomes in human cancer. Our studies show that TRß is a newly identified transcription regulator that acts to suppress CSC activity and that TRß could be considered as a molecular target for therapeutic intervention of ATC.

Distinct DNA Methylation Signatures in Neuroendocrine Tumors Specific for Primary Site and Inherited Predisposition.

PURPOSE: To compare the deoxyribonucleic acid (DNA) methylation signature of neuroendocrine tumors (NETs) by primary tumor site and inherited predisposition syndromes von Hippel-Lindau disease (VHL) and multiple endocrine neoplasia type 1 (MEN1). METHODS: Genome-wide DNA methylation (835 424 CpGs) of 96 NET samples. Principal components analysis (PCA) and unsupervised hierarchical clustering analyses were used to determine DNA methylome signatures. RESULTS: Hypomethylated CpGs were significantly more common in VHL-related versus sporadic and MEN1-related NETs (P < .001 for both comparisons). Small-intestinal NETs (SINETs) had the most differentially methylated CpGs, either hyper- or hypomethylated, followed by duodenal NETs (DNETs) and pancreatic NETs (PNETs, P < .001 for all comparisons). There was complete separation of SINETs on PCA, and 3 NETs of unknown origin clustered with the SINET samples. Sporadic, VHL-related, and MEN1-related PNETs formed distinct groups on PCA, and VHL clustered separately, showing pronounced DNA hypomethylation, while sporadic and MEN1-related NETs clustered together. MEN1-related PNETs, DNETs, and gastric NETs each had a distinct DNA methylome signature, with complete separation by PCA and unsupervised clustering. Finally, we identified 12 hypermethylated CpGs in the 1A promoter of the APC (adenomatous polyposis coli) gene, with higher methylation levels in MEN1-related NETs versus VHL-related and sporadic NETs (P < .001 for both comparisons). CONCLUSIONS: DNA CpG methylation profiles are unique in different primary NET types even when occurring in MEN1-related NETs. This tumor DNA methylome signature may be utilized for noninvasive molecular characterization of NETs, through DNA methylation profiling of biopsy samples or even circulating tumor DNA in the near future.

Author Correction: Recurrent PTPRT/JAK2 mutations in lung adenocarcinoma among African Americans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Generation of Novel Genetic Models to Dissect Resistance to Thyroid Hormone Receptor α in Zebrafish.

Background: Patients with mutations of the thyroid hormone receptor alpha (THRA) gene show resistance to thyroid hormone alpha (RTHα). No amendable mouse models are currently available to elucidate deleterious effects of TRα1 mutants during early development. Zebrafish with transient suppressed expression by morpholino knockdown and ectopic expression of TRα1 mutants in the embryos have been reported. However, zebrafish with germline transmittable mutations have not been reported. The stable expression of thra mutants from embryos to adulthood facilitated the study of molecular actions of TRα1 mutants during development. Methods: In contrast to human and mice, the thra gene is duplicated in zebrafish, thraa, and thrab. Using CRISPR/Cas9-mediated targeted mutagenesis, we created dominant negative mutations in the two duplicated thra genes. We comprehensively analyzed the molecular and phenotypic characteristics of mutant fish during development. Results: Adult and juvenile homozygous thrab 1-bp ins (m/m) mutants exhibited severe growth retardation, but adult homozygous thraa 8-bp ins (m/m) mutants had very mild growth impairment. Expression of the growth hormone (gh1) and insulin-like growth factor 1 was markedly suppressed in homozygous thrab 1-bp ins (m/m) mutants. Decreased messenger RNA and protein levels of triiodothyronine-regulated keratin genes and inhibited keratinocyte proliferation resulted in hypoplasia of the epidermis in adult and juvenile homozygous thrab 1-bp ins (m/m) mutants, but not homozygous thraa 8-bp ins (m/m) mutants. RNA-seq analysis showed that homozygous thrab 1-bp ins (m/m) mutation had global impact on the functions of the adult pituitary. However, no morphological defects nor any changes in the expression of gh1 and keratin genes were observed in the embryos and early larvae. Thus, mutations of either the thraa or thrab gene did not affect initiation of embryogenesis. But the mutation of the thrab gene, but not the thraa gene, is detrimental in postlarval growth and skin development. Conclusions: The thra duplicated genes are essential to control temporal coordination in postlarval growth and development in a tissue-specific manner. We uncovered novel functions of the duplicated thra genes in zebrafish in development. These mutant zebrafish could be used as a model for further analysis of TRα1 mutant actions and for rapid screening of therapeutics for RTHα.

ONCOGENE PANEL SEQUENCING ANALYSIS IDENTIFIES CANDIDATE ACTIONABLE GENES IN ADVANCED WELL-DIFFERENTIATED GASTROENTEROPANCREATIC NEUROENDOCRINE TUMORS.

Objective: To report the rate of candidate actionable somatic mutations in patients with locally advanced and metastatic gastro-enteropancreatic (GEP) neuroendocrine tumors (NET) and of other genetic alterations that may be associated with tumorigenesis. Methods: A phase II mutation targeted therapy trial was conducted in patients with advanced well-differentiated G1/G2 GEP-NET. Mutations found in the mTOR pathway-associated genes led to treatment with the mTOR inhibitor everolimus, and were defined as actionable. Tumor deoxyribonucleic acid (DNA) from GEP-NET were sequenced and compared with germline DNA, using the OncoVAR-NET assay, designed for hybrid capture sequencing of 500 tumor suppressor genes and oncogenes. Somatic variants were called and copy-number (CN) variant analysis was performed. Results: Thirty patients (14 small-intestine, 8 pancreatic, 3 unknown primary NET, and 5 of other primary sites) harbored 37 lesions (4 patients had DNA of multiple lesions sequenced). Only 2 patients with sporadic NET (n = 26) had an actionable mutation leading to treatment with everolimus. Driver somatic mutations were detected in 18 of 30 patients (21/37 lesions sequenced). In the remaining samples without a driver mutation, CN alterations were found in 11/16 tumors (10/12 patients), including CN loss of chromosome (Chr) 18 (P<.05), CN gain of Chr 5, and loss of Chr 13. CN losses in Chr 18 were more common in patients without driver mutations detected. Pronounced genetic heterogeneity was detected in patients with multiple lesions sequenced. Conclusion: Genome-wide DNA sequencing may identify candidate actionable genes and lead to the identification of novel target genes for advanced well-differentiated GEP-NET. Abbreviations: Chr = chromosome; CN = copy number; DNA = deoxyribonucleic acid; FDA = Food and Drug Administration; GEP = gastro-enteropancreatic; MEN-1 = multiple endocrine neoplasia syndrome type 1; mTOR = mammalian target of rapamycin; NET = neuroendocrine tumor; PFS = progression-free survival; PNET = pancreatic neuroendocrine tumors; SINET = small-intestine neuroendocrine tumor.

Aneuploidy, TP53 mutation, and amplification of MYC correlate with increased intratumor heterogeneity and poor prognosis of breast cancer patients.

The clinical course of breast cancer varies from one patient to another. Currently, the choice of therapy relies on clinical parameters and histological and molecular tumor features. Alas, these markers are informative in only a subset of patients. Therefore, additional predictors of disease outcome would be valuable for treatment stratification. Extensive studies showed that the degree of variation of the nuclear DNA content, i.e., aneuploidy, determines prognosis. Our aim was to further elucidate the molecular basis of aneuploidy. We analyzed five diploid and six aneuploid tumors with more than 20 years of follow-up. By performing FISH with a multiplexed panel of 10 probes to enumerate copy numbers in individual cells, and by sequencing 563 cancer-related genes, we analyzed how aneuploidy is linked to intratumor heterogeneity. In our cohort, none of the patients with diploid tumors died of breast cancer during follow-up in contrast to four of six patients with aneuploid tumors (mean survival 86.4 months). The FISH analysis showed markedly increased genomic instability and intratumor heterogeneity in aneuploid tumors. MYC gain was observed in only 20% of the diploid cancers, while all aneuploid cases showed a gain. The mutation burden was similar in diploid and aneuploid tumors, however, TP53 mutations were not observed in diploid tumors, but in all aneuploid tumors in our collective. We conclude that quantitative measurements of intratumor heterogeneity by multiplex FISH, detection of MYC amplification and TP53 mutation could augment prognostication in breast cancer patients.

Recurrent epimutation of SDHC in gastrointestinal stromal tumors.

Succinate dehydrogenase (SDH) is a conserved effector of cellular metabolism and energy production, and loss of SDH function is a driver mechanism in several cancers. SDH-deficient gastrointestinal stromal tumors (dSDH GISTs) collectively manifest similar phenotypes, including hypermethylated epigenomic signatures, tendency to occur in pediatric patients, and lack of KIT/PDGFRA mutations. dSDH GISTs often harbor deleterious mutations in SDH subunit genes (SDHA, SDHB, SDHC, and SDHD, termed SDHx), but some are SDHx wild type (WT). To further elucidate mechanisms of SDH deactivation in SDHx-WT GIST, we performed targeted exome sequencing on 59 dSDH GISTs to identify 43 SDHx-mutant and 16 SDHx-WT cases. Genome-wide DNA methylation and expression profiling exposed SDHC promoter-specific CpG island hypermethylation and gene silencing in SDHx-WT dSDH GISTs [15 of 16 cases (94%)]. Six of 15 SDHC-epimutant GISTs occurred in the setting of the multitumor syndrome Carney triad. We observed neither SDHB promoter hypermethylation nor large deletions on chromosome 1q in any SDHx-WT cases. Deep genome sequencing of a 130-kbp (kilo-base pair) window around SDHC revealed no recognizable sequence anomalies in SDHC-epimutant tumors. More than 2000 benign and tumor reference tissues, including stem cells and malignancies with a hypermethylator epigenotype, exhibit solely a non-epimutant SDHC promoter. Mosaic constitutional SDHC promoter hypermethylation in blood and saliva from patients with SDHC-epimutant GIST implicates a postzygotic mechanism in the establishment and maintenance of SDHC epimutation. The discovery of SDHC epimutation provides a unifying explanation for the pathogenesis of dSDH GIST, whereby loss of SDH function stems from either SDHx mutation or SDHC epimutation.

Immunohistochemical loss of succinate dehydrogenase subunit A (SDHA) in gastrointestinal stromal tumors (GISTs) signals SDHA germline mutation.

A subset (7% to 10%) of gastric gastrointestinal stromal tumors (GISTs) is notable for the immunohistochemical loss of succinate dehydrogenase (SDH) subunit B (SDHB), which signals the loss of function of the SDH complex consisting of mitochondrial inner membrane proteins. These SDH-deficient GISTs are known to be KIT/PDGFRA wild type, and most patients affected by this subset of GISTs are young. Some of these patients have germline mutations of SDH subunit genes SDHB, SDHC, or SDHD, known as Carney-Stratakis syndrome when combined with paraganglioma. More recently, germline mutations in SDH subunit A gene (SDHA) have also been reported in few patients with KIT/PDGFRA wild-type GISTs. In this study we immunohistochemically examined 127 SDHB-negative and 556 SDHB-positive gastric GISTs and 261 SDHB-positive intestinal GISTs for SDHA expression using a mouse monoclonal antibody 2E3 (Abcam). Cases with available DNA were tested for SDHA, SDHB, SDHC, and SDHD gene mutations using a hybridization-based custom capture next-generation sequencing assay. A total of 36 SDHA-negative GISTs (28%) were found among 127 SDHB-negative gastric GISTs. No SDHB-positive GIST was SDHA negative. Among 7 SDHA-negative tumors analyzed, there were 7 SDHA mutants, most germline. A second hit indicating biallelic inactivation of SDHA was present in 6 of those cases. These patients had no other SDH subunit gene mutations. Among the 25 SDHA-positive, SDHB-negative GISTs analyzed, we identified 3 SDHA mutations (1 germline), and 11 SDHB, SDHC, or SDHD mutations (mostly germline), and 11 patients with no SDH mutations. Compared with patients with SDHA-positive GISTs, those with SDHA-negative GISTs had an older median age (34 vs. 21 y), lower female to male ratio (1.8 vs. 3.1) but similar mitotic counts and median tumor sizes, with a slow course of disease in most cases, despite a slightly higher rate of liver metastases. SDHA-negative GISTs comprise approximately 30% of SDHB-negative/SDH-deficient GISTs, and SDHA loss generally correlates with SDHA mutations.

Two types of human malignant melanoma cell lines revealed by expression patterns of mitochondrial and survival-apoptosis genes: implications for malignant melanoma therapy.

Human malignant melanoma has poor prognosis because of resistance to apoptosis and therapy. We describe identification of the expression profile of 1,037 mitochondria-focused genes and 84 survival-apoptosis genes in 21 malignant melanoma cell lines and 3 normal melanocyte controls using recently developed hMitChip3 cDNA microarrays. Unsupervised hierarchical clustering analysis of 1,037 informative genes, and 84 survival-apoptosis genes, classified these malignant melanoma cell lines into type A (n = 12) and type B (n = 9). Three hundred fifty-five of 1,037 (34.2%) genes displayed significant (P ≤ 0.030; false discovery rate ≤ 3.68%) differences (± ≥ 2.0-fold) in average expression, with 197 genes higher and 158 genes lower in type A than in type B. Of 84 genes with known survival-apoptosis functions, 38 (45.2%) displayed the significant (P < 0.001; false discovery rate < 0.15%) difference. Antiapoptotic (BCL2, BCL2A1, PPARD, and RAF1), antioxidant (MT3, PRDX5, PRDX3, GPX4, GLRX2, and GSR), and proapoptotic (BAD, BNIP1, APAF1, BNIP3L, CASP7, CYCS, CASP1, and VDAC1) genes expressed at higher levels in type A than in type B, whereas the different set of antiapoptotic (PSEN1, PPP2CA, API5, PPP2R1B, PPP2R1A, and FIS1), antioxidant (HSPD1, GSS, SOD1, ATOX1, and CAT), and proapoptotic (ENDOG, BAK1, CASP2, CASP4, PDCD5, HTRA2, SEPT4, TNFSF10, and PRODH) genes expressed at lower levels in type A than in type B. Microarray data were validated by quantitative reverse transcription-PCR. These results showed the presence of two types of malignant melanoma, each with a specific set of dysregulated survival-apoptosis genes, which may prove useful for development of new molecular targets for therapeutic intervention and novel diagnostic biomarkers for treatment and prognosis of malignant melanoma.