Proteogenomic Characterization Reveals Estrogen Signaling as a Target for Never-Smoker Lung Adenocarcinoma Patients without EGFR or ALK Alterations.

Never-smoker lung adenocarcinoma (NSLA) is prevalent in Asian populations, particularly in women. EGFR mutations and anaplastic lymphoma kinase (ALK) fusions are major genetic alterations observed in NSLA, and NSLA with these alterations have been well studied and can be treated with targeted therapies. To provide insights into the molecular profile of NSLA without EGFR and ALK alterations (NENA), we selected 141 NSLA tissues and performed proteogenomic characterization, including whole genome sequencing (WGS), transcriptomic, methylation EPIC array, total proteomic, and phosphoproteomic analyses. Forty patients with NSLA harboring EGFR and ALK alterations and seven patients with NENA with microsatellite instability were excluded. Genome analysis revealed that TP53 (25%), KRAS (22%), and SETD2 (11%) mutations and ROS1 fusions (14%) were the most frequent genetic alterations in NENA patients. Proteogenomic impact analysis revealed that STK11 and ERBB2 somatic mutations had broad effects on cancer-associated genes in NENA. DNA copy number alteration analysis identified 22 prognostic proteins that influenced transcriptomic and proteomic changes. Gene set enrichment analysis revealed estrogen signaling as the key pathway activated in NENA. Increased estrogen signaling was associated with proteogenomic alterations, such as copy number deletions in chromosomes 14 and 21, STK11 mutation, and DNA hypomethylation of LLGL2 and ST14. Finally, saracatinib, an Src inhibitor, was identified as a potential drug for targeting activated estrogen signaling in NENA and was experimentally validated in vitro. Collectively, this study enhanced our understanding of NENA NSLA by elucidating the proteogenomic landscape and proposed saracatinib as a potential treatment for this patient population that lacks effective targeted therapies. SIGNIFICANCE: The proteogenomic landscape in never-smoker lung cancer without known driver mutations reveals prognostic proteins and enhanced estrogen signaling that can be targeted as a potential therapeutic strategy to improve patient outcomes.

Introduction of the Korea BioData Station (K-BDS) for sharing biological data.

A wave of new technologies has created opportunities for the cost-effective generation of high-throughput profiles of biological systems, foreshadowing a "data-driven science" era. The large variety of data available from biological research is also a rich resource that can be used for innovative endeavors. However, we are facing considerable challenges in big data deposition, integration, and translation due to the complexity of biological data and its production at unprecedented exponential rates. To address these problems, in 2020, the Korean government officially announced a national strategy to collect and manage the biological data produced through national R&D fund allocations and provide the collected data to researchers. To this end, the Korea Bioinformation Center (KOBIC) developed a new biological data repository, the Korea BioData Station (K-BDS), for sharing data from individual researchers and research programs to create a data-driven biological study environment. The K-BDS is dedicated to providing free open access to a suite of featured data resources in support of worldwide activities in both academia and industry.

Aberrant Methylation of Somatostatin Receptor 2 Gene Is Initiated in Aged Gastric Mucosa Infected with Helicobacter pylori and Consequential Gene Silencing Is Associated with Establishment of Inflammatory Microenvironment In Vitro Study.

The loss-of-function variants are thought to be associated with inflammation in the stomach. We here aimed to evaluate the extent and role of methylation at the SSTR2 promoter in inflammation and gastric tumor formation. A whole-genome bisulfite sequencing analysis revealed that the SSTR2 promoter was significantly hypermethylated in gastric tumors, dysplasia, and intestinal metaplasia compared to non-tumor tissues from patients with gastric cancer. Using public data, we confirmed SSTR2 promoter methylation in primary gastric tumors and intestinal metaplasia, and even aged gastric mucosae infected with Helicobacter pylori, suggesting that aberrant methylation is initiated in normal gastric mucosa. The loss-of-function of SSTR2 in SNU638 cell-induced cell proliferation in vitro, while stable transfection of SSTR2 in AGS and MKN74 cells inhibited cell proliferation and tumorigenesis in vitro and in vivo. As revealed by a comparison of target genes differentially expressed in these cells with hallmark molecular signatures, inflammation-related pathways were distinctly induced in SSTR2-KO SNU638 cell. By contrast, inflammation-related pathways were inhibited in AGS and MKN74 cells ectopically expressing SSTR2. Collectively, we propose that SSTR2 silencing upon promoter methylation is initiated in aged gastric mucosae infected with H. pylori and promotes the establishment of an inflammatory microenvironment via the intrinsic pathway. These findings provide novel insights into the initiation of gastric carcinogenesis.

Intermediate cells of in vitro cellular reprogramming and in vivo tissue regeneration require desmoplakin.

Amphibians and fish show considerable regeneration potential via dedifferentiation of somatic cells into blastemal cells. In terms of dedifferentiation, in vitro cellular reprogramming has been proposed to share common processes with in vivo tissue regeneration, although the details are elusive. Here, we identified the cytoskeletal linker protein desmoplakin (Dsp) as a common factor mediating both reprogramming and regeneration. Our analysis revealed that Dsp expression is elevated in distinct intermediate cells during in vitro reprogramming. Knockdown of Dsp impedes in vitro reprogramming into induced pluripotent stem cells and induced neural stem/progenitor cells as well as in vivo regeneration of zebrafish fins. Notably, reduced Dsp expression impairs formation of the intermediate cells during cellular reprogramming and tissue regeneration. These findings suggest that there is a Dsp-mediated evolutionary link between cellular reprogramming in mammals and tissue regeneration in lower vertebrates and that the intermediate cells may provide alternative approaches for mammalian regenerative therapy.

Glycogen Storage Disease Phenotypes Accompanying the Perturbation of the Methionine Cycle in NDRG3-Deficient Mouse Livers.

N-Myc downstream regulated gene 3 (NDRG3) is a unique pro-tumorigenic member among NDRG family genes, mediating growth signals. Here, we investigated the pathophysiological roles of NDRG3 in relation to cell metabolism by disrupting its functions in liver. Mice with liver-specific KO of NDRG3 (Ndrg3 LKO) exhibited glycogen storage disease (GSD) phenotypes including excessive hepatic glycogen accumulation, hypoglycemia, elevated liver triglyceride content, and several signs of liver injury. They suffered from impaired hepatic glucose homeostasis, due to the suppression of fasting-associated glycogenolysis and gluconeogenesis. Consistently, the expression of glycogen phosphorylase (PYGL) and glucose-6-phosphate transporter (G6PT) was significantly down-regulated in an Ndrg3 LKO-dependent manner. Transcriptomic and metabolomic analyses revealed that NDRG3 depletion significantly perturbed the methionine cycle, redirecting its flux towards branch pathways to upregulate several metabolites known to have hepatoprotective functions. Mechanistically, Ndrg3 LKO-dependent downregulation of glycine N-methyltransferase in the methionine cycle and the resultant elevation of the S-adenosylmethionine level appears to play a critical role in the restructuring of the methionine metabolism, eventually leading to the manifestation of GSD phenotypes in Ndrg3 LKO mice. Our results indicate that NDRG3 is required for the homeostasis of liver cell metabolism upstream of the glucose-glycogen flux and methionine cycle and suggest therapeutic values for regulating NDRG3 in disorders with malfunctions in these pathways.

Protein Arginine Methyltransferase 1 Is Essential for the Meiosis of Male Germ Cells.

Protein arginine methyltransferase 1 (PRMT1) is a major enzyme responsible for the formation of methylarginine in mammalian cells; however, its function in vivo is not well understood due to its early embryonic lethality in null mice exhibiting spontaneous DNA damage, cell cycle delays, and defects in check point activation. Here, we generated germ cell-specific Prmt1 knock-out (KO) mice to evaluate the function of PRMT1 in spermatogenesis. Our findings demonstrate that PRMT1 is vital for male fertility in mice. Spermatogenesis in Prmt1 KO mice was arrested at the zygotene-like stage of the first meiotic division due to an elevated number of DNA double-strand breaks (DSBs). There was a loss of methylation in meiotic recombination 11 (MRE11), the key endonuclease in MRE11/RAD50/NBS 1 (MRN) complex, resulting in the accumulation of SPO11 protein in DSBs. The ATM-mediated negative feedback control over SPO11 was lost and, consequently, the repair pathway of DSBs was highly affected in PRMT1 deficient male germ cells. Our findings provide a novel insight into the role of PRMT1-mediated asymmetric demethylation in mouse spermatogenesis.

STK31 upregulation is associated with chromatin remodeling in gastric cancer and induction of tumorigenicity in a xenograft mouse model.

Pathological changes in the epigenetic landscape of chromatin are hallmarks of cancer. Our previous study showed that global methylation of promoters may increase or decrease during the transition from gastric mucosa to intestinal metaplasia (IM) to gastric cancer (GC). Here, CpG hypomethylation of the serine/threonine kinase STK31 promoter in IM and GC was detected in a reduced representation bisulfite sequencing database. STK31 hypomethylation, which resulted in its upregulation in 120 cases of primary GC, was confirmed. Using public genome­wide histone modification data, upregulation of STK31 promoter activity was detected in primary GC but not in normal mucosae, suggesting that STK31 may be repressed in gastric mucosa but activated in GC as a consequence of hypomethylation­associated chromatin remodeling. STK31 knockdown suppressed the proliferation, colony formation and migration activities of GC cells in vitro, whereas stable overexpression of STK31 promoted the proliferation, colony formation, and migration activities of GC cells in vitro and tumorigenesis in nude mice. Patients with GC in which STK31 was upregulated exhibited significantly shorter survival times in a combined cohort. Thus, activation of STK31 by chromatin remodeling may be associated with gastric carcinogenesis and also may help predict GC prognosis.

BaSDAS: a web-based pooled CRISPR-Cas9 knockout screening data analysis system.

We developed the BaSDAS (Barcode-Seq Data Analysis System), a GUI-based pooled knockout screening data analysis system, to facilitate the analysis of pooled knockout screen data easily and effectively by researchers with limited bioinformatics skills. The BaSDAS supports the analysis of various pooled screening libraries, including yeast, human, and mouse libraries, and provides many useful statistical and visualization functions with a user-friendly web interface for convenience. We expect that BaSDAS will be a useful tool for the analysis of genome-wide screening data and will support the development of novel drugs based on functional genomics information.

PRMT1 Is Required for the Maintenance of Mature ß-Cell Identity.

Loss of functional ß-cell mass is an essential feature of type 2 diabetes, and maintaining mature ß-cell identity is important for preserving a functional ß-cell mass. However, it is unclear how ß-cells achieve and maintain their mature identity. Here we demonstrate a novel function of protein arginine methyltransferase 1 (PRMT1) in maintaining mature ß-cell identity. Prmt1 knockout in fetal and adult ß-cells induced diabetes, which was aggravated by high-fat diet-induced metabolic stress. Deletion of Prmt1 in adult ß-cells resulted in the immediate loss of histone H4 arginine 3 asymmetric dimethylation (H4R3me2a) and the subsequent loss of ß-cell identity. The expression levels of genes involved in mature ß-cell function and identity were robustly downregulated as soon as Prmt1 deletion was induced in adult ß-cells. Chromatin immunoprecipitation sequencing and assay for transposase-accessible chromatin sequencing analyses revealed that PRMT1-dependent H4R3me2a increases chromatin accessibility at the binding sites for CCCTC-binding factor (CTCF) and ß-cell transcription factors. In addition, PRMT1-dependent open chromatin regions may show an association with the risk of diabetes in humans. Together, our results indicate that PRMT1 plays an essential role in maintaining ß-cell identity by regulating chromatin accessibility.

GENT2: an updated gene expression database for normal and tumor tissues.

BACKGROUND: Gene Expression database of Normal and Tumor tissues 2 (GENT2) is an updated version of GENT, which has provided a user-friendly search platform for gene expression patterns across different normal and tumor tissues compiled from public gene expression data sets. RESULTS: We refactored GENT2 with recent technologies such as Apache Lucene indexing for fast search and Google Web Toolkit (GWT) framework for a user-friendly web interface. Now, GENT2 contains more than 68,000 samples and has several new useful functions. First, GENT2 now provides gene expression across 72 different tissues compared to 57 in GENT. Second, with increasing importance of tumor subtypes, GENT2 provides an option to study the differential expression and its prognostic significance based on tumor subtypes. Third, whenever available, GENT2 provides prognostic information of a gene of interest. Fourth, GENT2 provides a meta-analysis of survival information to provide users more reliable prognostic value of a gene of interest. CONCLUSIONS: In conclusion, with these significant improvements, GENT2 will continue to be a useful tool to a wide range of researchers. GENT2 is freely available at http://gent2.appex.kr .

Dynamic Transcriptome, DNA Methylome, and DNA Hydroxymethylome Networks During T-Cell Lineage Commitment.

The stepwise development of T cells from a multipotent precursor is guided by diverse mechanisms, including interactions among lineage-specific transcription factors (TFs) and epigenetic changes, such as DNA methylation and hydroxymethylation, which play crucial roles in mammalian development and lineage commitment. To elucidate the transcriptional networks and epigenetic mechanisms underlying T-cell lineage commitment, we investigated genome-wide changes in gene expression, DNA methylation and hydroxymethylation among populations representing five successive stages of T-cell development (DN3, DN4, DP, CD4+, and CD8+) by performing RNA-seq, MBD-seq and hMeDIP-seq, respectively. The most significant changes in the transcriptomes and epigenomes occurred during the DN4 to DP transition. During the DP stage, many genes involved in chromatin modification were up-regulated and exhibited dramatic changes in DNA hydroxymethylation. We also observed 436 alternative splicing events, and approximately 57% (252) of these events occurred during the DP stage. Many stage-specific, differentially methylated regions were observed near the stage-specific, differentially expressed genes. The dynamic changes in DNA methylation and hydroxymethylation were associated with the recruitment of stage-specific TFs. We elucidated interactive networks comprising TFs, chromatin modifiers, and DNA methylation and hope that this study provides a framework for the understanding of the molecular networks underlying T-cell lineage commitment.

Vitamin C Promotes Astrocyte Differentiation Through DNA Hydroxymethylation.

Previous studies have reported that vitamin C (VC) promotes neural stem/precursor cell (NSC) differentiation toward dopamine (DA) neurons via DNA hydroxymethylation-induced transcriptional activation of DA neuron-specific genes. To further understand the VC effects on NSC differentiation, we profiled the transcriptome and DNA methylome/hydroxymethylome using high-throughput sequencing. Interestingly, RNA sequencing analyses have shown that, in addition to DA neuronal genes, astrocytic genes Gfap, Slc1a3, and S100a16 were also upregulated in NSC cultures differentiated with VC treatment. Consistently, enhanced GFAP+ astrocytic yields were manifested in the differentiated cultures with VC treatment, collectively indicating that VC promotes astrocytic differentiation. In genome-wide hydroxymethylome analyses, VC treatment induces enrichment of DNA hydroxymethylation (5-hydroxymethyl cytosine; 5hmC) near the consensus binding motifs of nuclear factor I (NFI). Furthermore, we showed that VC significantly enhanced recruitment of NFI and STAT3, key transcription factors for astrogenesis, in the 5hmC-enriched regions of the astrocyte-specific genes. These findings suggest that VC play important roles in astrocytogenesis during brain development. Stem Cells 2018;36:1578-1588.

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.

Optimizing de novo transcriptome assembly and extending genomic resources for striped catfish (Pangasianodon hypophthalmus).

Striped catfish (Pangasianodon hypophthalmus) is a commercially important freshwater fish used in inland aquaculture in the Mekong Delta, Vietnam. The culture industry is facing a significant challenge however from saltwater intrusion into many low topographical coastal provinces across the Mekong Delta as a result of predicted climate change impacts. Developing genomic resources for this species can facilitate the production of improved culture lines that can withstand raised salinity conditions, and so we have applied high-throughput Ion Torrent sequencing of transcriptome libraries from six target osmoregulatory organs from striped catfish as a genomic resource for use in future selection strategies. We obtained 12,177,770 reads after trimming and processing with an average length of 97bp. De novo assemblies were generated using CLC Genomic Workbench, Trinity and Velvet/Oases with the best overall contig performance resulting from the CLC assembly. De novo assembly using CLC yielded 66,451 contigs with an average length of 478bp and N50 length of 506bp. A total of 37,969 contigs (57%) possessed significant similarity with proteins in the non-redundant database. Comparative analyses revealed that a significant number of contigs matched sequences reported in other teleost fishes, ranging in similarity from 45.2% with Atlantic cod to 52% with zebrafish. In addition, 28,879 simple sequence repeats (SSRs) and 55,721 single nucleotide polymorphisms (SNPs) were detected in the striped catfish transcriptome. The sequence collection generated in the current study represents the most comprehensive genomic resource for P. hypophthalmus available to date. Our results illustrate the utility of next-generation sequencing as an efficient tool for constructing a large genomic database for marker development in non-model species.

Bacterial genome mapper: A comparative bacterial genome mapping tool.

UNLABELLED: Recently, next generation sequencing (NGS) technologies have led to a revolutionary increase in sequencing speed and costefficacy. Consequently, a vast number of contigs from many recently sequenced bacterial genomes remain to be accurately mapped and annotated, requiring the development of more convenient bioinformatics programs. In this paper, we present a newly developed web-based bioinformatics program, Bacterial Genome Mapper, which is suitable for mapping and annotating contigs that have been assembled from bacterial genome sequence raw data. By constructing a multiple alignment map between target contig sequences and two reference bacterial genome sequences, this program also provides very useful comparative genomics analysis of draft bacterial genomes. AVAILABILITY: The database is available for free at http://mbgm.kribb.re.kr.

Genome analysis of the domestic dog (Korean Jindo) by massively parallel sequencing.

Although pioneering sequencing projects have shed light on the boxer and poodle genomes, a number of challenges need to be met before the sequencing and annotation of the dog genome can be considered complete. Here, we present the DNA sequence of the Jindo dog genome, sequenced to 45-fold average coverage using Illumina massively parallel sequencing technology. A comparison of the sequence to the reference boxer genome led to the identification of 4 675 437 single nucleotide polymorphisms (SNPs, including 3 346 058 novel SNPs), 71 642 indels and 8131 structural variations. Of these, 339 non-synonymous SNPs and 3 indels are located within coding sequences (CDS). In particular, 3 non-synonymous SNPs and a 26-bp deletion occur in the TCOF1 locus, implying that the difference observed in cranial facial morphology between Jindo and boxer dogs might be influenced by those variations. Through the annotation of the Jindo olfactory receptor gene family, we found 2 unique olfactory receptor genes and 236 olfactory receptor genes harbouring non-synonymous homozygous SNPs that are likely to affect smelling capability. In addition, we determined the DNA sequence of the Jindo dog mitochondrial genome and identified Jindo dog-specific mtDNA genotypes. This Jindo genome data upgrade our understanding of dog genomic architecture and will be a very valuable resource for investigating not only dog genetics and genomics but also human and dog disease genetics and comparative genomics.

Novel mechanism of conjoined gene formation in the human genome.

Recently, conjoined genes (CGs) have emerged as important genetic factors necessary for understanding the human genome. However, their formation mechanism and precise structures have remained mysterious. Based on a detailed structural analysis of 57 human CG transcript variants (CGTVs, discovered in this study) and all (833) known CGs in the human genome, we discovered that the poly(A) signal site from the upstream parent gene region is completely removed via the skipping or truncation of the final exon; consequently, CG transcription is terminated at the poly(A) signal site of the downstream parent gene. This result led us to propose a novel mechanism of CG formation: the complete removal of the poly(A) signal site from the upstream parent gene is a prerequisite for the CG transcriptional machinery to continue transcribing uninterrupted into the intergenic region and downstream parent gene. The removal of the poly(A) signal sequence from the upstream gene region appears to be caused by a deletion or truncation mutation in the human genome rather than post-transcriptional trans-splicing events. With respect to the characteristics of CG sequence structures, we found that intergenic regions are hot spots for novel exon creation during CGTV formation and that exons farther from the intergenic regions are more highly conserved in the CGTVs. Interestingly, many novel exons newly created within the intergenic and intragenic regions originated from transposable element sequences. Additionally, the CGTVs showed tumor tissue-biased expression. In conclusion, our study provides novel insights into the CG formation mechanism and expands the present concepts of the genetic structural landscape, gene regulation, and gene formation mechanisms in the human genome.