Je-Keun Rhee

Molecular Basis for the Recognition of Primary microRNAs by the Drosha-DGCR8 Complex

The Drosha-DGCR8 complex initiates microRNA maturation by precise cleavage of the stem loops that are embedded in primary transcripts (pri-miRNAs). Here we propose a model for this process that is based upon evidence from both computational and biochemical analyses. A typical metazoan pri-miRNA consists of a stem of approximately 33 bp, with a terminal loop and flanking segments. The terminal loop is unessential, whereas the flanking ssRNA segments are critical for processing. The cleavage site is determined mainly by the distance (approximately 11 bp) from the stem-ssRNA junction. Purified DGCR8, but not Drosha, interacts with pri-miRNAs both directly and specifically, and the flanking ssRNA segments are vital for this binding to occur. Thus, DGCR8 may function as the molecular anchor that measures the distance from the dsRNA-ssRNA junction. Our current study thus facilitates the prediction of novel microRNAs and will assist in the rational design of small hairpin RNAs for RNA interference.

Modulation of CaV2.3 Calcium Channel Currents by Eugenol

Eugenol, a natural congener of capsaicin, is a routine analgesic agent in dentistry. We have recently demonstrated the inhibition of CaV2.2 calcium channel and sodium channel currents to be molecular mechanisms underlying the analgesic effect of eugenol. We hypothesized that CaV2.3 channels are also modulated by eugenol and investigated its mode of action using the whole-cell patch-clamp technique in a heterologous expression system. Eugenol inhibited calcium currents in the E52 cell line, stably expressing the human CaV2.3 calcium channels, where TRPV1 is not endogenously expressed. The extent of current inhibition was not significantly different between naïve E52 cells and TRPV1-expressing E52 cells, suggesting no involvement of TRPV1. In contrast, TRPV1 activation is prerequisite for the inhibition of CaV2.3 calcium channels by capsaicin. The results indicate that eugenol has mechanisms distinct from those of capsaicin for modulating CaV2.3 channels. We suggest that inhibition of CaV2.3 channels by eugenol might contribute to its analgesic effect.

Survey of computational haplotype determination methods for single individual

Abstract Genome-wide association studies have expanded our understanding of the relationship between the human genome and disease. However, because of current technical limitations, it is still challenging to clearly resolve diploid sequences, that is, two copies for each chromosome. One copy of each chromosome is inherited from each parent and the genomic function is determined by the interplay between the alleles represented as genotypes in the diploid sequences. Thus, to understand the nature of genetic variation in biological processes, including disease, it is necessary to determine the complete genomic sequence of each haplotype. Although there are experimental approaches for haplotype sequencing that physically separate the chromosomes, these methods are expensive and laborious and require special equipment. Here, we review the computational approaches that can be used to determine the haplotype phase. Since 1990, many researchers have tried to reconstruct the haplotype phase using a variety of computational methods, and some researches have been successfully help to determine the haplotype phase. In this review, we investigate how the computational haplotype determination methods have been developed, and we present the remaining problems affecting the determination of the haplotype of single individual using next-generation sequencing methods.

Clonal origins and parallel evolution of regionally synchronous colorectal adenoma and carcinoma

Although the colorectal adenoma-to-carcinoma sequence represents a classical cancer progression model, the evolution of the mutational landscape underlying this model is not fully understood. In this study, we analyzed eight synchronous pairs of colorectal high-grade adenomas and carcinomas, four microsatellite-unstable (MSU) and four -stable (MSS) pairs, using whole-exome sequencing. In the MSU adenoma-carcinoma pairs, we observed no subclonal mutations in adenomas that became fixed in paired carcinomas, suggesting a ‘parallel’ evolution of synchronous adenoma-to-carcinoma, rather than a ‘stepwise’ evolution. The abundance of indel (in MSU and MSS pairs) and microsatellite instability (in MSU pairs) was noted in the later adenoma- or carcinoma-specific mutations, indicating that the mutational processes and functional constraints operative in early and late colorectal carcinogenesis are different. All MSU cases exhibited clonal, truncating mutations in ACVR2A, TGFBR2, and DNA mismatch repair genes, but none were present in APC or KRAS. In three MSS pairs, both APC and KRAS mutations were identified as both early and clonal events, often accompanying clonal copy number changes. An MSS case uniquely exhibited clonal ERBB2 amplification, followed by APC and TP53 mutations as carcinoma-specific events. Along with the previously unrecognized clonal origins of synchronous colorectal adenoma-carcinoma pairs, our study revealed that the preferred sequence of mutational events during colorectal carcinogenesis can be context-dependent.

Intraindividual genomic heterogeneity of high-grade serous carcinoma of the ovary and clinical utility of ascitic cancer cells for mutation profiling

Intraindividual tumoural heterogeneity (ITH) is a hallmark of solid tumours and impedes accurate genomic diagnosis and selection of proper therapy. The aim of this study was to identify ITH of ovarian high‐grade serous carcinomas (OSCs) and to determine the utility of ascitic cancer cells as a resource for mutation profiling in spite of ITH. We performed whole‐exome sequencing, copy number profiling and DNA methylation profiling of four OSC genomes by using multiregional biopsies from 13 intraovarian lesions, 12 extraovarian tumour lesions (omentum/peritoneum), and ascitic cells. We observed substantial levels of heterogeneity in mutations and copy number alterations (CNAs) of the OSCs. We categorized the mutations into ‘common’, ‘shared’ and ‘private’ according to the regional distribution. Six common, eight shared and 24 private mutations were observed in known cancer‐related genes. Common mutations had a higher mutant allele frequency, and included TP53 mutations in all four OSCs. Region‐specific chromosomal amplifications and deletions involving BRCA1, PIK3CA and RB1 were also identified. It is of note that the mutations detected in ascitic cancer cells represented 92.3–100% of overall somatic mutations in the given case. Phylogenetic analyses of ascitic genomes predicted a polyseeding origin of somatic mutations in ascitic cells. Our results demonstrate that, despite ITH, somatic mutations, CNAs and DNA methylations in both ‘common’ category and cancer‐related genes were highly conserved in ascitic cells of OSCs, highlighting the clinical relevance of genome analysis of ascitic cells. Ascitic tumour cells may serve as a potential resource for discovering somatic mutations of primary OSC with diagnostic and therapeutic relevance. Copyright

AESNB: Active Example Selection with Naïve Bayes Classifier for Learning from Imbalanced Biomedical Data

Various real-world biomedical classification tasks suffer from the imbalanced data problem which tends to make the prediction performance of some classes significantly decrease. In this paper, we present an active example selection method with naïve Bayes classifier (AESNB) as a solution for the imbalanced data problem. The proposed method starts with a small balanced subset of training examples. A naïve Bayes classifier is trained incrementally by actively selecting and adding informative examples regardless of the original class distribution. Informative examples are defined as examples that produce high error scores by the current classifier. We examined the performance of AESNB algorithm by using five imbalanced biomedical datasets. Our experimental results show that the naïve Bayes classifier with our active example selection method achieves a competitive classification performance compared to the classifier with sampling or cost-sensitive methods.

Identification of cell cycle-related regulatory motifs using a kernel canonical correlation analysis

BackgroundGene regulation is a key mechanism in higher eukaryotic cellular processes. One of the major challenges in gene regulation studies is to identify regulators affecting the expression of their target genes in specific biological processes. Despite their importance, regulators involved in diverse biological processes still remain largely unrevealed. In the present study, we propose a kernel-based approach to efficiently identify core regulatory elements involved in specific biological processes using gene expression profiles.ResultsWe developed a framework that can detect correlations between gene expression profiles and the upstream sequences on the basis of the kernel canonical correlation analysis (kernel CCA). Using a yeast cell cycle dataset, we demonstrated that upstream sequence patterns were closely related to gene expression profiles based on the canonical correlation scores obtained by measuring the correlation between them. Our results showed that the cell cycle-specific regulatory motifs could be found successfully based on the motif weights derived through kernel CCA. Furthermore, we identified co-regulatory motif pairs using the same framework.ConclusionGiven expression profiles, our method was able to identify regulatory motifs involved in specific biological processes. The method could be applied to the elucidation of the unknown regulatory mechanisms associated with complex gene regulatory processes.

Generation and molecular characterization of pancreatic cancer patient-derived xenografts reveals their heterologous nature

Pancreatic ductal adenocarcinoma (PDAC) is the most challenging type of cancer to treat, with a 5-year survival rate of <10%. Furthermore, because of the large portion of the inoperable cases, it is difficult to obtain specimens to study the biology of the tumors. Therefore, a patient-derived xenograft (PDX) model is an attractive option for preserving and expanding these tumors for translational research. Here we report the generation and characterization of 20 PDX models of PDAC. The success rate of the initial graft was 74% and most tumors were re-transplantable. Histological analysis of the PDXs and primary tumors revealed a conserved expression pattern of p53 and SMAD4; an exome single nucleotide polymorphism (SNP) array and Comprehensive Cancer Panel showed that PDXs retained over 94% of cancer-associated variants. In addition, Polyphen2 and the Sorting Intolerant from Tolerant (SIFT) prediction identified 623 variants among the functional SNPs, highlighting the heterologous nature of pancreatic PDXs; an analysis of 409 tumor suppressor genes and oncogenes in Comprehensive Cancer Panel revealed heterologous cancer gene mutation profiles for each PDX-primary tumor pair. Altogether, we expect these PDX models are a promising platform for screening novel therapeutic agents and diagnostic markers for the detection and eradication of PDAC.

Whole-exome sequencing identified mutational profiles of high-grade colon adenomas

Although gene-to-gene analyses identified genetic alterations such as APC, KRAS and TP53 mutations in colon adenomas, it is largely unknown whether there are any others in them. Mutational profiling of high-grade colon adenoma (HGCA) that just precedes colon carcinoma might identify not only novel adenoma-specific genes but also critical genes for its progression to carcinoma. For this, we performed whole-exome sequencing (WES) of 12 HGCAs and identified 11 non-hypermutated and one hypermutated (POLE-mutated) cases. We identified 22 genes including APC, KRAS, TP53, GNAS, NRAS, SMAD4, ARID2, and PIK3CA with non-silent mutations in the cancer Census Genes. Bi-allelic and mono-allelic APC alterations were found in nine and one HGCAs, respectively, while the other two harbored wild-type APC. Five HGCAs harbored either mono-allelic (four HGCAs) or bi-allelic (one HGCA) SMAD4 mutation or 18q loss that had been known as early carcinoma-specific changes. We identified MTOR, ACVR1B, GNAQ, ATM, CNOT1, EP300, ARID2, RET and MAP2K4 mutations for the first time in colon adenomas. Our WES data is largely matched with the earlier ‘adenoma-carcinoma model’ (APC, KRAS, NRAS and GNAS mutations), but there are newly identified SMAD4, MTOR, ACVR1B, GNAQ, ATM, CNOT1, EP300, ARID2, RET and MAP2K4 mutations in this study. Our findings provide resource for understanding colon premalignant lesions and for identifying genomic clues for differential diagnosis and therapy options for colon adenomas and carcinomas.

The chronological sequence of somatic mutations in early gastric carcinogenesis inferred from multiregion sequencing of gastric adenomas

Mutation profiles and intratumoral heterogeneity are not well understood for benign gastric adenomas, some of which progress into malignant gastric adenocarcinomas. In this study, we performed whole-exome sequencing of three microsatellite stable (MSS) and two microsatellite instability-high (MSI-H) gastric adenomas with three regional tumor biopsies per case. We observed that the mutation abundance of benign gastric adenomas was comparable to those of gastric adenocarcinomas, suggesting that the mutational makeup for gastric carcinogenesis may already be achieved in benign adenomas. The extent of intratumoral heterogeneity was more substantial for MSS genomes in that only 1% - 14% of somatic mutations were common across the regional biopsies or ‘public’, while 50% - 94% of mutations were public in MSI-H gastric adenomas. We observed biallelic, loss-of-functional events of APC with truncating mutations and/or 5q losses for all cases, mostly as public events. All MSS gastric adenomas also harbored ARID2 truncating mutations, often as multiple, region-specific ones indicative of convergent evolution. Hotspot missense mutations on known cancer genes such as ERBB2 and KRAS were largely observed as region-specific aberrations. These findings suggest that biallelic functional APC inactivation initiates the gastric carcinogenesis and is followed by mutations of histone modifiers and then activation of known cancer-related genes. As the first exome-wide multi-region mutational profiling of gastric adenomas, our study provides clues on the chronological sequence of somatic mutations and their clonal architectures in early gastric carcinogenesis.