Transcriptomic Analysis of Papillary Thyroid Cancer: A Focus on Immune-Subtyping, Oncogenic Fusion, and Recurrence

Objectives@#. Thyroid cancer is the most common endocrine tumor, with rapidly increasing incidence worldwide. However, its transcriptomic characteristics associated with immunological signatures, driver fusions, and recurrence markers remain unclear. We aimed to investigate the transcriptomic characteristics of advanced papillary thyroid cancer. @*Methods@#. This study included 282 papillary thyroid cancer tumor samples and 155 normal samples from Chungnam National University Hospital and Seoul National University Hospital. Transcriptomic quantification was determined by high-throughput RNA sequencing. We investigated the associations of clinical parameters and molecular signatures using RNA sequencing. We validated predictive biomarkers using the Cancer Genome Atlas database. @*Results@#. Through a comparison of differentially expressed genes, gene sets, and pathways in papillary thyroid cancer compared to normal tumor-adjacent tissue, we found increased immune signaling associated with cytokines or T cells and decreased thyroid hormone synthetic pathways. In addition, patients with recurrence presented increased CD8+ T-cell and Th1-cell signatures. Interestingly, we found differentially overexpressed genes related to immune-escape signaling such as CTLA4, IDO1, LAG3, and PDCD1 in advanced papillary thyroid cancer with a low thyroid differentiation score. Fusion analysis showed that the PI3K and mitogen-activated protein kinase (MAPK) signaling pathways were regulated differently according to the RET fusion partner genes (CCDC6 or NCOA4). Finally, we identified HOXD9 as a novel molecular biomarker that predicts the recurrence of thyroid cancer in addition to known risk factors (tumor size, lymph node metastasis, and extrathyroidal extension). @*Conclusion@#. We identified a high association with immune-escape signaling in the immune-hot group with aggressive clinical characteristics among Korean thyroid cancer patients. Moreover, RET fusion differentially regulated PI3K and MAPK signaling depending on the partner gene of RET, and HOXD9 was found to be a recurrence marker for advanced papillary thyroid cancer.

Generation of Whole-Genome Sequencing Data for Comparing Primary and Castration-Resistant Prostate Cancer

Because castration-resistant prostate cancer (CRPC) does not respond to androgen deprivation therapy and has a very poor prognosis, it is critical to identify a prognostic indicator for predicting high-risk patients who will develop CRPC. Here, we report a dataset of whole genomes from four pairs of primary prostate cancer (PC) and CRPC samples. The analysis of the paired PC and CRPC samples in the whole-genome data showed that the average number of somatic mutations per patients was 7,927 in CRPC tissues compared with primary PC tissues (range, 1,691 to 21,705). Our whole-genome sequencing data of primary PC and CRPC may be useful for understanding the genomic changes and molecular mechanisms that occur during the progression from PC to CRPC.

Bioinformatic identification of prognostic signature defined by copy number alteration and expression of CCNE1 in non-muscle invasive bladder cancer

Non-muscle invasive bladder cancer (NMIBC) patients frequently fail to respond to treatment and experience disease progression because of their clinical and biological diversity. In this study, we identify a prognostic molecular signature for predicting the heterogeneity of NMIBC by using an integrative analysis of copy number and gene expression data. We analyzed the copy number and gene expression profiles of 404 patients with bladder cancer obtained from The Cancer Genome Atlas (TCGA) consortium. Of the 14 molecules with significant copy number alterations that were previously reported, 13 were significantly correlated with copy number and expression changes. Prognostic gene sets based on the 13 genes were developed, and their prognostic values were verified in three independent patient cohorts (n=501). Among them, a signature of CCNE1 and its coexpressed genes was significantly associated with disease progression and validated in the independent cohorts. The CCNE1 signature was an independent risk factor based on the result of a multivariate analysis (hazard ratio=6.849, 95% confidence interval=1.613–29.092, P=0.009). Finally, gene network and upstream regulator analyses revealed that NMIBC progression is potentially mediated by CCND1-CCNE1-SP1 pathways. The prognostic molecular signature defined by copy number and expression changes of CCNE1 suggests a novel diagnostic tool for predicting the likelihood of NMIBC progression.

Use of Graph Database for the Integration of Heterogeneous Biological Data

Understanding complex relationships among heterogeneous biological data is one of the fundamental goals in biology. In most cases, diverse biological data are stored in relational databases, such as MySQL and Oracle, which store data in multiple tables and then infer relationships by multiple-join statements. Recently, a new type of database, called the graph-based database, was developed to natively represent various kinds of complex relationships, and it is widely used among computer science communities and IT industries. Here, we demonstrate the feasibility of using a graph-based database for complex biological relationships by comparing the performance between MySQL and Neo4j, one of the most widely used graph databases. We collected various biological data (protein-protein interaction, drug-target, gene-disease, etc.) from several existing sources, removed duplicate and redundant data, and finally constructed a graph database containing 114,550 nodes and 82,674,321 relationships. When we tested the query execution performance of MySQL versus Neo4j, we found that Neo4j outperformed MySQL in all cases. While Neo4j exhibited a very fast response for various queries, MySQL exhibited latent or unfinished responses for complex queries with multiple-join statements. These results show that using graph-based databases, such as Neo4j, is an efficient way to store complex biological relationships. Moreover, querying a graph database in diverse ways has the potential to reveal novel relationships among heterogeneous biological data.

Erratum: Increased Expression of Herpes Virus-Encoded hsv1-miR-H18 and hsv2-miR-H9-5p in Cancer-Containing Prostate Tissue Compared to That in Benign Prostate Hyperplasia Tissue / 대한배뇨장애요실금학회지

In this article, a part of fund and grant supports was omitted unintentionally.

Increased Expression of Herpes Virus-Encoded hsv1-miR-H18 and hsv2-miR-H9-5p in Cancer-Containing Prostate Tissue Compared to That in Benign Prostate Hyperplasia Tissue / 대한배뇨장애요실금학회지

PURPOSE: Previously, we reported the presence of virus-encoded microRNAs (miRNAs) in the urine of prostate cancer (CaP) patients. In this study, we investigated the expression of two herpes virus-encoded miRNAs in prostate tissue. METHODS: A total of 175 tissue samples from noncancerous benign prostatic hyperplasia (BPH), 248 tissue samples from patients with CaP and BPH, and 50 samples from noncancerous surrounding tissues from these same patients were analyzed for the expression of two herpes virus-encoded miRNAs by real-time polymerase chain reaction (PCR) and immunocytochemistry using nanoparticles as molecular beacons. RESULTS: Real-time reverse transcription-PCR results revealed significantly higher expression of hsv1-miR-H18 and hsv2-miRH9- 5p in surrounding noncancerous and CaP tissues than that in BPH tissue (each comparison, P<0.001). Of note, these miRNA were expressed equivalently in the CaP tissues and surrounding noncancerous tissues. Moreover, immunocytochemistry clearly demonstrated a significant enrichment of both hsv1-miR-H18 and hsv2-miR-H9 beacon-labeled cells in CaP and surrounding noncancerous tissue compared to that in BPH tissue (each comparison, P<0.05 for hsv1-miR-H18 and hsv2- miR-H9). CONCLUSIONS: These results suggest that increased expression of hsv1-miR-H18 and hsv2-miR-H95p might be associated with tumorigenesis in the prostate. Further studies will be required to elucidate the role of these miRNAs with respect to CaP and herpes viral infections.

Next generation sequencing and urologic cancer

No abstract available.

Urinary MicroRNAs of Prostate Cancer: Virus-Encoded hsv1-miRH18 and hsv2-miR-H9-5p Could Be Valuable Diagnostic Markers / 대한배뇨장애요실금학회지

PURPOSE: MicroRNAs (miRNAs) in biological fluids are potential biomarkers for the diagnosis and assessment of urological diseases such as benign prostatic hyperplasia (BPH) and prostate cancer (PCa). The aim of the study was to identify and validate urinary cell-free miRNAs that can segregate patients with PCa from those with BPH. METHODS: In total, 1,052 urine, 150 serum, and 150 prostate tissue samples from patients with PCa or BPH were used in the study. A urine-based miRNA microarray analysis suggested the presence of differentially expressed urinary miRNAs in patients with PCa, and these were further validated in three independent PCa cohorts, using a quantitative reverse transcriptionpolymerase chain reaction analysis. RESULTS: The expression levels of hsa-miR-615-3p, hsv1-miR-H18, hsv2-miR-H9-5p, and hsa-miR-4316 were significantly higher in urine samples of patients with PCa than in those of BPH controls. In particular, herpes simplex virus (hsv)-derived hsv1-miR-H18 and hsv2-miR-H9-5p showed better diagnostic performance than did the serum prostate-specific antigen (PSA) test for patients in the PSA gray zone. Furthermore, a combination of urinary hsv2-miR-H9-5p with serum PSA showed high sensitivity and specificity, providing a potential clinical benefit by reducing unnecessary biopsies. CONCLUSIONS: Our findings showed that hsv-encoded hsv1-miR-H18 and hsv2-miR-H9-5p are significantly associated with PCa and can facilitate early diagnosis of PCa for patients within the serum PSA gray zone.

A Database of Gene Expression Profiles of Korean Cancer Genome

Because there are clear molecular differences entailing different treatment effectiveness between Korean and non-Korean cancer patients, identifying distinct molecular characteristics of Korean cancers is profoundly important. Here, we report a web-based data repository, namely Korean Cancer Genome Database (KCGD), for searching gene signatures associated with Korean cancer patients. Currently, a total of 1,403 cancer genomics data were collected, processed and stored in our repository, an ever-growing database. We incorporated most widely used statistical survival analysis methods including the Cox proportional hazard model, log-rank test and Kaplan-Meier plot to provide instant significance estimation for searched molecules. As an initial repository with the aim of Korean-specific marker detection, KCGD would be a promising web application for users without bioinformatics expertise to identify significant factors associated with cancer in Korean.

A known expressed sequence tag, BM742401, is a potent lincRNA inhibiting cancer metastasis

Long intergenic non-coding RNAs (lincRNAs) have historically been ignored in cancer biology. However, thousands of lincRNAs have been identified in mammals using recently developed genomic tools, including microarray and high-throughput RNA sequencing (RNA-seq). Several of the lincRNAs identified have been well characterized for their functions in carcinogenesis. Here we performed RNA-seq experiments comparing gastric cancer with normal tissues to find differentially expressed transcripts in intergenic regions. By analyzing our own RNA-seq and public microarray data, we identified 31 transcripts, including a known expressed sequence tag, BM742401. BM742401 was downregulated in cancer, and its downregulation was associated with poor survival in gastric cancer patients. Ectopic overexpression of BM742401 inhibited metastasis-related phenotypes and decreased the concentration of extracellular MMP9. These results suggest that BM742401 is a potential lincRNA marker and therapeutic target.

Patome: Database of Patented Bio-sequences

We have built a database server called Patome which contains the annotation information for patented bio-sequences from the Korean Intellectual Property Office (KIPO). The aims of the Patome are to annotate Korean patent bio-sequences and to provide information on patent relationship of public database entries. The patent sequences were annotated with Reference Sequence (RefSeq) or NCBI's nr database. The raw patent data and the annotated data were stored in the database. Annotation information can be used to determine whether a particular RefSeq ID or NCBI's nr ID is related to Korean patent. Patome infrastructure consists of three components-the database itself, a sequence data loader, and an online database query interface. The database can be queried using submission number, organism, title, applicant name, or accession number. Patome can be accessed at http://www.patome.net. The information will be updated every two months.

Microvasculature of the cruciate ligaments and the surrounding structures: a microangiographic study of the human knee specimen

No abstract available.