Charny Park

lncRNAtor: a comprehensive resource for functional investigation of long non-coding RNAs

MOTIVATION A number of long non-coding RNAs (lncRNAs) have been identified by deep sequencing methods, but their molecular and cellular functions are known only for a limited number of lncRNAs. Current databases on lncRNAs are mostly for cataloging purpose without providing in-depth information required to infer functions. A comprehensive resource on lncRNA function is an immediate need. RESULTS We present a database for functional investigation of lncRNAs that encompasses annotation, sequence analysis, gene expression, protein binding and phylogenetic conservation. We have compiled lncRNAs for six species (human, mouse, zebrafish, fruit fly, worm and yeast) from ENSEMBL, HGNC, MGI and lncRNAdb. Each lncRNA was analyzed for coding potential and phylogenetic conservation in different lineages. Gene expression data of 208 RNA-Seq studies (4995 samples), collected from GEO, ENCODE, modENCODE and TCGA databases, were used to provide expression profiles in various tissues, diseases and developmental stages. Importantly, we analyzed RNA-Seq data to identify coexpressed mRNAs that would provide ample insights on lncRNA functions. The resulting gene list can be subject to enrichment analysis such as Gene Ontology or KEGG pathways. Furthermore, we compiled protein-lncRNA interactions by collecting and analyzing publicly available CLIP-seq or PAR-CLIP sequencing data. Finally, we explored evolutionarily conserved lncRNAs with correlated expression between human and six other organisms to identify functional lncRNAs. The whole contents are provided in a user-friendly web interface. AVAILABILITY AND IMPLEMENTATION lncRNAtor is available at http://lncrnator.ewha.ac.kr/. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

miRGator v2.0 : an integrated system for functional investigation of microRNAs

miRGator is an integrated database of microRNA (miRNA)-associated gene expression, target prediction, disease association and genomic annotation, which aims to facilitate functional investigation of miRNAs. The recent version of miRGator v2.0 contains information about (i) human miRNA expression profiles under various experimental conditions, (ii) paired expression profiles of both mRNAs and miRNAs, (iii) gene expression profiles under miRNA-perturbation (e.g. miRNA knockout and overexpression), (iv) known/predicted miRNA targets and (v) miRNA-disease associations. In total, >8000 miRNA expression profiles, ∼300 miRNA-perturbed gene expression profiles and ∼2000 mRNA expression profiles are compiled with manually curated annotations on disease, tissue type and perturbation. By integrating these data sets, a series of novel associations (miRNA–miRNA, miRNA–disease and miRNA–target) is extracted via shared features. For example, differentially expressed genes (DEGs) after miRNA knockout were systematically compared against miRNA targets. Likewise, differentially expressed miRNAs (DEmiRs) were compared with disease-associated miRNAs. Additionally, miRNA expression and disease-phenotype profiles revealed miRNA pairs whose expression was regulated in parallel in various experimental and disease conditions. Complex associations are readily accessible using an interactive network visualization interface. The miRGator v2.0 serves as a reference database to investigate miRNA expression and function (http://miRGator.kobic.re.kr).

Deregulation of Immune Response Genes in Patients With Epstein-Barr Virus-Associated Gastric Cancer and Outcomes

BACKGROUND & AIMS Patients with Epstein-Barr virus-associated gastric carcinoma (EBVaGC) have a better prognosis than those with gastric cancer not associated with EBV infection (EBVnGC). This is partly because EBV infection recruits lymphocytes, which infiltrate the tumor. A high degree of tumor heterogeneity is likely to be associated with poor response. We investigated differences in gene expression patterns between EBVaGC and EBVnGC. METHODS We used gene expression profile analysis to compare tumor and nontumor gastric tissues from 12 patients with EBVaGC and 14 patients with EBVnGC. Findings were validated by whole transcriptome RNAseq and real-time quantitative polymerase chain reaction analyses. CD3(+) primary T cells were isolated from human blood samples; migration of these cells and of Jurkat cells were measured in culture with EBV-infected and uninfected gastric cancer cells. RESULTS Based on Pearson correlation matrix analysis, EBVaGCs had a higher degree of homogeneity than EBVnGCs. Although 4550 genes were differentially expressed between tumor and nontumor gastric tissues of patients with EBVnGC, only 186 genes were differentially expressed between tumor and nontumor gastric tissues of patients with EBVaGC (P < .001). This finding supports the concept that EBVaGCs have fewer genetic and epigenetic alterations than EBVnGCs. Expression of major histocompatibility complex class II genes and genes that regulate chemokine activity were more often deregulated in EBVaGCs compared with nontumor tissues. In culture, more T cells migrated to EBV-infected gastric cancer cells than to uninfected cells; migration was blocked with a neutralizing antibody against CXCR3 (a receptor for many chemokines). CONCLUSIONS Fewer genes are deregulated in EBVaGC than in EBVnGC. Most changes in EBVaGCs occur in immune response genes. These changes might allow EBVaGC to recruit reactive immune cells; this might contribute to the better outcomes of these patients compared with those with EBVnGC.

Programming of Plant Leaf Senescence with Temporal and Inter-Organellar Coordination of Transcriptome in Arabidopsis.

RNA-seq analysis of total and small RNAs throughout the lifespan of Arabidopsis leaves revealed that leaf senescence proceeds with tight temporal and distinctive inter-organellar coordination of transcriptomes. Plant leaves, harvesting light energy and fixing CO2, are a major source of foods on the earth. Leaves undergo developmental and physiological shifts during their lifespan, ending with senescence and death. We characterized the key regulatory features of the leaf transcriptome during aging by analyzing total- and small-RNA transcriptomes throughout the lifespan of Arabidopsis (Arabidopsis thaliana) leaves at multidimensions, including age, RNA-type, and organelle. Intriguingly, senescing leaves showed more coordinated temporal changes in transcriptomes than growing leaves, with sophisticated regulatory networks comprising transcription factors and diverse small regulatory RNAs. The chloroplast transcriptome, but not the mitochondrial transcriptome, showed major changes during leaf aging, with a strongly shared expression pattern of nuclear transcripts encoding chloroplast-targeted proteins. Thus, unlike animal aging, leaf senescence proceeds with tight temporal and distinct interorganellar coordination of various transcriptomes that would be critical for the highly regulated degeneration and nutrient recycling contributing to plant fitness and productivity.

Gastrointestinal malignancies harbor actionable MET exon 14 deletions

Recently, MET exon 14 deletion (METex14del) has been postulated to be one potential mechanism for MET protein overexpression. We screened for the presence of METex14del transcript by multiplexed fusion transcript analysis using nCounter assay followed by confirmation with quantitative reverse transcription PCR with correlation to MET protein expression by immunohistochemistry (IHC) and MET amplification by fluorescence in situ hybridization (FISH). We extracted RNAs from 230 patients enrolled onto the prospective molecular profiling clinical trial (NEXT-1) (NCT02141152) between November 2013 and August 2014. Thirteen METex14del cases were identified including 3 gastric cancer, 4 colon cancer, 5 non-small cell lung cancer, and one adenocarcinoma of unknown primary. Of these 13 METex14del cases, 11 were MET IHC 3+ and 2 were 2+. Only one out of the 13 METex14del cases was MET amplified (MET/CEP ratio > 2.0). Growths of two (gastric, colon) METex14del+ patient tumor derived cell lines were profoundly inhibited by both MET tyrosine kinase inhibitors and a monoclonal antibody targeting MET. In conclusion, METex14del is a unique molecular aberration present in gastrointestinal (GI) malignancies corresponding with overexpression of MET protein but rarely with MET amplification. Substantial growth inhibition of METex14del+ patient tumor derived cell lines by several MET targeting drugs strongly suggests METex14del is a potential actionable driver mutation in GI malignancies.

Patient-derived cell models as preclinical tools for genome-directed targeted therapy

Background In this study, we established patient-derived tumor cell (PDC) models using tissues collected from patients with metastatic cancer and assessed whether these models could be used as a tool for genome-based cancer treatment. Methods PDCs were isolated and cultured from malignant effusions including ascites and pleural fluid. Pathological examination, immunohistochemical analysis, and genomic profiling were performed to compare the histological and genomic features of primary tumors, PDCs. An exploratory gene expression profiling assay was performed to further characterize PDCs. Results From January 2012 to May 2013, 176 samples from patients with metastatic cancer were collected. PDC models were successfully established in 130 (73.6%) samples. The median time from specimen collection to passage 1 (P1) was 3 weeks (range, 0.5–4 weeks), while that from P1 to P2 was 2.5 weeks (range, 0.5–5 weeks). Sixteen paired samples of genomic alterations were highly concordant between each primary tumor and progeny PDCs, with an average variant allele frequency (VAF) correlation of 0.878. We compared genomic profiles of the primary tumor (P0), P1 cells, P2 cells, and patient-derived xenografts (PDXs) derived from P2 cells and found that three samples (P0, P1, and P2 cells) were highly correlated (0.99–1.00). Moreover, PDXs showed more than 100 variants, with correlations of only 0.6–0.8 for the other samples. Drug responses of PDCs were reflective of the clinical response to targeted agents in selected patient PDC lines. Conclusion(s) Our results provided evidence that our PDC model was a promising model for preclinical experiments and closely resembled the patient tumor genome and clinical response.

Identification of the BRAF V600E mutation in gastroenteropancreatic neuroendocrine tumors

Genomic profiles of gastroenteropancreatic neuroendocrine tumors (GEP-NETs) are still insufficiently understood, and the genetic alterations associated with drug responses have not been studied. Here, we performed whole exome sequencing of 12 GEP-NETs from patients enrolled in a nonrandomized, open-labeled, single-center phase II study for pazopanib, and integrated our results with previously published results on pancreas (n = 12) and small intestine NETs (n = 50). The mean numbers of somatic mutations in each case varied widely from 20 to 4682. Among 12 GEP-NETs, eight showed mutations of more than one cancer-related gene, including TP53, CNBD1, RB1, APC, BCOR, BRAF, CTNNB1, EGFR, EP300, ERBB3, KDM6A, KRAS, MGA, MLL3, PTEN, RASA1, SMARCB1, SPEN, TBC1D12, and VHL. TP53 was recurrently mutated in three cases, whereas CNBD1 and RB1 mutations were identified in two cases. Three GEP-NET patients with TP53 mutations demonstrated a durable response and one small intestinal grade (G) 1 NET patient with BRAF V600E mutation showed progression after pazopanib treatment. We found BRAF V600E (G1 NET from rectum and two G3 NETs from colon) and BRAF G593S (G2 NET from pancreas) missense mutations (9.1%) in an independent cohort of 44 GEP-NETs from the rectum (n = 26), colon (n = 7), pancreas (n = 4), small intestine (n = 3), stomach (n = 3) and appendix (n = 1) by Sanger sequencing. All tumor specimens were obtained before chemotherapy. In conclusion, BRAF V600E mutation is likely to result in resistance to pazopanib but may be a potentianally actionable mutation in metastatic GEP-NETs patients.

Host immune response index in gastric cancer identified by comprehensive analyses of tumor immunity

ABSTRACT Tumor infiltrating lymphocytes (TIL) in Epstein-Barr virus (EBV)-associated/microsatellite-unstable (MSI) gastric carcinomas (GC) constitute immune-active principal cellular components of tumor microenvironment and contribute to better prognosis. With the remarkable success of cancer immunotherapies, there is an urgent need for a comprehensive understanding of tumor-immune interactions in patients with GC in the context of host immune response. To identify GC subtype-specific immune response gene set, we tested differentially expressed genes for MSI and EBV+ GC subtypes in randomly selected test set (n = 278) in merged ACRG-SMC microarray and TCGA RNA sequencing data set. We identified Host ImmunE Response index (HIERÏ) consisting of 29 immune genes classifying GC patients into robust 3 groups with prognostic significance. Immune-high cluster 1 was enriched with PD-L1High/EBV+/MSI/TILHigh with the best clinical outcome while immune-low cluster 3 displayed worst outcome and exemplified with PD-L1Low/EBV-/MSS. The results were validated in the same cohort (n = 279) and independent cohort (n = 181) with RNA from formalin-fixed paraffin-embedded (FFPE) tissue. Unexpectedly, nearly half of GC in cluster 1 were EBV-/MSS and 10% of cluster 3 GC were EBV+/MSI GC patients, suggesting that in addition to EBV+/MSI GC subtypes, EBV-/MSS subtype also constitutes almost half of high immune cluster and would be a good candidate for immune checkpoint inhibitor therapy. In contrary, almost 10% of EBV+/MSI GC patients may not respond to immune checkpoint inhibitor therapy. Thus, our HIERÏ gene signature demonstrates the potential to subclassify tumor immunity levels, predict prognosis and help immunotherapeutic decisions.

Cytosolic Hsp60 orchestrates the survival and inflammatory responses of vascular smooth muscle cells in injured aortic vessels

AIMS Pro-inflammatory response of vascular smooth muscle cells (VSMCs) is triggered by endothelial damage and a causative step for thrombosis and neointimal thickening in the injured arterial vessels. Therefore, we investigate a role of cytosolic Hsp60 as a novel pro-inflammatory mediator in VSMCs. METHODS AND RESULTS Hsp60 was detected in the cytosol of VSMCs. The selective depletion of cytosolic Hsp60 in VSMCs reduced the IκB kinase activation, repressed the induction of nuclear factor (NF)-κB-dependent survival genes (MnSOD and Bfl-1/A1), and enhanced apoptotic death in response to TNF-α. Moreover, a quantitative RNA sequencing revealed that the expression of 75 genes among the 774 TNF-α-inducible genes was significantly reduced by the depletion of cytosolic Hsp60. In particular, the expression of pro-inflammatory cytokines/chemokines, such as CCL2, CCL20, and IL-6, was regulated by the cytosolic Hsp60 in VSMCs. Finally, the depletion of cytosolic Hsp60 markedly inhibited the neointimal thickening in the balloon-injured arterial vessels by inducing apoptotic cell death and inhibiting chemokine production. CONCLUSIONS This study provides the first evidence that cytosolic Hsp60 could be a therapeutic target for preventing VSMC hyperplasia and inflammatory response in the injured vessels.

ASpedia: a comprehensive encyclopedia of human alternative splicing

Abstract Alternative splicing confers the human genome complexity by increasing the diversity of expressed mRNAs. Hundreds or thousands of splicing regions have been identified through differential alternative splicing analysis of high-throughput datasets. However, it is hard to explain the functional impact of each splicing event. Protein domain formation and nonsense-mediated decay are considered the main functional features of splicing. However, other functional features such as miRNA target sites, phosphorylation sites and single-nucleotide variations are directly affected by alternative splicing and affect downstream function. Hence, we established ASpedia: a comprehensive database for human alternative splicing annotation, which encompasses a range of functions, from genomic annotation to isoform-specific function (ASpedia, http://combio.snu.ac.kr/aspedia). The database provides three features: (i) genomic annotation extracted from DNA, RNA and proteins; (ii) transcription and regulation elements analyzed from next-generation sequencing datasets; and (iii) isoform-specific functions collected from known and published datasets. The ASpedia web application includes three components: an annotation database, a retrieval system and a browser specialized in the identification of human alternative splicing events. The retrieval system supports multiple AS event searches resulting from high-throughput analysis and the AS browser comprises genome tracks. Thus, ASpedia facilitates the systemic annotation of the functional impacts of multiple AS events.