key: cord-1019017-zvu9ww8j authors: Li, Pengfei; Zhai, Pei; Ye, Zengjie; Deng, Peng; Fan, Yueguang; Zeng, Yirong; Pang, Zhihui; Zeng, Jianchun; Li, Jie; Feng, Wenjun title: Differential expression of miR-195-5p in collapse of steroid-induced osteonecrosis of the femoral head date: 2017-04-21 journal: Oncotarget DOI: 10.18632/oncotarget.17333 sha: 91eedbad589b474aad1ac0b96957078885146537 doc_id: 1019017 cord_uid: zvu9ww8j BACKGROUND: Femoral head collapse is a key reference point for determining a treatment regimen of femoral head osteonecrosis. However, there are no effective preventive measures and the efficacy of hip-preserving surgery is unsatisfactory due to the unclear mechanism of collapse. This study aimed to identify and validate miRNAs differentially expressed in collapse and non-collapse areas of the osteonecrotic femoral head, and to predict the target genes and pathways of these miRNAs. RESULTS: Nine samples passed the quality control test. A total of 2085 differentially expressed miRNAs were detected, among which 433 miRNAs showed differential expression in the T1 group compared to the W1 group; 344 miRNAs showed differential expression in the T2 group compared to the W2 group; 107 miRNAs showed differential expression in the T3 group compared to the W3 group. After combining data from all three patients, 10 miRNAs showed differential expression in the collapse area (T1+T2+T3) compared to the non-collapse area (W1+W2+W3). Compared to the normal area, has-miR-195-5p showed the most significant downregulation. Expression results from RT-PCR revealed that the expression of hsa-miR-195-5p in the collapse area (T1+T2+T3) was significantly lower than that in the non-collapse area (W1+W2+W3) and normal area (Z1+Z2+Z3). 157 genes were perdicted as the target gene of hsa-miR-195-5p. MATERIALS AND METHODS: Femoral heads of three patients (2 males and 1 female) treated by total hip arthroplasty surgery for steroid-induced femoral head osteonecrosis were selected based on inclusion and exclusion criteria. Bone tissue samples were obtained from the collapse area (T), non-collapse area (W), and normal area (Z) according to the anatomical structure of osteonecrotic femoral heads. Total RNA was extracted from the samples and the microarray chip was scanned. miRNAs showing differential expressions of more than 1.5-fold were selected and was validated by RT-PCR. TargetScan, mirBase and miRanda bioinformatics software was used to predict target genes and identify possible pathways involving these genes. CONCLUSIONS: miR-195-5p showed the most significant difference in the collapse area of osteonecrotic femoral heads, suggesting that collapse may be related to the downregulation of miR-195-5p. Steroid-induced osteonecrosis of the femoral head is a disabling disease, particularly after the 2003 SARS(severe acute respiratory syndrome) outbreak in China, but the mechanism of whom is not completely understood [1] . Approximately 20,000-30,000 people are diagnosed with osteonecrosis in the United States each year [2] ; the morbidity rate of non-traumatic femoral head osteonecrosis in Japan was 1.91/100,000 per year [3] . The natural course of femoral head osteonecrosis is necrosis, followed by collapse and osteoarthritis. Once collapse occurs, osteoarthritis and joint replacement is inevitable [4, 5] . Therefore, collapse must be treated surgically for femoral head osteonecrosis, and preventing collapse has been the focus of non-joint replacement therapy [6] [7] [8] [9] [10] [11] [12] [13] [14] . Current studies of collapse are typically performed from the biological, biomechanical, or mechanobiological perspective [15] [16] [17] [18] [19] ; however, the mechanism of collapse is unclear, leading to insufficient effective preventive measures and unsatisfactory efficacy of hippreserving surgery. As the number of studies of miRNA in orthopedics has increased, specific miRNAs involved in femoral head osteonecrosis have been identified [20] [21] [22] [23] [24] . However, the developmental course of femoral head osteonecrosis is a dynamic process. Whether miRNA is differentially expressed in different states of osteonecrotic femoral head (non-collapse vs. collapse) remains unclear. Studies of differential miRNA expression may increase the understanding of the biological mechanisms underlying collapse. This study aimed to identify miRNAs differentially expressed in collapse and non-collapse areas of the osteonecrotic femoral head and validate these miRNAs using reverse transcriptionpolymerase chain reaction (RT-PCR) assays, and eventually predict the target genes and pathways of these miRNAs. The 9 selected samples passed the quality control test .A total of 2085 differentially expressed miRNAs were detected in the 9 samples. Among the detected miRNAs, 433 miRNAs showed more than 1.5-fold differential expression in the T1 group compared to the W1 group (247 miRNAs were up-regulated and 186 miRNA were downregulated); 344 miRNAs were differentially expressed by more than two-fold in the T2 group compared to the W2 group (225 miRNAs were up-regulated and 119 miRNAs were down-regulated); 107 miRNAs showed differential expression by more than two-fold in the T3 group compared to the W3 group (51 miRNAs were up-regulated and 56 miRNAs were down-regulated) ( Figure 1 ). Combining the data from all three patients, 10 miRNA were differentially expressed by more than twofold in the collapse area (T1+T2+T3) compared to the non-collapse area (W1+W2+W3) (8 miRNAs were upregulated: hsa-miR-4472, hsa-miR-4306. hsa-miR-4747-5p, hsa-miR-4441, hsa-miR-4709-3p, ebv-miR-BHRF1-2-3p, hsa-miR-585-3p, and hsa-miR-5572, while 2 miRNAs were down-regulated: hsa-miR-195-5p and hsa-miR-645) . Compared to the normal region, has-miR-195-5p showed the most significant down-regulation (Tables 1-3) . The following U6 primers were used for RT-PCR:F: 5′GCTTCGGCAGCACATAT ACTAAAAT3′, R: 5′CGCTTCACGAATTTGCGTGT CAT3′.The primer sequences of hsa-miR-195-5p were: GSP: 5′GGGGTAGCAGCACAGAAAT3′, R: 5′CAGTGCGTGTCGTGGAGT3′(GSP is a specific primer for miRNA, while R is a primer that matches the RT primer). The RT-PCR expression level of hsa-miR-195-5p was significantly lower in the collapse area (T1 T2 T3)than in the non-collapse area (W1 W2 W3) and normal area (Z1 Z2 Z3) . (Table 4 ) Comparison of the results from RT-PCR and the microarray data also indicated that hsa-miR-195-5p expression was downregulated.Amplification and melting curves were obtained ( Figure 2 ). 157 genes were perdicted as the target gene of hsa-miR-195-5p (Figures 3, 4) . GO analysis and KEGG Pathway analysis were showed that pathway in cancer is significant (Figures 5, 6 ). Previous studies have suggested that collapse occurs because of an imbalance between bone resorption and bone formation after femoral head osteonecrosis, resulting in femoral heads with reduced structural strength that are unable to bear external forces. The studies have also suggested that mechanical support and adjustment of bone cell metabolism can prevent collapse [25] . To comprehensively understand the biological behavior of osteocytes, an increasing number of studies have focused on the differential expression of miRNAs associated with osteonecrosis of the femoral head [23] [24] [26] [27] [28] . We evaluated whether the decrease in structural strength was associated with the differential expression of relevant miRNAs, which may lead to faster apoptosis and an increased imbalance with bone metabolism. We found that in the osteonecrotic region, hsa-miR-195-5p was down-regulated in the collapse area compared to in the non-collapsed area. This is the first study examining the relationship between miRNA and femoral head collapse. www.impactjournals.com/oncotarget The anatomical structure of osteonecrotic femoral head has four layers: cartilage, osteonecrotic area, sclerotic area, and normal area. Recent studies have proposed that the decrease in structural strength is correlated with sclerotic area and osteonecrotic edge, enabling collapse prediction [16, 19] . During the pathological process, collapse of the osteonecrotic femoral head occurs because of continuous spreading of trabecular fracture, which is not a sudden process, and thus the osteonecrotic area includes the collapse area and non-collapse area. Femoral head collapse originates from subchondral fracture. Therefore, in this study, subchondral fracture in the osteonecrotic Oncotarget 42643 www.impactjournals.com/oncotarget femoral head was defined as the collapse area, while the area in the osteonecrotic femoral head without subchondral fracture was defined as the non-collapse area. Non-collapse and collapse are different pathological stages of osteonecrotic femoral head. We not only developed a method for predicting the risk of collapse, but also provide a new insight into preventing collapse using gene editing techniques. It has been reported that miR-195-5p is a tumorsuppressor gene that mainly suppresses cell proliferation and invasion and is down regulated in cancer patients [29] [30] [31] [32] [33] [34] . This gene inhibits hepatocellular carcinoma by regulating PHF19 [30] , regulates rectal cancer through CDK8, regulates osteosarcoma by inhibiting NKD1 [33] , regulates prostate cancer through Fra-1 or RPS6KB1 [34] . Based on these results, we hypothesized that the mechanism of collapse involves down-regulation of miR-195-5p, which disrupts the proliferation of normal osteoblasts and accelerates cell apoptosis. This leads to a further decrease in structural strength within the femoral head and failure to support external forces, eventually resulting in collapse. Future studies should examine methods of increasing the expression of miR-195-5p to treat cancer and the collapse of femoral head osteonecrosis. The limitations of this study were as follows: firstly, small sample size, but with reproducibility; secondly, target genes and pathways were only predicted and require further validation; In addition, whether collapse leads to miRNA changes or miRNA changes lead to collapse remains unknown. These factors exhibit a mutually causal relationship. In future studies, we will evaluate the genetic pathways of the collapse-associated differential miRNAs and gene editing technology to prevent and treat collapse. Three patients were randomly selected from the hospital system in our department. Inclusion criteria were as follows: 1. diagnosed with steroid-induced femoral head osteonecrosis (a clear history of hormone use); 2. diagnosed with ACRO III/IV stage based on X-ray by 3 doctors; and 3. the surgeon and patient decided to perform the total hip joint replacement surgery. Exclusion criteria were as follows: 1. The cause of femoral head osteonecrosis was not clear or involved multiple factors; 2. Severe liver and kidney failure in which the patient could not tolerate surgery; And 3. Surgical strategy in which patient consent was not obtained. Informed consent was obtained from all patients and the ethical review committee of our hospital. According to the inclusion and exclusion criteria, three patients were selected. There were 2 males, one of whom at the age of 30 years used hormones to treat thrombocytopenic purpura, while the other, at the age of 46 years, used hormones to treat skin disease. At the age of 20 years, the female patient used hormones to treat systemic lupus erythematosus. After surgery, femoral heads were removed and cleaned with saline. Three bone tissues (normal area, collapse area, and non-collapse area) were taken from femoral heads (Figure 7) . Two doctors identified these tissues, and a total of 9 samples were obtained (T1 T2 T3, W1 W2 W3, and Z1 Z2 Z3). Bone tissues were immediately stored in liquid nitrogen and ground into powder (≥ 100 mg) in liquid nitrogen. The powder was transferred to an microcentrifuge tube, mixed with 1 mL Trizol, and stored at −70°C until use. Total RNA extraction, RNA labeling, hybridization and microarray chip scanning [23] Total RNA extraction and quality control Total RNA was extracted from tissue cells using Trizol Reagent (Invitrogen life technologies) according to Oncotarget 42645 www.impactjournals.com/oncotarget the manufacturer's instructions; total RNA was extracted from whole blood using TRI reagent BD (MRCgene, TB-126). Total RNA was extracted from plasma and serum exosomes using TRIzol LS. RNAsey Mini Kit (QIAGEN) was used for RNA purification; NanoDrop ND-1000 was employed for measurement of RNA concentration after purification. RNA integrity was detected by electrophoresis. After extraction, RNA samples underwent quality control and miRNA was labeled using miRCURYTM Array Power Labeling kit (Cat #208032-A, Exiqon). The detailed steps are described as follows:One microgram RNA sample was added to 2 μL water, and then 1 μL CIP buffer and CIP enzyme (Exiqon) were added to the RNA sample solution, and the mixture solution was incubated at 37°C for 30 min. The reaction in the sample solution was terminated at 95°C for 5 min, and then 3 μL labeling buffer, 1.5μL fluorescent label (Hy3 ™ ), 2.0 μL DMSO, and 2.0 μL labeling enzyme were added and the solution was incubated at 16°C for 1 h.The reaction in the samples was terminated at 65°C for 15 min. After labeling, the sample was hybridized with miRCURYTM LNA Array (v.18.0) (Exiqon), and the subsequent procedures were conducted according to Exiqon's experimental methods.the hybridization system used was the Nimblegen system (Nimblegen Systems Inc., Addison, WI, USA).The microarray chip was rinsed with Wash buffer kit (Exiqon) after hybridization. GenePix Pro V6.0 was utilized to read the microarray scan images and the probe signals were selected. The median was taken and combined for the same probes. Probes >= 30.0 for the samples were kept and the median was standardized in all microarray chips. (Figure 8 ) miRNAs showing differential expressions of more than 1.5-fold were selected [24] . cDNA synthesis of RNA samples [23] . 20 μL of RT reaction mixture was prepared, then all cDNA samples were separately set up in the RT-PCR System (Applied Biosystems); The solution was mixed up to 8 μL and centrifuged at 5000 rpm for a short period. The mixture was added to each well of a 384-well PCR plate, and then 2 µL corresponding cDNA was added to each well. The plate was sealed with sealing film and centrifuged. The 384-well PCR plate was placed in the RT-PCR amplification instrument to carry out the PCR reaction. U6 (internal reference) and all parameters were conducted. The RT-PCR reaction was performed for the target miRNA in all samples and the internal reference gene (U6) separately. In Table 1 , data analysis was performed using 2 -∆∆CT method. SPSS 20.0 (SPSS Inc., Chicago, Ill., USA) was employed for statistical analysis of the data. TargetScan, mirBase and miRanda bioinformatics software were used to predict target genes. Gene Ontology(GO, http://www.geneontology.org/) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG, www.genome.jp/kegg/) pathway analysis were done respectively. The GO project covers three domains: Biological Process, Cellular Component and Molecular Function. Fisher's exact test is used to find if there is more overlap between the differential expression (DE) list and the GO annotation list than would be expected by chance. The p-value denotes the significance of GO terms enrichment in the DE genes. The lower the p-value, the more significant the GO Term (p-value <= 0.05 ).Pathway analysis is a functional analysis mapping genes to KEGG pathways. The p-value denotes the significance of the Pathway correlated to the conditions. Lower the p-value, more significant is the Pathway(p-value cut-off is 0.05). None. The influence of age, gender and treatment with steroids on the incidence of osteonecrosis of the femoral head during the management of severe acute respiratory syndrome: a retrospective study Osteonecrosis of the hip: management in the 21st century Epidemiology of nontraumatic osteonecrosis of the femoral head in Japan The natural history of untreated asymptomatic osteonecrosis of the femoral head: a systematic literature review Current management options for osteonecrosis of the femoral head: part 1, diagnosis and nonoperative management Which factors influence preservation of the osteonecrotic femoral head Osteonecrosis of the femoral head: evaluation and treatment Osteonecrosis: avoiding total hip arthroplasty Fate of untreated asymptomatic osteonecrosis of the femoral head Treatment of early stage osteonecrosis of the femoral head Current concepts on osteonecrosis of the femoral head Management of avascular necrosis of femoral head at pre-collapse stage Management of femoral head osteonecrosis: Current concepts Relationship between preservation of the lateral pillar and collapse of the femoral head in patients with osteonecrosis The role of sclerotic changes in the starting mechanisms of collapse: A histomorphometric and FEM study on the femoral head of osteonecrosis Does Zoledronate Prevent Femoral Head Collapse from Osteonecrosis? A Prospective, Randomized, Open-Label, Multicenter Study Prediction of femoral head collapse in osteonecrosis using 18F-fluoride positron emission tomography Prediction of osteonecrosis collapse of the femoral head based on the proportion of the proximal sclerotic rim Angiogenic microRNA-210 is present in cells surrounding osteonecrosis Involvement of MicroRNA-210 Demethylation in Steroid-associated Osteonecrosis of the Femoral Head MiR-708 promotes steroid-induced osteonecrosis of femoral head, suppresses osteogenic differentiation by targeting SMAD3 Differential expression of miR-672-5p and miR-146a-5p in osteoblasts in rats after steroid intervention Identification of differentially expressed microRNAs involved in non-traumatic osteonecrosis through microRNA expression profiling Key pathway to prevent the collapse of femoral head in osteonecrosis A single-nucleotide polymorphism in MMP9 is associated with decreased risk of steroid-induced osteonecrosis of the femoral head Icariside II activates EGFR-Akt-Nrf2 signaling and protects osteoblasts from dexamethasone miR-27a attenuates adipogenesis and promotes osteogenesis in steroid-induced rat BMSCs by targeting PPARγ and GREM1 MiR-497~195 cluster microRNAs regulate osteoblast differentiation by targeting BMP signaling MicroRNA-195-5p acts as an anti-oncogene by targeting PHF19 in hepatocellular carcinoma MicroRNA-195-5p suppresses glucose uptake and proliferation of human bladder cancer T24 cells by regulating GLUT3 expression Circulating Oncotarget 42647 www.impactjournals.com/oncotarget microRNAs as prognostic therapy biomarkers in head and neck cancer patients Identification of a plasma four-microRNA panel as potential noninvasive biomarker for osteosarcoma MicroRNA-195-5p, a new regulator of Fra-1, suppresses the migration and invasion of prostate cancer cells This research was supported by grants from the National Natural Science Foundation of China (No. 81373652).