Jong Hwan Kim

Identification and evaluation of age-correlated DNA methylation markers for forensic use

In forensics, age prediction is useful to narrow down the number of potential suspects because it can provide some general characteristics for predicting appearance. Previous genome-wide studies based on DNA methylation have reported age prediction algorithms using a penalized multivariate regression method known as elastic net and a few dozen to hundreds of CpG sites. Although more CpG sites may provide better accuracy than fewer CpG sites, this approach is not applicable to forensics because the amounts of crime-scene DNA are usually limited. In this study, we selected three age-correlated CpG sites, namely cg16867657 (ELOVL2), which is known to be an excellent age predictor, cg04208403 (ZNF423), and cg19283806 (CCDC102B), from HumanMethylation450 BeadChip datasets of 1415 individuals. Furthermore, we evaluated these markers in a 535-sample training set and a 230-sample validation set from Korean individuals using a pyrosequencing platform. From the training set, an age prediction model using the multiple linear regression method explained 91.44% of age-correlated variation in DNA methylation patterns. The standard error of estimate and mean absolute deviation were 6.320 and 3.156 years, respectively. In the validation set, the standard error of estimate and mean absolute deviation were estimated as 6.853 and 3.346 years, respectively. For the validation set, the model explained 91.08% of the variation in methylation and predicted age (±6years) with accuracy of 77.30% in the <60years age group and 57.30% in the older group (≥60 years). These results suggest that our three DNA methylation markers may be useful for age prediction in samples from Asian populations.

Infection-specific phosphorylation of glutamyl-prolyl tRNA synthetase induces antiviral immunity

The mammalian cytoplasmic multi-tRNA synthetase complex (MSC) is a depot system that regulates non-translational cellular functions. Here we found that the MSC component glutamyl-prolyl-tRNA synthetase (EPRS) switched its function following viral infection and exhibited potent antiviral activity. Infection-specific phosphorylation of EPRS at Ser990 induced its dissociation from the MSC, after which it was guided to the antiviral signaling pathway, where it interacted with PCBP2, a negative regulator of mitochondrial antiviral signaling protein (MAVS) that is critical for antiviral immunity. This interaction blocked PCBP2-mediated ubiquitination of MAVS and ultimately suppressed viral replication. EPRS-haploid (Eprs+/−) mice showed enhanced viremia and inflammation and delayed viral clearance. This stimulus-inducible activation of MAVS by EPRS suggests an unexpected role for the MSC as a regulator of immune responses to viral infection.

Genomic and epigenomic heterogeneity in molecular subtypes of gastric cancer

Gastric cancer is a complex disease that is affected by multiple genetic and environmental factors. For the precise diagnosis and effective treatment of gastric cancer, the heterogeneity of the disease must be simplified; one way to achieve this is by dividing the disease into subgroups. Toward this effort, recent advances in high-throughput sequencing technology have revealed four molecular subtypes of gastric cancer, which are classified as Epstein-Barr virus-positive, microsatellite instability, genomically stable, and chromosomal instability subtypes. We anticipate that this molecular subtyping will help to extend our knowledge for basic research purposes and will be valuable for clinical use. Here, we review the genomic and epigenomic heterogeneity of the four molecular subtypes of gastric cancer. We also describe a mutational meta-analysis and a reanalysis of DNA methylation that were performed using previously reported gastric cancer datasets.

TWEAK/Fn14 signaling mediates gastric cancer cell resistance to 5-fluorouracil via NF-κB activation

Chemoresistance is one of the most serious problems in the treatment of cancer. In the present study, we show that Fn14 promotes resistance to 5-fluorouracil (5-FU) in gastric cancer (GC). We found that 5-FU treatment upregulated Fn14 expression in various cancer cell lines, including GC cell lines, and that knockdown of Fn14 using shRNA accelerated 5-FU sensitivity. In contrast, Fn14 overexpression or TWEAK treatment promoted resistance to 5-FU. Furthermore, we investigated the mechanisms underlying Fn14-mediated chemoresistance. We first revealed that 5-FU-mediated upregulation of Fn14 occurred as a result of NF-κB activation, indicating that 5-FU-mediated NF-κB activation was the principal event underlying Fn14 upregulation and 5-FU resistance in GC. Taken together, our results suggest that Fn14 is a novel therapeutic target and that inhibition of Fn14 combined with 5-FU treatment may be an effective molecular therapeutic strategy to treat 5-FU-resistant gastric cancers.

Genome-wide DNA methylation profiles altered by Helicobacter pylori in gastric mucosa and blood leukocyte DNA

Purpose To investigate Helicobacter pylori (H.pylori) associated genome-wide aberrant methylation patterns in gastric mucosa and blood leukocyte DNA, a population-based study was conducted in Linqu County. Results A total of 3000 and 386 CpGs were differentially methylated after successful H.pylori eradication in gastric mucosa and blood leukocyte DNA respectively, and 17 were the same alteration trend in the both tissues. The differentially methylated CpGs were located more frequently in promoters or CpG islands for gastric mucosa and gene body or open sea for blood leukocyte DNA. In eradicated gastric mucosa, the hypermethylated CpGs were enriched across inflammatory pathways, while the hypomethylated CpGs in tube morphogenesis, development and so on. The final validation found lower SPI1, PRIC285 and S1PR4 methylation levels in H.pylori positive subjects by case-control comparison, and increased methylation levels in H.pylori eradicated gastric mucosa by self-comparison. The Cancer Genome Atlas (TCGA) database analysis suggested that the up-regulation of the three genes by hypomethylation might be associated with gastric carcinogenesis. Experimental Design Infinium HumanMethylation 450K BeadChip was used to compare methylation profiles prior to and after eradication treatment. The methylation levels of identified candidate differentially methylated genes before and after H.pylori eradication were further validated by two stages (Stage I: self-comparison of 16 subjects before and after anti-H.pylori treatment; Stage II: case-control comparison of 25 H.pylori positive and 25 negative subjects and self-comparison of 50 anti-H.pylori treated subjects). Conclusions Novel H.pylori associated aberrant methylated genes were identified across the whole genome both in gastric mucosa and blood leukocyte DNA.

APPEX: analysis platform for the identification of prognostic gene expression signatures in cancer

SUMMARY Because cancer has heterogeneous clinical behaviors due to the progressive accumulation of multiple genetic and epigenetic alterations, the identification of robust molecular signatures for predicting cancer outcome is profoundly important. Here, we introduce the APPEX Web-based analysis platform as a versatile tool for identifying prognostic molecular signatures that predict cancer diversity. We incorporated most of statistical methods for survival analysis and implemented seven survival analysis workflows, including CoxSingle, CoxMulti, IntransSingle, IntransMulti, SuperPC, TimeRoc and multivariate. A total of 236 publicly available datasets were collected, processed and stored to support easy independent validation of prognostic signatures. Two case studies including disease recurrence and bladder cancer progression were described using different combinations of the seven workflows. AVAILABILITY AND IMPLEMENTATION APPEX is freely available at http://www.appex.kr. CONTACT kimsy@kribb.re.kr SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Development of zebrafish medulloblastoma-like PNET model by TALEN-mediated somatic gene inactivation

Genetically engineered animal tumor models have traditionally been generated by the gain of single or multiple oncogenes or the loss of tumor suppressor genes; however, the development of live animal models has been difficult given that cancer phenotypes are generally induced by somatic mutation rather than by germline genetic inactivation. In this study, we developed somatically mutated tumor models using TALEN-mediated somatic gene inactivation of cdkn2a/b or rb1 tumor suppressor genes in zebrafish. One-cell stage injection of cdkn2a/b-TALEN mRNA resulted in malignant peripheral nerve sheath tumors with high frequency (about 39%) and early onset (about 35 weeks of age) in F0 tp53e7/e7 mutant zebrafish. Injection of rb1-TALEN mRNA also led to the formation of brain tumors at high frequency (58%, 31 weeks of age) in F0 tp53e7/e7 mutant zebrafish. Analysis of each tumor induced by somatic inactivation showed that the targeted genes had bi-allelic mutations. Tumors induced by rb1 somatic inactivation were characterized as medulloblastoma-like primitive neuroectodermal tumors based on incidence location, histopathological features, and immunohistochemical tests. In addition, 3′ mRNA Quanti-Seq analysis showed differential activation of genes involved in cell cycle, DNA replication, and protein synthesis; especially, genes involved in neuronal development were up-regulated.

A ChIP-Seq Data Analysis Pipeline Based on Bioconductor Packages

Nowadays, huge volumes of chromatin immunoprecipitation-sequencing (ChIP-Seq) data are generated to increase the knowledge on DNA-protein interactions in the cell, and accordingly, many tools have been developed for ChIP-Seq analysis. Here, we provide an example of a streamlined workflow for ChIP-Seq data analysis composed of only four packages in Bioconductor: dada2, QuasR, mosaics, and ChIPseeker. ‘dada2’ performs trimming of the high-throughput sequencing data. ‘QuasR’ and ‘mosaics’ perform quality control and mapping of the input reads to the reference genome and peak calling, respectively. Finally, ‘ChIPseeker’ performs annotation and visualization of the called peaks. This workflow runs well independently of operating systems (e.g., Windows, Mac, or Linux) and processes the input fastq files into various results in one run. R code is available at github: https://github.com/ddhb/Workflow_of_Chipseq.git.

Non-negligible Occurrence of Errors in Gender Description in Public Data Sets

Due to advances in omics technologies, numerous genome-wide studies on human samples have been published, and most of the omics data with the associated clinical information are available in public repositories, such as Gene Expression Omnibus and ArrayExpress. While analyzing several public datasets, we observed that errors in gender information occur quite often in public datasets. When we analyzed the gender description and the methylation patterns of gender-specific probes (glucose-6-phosphate dehydrogenase [G6PD], ephrin-B1 [EFNB1], and testis specific protein, Y-linked 2 [TSPY2]) in 5,611 samples produced using Infinium 450K HumanMethylation arrays, we found that 19 samples from 7 datasets were erroneously described. We also analyzed 1,819 samples produced using the Affymetrix U133Plus2 array using several gender-specific genes (X (inactive)-specific transcript [XIST], eukaryotic translation initiation factor 1A, Y-linked [EIF1AY], and DEAD [Asp-Glu-Ala-Asp] box polypeptide 3, Y-linked [DDDX3Y]) and found that 40 samples from 3 datasets were erroneously described. We suggest that the users of public datasets should not expect that the data are error-free and, whenever possible, that they should check the consistency of the data.

Vascular defects of DYRK1A knockouts are ameliorated by modulating calcium signaling in zebrafish

DYRK1A is a major causative gene in Down syndrome (DS). Reduced incidence of solid tumors and vascular anomalies in DS patients suggests a role of DYRK1A in angiogenesis, but in vivo evidence is lacking. Here, we used zebrafish dyrk1aa mutant embryos to understand DYRK1A function in the cerebral vasculature formation. Zebrafish dyrk1aa mutants exhibited cerebral hemorrhage and defects in angiogenesis of central arteries in the developing hindbrain. Such phenotypes were rescued by wild-type dyrk1aa mRNA, but not by a kinase-dead form, indicating the importance of DYRK1A kinase activity. Chemical screening using a bioactive small molecule library identified a calcium chelator, EGTA, as one of the hits that most robustly rescued the hemorrhage. Vascular defects of mutants were also rescued by independent modulation of calcium signaling by FK506. Furthermore, the transcriptomic analyses supported the alterations of calcium signaling networks in dyrk1aa mutants. Together, our results suggest that dyrk1aa plays an essential role in angiogenesis and in maintenance of the developing cerebral vasculature via regulation of calcium signaling, which may have therapeutic potential for DYRK1A-related vascular diseases.