Differential gene expression profile by RNA sequencing study of elderly osteoporotic hip fracture patients with sarcopenia

Abstract

Background: The purpose of this study was to report the RNA sequencing profile according to the presence or absence of sarcopenia in elderly patients with osteoporotic hip fracture. Therefore, an important genetic factor candidate for sarcopenia causing hip fracture in elderly with osteoporosis has been identified. Methods: The patient group involved subjects over 65 years who had undergone hip fracture surgery. Among 323 hip fracture (HF) patients identified from May 2017 to December 2019, 162 HF patients (90 non-sarcopenia and 72 sarcopenia groups), excluding subjects with high energy trauma and non-osteoporosis, were finally included in the analysis. For RNA sequencing, each patient with hand grip strength (HGS) values in the top 10% were enrolled in the control group and with the bottom 10% in the patient group. After excluding patients with poor tissue quality, 6 patients and 5 patients were selected for sarcopenia and non-sarcopenia groups, respectively. For qPCR validation, each patient with HGS values in the top 20% and bottom 20% was enrolled in the control and patient groups, respectively. After excluding patients with poor tissue quality, 12 patients and 12 patients were enrolled in the sarcopenia and non-sarcopenia groups, respectively. Sarcopenia was defined according to the Asia Working Group for Sarcopenia (AWGS) criteria for low muscle strength (hand grip strength below 18 kg in women and 28 kg in men) and low muscle mass (SMI below 5.4 kg/m(2) in women and 7.0 kg/m(2) in men). The libraries were prepared for 100 bp paired-end sequencing using TruSeq Stranded mRNA Sample Preparation Kit (Illumina, CA, USA). The gene expression counts were supplied to Deseq(2) to extract possible gene sets as differentially expressed genes (DEG) that discriminate between sarcopenia and non-sarcopenia groups that were carefully assigned by clinical observation. For the classification of the candidate genes from DEG analysis, we used the public databases; gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Quantitative real-time PCR was performed for validation. Results: Samples collected were subjected to RNAseq using the Illumina platform. A total of 11 samples from both sarcopenia and non-sarcopenia groups were sequenced. Fifteen genes (RUNX 1, NGFR, CH3L1, BCL3, PLA2G2A, MYBPH, TEP1, SEMA6B, CSPG4, ACSL5, SLC25A3, NDUFB5, CYC1, ACAT1, and TCAP) were identified as differentially expressed genes (DEG) in both the groups. In the qPCR results, the expression levels of SLC25A3 and TCAP gene in the OS group were significantly lower than in the non-OS groups whereas an increase in RUNX1 mRNA level was observed in the OS samples (p < 0.05). Conclusions: In summary, this study detected gene expression difference according to the presence or absence of sarcopenia in elderly osteoporosis female patients with hip fracture. We have also identified 15 important genes (RUNX 1, NGFR, CH3L1, BCL3, PLA2G2A, MYBPH, TEP1, SEMA6B, CSPG4, ACSL5, SLC25A3, NDUFB5, CYC1, ACAT1, TCAP), a few GO categories and biological pathways that may be associated with the osteosarcopenia. Our study may provide effective means for the prevention, diagnosis and treatment sarcopenia in elderly osteoporosis female patients. The Translational potential of this article: These findings provide a novel insight into the effects of aging on the response in women with postmenopausal osteoporosis. Further studies are underway to identify the specific signalling pathways involved. These results reveal potential therapeutic targets that could aid the regenerative capacity of aging skeletal muscle.