Hwanseok Rhee

The transcriptional landscape and mutational profile of lung adenocarcinoma

All cancers harbor molecular alterations in their genomes. The transcriptional consequences of these somatic mutations have not yet been comprehensively explored in lung cancer. Here we present the first large scale RNA sequencing study of lung adenocarcinoma, demonstrating its power to identify somatic point mutations as well as transcriptional variants such as gene fusions, alternative splicing events, and expression outliers. Our results reveal the genetic basis of 200 lung adenocarcinomas in Koreans including deep characterization of 87 surgical specimens by transcriptome sequencing. We identified driver somatic mutations in cancer genes including EGFR, KRAS, NRAS, BRAF, PIK3CA, MET, and CTNNB1. Candidates for novel driver mutations were also identified in genes newly implicated in lung adenocarcinoma such as LMTK2, ARID1A, NOTCH2, and SMARCA4. We found 45 fusion genes, eight of which were chimeric tyrosine kinases involving ALK, RET, ROS1, FGFR2, AXL, and PDGFRA. Among 17 recurrent alternative splicing events, we identified exon 14 skipping in the proto-oncogene MET as highly likely to be a cancer driver. The number of somatic mutations and expression outliers varied markedly between individual cancers and was strongly correlated with smoking history of patients. We identified genomic blocks within which gene expression levels were consistently increased or decreased that could be explained by copy number alterations in samples. We also found an association between lymph node metastasis and somatic mutations in TP53. These findings broaden our understanding of lung adenocarcinoma and may also lead to new diagnostic and therapeutic approaches.

An epigenomic roadmap to induced pluripotency reveals DNA methylation as a reprogramming modulator

Reprogramming of somatic cells to induced pluripotent stem cells involves a dynamic rearrangement of the epigenetic landscape. To characterize this epigenomic roadmap, we have performed MethylC-seq, ChIP-seq (H3K4/K27/K36me3) and RNA-Seq on samples taken at several time points during murine secondary reprogramming as part of Project Grandiose. We find that DNA methylation gain during reprogramming occurs gradually, while loss is achieved only at the ESC-like state. Binding sites of activated factors exhibit focal demethylation during reprogramming, while ESC-like pluripotent cells are distinguished by extension of demethylation to the wider neighbourhood. We observed that genes with CpG-rich promoters demonstrate stable low methylation and strong engagement of histone marks, whereas genes with CpG-poor promoters are safeguarded by methylation. Such DNA methylation-driven control is the key to the regulation of ESC-pluripotency genes, including Dppa4, Dppa5a and Esrrb. These results reveal the crucial role that DNA methylation plays as an epigenetic switch driving somatic cells to pluripotency.

Gene expression defines natural changes in mammalian lifespan

Mammals differ more than 100‐fold in maximum lifespan, which can be altered in either direction during evolution, but the molecular basis for natural changes in longevity is not understood. Divergent evolution of mammals also led to extensive changes in gene expression within and between lineages. To understand the relationship between lifespan and variation in gene expression, we carried out RNA‐seq‐based gene expression analyses of liver, kidney, and brain of 33 diverse species of mammals. Our analysis uncovered parallel evolution of gene expression and lifespan, as well as the associated life‐history traits, and identified the processes and pathways involved. These findings provide direct insights into how nature reversibly adjusts lifespan and other traits during adaptive radiation of lineages.

Heterozygous mutations in cyclic AMP phosphodiesterase-4D (PDE4D) and protein kinase A (PKA) provide new insights into the molecular pathology of acrodysostosis

Acrodysostosis without hormone resistance is a rare skeletal disorder characterized by brachydactyly, nasal hypoplasia, mental retardation and occasionally developmental delay. Recently, loss-of-function mutations in the gene encoding cAMP-hydrolyzing phosphodiesterase-4D (PDE4D) have been reported to cause this rare condition but the pathomechanism has not been fully elucidated. To understand the pathogenetic mechanism of PDE4D mutations, we conducted 3D modeling studies to predict changes in the binding efficacy of cAMP to the catalytic pocket in PDE4D mutants. Our results indicated diminished enzyme activity in the two mutants we analyzed (Gly673Asp and Ile678Thr; based on PDE4D4 residue numbering). Ectopic expression of PDE4D mutants in HEK293 cells demonstrated this reduction in activity, which was identified by increased cAMP levels. However, the cells from an acrodysostosis patient showed low cAMP accumulation, which resulted in a decrease in the phosphorylated cAMP Response Element-Binding Protein (pCREB)/CREB ratio. The reason for this discrepancy was due to a compensatory increase in expression levels of PDE4A and PDE4B isoforms, which accounted for the paradoxical decrease in cAMP levels in the patient cells expressing mutant isoforms with a lowered PDE4D activity. Skeletal radiographs of 10-week-old knockout (KO) rats showed that the distal part of the forelimb was shorter than in wild-type (WT) rats and that all the metacarpals and phalanges were also shorter in KO, as the name acrodysostosis implies. Like the G-protein α-stimulatory subunit and PRKAR1A, PDE4D critically regulates the cAMP signal transduction pathway and influences bone formation in a way that activity-compromising PDE4D mutations can result in skeletal dysplasia. We propose that specific inhibitory PDE4D mutations can lead to the molecular pathology of acrodysostosis without hormone resistance but that the pathological phenotype may well be dependent on an over-compensatory induction of other PDE4 isoforms that can be expected to be targeted to different signaling complexes and exert distinct effects on compartmentalized cAMP signaling.

Targeted exon sequencing fails to identify rare coding variants with large effect in rheumatoid arthritis.

IntroductionAlthough it has been suggested that rare coding variants could explain the substantial missing heritability, very few sequencing studies have been performed in rheumatoid arthritis (RA). We aimed to identify novel functional variants with rare to low frequency using targeted exon sequencing of RA in Korea.MethodsWe analyzed targeted exon sequencing data of 398 genes selected from a multifaceted approach in Korean RA patients (n = 1,217) and controls (n = 717). We conducted a single-marker association test and a gene-based analysis of rare variants. For meta-analysis or enrichment tests, we also used ethnically matched independent samples of Korean genome-wide association studies (GWAS) (n = 4,799) or immunochip data (n = 4,722).ResultsAfter stringent quality control, we analyzed 10,588 variants of 398 genes from 1,934 Korean RA case controls. We identified 13 nonsynonymous variants with nominal association in single-variant association tests. In a meta-analysis, we did not find any novel variant with genome-wide significance for RA risk. Using a gene-based approach, we identified 17 genes with nominal burden signals. Among them, VSTM1 showed the greatest association with RA (P = 7.80 × 10-4). In the enrichment test using Korean GWAS, although the significant signal appeared to be driven by total genic variants, we found no evidence for enriched association of coding variants only with RA.ConclusionsWe were unable to identify rare coding variants with large effect to explain the missing heritability for RA in the current targeted resequencing study. Our study raises skepticism about exon sequencing of targeted genes for complex diseases like RA.

miR‐204 downregulates EphB2 in aging mouse hippocampal neurons

Hippocampal synaptic function and plasticity deteriorate with age, often resulting in learning and memory deficits. As MicroRNAs (miRNAs) are important regulators of neuronal protein expression, we examined whether miRNAs may contribute to this age‐associated decline in hippocampal function. We first compared the small RNA transcriptome of hippocampal tissues from young and old mice. Among 269 hippocampal miRNAs, 80 were differentially expressed (≥ twofold) among the age groups. We focused on 36 miRNAs upregulated in the old mice compared with those in the young mice. The potential targets of these 36 miRNAs included 11 critical Eph/Ephrin synaptic signaling components. The expression levels of several genes in the Eph/Ephrin pathway, including EphB2, were significantly downregulated in the aged hippocampus. EphB2 is a known regulator of synaptic plasticity in hippocampal neurons, in part by regulating the surface expression of the NMDA receptor NR1 subunit. We found that EphB2 is a direct target of miR‐204 among miRNAs that were upregulated with age. The transfection of primary hippocampal neurons with a miR‐204 mimic suppressed both EphB2 mRNA and protein expression and reduced the surface expression of NR1. Transfection of miR‐204 also decreased the total expression of NR1. miR‐204 induces senescence‐like phenotype in fully matured neurons as evidenced by an increase in p16‐positive cells. We suggest that aging is accompanied by the upregulation of miR‐204 in the hippocampus, which downregulates EphB2 and results in reduced surface and total NR1 expression. This mechanism may contribute to age‐associated decline in hippocampal synaptic plasticity and the related cognitive functions.

Recurrent fusion transcripts detected by whole-transcriptome sequencing of 120 primary breast cancer samples

Relatively few recurrent gene fusion events have been associated with breast cancer to date. In an effort to uncover novel fusion transcripts, we performed whole‐transcriptome sequencing of 120 fresh‐frozen primary breast cancer samples and five adjacent normal breast tissues using the Illumina HiSeq2000 platform. Three different fusion‐detecting tools (deFuse, Chimerascan, and TopHatFusion) were used, and the results were compared. These tools detected 3,831, 6,630 and 516 fusion transcripts (FTs) overall. We primarily focused on the results obtained using the deFuse software. More FTs were identified from HER2 subtype breast cancer samples than from the luminal or triple‐negative subtypes (P < 0.05). Seventy fusion candidates were selected for validation, and 32 (45.7%) were confirmed by RT‐PCR and Sanger sequencing. Of the validated fusions, six were recurrent (found in 2 or more samples), three were in‐frame (PRDX1‐AKR1A1, TACSTD2‐OMA1, and C2CD2‐TFF1) and three were off‐frame (CEACAM7‐CEACAM6, CYP4X1‐CYP4Z2P, and EEF1DP3‐FRY). Notably, the novel read‐through fusion, EEF1DP3‐FRY, was identified and validated in 6.7% (8/120) of the breast cancer samples. This off‐frame fusion results in early truncation of the FRY gene, which plays a key role in the structural integrity during mitosis. Three previously reported fusions, PPP1R1B‐STARD3, MFGE8‐HAPL, and ETV6‐NTRK3, were detected in 8.3, 3.3, and 0.8% of the 120 samples, respectively, by both deFuse and Chimerascan. The recently reported MAGI3‐AKT3 fusion was not detected in our analysis. Although future work will be needed to examine the biological significance of our new findings, we identified a number of novel fusions and confirmed some previously reported fusions.