Whole-Transcriptome Sequencing of Knee Joint Cartilage from Kashin–Beck Disease and Osteoarthritis Patients

Lixin Han,Bolun Cheng,Wenming Wei,Li Liu,Shiqiang Cheng,Huan Liu,Yumeng Jia,Yan Wen,Feng Zhang
DOI: https://doi.org/10.3390/ijms25084348
IF: 5.6
2024-04-16
International Journal of Molecular Sciences
Abstract:The aim of this study was to provide a comprehensive understanding of similarities and differences in mRNAs, lncRNAs, and circRNAs within cartilage for Kashin–Beck disease (KBD) compared to osteoarthritis (OA). We conducted a comparison of the expression profiles of mRNAs, lncRNAs, and circRNAs via whole-transcriptome sequencing in eight KBD and ten OA individuals. To facilitate functional annotation-enriched analysis for differentially expressed (DE) genes, DE lncRNAs, and DE circRNAs, we employed bioinformatic analysis utilizing Gene Ontology (GO) and KEGG. Additionally, using quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), we validated the expression levels of four cartilage-related genes in chondrocytes. We identified a total of 43 DE mRNAs, 1451 DE lncRNAs, and 305 DE circRNAs in KBD cartilage tissue compared to OA (q value 1). We also performed competing endogenous RNA network analysis, which identified a total of 65 lncRNA-mRNA interactions and 4714 miRNA-circRNA interactions. In particular, we observed that circRNA12218 had binding sites for three miRNAs targeting ACAN, while circRNA12487 had binding sites for seven miRNAs targeting COL2A1. Our results add a novel set of genes and non-coding RNAs that could potentially serve as candidate diagnostic biomarkers or therapeutic targets for KBD patients.
biochemistry & molecular biology,chemistry, multidisciplinary
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to comprehensively understand the similarities and differences of mRNA, long - non - coding RNA (lncRNA) and circular RNA (circRNA) in the knee joint cartilage of patients with Kashin - Beck disease (KBD) and osteoarthritis (OA) through whole - transcriptome sequencing. Specifically, the research aims to: 1. **Compare the gene expression patterns between the two diseases**: By performing whole - transcriptome sequencing on knee joint cartilage samples from patients with KBD and OA, identify the mRNA, lncRNA and circRNA that are differentially expressed between KBD and OA. 2. **Reveal potential biological functions and signaling pathways**: Utilize Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to determine the key biological functions and signaling pathways in which these differentially - expressed genes and non - coding RNAs are involved. 3. **Construct a competing endogenous RNA network**: By analyzing the interactions between lncRNA - mRNA and circRNA - miRNA, construct a competing endogenous RNA (ceRNA) network to further understand the regulatory mechanisms of these non - coding RNAs in KBD. 4. **Verify the expression levels of candidate genes**: Use quantitative reverse transcriptase polymerase chain reaction (qRT - PCR) to verify the expression levels of four cartilage - related genes in chondrocytes to confirm the accuracy of the sequencing results. Through these studies, the authors hope to discover new genes and non - coding RNAs, which may serve as diagnostic markers or therapeutic targets for KBD, thereby providing new insights into the pathological mechanisms and treatment methods of KBD. ### Key formulas In this paper, the formulas related to the statistical significance threshold and fold - change of differential expression are as follows: - Screening conditions for differentially - expressed genes: \[ q\text{value} < 0.05 \] - Logarithmic fold - change: \[ |\log_2 \text{FC}| > 1 \] These conditions are used to screen out the mRNA, lncRNA and circRNA that are significantly differentially expressed between KBD and OA.