Global RNA Expression and DNA Methylation Patterns in Primary Anaplastic Thyroid Cancer

Prediction of the Aggressive Clinical Course of Papillary Thyroid Carcinoma Based on Fine Needle Aspiration Biopsy Molecular Testing

Molecular genetic events are among the numerous factors affecting the clinical course of papillary thyroid carcinoma (PTC). Recent studies have demonstrated that aberrant expression of miRNA, as well as different thyroid-related genes, correlate with the aggressive clinical course of PTC and unfavorable treatment outcomes, which opens up new avenues for using them in the personalization of the treatment strategy for patients with PTC. In the present work, our goal was to assess the applicability of molecular markers in the preoperative diagnosis of aggressive variants of papillary thyroid cancer. The molecular genetic profile (expression levels of 34 different markers and BRAF mutations) was studied for 108 cytology specimens collected by fine-needle aspiration biopsy in patients with PTC having different clinical manifestations. Statistically significant differences with adjustment for multiple comparisons (p < 0.0015) for clinically aggressive variants of PTC were obtained for four markers: miRNA-146b, miRNA-221, fibronectin 1 (FN1), and cyclin-dependent kinase inhibitor 2A (CDKN2A) genes. A weak statistical correlation (0.0015 < p < 0.05) was observed for miRNA-31, -375, -551b, -148b, -125b, mtDNA, CITED1, TPO, HMGA2, CLU, NIS, SERPINA1, TFF3, and TMPRSS4. The recurrence risk of papillary thyroid carcinoma can be preoperatively predicted using miRNA-221, FN1, and CDKN2A genes.

Bioinformatics Analysis of Human Papillomavirus 16 Integration in Cervical Cancer: Changes in MAGI-1 Expression in Premalignant Lesions and Invasive Carcinoma

HPV 16 integration is crucial for the onset and progression of premalignant lesions to invasive squamous cell carcinoma (ISCC) because it promotes the amplification of proto-oncogenes and the silencing of tumor suppressor genes; some of these are proteins with PDZ domains involved in homeostasis and cell polarity. Through a bioinformatics approach based on interaction networks, a group of proteins associated with HPV 16 infection, PDZ domains, and direct physical interaction with E6 and related to different hallmarks of cancer were identified. MAGI-1 was selected to evaluate the expression profile and subcellular localization changes in premalignant lesions and ISCC with HPV 16 in an integrated state in cervical cytology; the profile expression of MAGI-1 diminished according to lesion grade. Surprisingly, in cell lines CaSki and SiHa, the protein localization was cytoplasmic and nuclear. In contrast, in histological samples, a change in subcellular localization from the cytoplasm in low-grade squamous intraepithelial lesions (LSIL) to the nucleus in the high-grade squamous intraepithelial lesion (HSIL) was observed; in in situ carcinomas and ISCC, MAGI-1 expression was absent. In conclusion, MAGI-1 expression could be a potential biomarker for distinguishing those cells with normal morphology but with HPV 16 integrated from those showing morphology-related uterine cervical lesions associated with tumor progression.

Targeting the inward rectifier potassium channel 5.1 in thyroid cancer: artificial intelligence-facilitated molecular docking for drug discovery

BACKGROUND

Recurrent and metastatic thyroid cancer is more invasive and can transform to dedifferentiated thyroid cancer, thus leading to a severe decline in the 10-year survival. The thyroid-stimulating hormone receptor (TSHR) plays an important role in differentiation process. We aim to find a therapeutic target in redifferentiation strategies for thyroid cancer.

METHODS

Our study integrated the differentially expressed genes acquired from the Gene Expression Omnibus database by comparing TSHR expression levels in the Cancer Genome Atlas database. We conducted functional enrichment analysis and verified the expression of these genes by RT-PCR in 68 pairs of thyroid tumor and paratumor tissues. Artificial intelligence-enabled virtual screening was combined with the VirtualFlow platform for deep docking.

RESULTS

We identified five genes (KCNJ16, SLC26A4, TG, TPO, and SYT1) as potential cancer treatment targets. TSHR and KCNJ16 were downregulated in the thyroid tumor tissues, compared with paired normal tissues. In addition, KCNJ16 was lower in the vascular/capsular invasion group. Enrichment analyses revealed that KCNJ16 may play a significant role in cell growth and differentiation. The inward rectifier potassium channel 5.1 (Kir5.1, encoded by KCNJ16) emerged as an interesting target in thyroid cancer. Artificial intelligence-facilitated molecular docking identified Z2087256678_2, Z2211139111_1, Z2211139111_2, and PV-000592319198_1 (-7.3 kcal/mol) as the most potent commercially available molecular targeting Kir5.1.

CONCLUSION

This study may provide greater insights into the differentiation features associated with TSHR expression in thyroid cancer, and Kir5.1 may be a potential therapeutic target in the redifferentiation strategies for recurrent and metastatic thyroid cancer.

Genomic and epigenomic profile of thyroid cancer

Thyroid cancer is the most common malignancy of the endocrine system, and its incidence has been steadily increasing. Advances in sequencing have allowed analysis of the entire cancer genome, and has provided new information on the genetic lesions and modifications responsible for the onset, progression, dedifferentiation and metastasis of thyroid carcinomas. Moreover, integrated genomics has advanced our understanding of the development of cancer and its behavior, and has facilitated the identification of new genetic mutations and molecular pathways. The functional analysis of epigenetic modifications, such as DNA methylation, histone acetylation and non-coding RNAs, have contributed to define new regulatory mechanisms that control cell malignancy in thyroid cancer, especially aggressive forms. Here we review the most recent advances in genomics and epigenomics of thyroid cancer, which have resulted in a new classification and interpretation of the initiation and progression of thyroid tumors, providing new tools and opportunities for further investigation and for the clinical development of new treatment strategies.

Analyzing integrated network of methylation and gene expression profiles in lung squamous cell carcinoma

Gene expression, DNA methylation, and their organizational relationships are commonly altered in lung squamous cell carcinoma (LUSC). To elucidate these complex interactions, we reconstructed a differentially expressed gene network and a differentially methylated cytosine (DMC) network by partial information decomposition and an inverse correlation algorithm, respectively. Then, we performed graph union to integrate the networks. Community detection and enrichment analysis of the integrated network revealed close interactions between the cell cycle, keratinization, immune system, and xenobiotic metabolism gene sets in LUSC. DMC analysis showed that hypomethylation targeted the gene sets responsible for cell cycle, keratinization, and NRF2 pathways. On the other hand, hypermethylated genes affected circulatory system development, the immune system, extracellular matrix organization, and cilium organization. By centrality measurement, we identified NCAPG2, PSMG3, and FADD as hub genes that were highly connected to other nodes and might play important roles in LUSC gene dysregulation. We also found that the genes with high betweenness centrality are more likely to affect patients’ survival than those with low betweenness centrality. These results showed that the integrated network analysis enabled us to obtain a global view of the interactions and regulations in LUSC.

The magic of MAGI-1: A scaffolding protein with multi signalosomes and functional plasticity

MAGI-1 is a critical cellular scaffolding protein with over 110 different cellular and microbial protein interactors. Since the discovery of MAGI-1 in 1997, MAGI-1 has been implicated in diverse cellular functions such as polarity, cell-cell communication, neurological processes, kidney function, and a host of diseases including cancer and microbial infection. Additionally, MAGI-1 has undergone nomenclature changes in response to the discovery of an additional PDZ domain, leading to lack of continuity in the literature. We address the nomenclature of MAGI-1 as well as summarize many of the critical functions of the known interactions. Given the importance of many of the interactors, such as human papillomavirus E6, the Coxsackievirus and adenovirus receptor (CAR), and PTEN, the enhancement or disruption of MAGI-based interactions has the potential to affect cellular functions that can potentially be harnessed as a therapeutic strategy for a variety of diseases.

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.

Characterization of the CpG island methylator phenotype subclass in papillary thyroid carcinoma

CpG island methylator phenotype (CIMP), characterized by the concurrent and widespread hypermethylation of a cluster of CpGs, has been reported to play an important role in carcinogenesis. Limited studies have explored the role of CIMP in papillary thyroid carcinomas (PTCs). Here, in genome-wide DNA methylation analysis of 350 primary PTCs from the Cancer Genome Atlas database that were assessed using the Illumina HumanMethylation450K platform, our study helps to identify two subtypes displayed markedly distinct DNA methylation levels, termed CIMP (high levels of DNA methylation) and nCIMP subgroup (low levels of DNA methylation). Interestingly, PTCs with CIMP tend to have a higher degree of malignancy, since this subtype was tightly associated with older age, advanced pathological stage, and lymph node metastasis (all P < 0.05). Differential methylation analysis showed a broad methylation gain in CIMP and subsequent generalized gene set testing analysis based on the significantly methylated probes in CIMP showed remarkable enrichment in epithelial mesenchymal transition and angiogenesis hallmark pathways, confirming that the CIMP phenotype may promote the tumor progression from another perspective. Analysis of tumor microenvironment showed that CIMP PTCs are in an immune-depletion status, which may affect the effectiveness of immunotherapy. Genetically, the significantly higher tumor mutation burden and copy number alteration both at the genome and focal level confirmed the genomic heterogeneity and chromosomal instability of CIMP. tumor Corresponding to the above findings, PTC patients with CIMP showed remarkable poor clinical outcome as compared to nCIMP regarding overall survival and progression-free survival. More importantly, CIMP was associated with worse survival independent of known prognostic factors.

Pan-Genomic Sequencing Reveals Actionable CDKN2A/2B Deletions and Kataegis in Anaplastic Thyroid Carcinoma

Anaplastic thyroid carcinoma (ATC) is a lethal malignancy characterized by poor response to conventional therapies. Whole-genome sequencing (WGS) analyses of this tumor type are limited, and we therefore interrogated eight ATCs using WGS and RNA sequencing. Five out of eight cases (63%) displayed cyclin-dependent kinase inhibitor 2A (CDKN2A) abnormalities, either copy number loss (n = 4) or truncating mutations (n = 1). All four cases with loss of the CDKN2A locus (encoding p16 and p14arf) also exhibited loss of the neighboring CDKN2B gene (encoding p15ink4b), and displayed reduced CDKN2A/2B mRNA levels. Mutations in established ATC-related genes were observed, including TP53, BRAF, ARID1A, and RB1, and overrepresentation of mutations were also noted in 13 additional cancer genes. One of the more predominant mutational signatures was intimately coupled to the activity of Apolipoprotein B mRNA-editing enzyme, the catalytic polypeptide-like (APOBEC) family of cytidine deaminases implied in kataegis, a focal hypermutation phenotype, which was observed in 4/8 (50%) cases. We corroborate the roles of CDKN2A/2B in ATC development and identify kataegis as a recurrent phenomenon. Our findings pinpoint clinically relevant alterations, which may indicate response to CDK inhibitors, and focal hypermutational phenotypes that may be coupled to improved responses using immune checkpoint inhibitors.

Alterations in SLC4A2, SLC26A7 and SLC26A9 Drive Acid-Base Imbalance in Gastric Neuroendocrine Tumors and Uncover a Novel Mechanism for a Co-Occurring Polyautoimmune Scenario

Autoimmune polyendocrine syndrome (APS) is assumed to involve an immune system malfunction and entails several autoimmune diseases co-occurring in different tissues of the same patient; however, they are orphans of its accurate diagnosis, as its genetic basis and pathogenic mechanism are not understood. Our previous studies uncovered alterations in the ATPase H+/K+ Transporting Subunit Alpha (ATP4A) proton pump that triggered an internal cell acid-base imbalance, offering an autoimmune scenario for atrophic gastritis and gastric neuroendocrine tumors with secondary autoimmune pathologies. Here, we propose the genetic exploration of APS involving gastric disease to understand the underlying pathogenic mechanism of the polyautoimmune scenario. The whole exome sequencing (WES) study of five autoimmune thyrogastric families uncovered different pathogenic variants in SLC4A2, SLC26A7 and SLC26A9, which cotransport together with ATP4A. Exploratory in vitro studies suggested that the uncovered genes were involved in a pathogenic mechanism based on the alteration of the acid-base balance. Thus, we built a custom gene panel with 12 genes based on the suggested mechanism to evaluate a new series of 69 APS patients. In total, 64 filtered putatively damaging variants in the 12 genes of the panel were found in 54.17% of the studied patients and none of the healthy controls. Our studies reveal a constellation of solute carriers that co-express in the tissues affected with different autoimmune diseases, proposing a unique genetic origin for co-occurring pathologies. These results settle a new-fangled genetics-based mechanism for polyautoimmunity that explains not only gastric disease, but also thyrogastric pathology and disease co-occurrence in APS that are different from clinical incidental findings. This opens a new window leading to the prediction and diagnosis of co-occurring autoimmune diseases and clinical management of patients.

Identification of 4-methylation driven genes based prognostic signature in thyroid cancer: an integrative analysis based on the methylmix algorithm

Thyroid cancer (TC) is known with a high rate of persistence and recurrence. We aimed to develop a prognostic signature to monitor and assess the survival of TC patients. mRNA expression and methylation data were downloaded from the TCGA database. Then, R package methylmix was applied to construct a mixed model was used to identify methylation-driven genes (MDGs) according to the methylation levels. Furthermore, an MDGs based prognostic signature and predictive nomogram were constructed according to the analysis of univariate and multivariate Cox regression. Totally 62 methylation-driven genes that were mainly enriched in substrate-dependent cell migration, cellular response to mechanical stimulus, et al. were found in TC tissues. aldolase C (AldoC), C14orf62, dishevelled 1 (DVL1), and protein tyrosine phosphatase receptor type C (PTPRC) were identified to be significantly related to patients' survival, and may serve as independent prognostic biomarkers for TC. Additionally, the prognostic methylation signature and a novel prognostic, predictive nomogram was established based on the methylation level of 4 MDGs. In this study, we developed a 4-MDGs based prognostic model, which might be the potential predictors for the survival rate of TC patients, and this findings might provide a novel sight for accurate monitoring and prognosis assessment.

A New Story of the Three Magi: Scaffolding Proteins and lncRNA Suppressors of Cancer

Scaffolding molecules exert a critical role in orchestrating cellular response through the spatiotemporal assembly of effector proteins as signalosomes. By increasing the efficiency and selectivity of intracellular signaling, these molecules can exert (anti/pro)oncogenic activities. As an archetype of scaffolding proteins with tumor suppressor property, the present review focuses on MAGI1, 2, and 3 (membrane-associated guanylate kinase inverted), a subgroup of the MAGUK protein family, that mediate networks involving receptors, junctional complexes, signaling molecules, and the cytoskeleton. MAGI1, 2, and 3 are comprised of 6 PDZ domains, 2 WW domains, and 1 GUK domain. These 9 protein binding modules allow selective interactions with a wide range of effectors, including the PTEN tumor suppressor, the β-catenin and YAP1 proto-oncogenes, and the regulation of the PI3K/AKT, the Wnt, and the Hippo signaling pathways. The frequent downmodulation of MAGIs in various human malignancies makes these scaffolding molecules and their ligands putative therapeutic targets. Interestingly, MAGI1 and MAGI2 genetic loci generate a series of long non-coding RNAs that act as a tumor promoter or suppressor in a tissue-dependent manner, by selectively sponging some miRNAs or by regulating epigenetic processes. Here, we discuss the different paths followed by the three MAGIs to control carcinogenesis.

MAGI1, a Scaffold Protein with Tumor Suppressive and Vascular Functions

MAGI1 is a cytoplasmic scaffolding protein initially identified as a component of cell-to-cell contacts stabilizing cadherin-mediated cell–cell adhesion in epithelial and endothelial cells. Clinical-pathological and experimental evidence indicates that MAGI1 expression is decreased in some inflammatory diseases, and also in several cancers, including hepatocellular carcinoma, colorectal, cervical, breast, brain, and gastric cancers and appears to act as a tumor suppressor, modulating the activity of oncogenic pathways such as the PI3K/AKT and the Wnt/β-catenin pathways. Genomic mutations and other mechanisms such as mechanical stress or inflammation have been described to regulate MAGI1 expression. Intriguingly, in breast and colorectal cancers, MAGI1 expression is induced by non-steroidal anti-inflammatory drugs (NSAIDs), suggesting a role in mediating the tumor suppressive activity of NSAIDs. More recently, MAGI1 was found to localize at mature focal adhesion and to regulate integrin-mediated adhesion and signaling in endothelial cells. Here, we review MAGI1′s role as scaffolding protein, recent developments in the understanding of MAGI1 function as tumor suppressor gene, its role in endothelial cells and its implication in cancer and vascular biology. We also discuss outstanding questions about its regulation and potential translational implications in oncology.

Computational methods for the prediction of chromatin interaction and organization using sequence and epigenomic profiles

The exploration of three-dimensional chromatin interaction and organization provides insight into mechanisms underlying gene regulation, cell differentiation and disease development. Advances in chromosome conformation capture technologies, such as high-throughput chromosome conformation capture (Hi-C) and chromatin interaction analysis by paired-end tag (ChIA-PET), have enabled the exploration of chromatin interaction and organization. However, high-resolution Hi-C and ChIA-PET data are only available for a limited number of cell lines, and their acquisition is costly, time consuming, laborious and affected by theoretical limitations. Increasing evidence shows that DNA sequence and epigenomic features are informative predictors of regulatory interaction and chromatin architecture. Based on these features, numerous computational methods have been developed for the prediction of chromatin interaction and organization, whereas they are not extensively applied in biomedical study. A systematical study to summarize and evaluate such methods is still needed to facilitate their application. Here, we summarize 48 computational methods for the prediction of chromatin interaction and organization using sequence and epigenomic profiles, categorize them and compare their performance. Besides, we provide a comprehensive guideline for the selection of suitable methods to predict chromatin interaction and organization based on available data and biological question of interest.