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

Comparative analysis of commonly used peak calling programs for ChIP-Seq analysis

Chromatin immunoprecipitation coupled with high-throughput DNA sequencing (ChIP-Seq) is a powerful technology to profile the location of proteins of interest on a whole-genome scale. To identify the enrichment location of proteins, many programs and algorithms have been proposed. However, none of the commonly used peak calling programs could accurately explain the binding features of target proteins detected by ChIP-Seq. Here, publicly available data on 12 histone modifications, including H3K4ac/me1/me2/me3, H3K9ac/me3, H3K27ac/me3, H3K36me3, H3K56ac, and H3K79me1/me2, generated from a human embryonic stem cell line (H1), were profiled with five peak callers (CisGenome, MACS1, MACS2, PeakSeq, and SISSRs). The performance of the peak calling programs was compared in terms of reproducibility between replicates, examination of enriched regions to variable sequencing depths, the specificity-to-noise signal, and sensitivity of peak prediction. There were no major differences among peak callers when analyzing point source histone modifications. The peak calling results from histone modifications with low fidelity, such as H3K4ac, H3K56ac, and H3K79me1/me2, showed low performance in all parameters, which indicates that their peak positions might not be located accurately. Our comparative results could provide a helpful guide to choose a suitable peak calling program for specific histone modifications.

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

Comparative analysis of commonly used peak calling programs for ChIP-Seq analysis

Chromatin immunoprecipitation coupled with high-throughput DNA sequencing (ChIP-Seq) is a powerful technology to profile the location of proteins of interest on a whole-genome scale. To identify the enrichment location of proteins, many programs and algorithms have been proposed. However, none of the commonly used peak calling programs could accurately explain the binding features of target proteins detected by ChIP-Seq. Here, publicly available data on 12 histone modifications, including H3K4ac/me1/me2/me3, H3K9ac/me3, H3K27ac/me3, H3K36me3, H3K56ac, and H3K79me1/me2, generated from a human embryonic stem cell line (H1), were profiled with five peak callers (CisGenome, MACS1, MACS2, PeakSeq, and SISSRs). The performance of the peak calling programs was compared in terms of reproducibility between replicates, examination of enriched regions to variable sequencing depths, the specificity-to-noise signal, and sensitivity of peak prediction. There were no major differences among peak callers when analyzing point source histone modifications. The peak calling results from histone modifications with low fidelity, such as H3K4ac, H3K56ac, and H3K79me1/me2, showed low performance in all parameters, which indicates that their peak positions might not be located accurately. Our comparative results could provide a helpful guide to choose a suitable peak calling program for specific histone modifications.

Evaluation and Interpretation of Transcriptome Data Underlying Heterogeneous Chronic Obstructive Pulmonary Disease

Chronic obstructive pulmonary disease (COPD) is a type of progressive lung disease, featured by airflow obstruction. Recently, a comprehensive analysis of the transcriptome in lung tissue of COPD patients was performed, but the heterogeneity of the sample was not seriously considered in characterizing the mechanistic dysregulation of COPD. Here, we established a new transcriptome analysis pipeline using a deconvolution process to reduce the heterogeneity and clearly identified that these transcriptome data originated from the mild or moderate stage of COPD patients. Differentially expressed or co-expressed genes in the protein interaction subnetworks were linked with mitochondrial dysfunction and the immune response, as expected. Computational protein localization prediction revealed that 19 proteins showing changes in subcellular localization were mostly related to mitochondria, suggesting that mislocalization of mitochondria-targeting proteins plays an important role in COPD pathology. Our extensive evaluation of COPD transcriptome data could provide guidelines for analyzing heterogeneous gene expression profiles and classifying potential candidate genes that are responsible for the pathogenesis of COPD.

A Follow-up Association Study of Genetic Variants for Bone Mineral Density in a Korean Population

Bone mineral density (BMD) is one of the quantitative traits that are genetically inherited and affected by various factors. Over the past years, genome-wide association studies (GWASs) have searched for many genetic loci that influence BMD. A recent meta-analysis of 17 GWASs for BMD of the femoral neck and lumbar spine is the largest GWAS for BMD to date and offers 64 single-nucleotide polymorphisms (SNPs) in 56 associated loci. We investigated these BMD loci in a Korean population called Korea Association REsource (KARE) to identify their validity in an independent study. The KARE population contains genotypes from 8,842 individuals, and their BMD levels were measured at the distal radius (BMD-RT) and midshaft tibia (BMD-TT). Thirteen genomic loci among 56 loci were significantly associated with BMD variations, and 3 loci were involved in known biological pathways related to BMD. In order to find putative functional variants, nearby SNPs in relation to linkage equilibrium were annotated, and their possible functional effects were predicted. These findings reveal that tens of variants, not a single factor, may contribute to the genetic architecture of BMD; have an important role regardless of ethnic group; and may highlight the importance of a replication study in GWASs to validate genuine loci for BMD variation.

Characterization of Chromatin Structure-associated Histone Modifications in Breast Cancer Cells

Chromatin structure and dynamics that are influenced by epigenetic marks, such as histone modification and DNA methylation, play a crucial role in modulating gene transcription. To understand the relationship between histone modifications and regulatory elements in breast cancer cells, we compared our chromatin immunoprecipitation sequencing (ChIP-Seq) histone modification patterns for histone H3K4me1, H3K4me3, H3K9/16ac, and H3K27me3 in MCF-7 cells with publicly available formaldehyde-assisted isolation of regulatory elements (FAIRE)-chip signals in human chromosomes 8, 11, and 12, identified by a method called FAIRE. Active regulatory elements defined by FAIRE were highly associated with active histone modifications, like H3K4me3 and H3K9/16ac, especially near transcription start sites. The H3K9/16ac-enriched genes that overlapped with FAIRE signals (FAIRE-H3K9/14ac) were moderately correlated with gene expression levels. We also identified functional sequence motifs at H3K4me1-enriched FAIRE sites upstream of putative promoters, suggesting that regulatory elements could be associated with H3K4me1 to be regarded as distal regulatory elements. Our results might provide an insight into epigenetic regulatory mechanisms explaining the association of histone modifications with open chromatin structure in breast cancer cells.