KoNA: Korean Nucleotide Archive as A New Data Repository for Nucleotide Sequence Data.

During the last decade, the generation and accumulation of petabase-scale high-throughput sequencing data have resulted in great challenges, including access to human data, as well as transfer, storage, and sharing of enormous amounts of data. To promote data-driven biological research, the Korean government announced that all biological data generated from government-funded research projects should be deposited at the Korea BioData Station (K-BDS), which consists of multiple databases for individual data types. Here, we introduce the Korean Nucleotide Archive (KoNA), a repository of nucleotide sequence data. As of July 2022, the Korean Read Archive in KoNA has collected over 477 TB of raw next-generation sequencing data from national genome projects. To ensure data quality and prepare for international alignment, a standard operating procedure was adopted, which is similar to that of the International Nucleotide Sequence Database Collaboration. The standard operating procedure includes quality control processes for submitted data and metadata using an automated pipeline, followed by manual examination. To ensure fast and stable data transfer, a high-speed transmission system called GBox is used in KoNA. Furthermore, the data uploaded to or downloaded from KoNA through GBox can be readily processed using a cloud computing service called Bio-Express. This seamless coupling of KoNA, GBox, and Bio-Express enhances the data experience, including submission, access, and analysis of raw nucleotide sequences. KoNA not only satisfies the unmet needs for a national sequence repository in Korea but also provides datasets to researchers globally and contributes to advances in genomics. The KoNA is available at https://www.kobic.re.kr/kona/.

Visualizing a Cold Stress-Specific Pulse Wave in Traditional Pulse Diagnosis ('Tight Pulse') Correlated with Vascular Changes in the Radial Artery Induced by a Cold Pressor Trial.

Radial pulse diagnosis is the most common method to examine the human health state in Traditional East Asian Medicine (TEAM). A cold stress-related suboptimal health state (subhealth) is often undetectable during routine medical examinations, however, it can be detected through the palpation of specific pulse waves, particularly a 'tight pulse', in TEAM. Therefore, this study examined a correlation between 'tight pulse' and vascular changes in the radial artery (RA) induced by a cold pressor trial (CPT). Twenty healthy subjects underwent sequentially control trial and CPT with room-temperature and ice-cold water, respectively, on the right forearm. The radial pulse and vascular changes were then examined on the left arm. The radial pulse scores for frequencies of 'tight pulse' with strong arterial tension increased after the CPT compared with the control trial. The pulse scores were reversely correlated with the RA thickness and volumes in ultrasonography, but not with changes in the systolic/diastolic blood pressure. The RA thickness-based vascular surface and three-dimensional images visualized a 'tight pulse' showing the vasoconstriction and bumpy-/rope-shaped vascular changes in the radial pulse diagnostic region after the CPT. These findings provide valuable insights into the potential integration of clinical radial pulse diagnosis with ultrasonography for cold-related subhealth.

CanISO: a database of genomic and transcriptomic variations in domestic dog (Canis lupus familiaris).

BACKGROUND: The domestic dog, Canis lupus familiaris, is a companion animal for humans as well as an animal model in cancer research due to similar spontaneous occurrence of cancers as humans. Despite the social and biological importance of dogs, the catalogue of genomic variations and transcripts for dogs is relatively incomplete. RESULTS: We developed CanISO, a new database to hold a large collection of transcriptome profiles and genomic variations for domestic dogs. CanISO provides 87,692 novel transcript isoforms and 60,992 known isoforms from whole transcriptome sequencing of canine tumors (N = 157) and their matched normal tissues (N = 64). CanISO also provides genomic variation information for 210,444 unique germline single nucleotide polymorphisms (SNPs) from the whole exome sequencing of 183 dogs, with a query system that searches gene- and transcript-level information as well as covered SNPs. Transcriptome profiles can be compared with corresponding human transcript isoforms at a tissue level, or between sample groups to identify tumor-specific gene expression and alternative splicing patterns. CONCLUSIONS: CanISO is expected to increase understanding of the dog genome and transcriptome, as well as its functional associations with humans, such as shared/distinct mechanisms of cancer. CanISO is publicly available at https://www.kobic.re.kr/caniso/ .

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.

A database of 5305 healthy Korean individuals reveals genetic and clinical implications for an East Asian population.

Despite substantial advances in disease genetics, studies to date have largely focused on individuals of European descent. This limits further discoveries of novel functional genetic variants in other ethnic groups. To alleviate the paucity of East Asian population genome resources, we established the Korean Variant Archive 2 (KOVA 2), which is composed of 1896 whole-genome sequences and 3409 whole-exome sequences from healthy individuals of Korean ethnicity. This is the largest genome database from the ethnic Korean population to date, surpassing the 1909 Korean individuals deposited in gnomAD. The variants in KOVA 2 displayed all the known genetic features of those from previous genome databases, and we compiled data from Korean-specific runs of homozygosity, positively selected intervals, and structural variants. In doing so, we found loci, such as the loci of ADH1A/1B and UHRF1BP1, that are strongly selected in the Korean population relative to other East Asian populations. Our analysis of allele ages revealed a correlation between variant functionality and evolutionary age. The data can be browsed and downloaded from a public website ( https://www.kobic.re.kr/kova/ ). We anticipate that KOVA 2 will serve as a valuable resource for genetic studies involving East Asian populations.

Temperature Characteristics of Traditional Indirect Moxibustion and Electronic Moxibustion.

Background: Electronic moxibustion (EM) was developed to minimize the side effects of traditional moxibustion, such as burns, and to overcome therapeutic compliances such as smoke or smell. Objectives: To investigate distributions and thermal stimulation of EM at various depths using silicon phantom and to compare this methodology to traditional indirect moxibustion (TIM). Methods: A silicon phantom composed of polydimethylsiloxane was heated and immersed in a hot plate containing warm water to set the phantom's temperature to that of biological tissue. K-type thermocouples were inserted into the phantom at depths of 0, 2, 5, 7, and 10 mm to measure temperature changes with thermal stimulation of EM or TIM placed on top of the phantom. Results: At the surface of the phantom, the peak temperature after applying TIM (55.04 ± 0.92℃ [Δ23.79 ± 0.96℃]) was significantly higher than after EM (43.25 ± 1.95℃ [Δ13.00 ± 2.23℃]), with both interventions reaching the highest temperature after 2 minutes. The temperature increase for TIM was also statistically significant compared to EM when measured at a depth of 2 mm. For the experimental setting with TIM, after reaching peak surface temperature, a rapid decrease was observed at the surface and 2 mm while EM showed a much more gradual decline. There was no significant difference in temperature change between the groups at depths of 5, 7, and 10 mm. Conclusion: TIM resulted in a higher temperature rise compared to EM at the surface and at a 2 mm depth reaching over 50℃, which creates risk of burns. Thermal stimulation with EM had a lower risk of burns with temperature increment not being statistically different from TIM below the depth of 5 mm.

Thermal Properties of Warm- versus Heated-Needle Acupuncture.

Background: Warm-needle acupuncture (WA) and fire-needle acupuncture are treatment techniques that use the combination of acupuncture and thermal stimulation. In clinical practice, a new method of fire-needle acupuncture called "heated-needle acupuncture (HA)" has been proposed, wherein the needle is directly heated after insertion. WA and HA share similarities in their methods, and no previous study has sought to assess whether their thermal outcomes are also similar. Methods: We controlled environmental variables and measured the maximum temperatures and temperature changes of a silicon phantom in which K-type thermocouples were embedded at depths of 0, 2, 5, 7, and 10 mm. WA and HA were also performed with acupuncture needles of various thicknesses (0.30 × 40 mm, 0.40 × 40 mm, and 0.50 × 40 mm). Results: Different time-dependent temperature distributions were observed between the two acupuncture methods: HA yielded a higher maximum temperature and temperature change on the surface, whereas WA yielded higher temperatures at the other tested depths. The thermal patterns were similar among the needles of different thicknesses for each method, with the following exception: while the temperature change and maximum temperature did not differ significantly by needle thickness for WA, these parameters increased significantly with needle thickness for HA. Conclusion: The two acupuncture procedures yielded different thermal patterns in a controlled environment. Further studies are necessary to reflect the effect of external environment variables occurring in reality.

Development of a quantitative prediction algorithm for target organ-specific similarity of human pluripotent stem cell-derived organoids and cells.

Human pluripotent stem cell (hPSC)-derived organoids and cells have similar characteristics to human organs and tissues. Thus, in vitro human organoids and cells serve as a superior alternative to conventional cell lines and animal models in drug development and regenerative medicine. For a simple and reproducible analysis of the quality of organoids and cells to compensate for the shortcomings of existing experimental validation studies, a quantitative evaluation method should be developed. Here, using the GTEx database, we construct a quantitative calculation system to assess similarity to the human organs. To evaluate our system, we generate hPSC-derived organoids and cells, and detected organ similarity. To facilitate the access of our system by researchers, we develop a web-based user interface presenting similarity to the appropriate organs as percentages. Thus, this program could provide valuable information for the generation of high-quality organoids and cells and a strategy to guide proper lineage-oriented differentiation.

Corrigendum: Characteristics of Genetic Variations Associated With Lennox-Gastaut Syndrome in Korean Families.

[This corrects the article DOI: 10.3389/fgene.2020.590924.].

3DIV update for 2021: a comprehensive resource of 3D genome and 3D cancer genome.

Three-dimensional (3D) genome organization is tightly coupled with gene regulation in various biological processes and diseases. In cancer, various types of large-scale genomic rearrangements can disrupt the 3D genome, leading to oncogenic gene expression. However, unraveling the pathogenicity of the 3D cancer genome remains a challenge since closer examinations have been greatly limited due to the lack of appropriate tools specialized for disorganized higher-order chromatin structure. Here, we updated a 3D-genome Interaction Viewer and database named 3DIV by uniformly processing ∼230 billion raw Hi-C reads to expand our contents to the 3D cancer genome. The updates of 3DIV are listed as follows: (i) the collection of 401 samples including 220 cancer cell line/tumor Hi-C data, 153 normal cell line/tissue Hi-C data, and 28 promoter capture Hi-C data, (ii) the live interactive manipulation of the 3D cancer genome to simulate the impact of structural variations and (iii) the reconstruction of Hi-C contact maps by user-defined chromosome order to investigate the 3D genome of the complex genomic rearrangement. In summary, the updated 3DIV will be the most comprehensive resource to explore the gene regulatory effects of both the normal and cancer 3D genome. '3DIV' is freely available at http://3div.kr.

iCSDB: an integrated database of CRISPR screens.

High-throughput screening based on CRISPR-Cas9 libraries has become an attractive and powerful technique to identify target genes for functional studies. However, accessibility of public data is limited due to the lack of user-friendly utilities and up-to-date resources covering experiments from third parties. Here, we describe iCSDB, an integrated database of CRISPR screening experiments using human cell lines. We compiled two major sources of CRISPR-Cas9 screening: the DepMap portal and BioGRID ORCS. DepMap portal itself is an integrated database that includes three large-scale projects of CRISPR screening. We additionally aggregated CRISPR screens from BioGRID ORCS that is a collection of screening results from PubMed articles. Currently, iCSDB contains 1375 genome-wide screens across 976 human cell lines, covering 28 tissues and 70 cancer types. Importantly, the batch effects from different CRISPR libraries were removed and the screening scores were converted into a single metric to estimate the knockout efficiency. Clinical and molecular information were also integrated to help users to select cell lines of interest readily. Furthermore, we have implemented various interactive tools and viewers to facilitate users to choose, examine and compare the screen results both at the gene and guide RNA levels. iCSDB is available at https://www.kobic.re.kr/icsdb/.

Bioinformatics services for analyzing massive genomic datasets.

The explosive growth of next-generation sequencing data has resulted in ultra-large-scale datasets and ensuing computational problems. In Korea, the amount of genomic data has been increasing rapidly in the recent years. Leveraging these big data requires researchers to use large-scale computational resources and analysis pipelines. A promising solution for addressing this computational challenge is cloud computing, where CPUs, memory, storage, and programs are accessible in the form of virtual machines. Here, we present a cloud computing-based system, Bio-Express, that provides user-friendly, cost-effective analysis of massive genomic datasets. Bio-Express is loaded with predefined multi-omics data analysis pipelines, which are divided into genome, transcriptome, epigenome, and metagenome pipelines. Users can employ predefined pipelines or create a new pipeline for analyzing their own omics data. We also developed several web-based services for facilitating downstream analysis of genome data. Bio-Express web service is freely available at https://www.bioexpress.re.kr/.

Characteristics of Genetic Variations Associated With Lennox-Gastaut Syndrome in Korean Families.

Lennox-Gastaut syndrome (LGS) is a severe type of childhood-onset epilepsy characterized by multiple types of seizures, specific discharges on electroencephalography, and intellectual disability. Most patients with LGS do not respond well to drug treatment and show poor long-term prognosis. Approximately 30% of patients without brain abnormalities have unidentifiable causes. Therefore, accurate diagnosis and treatment of LGS remain challenging. To identify causative mutations of LGS, we analyzed the whole-exome sequencing data of 17 unrelated Korean families, including patients with LGS and LGS-like epilepsy without brain abnormalities, using the Genome Analysis Toolkit. We identified 14 mutations in 14 genes as causes of LGS or LGS-like epilepsy. 64 percent of the identified genes were reported as LGS or epilepsy-related genes. Many of these variations were novel and considered as pathogenic or likely pathogenic. Network analysis was performed to classify the identified genes into two network clusters: neuronal signal transmission or neuronal development. Additionally, knockdown of two candidate genes with insufficient evidence of neuronal functions, SLC25A39 and TBC1D8, decreased neurite outgrowth and the expression level of MAP2, a neuronal marker. These results expand the spectrum of genetic variations and may aid the diagnosis and management of individuals with LGS.

ChimerDB 4.0: an updated and expanded database of fusion genes.

Fusion genes represent an important class of biomarkers and therapeutic targets in cancer. ChimerDB is a comprehensive database of fusion genes encompassing analysis of deep sequencing data (ChimerSeq) and text mining of publications (ChimerPub) with extensive manual annotations (ChimerKB). In this update, we present all three modules substantially enhanced by incorporating the recent flood of deep sequencing data and related publications. ChimerSeq now covers all 10 565 patients in the TCGA project, with compilation of computational results from two reliable programs of STAR-Fusion and FusionScan with several public resources. In sum, ChimerSeq includes 65 945 fusion candidates, 21 106 of which were predicted by multiple programs (ChimerSeq-Plus). ChimerPub has been upgraded by applying a deep learning method for text mining followed by extensive manual curation, which yielded 1257 fusion genes including 777 cases with experimental supports (ChimerPub-Plus). ChimerKB includes 1597 fusion genes with publication support, experimental evidences and breakpoint information. Importantly, we implemented several new features to aid estimation of functional significance, including the fusion structure viewer with domain information, gene expression plot of fusion positive versus negative patients and a STRING network viewer. The user interface also was greatly enhanced by applying responsive web design. ChimerDB 4.0 is available at http://www.kobic.re.kr/chimerdb/.

BiomeNet: a database for construction and analysis of functional interaction networks for any species with a sequenced genome.

MOTIVATION: Owing to advanced DNA sequencing and genome assembly technology, the number of species with sequenced genomes is rapidly increasing. The aim of the recently launched Earth BioGenome Project is to sequence genomes of all eukaryotic species on Earth over the next 10 years, making it feasible to obtain genomic blueprints of the majority of animal and plant species by this time. Genetic models of the sequenced species will later be subject to functional annotation, and a comprehensive molecular network should facilitate functional analysis of individual genes and pathways. However, network databases are lagging behind genome sequencing projects as even the largest network database provides gene networks for less than 10% of sequenced eukaryotic genomes, and the knowledge gap between genomes and interactomes continues to widen. RESULTS: We present BiomeNet, a database of 95 scored networks comprising over 8 million co-functional links, which can build and analyze gene networks for any species with the sequenced genome. BiomeNet transfers functional interactions between orthologous proteins from source networks to the target species within minutes and automatically constructs gene networks with the quality comparable to that of existing networks. BiomeNet enables assembly of the first-in-species gene networks not available through other databases, which are highly predictive of diverse biological processes and can also provide network analysis by extracting subnetworks for individual biological processes and network-based gene prioritizations. These data indicate that BiomeNet could enhance the benefits of decoding the genomes of various species, thus improving our understanding of the Earth' biodiversity. AVAILABILITY AND IMPLEMENTATION: The BiomeNet is freely available at http://kobic.re.kr/biomenet/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CaPSSA: visual evaluation of cancer biomarker genes for patient stratification and survival analysis using mutation and expression data.

SUMMARY: Predictive biomarkers for patient stratification play critical roles in realizing the paradigm of precision medicine. Molecular characteristics such as somatic mutations and expression signatures represent the primary source of putative biomarker genes for patient stratification. However, evaluation of such candidate biomarkers is still cumbersome and requires multistep procedures especially when using massive public omics data. Here, we present an interactive web application that divides patients from large cohorts (e.g. The Cancer Genome Atlas, TCGA) dynamically into two groups according to the mutation, copy number variation or gene expression of query genes. It further supports users to examine the prognostic value of resulting patient groups based on survival analysis and their association with the clinical features as well as the previously annotated molecular subtypes, facilitated with a rich and interactive visualization. Importantly, we also support custom omics data with clinical information. AVAILABILITY AND IMPLEMENTATION: CaPSSA (Cancer Patient Stratification and Survival Analysis) runs on a web-browser and is freely available without restrictions at http://www.kobic.re.kr/capssa/. The source code is available on https://github.com/yjjang/capssa. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

GEMiCCL: mining genotype and expression data of cancer cell lines with elaborate visualization.

Database URL: GEMiCCL is available at https://www.kobic.kr/GEMICCL/.

Closha: bioinformatics workflow system for the analysis of massive sequencing data.

BACKGROUND: While next-generation sequencing (NGS) costs have fallen in recent years, the cost and complexity of computation remain substantial obstacles to the use of NGS in bio-medical care and genomic research. The rapidly increasing amounts of data available from the new high-throughput methods have made data processing infeasible without automated pipelines. The integration of data and analytic resources into workflow systems provides a solution to the problem by simplifying the task of data analysis. RESULTS: To address this challenge, we developed a cloud-based workflow management system, Closha, to provide fast and cost-effective analysis of massive genomic data. We implemented complex workflows making optimal use of high-performance computing clusters. Closha allows users to create multi-step analyses using drag and drop functionality and to modify the parameters of pipeline tools. Users can also import the Galaxy pipelines into Closha. Closha is a hybrid system that enables users to use both analysis programs providing traditional tools and MapReduce-based big data analysis programs simultaneously in a single pipeline. Thus, the execution of analytics algorithms can be parallelized, speeding up the whole process. We also developed a high-speed data transmission solution, KoDS, to transmit a large amount of data at a fast rate. KoDS has a file transfer speed of up to 10 times that of normal FTP and HTTP. The computer hardware for Closha is 660 CPU cores and 800 TB of disk storage, enabling 500 jobs to run at the same time. CONCLUSIONS: Closha is a scalable, cost-effective, and publicly available web service for large-scale genomic data analysis. Closha supports the reliable and highly scalable execution of sequencing analysis workflows in a fully automated manner. Closha provides a user-friendly interface to all genomic scientists to try to derive accurate results from NGS platform data. The Closha cloud server is freely available for use from http://closha.kobic.re.kr/ .

3DIV: A 3D-genome Interaction Viewer and database.

Three-dimensional (3D) chromatin structure is an emerging paradigm for understanding gene regulation mechanisms. Hi-C (high-throughput chromatin conformation capture), a method to detect long-range chromatin interactions, allows extensive genome-wide investigation of 3D chromatin structure. However, broad application of Hi-C data have been hindered by the level of complexity in processing Hi-C data and the large size of raw sequencing data. In order to overcome these limitations, we constructed a database named 3DIV (a 3D-genome Interaction Viewer and database) that provides a list of long-range chromatin interaction partners for the queried locus with genomic and epigenomic annotations. 3DIV is the first of its kind to collect all publicly available human Hi-C data to provide 66 billion uniformly processed raw Hi-C read pairs obtained from 80 different human cell/tissue types. In contrast to other databases, 3DIV uniquely provides normalized chromatin interaction frequencies against genomic distance dependent background signals and a dynamic browsing visualization tool for the listed interactions, which could greatly advance the interpretation of chromatin interactions. '3DIV' is available at http://kobic.kr/3div.

Analysis of Herbal Medicine Prescriptions for Patients in An Academic Korean Medical Hospital: A Cross Sectional Study of Electronic Medical Records (2010-2013).

OBJECTIVE: To obtain fundamental information for the standardization of herbal medicine in Korea. METHODS: We analyzed the herbal medicine prescription data of patients at the Pusan National University Korean Medicine Hospital from March 2010 to February 2013. We used the Dongui-Bogam (Dong Yi Bao Jian) to classify prescribed herbal medicines. RESULTS: The study revealed that the most frequently prescribed herbal medicine was 'Liuwei Dihuang Pill (LWDHP, )' which was used for invigorating 'Shen (Kidndy)-yin'. 'LWDHP' was most frequently prescribed to male patients aged 50-59, 60-69, 70-79 and 80-89 years, and 'Xionggui Tiaoxue Decoction (XGTXD, )' was most frequently prescribed to female patients aged 30-39 and 40-49 years. According to the International Classification of Diseases (ICD) codes, 'Diseases of the musculoskeletal system and connective tissue' showed the highest prevalence. 'LWDHP' and 'XGTXD' was the most frequently prescribed in categories 5 and 3, respectively. Based on the percentage of prescriptions for each sex, 'Ziyin Jianghuo Decoction ()' was prescribed to mainly male patients, and 'XGTXD' with 'Guima Geban Decoction ()' were prescribed to mainly female patients. CONCLUSION: This study analysis successfully determined the frequency of a variety of herbal medicines, and many restorative herbal medicines were identified and frequently administered.