• The KIPS Transactions:PartD
• /
• v.16D no.3
• /
• pp.327-338
• /
• 2009

A microarray dataset is 2-dimensional dataset with a set of genes and a set of conditions. A bicluster is a subset of genes that show similar behavior within a subset of conditions. Genes that show similar behavior can be considered to have same cellular functions. Thus, biclustering algorithm is a useful tool to uncover groups of genes involved in the same cellular process and groups of conditions which take place in this process. We are proposing a polynomial time algorithm to identify functionally highly correlated biclusters. Our algorithm identifies 1) the gene set that has hidden patterns even if the level of noise is high, 2) the multiple, possibly overlapped, and diverse gene sets, 3) gene sets whose functional association is strongly high, and 4) deterministic biclustering results. We validated the level of functional association of our method, and compared with current methods using GO.

• Journal of Microbiology and Biotechnology
• /
• v.32 no.3
• /
• pp.365-377
• /
• 2022

Mammalian target of rapamycin (mTOR) is a serine-threonine kinase member of the cellular phosphatidylinositol 3-kinase (PI3K) pathway, which is involved in multiple biological functions by transcriptional and translational control. mTOR is a downstream mediator in the PI3K/Akt signaling pathway and plays a critical role in cell survival. In cancer, this pathway can be activated by membrane receptors, including the HER (or ErbB) family of growth factor receptors, the insulin-like growth factor receptor, and the estrogen receptor. In the present work, we congregated an electronic network of mTORC1 built on an assembly of data using natural language processing, consisting of 470 edges (activations/interactions and/or inhibitions) and 206 nodes representing genes/proteins, using the Cytoscape 3.6.0 editor and its plugins for analysis. The experimental design included the extraction of gene expression data related to five distinct types of cancers, namely, pancreatic ductal adenocarcinoma, hepatic cirrhosis, cervical cancer, glioblastoma, and anaplastic thyroid cancer from Gene Expression Omnibus (NCBI GEO) followed by pre-processing and normalization of the data using R & Bioconductor. ExprEssence plugin was used for network condensation to identify differentially expressed genes across the gene expression samples. Gene Ontology (GO) analysis was performed to find out the over-represented GO terms in the network. In addition, pathway enrichment and functional module analysis of the protein-protein interaction (PPI) network were also conducted. Our results indicated NOTCH1, NOTCH3, FLCN, SOD1, SOD2, NF1, and TLR4 as upregulated proteins in different cancer types highlighting their role in cancer progression. The MCODE analysis identified gene clusters for each cancer type with MYC, PCNA, PARP1, IDH1, FGF10, PTEN, and CCND1 as hub genes with high connectivity. MYC for cervical cancer, IDH1 for hepatic cirrhosis, MGMT for glioblastoma and CCND1 for anaplastic thyroid cancer were identified as genes with prognostic importance using survival analysis.

• Korean Journal of Ichthyology
• /
• v.32 no.3
• /
• pp.117-129
• /
• 2020

The body color pattern in fish is a distinctive feature for species identification. The blass bloched rockfish Sebastes pachycephalus is a commercially important marine fish species, distributed in the central and southern parts of Korea and south Hokkaido of Japan. It has a morphological feature divided into four subspecies according to with or lacking distinct spots on the body surface, and to the location of markings on the body surface. However, the genetic basis of body color pattern of S. pachycephalus is still unknown. Thus we analyzed the transcriptome of S. pachycephalus skin samples using RNA-seq analysis to investigate functional genes related to body color patterns. The experimental skin samples were prepared by classified into 'Wild type' (lacking distinct spots and markings) and 'Color type' (with distinct spots and marking). Two skin sample transcriptomes were compared pairwise and the results revealed that were 164 differentially expressed unigenes in the skin samples of 'Wild type' and 'Color type'. Gene Ontology analysis of 164 differentially expressed unigenes showed that these genes were included in the functional group of molecular function (2 genes), biological process (46 genes), and cellular component (6 genes). There were several genes that body color type skin specific expression and the genes were CTL (Galactose-specific lectin nattectin), CUL1 (Cullin-1), CMAS (N-acylneuraminate cytidylyltransferase), NMRK2 (Nicotinamide riboside kinase 2), ALOXE3 (Hydroperoxide isomerase ALOXE3), SLC4A7 (sodium bicarbonate cotransporter 3). Our study is the first attempt to search for functional genes involved in the formation of body color patterns in S. pachycephalus. The differentially expressed unigenes obtained in this study can be used as candidate genes for functional gene study related to body coloration of fish.

• The Korean Journal of Physiology and Pharmacology
• /
• v.23 no.5
• /
• pp.367-379
• /
• 2019

Although atopic dermatitis (AD) is known to be a representative skin disorder, it also affects the systemic immune response. In a recent study, myoblasts were shown to be involved in the immune regulation, but the roles of muscle cells in AD are poorly understood. We aimed to identify the relationship between mitochondria and atopy by genome-wide analysis of skeletal muscles in mice. We induced AD-like symptoms using house dust mite (HDM) extract in NC/Nga mice. The transcriptional profiles of the untreated group and HDM-induced AD-like group were analyzed and compared using microarray, differentially expressed gene and functional pathway analyses, and protein interaction network construction. Our microarray analysis demonstrated that immune response-, calcium handling-, and mitochondrial metabolism-related genes were differentially expressed. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology pathway analyses, immune response pathways involved in cytokine interaction, nuclear factor-kappa B, and T-cell receptor signaling, calcium handling pathways, and mitochondria metabolism pathways involved in the citrate cycle were significantly upregulated. In protein interaction network analysis, chemokine family-, muscle contraction process-, and immune response-related genes were identified as hub genes with many interactions. In addition, mitochondrial pathways involved in calcium signaling, cardiac muscle contraction, tricarboxylic acid cycle, oxidation-reduction process, and calcium-mediated signaling were significantly stimulated in KEGG and Gene Ontology analyses. Our results provide a comprehensive understanding of the genome-wide transcriptional changes of HDM-induced AD-like symptoms and the indicated genes that could be used as AD clinical biomarkers.

• The Korean Journal of Physiology and Pharmacology
• /
• v.19 no.2
• /
• pp.141-149
• /
• 2015

"G protein-coupled receptor 40" (GPR40), a receptor for long-chain fatty acids, mediates the stimulation of glucose-induced insulin secretion. We examined the profiles of differential gene expression in GPR40-activated cells treated with linoleic acid, and finally predicted the integral pathways of the cellular mechanism of GPR40-mediated insulinotropic effects. After constructing a GPR40-overexpressing stable cell line (RIN-40) from the rat pancreatic ${\beta}$-cell line RIN-5f, we determined the gene expression profiles of RIN-5f and RIN-40. In total, 1004 genes, the expression of which was altered at least twofold, were selected in RIN-5f versus RIN-40. Moreover, the differential genetic profiles were investigated in RIN-40 cells treated with $30{\mu}M$ linoleic acid, which resulted in selection of 93 genes in RIN-40 versus RIN-40 treated with linoleic acid. Based on the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG, http://www.genome.jp/kegg/), sets of genes induced differentially by treatment with linoleic acid in RIN-40 cells were found to be related to mitogen-activated protein (MAP) kinase- and neuroactive ligand-receptor interaction pathways. A gene ontology (GO) study revealed that more than 30% of the genes were associated with signal transduction and cell proliferation. Thus, this study elucidated a gene expression pattern relevant to the signal pathways that are regulated by GPR40 activation during the acute period. Together, these findings increase our mechanistic understanding of endogenous molecules associated with GPR40 function, and provide information useful for identification of a target for the management of type 2 diabetes mellitus.

• Asian-Australasian Journal of Animal Sciences
• /
• v.22 no.6
• /
• pp.871-884
• /
• 2009

Gene expression profiling is a useful tool for identifying critical genes and pathways in metabolism. The objective of this study was to determine the major differences in the expression of genes associated with metabolism and metabolic regulation in liver and mammary tissues of lactating cows. We used the Michigan State University bovine metabolism (BMET) microarray; previously, we have designed a bovine metabolism-focused microarray containing known genes of metabolic interest using publicly available genomic internet database resources. This is a high-density array of 70mer oligonucleotides representing 2,349 bovine genes. The expression of 922 genes was different at p<0.05, and 398 genes (17%) were differentially expressed by two-fold or more with 222 higher in liver and 176 higher in mammary tissue. Gene ontology categories with a high percentage of genes more highly expressed in liver than mammary tissues included carbohydrate metabolism (glycolysis, glucoenogenesis, propanoate metabolism, butanoate metabolism, electron carrier and donor activity), lipid metabolism (fatty acid oxidation, chylomicron/lipid transport, bile acid metabolism, cholesterol metabolism, steroid metabolism, ketone body formation), and amino acid/nitrogen metabolism (amino acid biosynthetic process, amino acid catabolic process, urea cycle, and glutathione metabolic process). Categories with more genes highly expressed in mammary than liver tissue included amino acid and sugar transporters and MAPK, Wnt, and JAK-STAT signaling pathways. Real-time PCR analysis showed consistent results with those of microarray analysis for all 12 genes tested. In conclusion, microarray analyses clearly identified differential gene expression profiles between hepatic and mammary tissues that are consistent with the differences in metabolism of these two tissues. This study enables understanding of the molecular basis of metabolic adaptation of the liver and mammary gland during lactation in bovine species.

• Asian Pacific Journal of Cancer Prevention
• /
• v.16 no.7
• /
• pp.2793-2800
• /
• 2015

In molecular-targeted cancer therapy, acquired resistance to gemcitabine is a major clinical problem that reduces its effectiveness, resulting in recurrence and metastasis of cancers. In spite of great efforts to reveal the overall mechanism of acquired gemcitabine resistance, no definitive genetic factors have been identified that are absolutely responsible for the resistance process. Therefore, we performed a cross-platform meta-analysis of three publically available microarray datasets for cancer cell lines with acquired gemcitabine resistance, using the R-based RankProd algorithm, and were able to identify a total of 158 differentially expressed genes (DEGs; 76 up- and 82 down-regulated) that are potentially involved in acquired resistance to gemcitabine. Indeed, the top 20 up- and down-regulated DEGs are largely associated with a common process of carcinogenesis in many cells. For the top 50 up- and down-regulated DEGs, we conducted integrated analyses of a gene regulatory network, a gene co-expression network, and a protein-protein interaction network. The identified DEGs were functionally enriched via Gene Ontology hierarchy and Kyoto Encyclopedia of Genes and Genomes pathway analyses. By systemic combinational analysis of the three molecular networks, we could condense the total number of DEGs to final seven genes. Notably, GJA1, LEF1, and CCND2 were contained within the lists of the top 20 up- or down-regulated DEGs. Our study represents a comprehensive overview of the gene expression patterns associated with acquired gemcitabine resistance and theoretical support for further clinical therapeutic studies.

• Asian Pacific Journal of Cancer Prevention
• /
• v.15 no.10
• /
• pp.4217-4224
• /
• 2014

Cholangiocarcinoma (CCA), a slow growing but highly metastatic tumor, is highly prevalent in Northeast Thailand. Specific tests that predict prognosis of CCA remain elusive. The present study was designed to investigate whether peripheral blood leukocyte (PBL) transcriptional profiles might be of use as a prognostic test in CCA patients. Gene expression profiles of PBLs from 9 CCA and 8 healthy subjects were conducted using the Affymetrix HG_U133 Plus 2.0 GeneChip. We indentified informative PBLs gene expression profiles that could reliably distinguish CCA patients from healthy subjects. Of these CCA specific genes, 117 genes were up regulated and 60 were down regulated. The molecular and cellular functions predicted for these CCA specific genes according to the Gene Ontology database indicated differential PBL expression of host immune response and tumor progression genes (EREG, TGF ${\beta}1$, CXCL2, CXCL3, IL-8, and VEGFA). The expression levels of 9 differentially expressed genes were verified in 36 CCA vs 20 healthy subjects. A set of three tumor invasion related genes (PLAU, CTSL and SERPINB2) computed as "prognostic index" was found to be an independent and statistically significant predictor for CCA patient survival. The present study shows that CCA PBLs may serve as disease predictive clinically accessible surrogates for indentifying expressed genes reflective of CCA disease severity.

• Asian Pacific Journal of Cancer Prevention
• /
• v.14 no.11
• /
• pp.6347-6350
• /
• 2013

Post-transplant lymphoproliferative disorder (PTLD) is a common complication of therapeutic immunosuppression after organ transplantation. Gene expression profile facilitates the identification of biological difference between Epstein-Barr virus (EBV) positive and negative PTLDs. Previous studies mainly implemented variance/regression analysis without considering unaccounted array specific factors. The aim of this study is to investigate the gene expression difference between EBV positive and negative PTLDs through partial least squares (PLS) based analysis. With a microarray data set from the Gene Expression Omnibus database, we performed PLS based analysis. We acquired 1188 differentially expressed genes. Pathway and Gene Ontology enrichment analysis identified significantly over-representation of dysregulated genes in immune response and cancer related biological processes. Network analysis identified three hub genes with degrees higher than 15, including CREBBP, ATXN1, and PML. Proteins encoded by CREBBP and PML have been reported to be interact with EBV before. Our findings shed light on expression distinction of EBV positive and negative PTLDs with the hope to offer theoretical support for future therapeutic study.

• Journal of Animal Science and Technology
• /
• v.58 no.4
• /
• pp.17.1-17.7
• /
• 2016

Background: The Sapsaree (Canis familiaris) is a Korean native dog that is very friendly, protective, and loyal to its owner, and is registered as a natural monument in Korea (number: 368). To investigate large-scale gene expression profiles and identify the genes related to exercise-induced stress in the Sapsaree, we performed whole-transcriptome RNA sequencing and analyzed gene expression patterns before and after exercise performance. Results: We identified 525 differentially expressed genes in ten dogs before and after exercise. Gene Ontology classification and KEGG pathway analysis revealed that the genes were mainly involved in metabolic processes, such as programmed cell death, protein metabolic process, phosphatidylinositol signaling system, and cation binding in cytoplasm. The ten Sapsarees could be divided into two groups based on the gene expression patterns before and after exercise. The two groups were significantly different in terms of their basic body type ($p{\leq}0.05$). Seven representative genes with significantly different expression patterns before and after exercise between the two groups were chosen and characterized. Conclusions: Body type had a significant effect on the patterns of differential gene expression induced by exercise. Whole-transcriptome sequencing is a useful method for investigating the biological characteristics of the Sapsaree and the large-scale genomic differences of canines in general.