The Roles of the Interaction of BCL2-Antagonist/Killer 1, Apoptotic Peptidase Activating Factor 1 and Selenium in the Pathogenesis of Kashin-Beck Disease

Pathway-based network analyses and candidate genes associated with Kashin-Beck disease

To perform a comprehensive analysis focusing on the biological functions and interactions of Kashin-Beck disease (KBD)-related genes to provide information towards understanding the pathogenesis of KBD.A retrospective, integrated bioinformatics analysis was designed and conducted. First, by reviewing the literature deposited in PubMed, we identified 922 genes genetically associated with KBD. Then, biological function and network analyses were conducted with Cytoscape software. Moreover, KBD specific molecular network analysis was conducted by Cytocluster using the Molecular Complex Detection Algorithm (MCODE).The biological function enrichment analysis suggested that collagen catabolic process, protein activation cascade, cellular response to growth factor stimulus, skeletal system development, and extrinsic apoptosis played important roles in KBD development. The apoptosis pathway, NF-kappa B signaling pathway, and the glutathione metabolism pathway were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway network, suggesting that these pathways may play key roles in KBD occurrence and development. MCODE clusters showed that in top 3 clusters, 54 of KBD-related genes were included in the network and 110 candidate genes were discovered might be potentially related to KBD.The 110 candidate genes discovered in the current study may be related to the development of KBD. The expression changes of apoptosis and oxidative stress-related genes might serve as biomarkers for early diagnosis and treatment of KBD.

The Importance of Se-Related Genes in the Chondrocyte of Kashin-Beck Disease Revealed by Whole Genomic Microarray and Network Analysis

Kashin-Beck disease (KBD) is an endemic, chronic, and degenerative osteoarthropathy. Selenium (Se) deficiency plays important role in the pathogenesis of KBD. We aimed to screen Se-related gene from chondrocytes of patients with KBD. Whole-genome oligonucleotide microarrays were used to detect differentially expressed genes. qRT-PCR was used to confirm the microarray results. Comparative Toxicogenomics Database (CTD) was used to screen Se-related genes from differentially expressed genes. Gene Ontology (GO) classifications and network analysis of Se-related genes were constituted by STRING online system. Three hundred ninety-nine differentially expressed genes were obtained from microarray. Among them, 54 Se-related genes were identified by CTD. The qRT-PCR validation showed that four genes expressed similarly with the ones in the microarray transcriptional profiles. The Se-related genes were categorized into 6 cellular components, 8 molecular functions, 44 biological processes, 10 pathways, and 1 network by STRING. The Se-related gene insulin-like growth factor binding protein 2 (IGFBP2), insulin-like growth factor binding protein 3 (IGFBP3), interleukin 6 (IL6), BCL2, apoptosis regulator (BCL2), and BCL2-associated X, apoptosis regulator (BAX), which involved in many molecular functions, biological processes, and apoptosis pathway may play important roles in the pathogenesis of KBD.