Anaplastic Thyroid Carcinoma: A ceRNA Analysis Pointed to a Crosstalk between SOX2, TP53, and microRNA Biogenesis

Risk stratification for radioactive iodine refractoriness using molecular alterations in distant metastatic differentiated thyroid cancer

Objective

Patients with radioactive iodine-refractory differentiated thyroid cancer (RAIR-DTC) are often diagnosed with delay and constrained to limited treatment options. The correlation between RAI refractoriness and the underlying genetic characteristics has not been extensively studied.

Methods

Adult patients with distant metastatic DTC were enrolled and assigned to undergo next-generation sequencing of a customized 26-gene panel (ThyroLead). Patients were classified into RAIR-DTC or non-RAIR groups to determine the differences in clinicopathological and molecular characteristics. Molecular risk stratification (MRS) was constructed based on the association between molecular alterations identified and RAI refractoriness, and the results were classified as high, intermediate or low MRS.

Results

A total of 220 patients with distant metastases were included, 63.2% of whom were identified as RAIR-DTC. Genetic alterations were identified in 90% of all the patients, with BRAF (59.7% vs. 17.3%), TERT promoter (43.9% vs. 7.4%), and TP53 mutations (11.5% vs. 3.7%) being more prevalent in the RAIR-DTC group than in the non-RAIR group, except for RET fusions (15.8% vs. 39.5%), which had the opposite pattern. BRAF and TERT promoter are independent predictors of RAIR-DTC, accounting for 67.6% of patients with RAIR-DTC. MRS was strongly associated with RAI refractoriness (P<0.001), with an odds ratio (OR) of high to low MRS of 7.52 [95% confidence interval (95% CI), 3.96-14.28; P<0.001] and an OR of intermediate to low MRS of 3.20 (95% CI, 1.01-10.14; P=0.041).

Conclusions

Molecular alterations were associated with RAI refractoriness, with BRAF and TERT promoter mutations being the predominant contributors, followed by TP53 and DICER1 mutations. MRS might serve as a valuable tool for both prognosticating clinical outcomes and directing precision-based therapeutic interventions.

MBS: a genome browser annotation track for high-confident microRNA binding sites in whole human transcriptome

MicroRNAs (miRNAs) are small non-coding ribonucleic acids (RNAs) that play a role in many regulatory pathways in eukaryotes. They usually exert their functions by binding mature messenger RNAs. The prediction of the binding targets of the endogenous miRNAs is crucial to unravel the processes they are involved in. In this work, we performed an extensive miRNA binding sites (MBS) prediction over all the annotated transcript sequences and made them available through an UCSC track. MBS annotation track allows to study and visualize the human miRNA binding sites transcriptome-wide in a genome browser, together with any other available information the user is interested in. In the creation of the database that underlies the MBS track, three consolidated algorithms of miRNA binding prediction have been used: PITA, miRanda and TargetScan, and information about the binding sites predicted by all of them has been collected. MBS track displays high-confident miRNA binding sites for the whole length of each human transcript, both coding and non-coding ones. Each annotation can redirect to a web page with the details of the miRNA binding and the involved transcripts. MBS can be easily applied to retrieve specific information such as the effects of alternative splicing on miRNA binding or when a specific miRNA binds an exon-exon junction in the mature RNA. Overall, MBS will be of great help for studying and visualizing, in a user-friendly mode, the predicted miRNA binding sites on all the transcripts arising from a gene or a region of interest. Database URL https://datasharingada.fondazionerimed.com:8080/MBS.

An extended catalogue of ncRNAs in Streptomyces coelicolor reporting abundant tmRNA, RNase-P RNA and RNA fragments derived from pre-ribosomal RNA leader sequences

Streptomyces coelicolor is a model organism for studying streptomycetes. This genus possesses relevant medical and economical roles, because it produces many biologically active metabolites of pharmaceutical interest, including the majority of commercialized antibiotics. In this bioinformatic study, the transcriptome of S. coelicolor has been analyzed to identify novel RNA species and quantify the expression of both annotated and novel transcripts in solid and liquid growth medium cultures at different times. The major characteristics disclosed in this study are: (i) the diffuse antisense transcription; (ii) the great abundance of transfer-messenger RNAs (tmRNA); (iii) the abundance of rnpB transcripts, paramount for the RNase-P complex; and (iv) the presence of abundant fragments derived from pre-ribosomal RNA leader sequences of unknown biological function. Overall, this study extends the catalogue of ncRNAs in S. coelicolor and suggests an important role of non-coding transcription in the regulation of biologically active molecule production.

Repetitive Sequence Transcription in Breast Cancer

Repetitive sequences represent about half of the human genome. They are actively transcribed and play a role during development and in epigenetic regulation. The altered activity of repetitive sequences can lead to genomic instability and they can contribute to the establishment or the progression of degenerative diseases and cancer transformation. In this work, we analyzed the expression profiles of DNA repetitive sequences in the breast cancer specimens of the HMUCC cohort. Satellite expression is generally upregulated in breast cancers, with specific families upregulated per histotype: in HER2-enriched cancers, they are the human satellite II (HSATII), in luminal A and B, they are part of the ALR family and in triple-negative, they are part of SAR and GSAT families, together with a perturbation in the transcription from endogenous retroviruses and their LTR sequences. We report that the background expression of repetitive sequences in healthy tissues of cancer patients differs from the tissues of non-cancerous controls. To conclude, peculiar patterns of expression of repetitive sequences are reported in each specimen, especially in the case of transcripts arising from satellite repeats.

ceRNA analysis of SARS-CoV-2

Viral RNAs can perturb the miRNA regulatory network, competing with host RNAs as part of their infective process. An in silico competing endogenous RNA (ceRNA) analysis has been carried on SARS-CoV-2. The results suggest that, in humans, the decrease of microRNA activity caused by viral RNAs can lead to a perturbation of vesicle trafficking and the inflammatory response, in particular by enhancing KLF10 activity. The results suggest also that, during the study of the mechanics of viral infections, it could be of general interest to investigate the competition of viral RNA with cellular transcripts for shared microRNAs.

lncRNA CERS6-AS1 as ceRNA promote cell proliferation of breast cancer by sponging miR-125a-5p to upregulate BAP1 expression

Long noncoding RNAs (lncRNAs) can act as oncogene and tumor suppressor genes in many types of cancers including breast cancer (BC). Our previous study has indicated microRNA (miR)-125a-5p was downregulated and function as a tumor suppressor in BC. However, its upstream regulation mechanism is still unclear. In this study, we used bioinformatics algorithms, RNA pulldown assay, and dual-luciferase reports assay to predict and confirm lncRNA CERS6-AS1 interacted with miR-125a-5p. Then we found CERS6-AS1 was upregulated in BC tissues. Experimental results of tumor growth in nude mice show that CERS6-AS1 promotes tumor growth. Furthermore, CERS6-AS1 regulated BC susceptibility gene 1-associated protein 1 (BAP1) expression via sponging miR-125a-5p via Western blot analysis and quantitative polymerase chain reaction arrays. Finally, we showed that miR-125a-5p had opposing effects to those of CERS6-AS1 on BC cells, demonstrating that CERS6-AS1 may promote cell proliferation and inhibit cell apoptosis via sponging miR-125a-5p. Our results indicated CERS6-AS1 promote BC cell proliferation and inhibit cell apoptosis via sponging miR-125a-5p to upregulate BAP1 expression.

Abemaciclib (CDK4/6 Inhibitor) Blockade Induces Cytotoxicity in Human Anaplastic Thyroid Carcinoma Cells

BACKGROUND

Thyroid cancer is the most prevalent endocrine malignancies globally. Anaplastic thyroid carcinoma (ATC) accounts for 1-3% of all Thyroid cancer. The evidence showed that ATC is a highly invasive solid tumor with poor prognosis. Despite conventional chemotherapy treatments, a considerable number of patients show developing resistance to therapeutic agents and tumor relapse. The aim of this study was the investigation anti-tumor effect of Abemaciclib (novel targeted cancer therapy drug) on Anaplastic Thyroid carcinoma SW1736 and C643 cell lines.

METHODS

SW1736 and C643 cell lines were treated by desire concentrations of Abemaciclib (0, 1, 2.5, 5, 10, and 20 μM) and cell viability was measured by MTT assay. Also, Anoikis resistance assay was conducted for non-adherent the cells in the exposure of Abemaciclib. The gene expression of apoptotic and anti-apoptotic genes was conducted by quantitative Real-time PCR.

RESULTS

Abemaciclib at the concentration of 10 and 20 μM effectively reduced cell proliferation and growth of the ATC cells compared to the control (p=0.000). Furthermore, we showed that 10 and 20 μM doses of the Abemaciclib inhibited the non-adherent ATC cells which were resistant to Anoikis death significantly (p=0.001). Moreover, we demonstrated this targeted therapy significantly reduced anti-apoptotic gene expression levels (BCL2 and CMYC) (p<0.05) and increased apoptotic gene expressions such as P21 and BAX (p<0.05).

CONCLUSION

Our data suggested that Abemaciclib can be utilized as a novel therapeutic agent in ATC cancer. Further in vivo and in vitro investigations are needed to evaluate molecular and clinical mechanisms of Abemaciclib.

Construction and investigation of lncRNA-associated ceRNA regulatory network in papillary thyroid cancer

Increasing evidence has experimentally proved the competitive endogenous RNA (ceRNA) hypothesis that long non-coding RNA (lncRNA) can affect the expression of RNA targets by competitively combining microRNA (miRNA) via miRNA response elements. However, an extensive ceRNA network of thyroid carcinoma in a large cohort has not been evaluated. We analyzed the RNAseq and miRNAseq data of 348 cases of primary papillary thyroid cancer (PTC) patients with clinical information downloaded from The Cancer Genome Atlas (TCGA) project to search for potential biomarkers or therapeutic targets. A computational approach was applied to build an lncRNA-miRNA-mRNA regulatory network of PTC. In total, 780 lncRNAs were detected as collectively dysregulated lncRNAs in all 3 PTC variants compared with normal tissues (fold change >2 and false discovery rate <0.05). The interactions among 45 lncRNAs, 13 miRNAs and 86 mRNAs constituted a ceRNA network of PTC. Nine out of the 45 aberrantly expressed lncRNAs were related to the clinical features of PTC patients. However, the expression levels of 3 lncRNAs (LINC00284, RBMS3-AS1 and ZFX-AS1) were identified to be tightly correlated with the patients overall survival (log-rank, P<0.05). The present study identified a list of specific lncRNAs associated with PTC progression and prognosis. This complex ceRNA interaction network in PTC may provide guidance for better understanding the molecular mechanisms underlying PTC.

Network analysis of pseudogene-gene relationships: from pseudogene evolution to their functional potentials

Pseudogenes are fossil relatives of genes. Pseudogenes have long been thought of as "junk DNAs", since they do not code proteins in normal tissues. Although most of the human pseudogenes do not have noticeable functions, ∼20% of them exhibit transcriptional activity. There has been evidence showing that some pseudogenes adopted functions as lncRNAs and work as regulators of gene expression. Furthermore, pseudogenes can even be "reactivated" in some conditions, such as cancer initiation. Some pseudogenes are transcribed in specific cancer types, and some are even translated into proteins as observed in several cancer cell lines. All the above have shown that pseudogenes could have functional roles or potentials in the genome. Evaluating the relationships between pseudogenes and their gene counterparts could help us reveal the evolutionary path of pseudogenes and associate pseudogenes with functional potentials. It also provides an insight into the regulatory networks involving pseudogenes with transcriptional and even translational activities.In this study, we develop a novel approach integrating graph analysis, sequence alignment and functional analysis to evaluate pseudogene-gene relationships, and apply it to human gene homologs and pseudogenes. We generated a comprehensive set of 445 pseudogene-gene (PGG) families from the original 3,281 gene families (13.56%). Of these 438 (98.4% PGG, 13.3% total) were non-trivial (containing more than one pseudogene). Each PGG family contains multiple genes and pseudogenes with high sequence similarity. For each family, we generate a sequence alignment network and phylogenetic trees recapitulating the evolutionary paths. We find evidence supporting the evolution history of olfactory family (both genes and pseudogenes) in human, which also supports the validity of our analysis method. Next, we evaluate these networks in respect to the gene ontology from which we identify functions enriched in these pseudogene-gene families and infer functional impact of pseudogenes involved in the networks. This demonstrates the application of our PGG network database in the study of pseudogene function in disease context.

Identification of competing endogenous RNAs of the tumor suppressor gene PTEN: A probabilistic approach

Regulation by microRNAs (miRNAs) and modulation of miRNA activity are critical components of diverse cellular processes. Recent research has shown that miRNA-based regulation of the tumor suppressor gene PTEN can be modulated by the expression of other miRNA targets acting as competing endogenous RNAs (ceRNAs). However, the key sequence-based features enabling a transcript to act as an effective ceRNA are not well understood and a quantitative model associating statistical significance to such features is currently lacking. To identify and assess features characterizing target recognition by PTEN-regulating miRNAs, we analyze multiple datasets from PAR-CLIP experiments in conjunction with RNA-Seq data. We consider a set of miRNAs known to regulate PTEN and identify high-confidence binding sites for these miRNAs on the 3′ UTR of protein coding genes. Based on the number and spatial distribution of these binding sites, we calculate a set of probabilistic features that are used to make predictions for novel ceRNAs of PTEN. Using a series of experiments in human prostate cancer cell lines, we validate the highest ranking prediction (TNRC6B) as a ceRNA of PTEN. The approach developed can be applied to map ceRNA networks of critical cellular regulators and to develop novel insights into crosstalk between different pathways involved in cancer.

Expression of MicroRNAs in Thyroid Carcinoma

MicroRNA (miRNA) are negative regulators of gene expression and subsequent protein production. This method of action translates into regulatory control over cellular processes, including development, signaling, metabolism, and apoptosis. A broad range of miRNA are shown to have abnormal expressions in thyroid cancers which could explain the pathology of tumor oncogenesis and disease progression. A review is conducted of the current research on miRNA dysregulation in thyroid cancers, including papillary thyroid carcinoma (PTC), follicular thyroid carcinoma (FTC), anaplastic thyroid cancer (ATC), and medullary thyroid carcinoma (MTC). Dysregulated miRNA and their associated regulatory pathways are identified and their oncogenic and pathological significance are discussed.

Extensive ceRNA-ceRNA interaction networks mediated by miRNAs regulate development in multiple rhesus tissues

Crosstalk between RNAs mediated by shared microRNAs (miRNAs) represents a novel layer of gene regulation, which plays important roles in development. In this study, we analyzed time series expression data for coding genes and long non-coding RNAs (lncRNAs) to identify thousands of interactions among competitive endogenous RNAs (ceRNAs) in four rhesus tissues. The ceRNAs exhibited dynamic expression and regulatory patterns during each tissue development process, which suggests that ceRNAs might work synergistically during different developmental stages or tissues to control specific functions. In addition, lncRNAs exhibit higher specificity as ceRNAs than coding-genes and their functions were predicted based on their competitive coding-gene partners to discover their important developmental roles. In addition to the specificity of tissue development, functional analyses demonstrated that the combined effects of multiple ceRNAs can have major impacts on general developmental and metabolic processes in multiple tissues, especially transcription-related functions where competitive interactions. Moreover, ceRNA interactions could sequentially and/or synergistically mediate the crosstalk among different signaling pathways during brain development. Analyzing ceRNA interactions during the development of multiple tissues will provideinsights in the regulation of normal development and the dysregulation of key mechanisms during pathogenesis.

A ceRNA approach may unveil unexpected contributors to deletion syndromes, the model of 5q- syndrome

In genomic deletions, gene haploinsufficiency might directly configure a specific disease phenotype. Nevertheless, in some cases no functional association can be identified between haploinsufficient genes and the deletion-associated phenotype. Transcripts can act as microRNA sponges. The reduction of transcripts from the hemizygous region may increase the availability of specific microRNAs, which in turn may exert in-trans regulation of target genes outside the deleted region, eventually contributing to the phenotype. Here we prospect a competing endogenous RNA (ceRNA) approach for the identification of candidate genes target of epigenetic regulation in deletion syndromes. As a model, we analyzed the 5q- myelodysplastic syndrome. Genes in haploinsufficiency within the common 5q deleted region in CD34+ blasts were identified in silico. Using the miRWalk 2.0 platform, we predicted microRNAs whose availability, and thus activity, could be enhanced by the deletion, and performed a genomewide analysis of the genes outside the 5q deleted region that could be targeted by the predicted miRNAs. The analysis pointed to two genes with altered expression in 5q- transcriptome, which have never been related with 5q- before. The prospected approach allows investigating the global transcriptional effect of genomic deletions, possibly prompting discovery of unsuspected contributors in the deletion-associated phenotype. Moreover, it may help in functionally characterizing previously reported unexpected interactions.