Comparison of structural variant callers for massive whole-genome sequence data.

BACKGROUND: Detecting structural variations (SVs) at the population level using next-generation sequencing (NGS) requires substantial computational resources and processing time. Here, we compared the performances of 11 SV callers: Delly, Manta, GridSS, Wham, Sniffles, Lumpy, SvABA, Canvas, CNVnator, MELT, and INSurVeyor. These SV callers have been recently published and have been widely employed for processing massive whole-genome sequencing datasets. We evaluated the accuracy, sequence depth, running time, and memory usage of the SV callers. RESULTS: Notably, several callers exhibited better calling performance for deletions than for duplications, inversions, and insertions. Among the SV callers, Manta identified deletion SVs with better performance and efficient computing resources, and both Manta and MELT demonstrated relatively good precision regarding calling insertions. We confirmed that the copy number variation callers, Canvas and CNVnator, exhibited better performance in identifying long duplications as they employ the read-depth approach. Finally, we also verified the genotypes inferred from each SV caller using a phased long-read assembly dataset, and Manta showed the highest concordance in terms of the deletions and insertions. CONCLUSIONS: Our findings provide a comprehensive understanding of the accuracy and computational efficiency of SV callers, thereby facilitating integrative analysis of SV profiles in diverse large-scale genomic datasets.

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/ .

Reduction of Squalene Epoxidase by Cholesterol Accumulation Accelerates Colorectal Cancer Progression and Metastasis.

BACKGROUND & AIMS: Squalene epoxidase (SQLE), a rate-limiting enzyme in cholesterol biosynthesis, is suggested as a proto-oncogene. Paradoxically, SQLE is degraded by excess cholesterol, and low SQLE is associated with aggressive colorectal cancer (CRC). Therefore, we studied the functional consequences of SQLE reduction in CRC progression. METHODS: Gene and protein expression data and clinical features of CRCs were obtained from public databases and 293 human tissues, analyzed by immunohistochemistry. In vitro studies showed underlying mechanisms of CRC progression mediated by SQLE reduction. Mice were fed a 2% high-cholesterol or a control diet before and after cecum implantation of SQLE genetic knockdown/control CRC cells. Metastatic dissemination and circulating cancer stem cells were demonstrated by in vivo tracking and flow cytometry analysis, respectively. RESULTS: In vitro studies showed that SQLE reduction helped cancer cells overcome constraints by inducing the epithelial-mesenchymal transition required to generate cancer stem cells. Surprisingly, SQLE interacted with GSK3ß and p53. Active GSK3ß contributes to the stability of SQLE, thereby increasing cell cholesterol content, whereas SQLE depletion disrupted the GSK3ß/p53 complex, resulting in a metastatic phenotype. This was confirmed in a spontaneous CRC metastasis mice model, where SQLE reduction, by a high-cholesterol regimen or genetic knockdown, strikingly promoted CRC aggressiveness through the production of migratory cancer stem cells. CONCLUSIONS: We showed that SQLE reduction caused by cholesterol accumulation aggravates CRC progression via the activation of the ß-catenin oncogenic pathway and deactivation of the p53 tumor suppressor pathway. Our findings provide new insights into the link between cholesterol and CRC, identifying SQLE as a key regulator in CRC aggressiveness and a prognostic biomarker.

The emerging genetic diversity of hereditary spastic paraplegia in Korean patients.

Hereditary Spastic Paraplegias (HSP) are a group of rare inherited neurological disorders characterized by progressive loss of corticospinal motor-tract function. Numerous patients with HSP remain undiagnosed despite screening for known genetic causes of HSP. Therefore, identification of novel genetic variations related to HSP is needed. In this study, we identified 88 genetic variants in 54 genes from whole-exome data of 82 clinically well-defined Korean HSP families. Fifty-six percent were known HSP genes, and 44% were composed of putative candidate HSP genes involved in the HSPome and originally reported neuron-related genes, not previously diagnosed in HSP patients. Their inheritance modes were 39, de novo; 33, autosomal dominant; and 10, autosomal recessive. Notably, ALDH18A1 showed the second highest frequency. Fourteen known HSP genes were firstly reported in Koreans, with some of their variants being predictive of HSP-causing protein malfunction. SPAST and REEP1 mutants with unknown function induced neurite abnormality. Further, 54 HSP-related genes were closely linked to the HSP progression-related network. Additionally, the genetic spectrum and variation of known HSP genes differed across ethnic groups. These results expand the genetic spectrum for HSP and may contribute to the accurate diagnosis and treatment for rare HSP.

Programming of Plant Leaf Senescence with Temporal and Inter-Organellar Coordination of Transcriptome in Arabidopsis.

Plant leaves, harvesting light energy and fixing CO2, are a major source of foods on the earth. Leaves undergo developmental and physiological shifts during their lifespan, ending with senescence and death. We characterized the key regulatory features of the leaf transcriptome during aging by analyzing total- and small-RNA transcriptomes throughout the lifespan of Arabidopsis (Arabidopsis thaliana) leaves at multidimensions, including age, RNA-type, and organelle. Intriguingly, senescing leaves showed more coordinated temporal changes in transcriptomes than growing leaves, with sophisticated regulatory networks comprising transcription factors and diverse small regulatory RNAs. The chloroplast transcriptome, but not the mitochondrial transcriptome, showed major changes during leaf aging, with a strongly shared expression pattern of nuclear transcripts encoding chloroplast-targeted proteins. Thus, unlike animal aging, leaf senescence proceeds with tight temporal and distinct interorganellar coordination of various transcriptomes that would be critical for the highly regulated degeneration and nutrient recycling contributing to plant fitness and productivity.

Global transcription network incorporating distal regulator binding reveals selective cooperation of cancer drivers and risk genes.

Global network modeling of distal regulatory interactions is essential in understanding the overall architecture of gene expression programs. Here, we developed a Bayesian probabilistic model and computational method for global causal network construction with breast cancer as a model. Whereas physical regulator binding was well supported by gene expression causality in general, distal elements in intragenic regions or loci distant from the target gene exhibited particularly strong functional effects. Modeling the action of long-range enhancers was critical in recovering true biological interactions with increased coverage and specificity overall and unraveling regulatory complexity underlying tumor subclasses and drug responses in particular. Transcriptional cancer drivers and risk genes were discovered based on the network analysis of somatic and genetic cancer-related DNA variants. Notably, we observed that the risk genes were functionally downstream of the cancer drivers and were selectively susceptible to network perturbation by tumorigenic changes in their upstream drivers. Furthermore, cancer risk alleles tended to increase the susceptibility of the transcription of their associated genes. These findings suggest that transcriptional cancer drivers selectively induce a combinatorial misregulation of downstream risk genes, and that genetic risk factors, mostly residing in distal regulatory regions, increase transcriptional susceptibility to upstream cancer-driving somatic changes.

miRseqViewer: multi-panel visualization of sequence, structure and expression for analysis of microRNA sequencing data.

SUMMARY: Deep sequencing of small RNAs has become a routine process in recent years, but no dedicated viewer is as yet available to explore the sequence features simultaneously along with secondary structure and gene expression of microRNA (miRNA). We present a highly interactive application that visualizes the sequence alignment, secondary structure and normalized read counts in synchronous multipanel windows. This helps users to easily examine the relationships between the structure of precursor and the sequences and abundance of final products and thereby will facilitate the studies on miRNA biogenesis and regulation. The project manager handles multiple samples of multiple groups. The read alignment is imported in BAM file format. Implemented features comprise sorting, zooming, highlighting, editing, filtering, saving, exporting, etc. Currently, miRseqViewer supports 84 organisms whose annotation is available at miRBase. AVAILABILITY AND IMPLEMENTATION: miRseqViewer, implemented in Java, is available at https://github.com/insoo078/mirseqviewer or at http://msv.kobic.re.kr. CONTACT: sanghyuk@ewha.ac.kr.

Exploring the DNA methylome of Korean patients with colorectal cancer consolidates the clinical implications of cancer-associated methylation markers.

Aberrant DNA methylation plays a critical role in the development and progression of colorectal cancer (CRC), which has high incidence and mortality rates in Korea. Various CRC-associated methylation markers for cancer diagnosis and prognosis have been developed; however, they have not been validated for Korean patients owing to the lack of comprehensive clinical and methylome data. Here, we obtained reliable methylation profiles for 228 tumor, 103 adjacent normal, and two unmatched normal colon tissues from Korean patients with CRC using an Illumina Infinium EPIC array; the data were corrected for biological and experiment biases. A comparative methylome analysis confirmed the previous findings that hypermethylated positions in the tumor were highly enriched in CpG island and promoter, 5' untranslated, and first exon regions. However, hypomethylated positions were enriched in the open-sea regions considerably distant from CpG islands. After applying a CpG island methylator phenotype (CIMP) to the methylome data of tumor samples to stratify the CRC patients, we consolidated the previously established clinicopathological findings that the tumors with high CIMP signatures were significantly enriched in the right colon. The results showed a higher prevalence of microsatellite instability status and MLH1 methylation in tumors with high CMP signatures than in those with low or non-CIMP signatures. Therefore, our methylome analysis and dataset provide insights into applying CRC-associated methylation markers for Korean patients regarding cancer diagnosis and prognosis. [BMB Reports 2024; 57(3): 161-166].

Comprehensive RNA-sequencing analysis of colorectal cancer in a Korean cohort.

Considering the recent increase in the number of colorectal cancer (CRC) cases in South Korea, we aimed to clarify the molecular characteristics of CRC unique to the Korean population. To gain insights into the complexities of CRC and promote the exchange of critical data, RNA-sequencing analysis was performed to reveal the molecular mechanisms that drive the development and progression of CRC; this analysis is critical for developing effective treatment strategies. We performed RNA-sequencing analysis of CRC and adjacent normal tissue samples from 214 Korean participants (comprising a total of 381 including 169 normal and 212 tumor samples) to investigate differential gene expression between the groups. We identified 19,575 genes expressed in CRC and normal tissues, with 3,830 differentially expressed genes (DEGs) between the groups. Functional annotation analysis revealed that the upregulated DEGs were significantly enriched in pathways related to the cell cycle, DNA replication, and IL-17, whereas the downregulated DEGs were enriched in metabolic pathways. We also analyzed the relationship between clinical information and subtypes using the Consensus Molecular Subtype (CMS) classification. Furthermore, we compared groups clustered within our dataset to CMS groups and performed additional analysis of the methylation data between DEGs and CMS groups to provide comprehensive biological insights from various perspectives. Our study provides valuable insights into the molecular mechanisms underlying CRC in Korean patients and serves as a platform for identifying potential target genes for this disease. The raw data and processed results have been deposited in a public repository for further analysis and exploration.

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/.

Comprehensive profiling of DNA methylation in Korean patients with colorectal cancer.

Alterations in DNA methylation play an important pathophysiological role in the development and progression of colorectal cancer. We comprehensively profiled DNA methylation alterations in 165 Korean patients with colorectal cancer (CRC), and conducted an in-depth investigation of cancer-specific methylation patterns. Our analysis of the tumor samples revealed a significant presence of hypomethylated probes, primarily within the gene body regions; few hypermethylated sites were observed, which were mostly enriched in promoter-like and CpG island regions. The CpG Island Methylator PhenotypeHigh (CIMP-H) exhibited notable enrichment of microsatellite instability-high (MSI-H). Additionally, our findings indicated a significant correlation between methylation of the MLH1 gene and MSI-H status. Furthermore, we found that the CIMP-H had a higher tendency to affect the right-side of the colon tissues and was slightly more prevalent among older patients. Through our methylome profile analysis, we successfully verified the thylation patterns and clinical characteristics of Korean patients with CRC. This valuable dataset lays a strong foundation for exploring novel molecular insights and potential therapeutic targets for the treatment of CRC. [BMB Reports 2024; 57(2): 110-115].

VnD: a structure-centric database of disease-related SNPs and drugs.

Numerous genetic variations have been found to be related to human diseases. Significant portion of those affect the drug response as well by changing the protein structure and function. Therefore, it is crucial to understand the trilateral relationship among genomic variations, diseases and drugs. We present the variations and drugs (VnD), a consolidated database containing information on diseases, related genes and genetic variations, protein structures and drug information. VnD was built in three steps. First, we integrated various resources systematically to deduce catalogs of disease-related genes, single nucleotide polymorphisms (SNPs), protein mutations and relevant drugs. VnD contains 137,195 disease-related gene records (13,940 distinct genes) and 16,586 genetic variation records (1790 distinct variations). Next, we carried out structure modeling and docking simulation for wild-type and mutant proteins to examine the structural and functional consequences of non-synonymous SNPs in the drug-related genes. Conformational changes in 590 wild-type and 4437 mutant proteins from drug-related genes were included in our database. Finally, we investigated the structural and biochemical properties relevant to drug binding such as the distribution of SNPs in proximal protein pockets, thermo-chemical stability, interactions with drugs and physico-chemical properties. The VnD database, available at http://vnd.kobic.re.kr:8080/VnD/ or vandd.org, would be a useful platform for researchers studying the underlying mechanism for association among genetic variations, diseases and drugs.

Accurate quantification of transcriptome from RNA-Seq data by effective length normalization.

We propose a novel, efficient and intuitive approach of estimating mRNA abundances from the whole transcriptome shotgun sequencing (RNA-Seq) data. Our method, NEUMA (Normalization by Expected Uniquely Mappable Area), is based on effective length normalization using uniquely mappable areas of gene and mRNA isoform models. Using the known transcriptome sequence model such as RefSeq, NEUMA pre-computes the numbers of all possible gene-wise and isoform-wise informative reads: the former being sequences mapped to all mRNA isoforms of a single gene exclusively and the latter uniquely mapped to a single mRNA isoform. The results are used to estimate the effective length of genes and transcripts, taking experimental distributions of fragment size into consideration. Quantitative RT-PCR based on 27 randomly selected genes in two human cell lines and computer simulation experiments demonstrated superior accuracy of NEUMA over other recently developed methods. NEUMA covers a large proportion of genes and mRNA isoforms and offers a measure of consistency ('consistency coefficient') for each gene between an independently measured gene-wise level and the sum of the isoform levels. NEUMA is applicable to both paired-end and single-end RNA-Seq data. We propose that NEUMA could make a standard method in quantifying gene transcript levels from RNA-Seq data.

MHC II immunogenicity shapes the neoepitope landscape in human tumors.

Despite advances in predicting physical peptide-major histocompatibility complex I (pMHC I) binding, it remains challenging to identify functionally immunogenic neoepitopes, especially for MHC II. By using the results of >36,000 immunogenicity assay, we developed a method to identify pMHC whose structural alignment facilitates T cell reaction. Our method predicted neoepitopes for MHC II and MHC I that were responsive to checkpoint blockade when applied to >1,200 samples of various tumor types. To investigate selection by spontaneous immunity at the single epitope level, we analyzed the frequency spectrum of >25 million mutations in >9,000 treatment-naive tumors with >100 immune phenotypes. MHC II immunogenicity specifically lowered variant frequencies in tumors under high immune pressure, particularly with high TCR clonality and MHC II expression. A similar trend was shown for MHC I neoepitopes, but only in particular tissue types. In summary, we report immune selection imposed by MHC II-restricted natural or therapeutic T cell reactivity.

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.

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

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

A Novel X-Linked Variant of IQSEC2 is Associated with Lennox-Gastaut Syndrome and Mild Intellectual Disability in Three Generations of a Korean Family.

Aim: Lennox-Gastaut syndrome (LGS) is a severe type of childhood-onset epilepsy with multiple types of seizures, specific discharges on electroencephalography, and intellectual disability. However, LGS-related genes are largely unknown. To identify causative genes related to LGS, we collected and analyzed data from a three-generation Korean family in which one member had LGS and two had intellectual disability. Methods: Genomic DNAs were extracted from blood samples of all participants and used in whole-exome sequencing (WES). Genetic variants were detected by the Genome Analysis Toolkit and confirmed by Sanger sequencing. Variant pathogenicity was evaluated by prediction programs and the American College of Medical Genetics criteria. The LGS patient had generalized slow spike-and-wave discharges, multiple types of seizures, and developmental delay. Results: Analyses of the WES data from the family revealed a novel variant (c.1048G>A, p.Ala350Thr) in the IQ motif and Sec7 domain 2 (IQSEC2). This variant is within a highly evolutionarily conserved IQ-like motif, indicating a decrease in the calmodulin-binding capacity or α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid transmission. The hemizygous variant in the male with LGS was a maternally inherited X-linked variant from the heterozygous maternal grandmother and mother, both of whom had intellectual disability. Conclusion: These findings indicate that the variant of IQSEC2 triggered both LGS and intellectual disability dependent on sex in this family. We report a novel X-linked inherited IQSEC2 variant for LGS and intellectual disability, which enhances the spectrum of variants in the IQ-like motif of IQSEC2.

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.

ssSNPTarget: genome-wide splice-site Single Nucleotide Polymorphism database.

Deep sequencing has shown that over 90% of human genes undergo alternative splicing. The splicing process requires exon-intron boundary recognition. SNPs located in the boundaries (splice sites) influence exon configuration. Also, splice site SNPs (ssSNPs) alter translation efficiency of the mRNA and lead to important changes in disease susceptibility. We developed the ssSNPTarget database to provide ssSNPs on human and mouse genes. It includes: 1) ssSNP distribution information in human and mouse genes; 2) effects of SNPs in splice sites: junction strength change, protein domain change, and alternative splicing events (exon skipping, 5'- or 3'-exon extension); 3) splice site conservation in eukaryotes; and 4) associated disease information derived from OMIM, GAD, and HGMD. ssSNPTarget contains 1,576 human ssSNPs associated with 1,193 genes and 538 mouse ssSNPs associated with 281 genes. Users can query ssSNPTarget with several types of search terms (gene symbol, SNP rs number, transcript ID, or genomic position), and the information can be accessed at http://variome.kobic.re.kr/ssSNPTarget/ or http://ssSNPTarget.org.

PDbase: a database of Parkinson's disease-related genes and genetic variation using substantia nigra ESTs.

BACKGROUND: Parkinson's disease (PD) is one of the most common neurodegenerative disorders, clinically characterized by impaired motor function. Since the etiology of PD is diverse and complex, many researchers have created PD-related research resources. However, resources for brain and PD studies are still lacking. Therefore, we have constructed a database of PD-related gene and genetic variations using the substantia nigra (SN) in PD and normal tissues. In addition, we integrated PD-related information from several resources. RESULTS: We collected the 6,130 SN expressed sequenced tags (ESTs) from brain SN normal tissues and PD patients SN tissues using full-cDNA library and normalized cDNA library construction methods from our previous study. The SN ESTs were clustered in 2,951 unigene clusters and assigned in 2,678 genes. We then found up-regulated 57 genes and down-regulated 48 genes by comparing normal and PD SN ESTs frequencies with over 0.9 cut-off probability of differential expression based on the Audic and Claverie method. In addition, we integrated disease-related information from public resources. To examine the characteristics of these PD-related genes, we analyzed alternative splicing events, single nucleotide polymorphism (SNP) markers located in the gene regions, repeat elements, gene regulation elements, and pathways and protein-protein interaction networks. CONCLUSION: We constructed the PDbase database to capture the PD-related gene, genetic variation, and functional elements. This database contains 2,698 PD-related genes through ESTs discovered from human normal and PD patients SN tissues, and through integrating several public resources. PDbase provides the mitochondrion proteins, microRNA gene regulation elements, single nucleotide polymorphisms (SNPs) markers within PD-related gene structures, repeat elements, and pathways and networks with protein-protein interaction information. The PDbase information can aid in understanding the causation of PD. It is available at http://bioportal.kobic.re.kr/PDbase/. Supplementary data is available at http://bioportal.kobic.re.kr/PDbase/suppl.jsp.