Animals and experimental design
New Zealand white rabbits (male, 6 months old, 2.0–2.5 kg) came from Laifu Breeding Farm (Nanjing, China). Eighteen rabbits were used for transcriptomic analysis, and another 60 were used for mechanistic studies. Housing conditions included 20–24 °C, 45–55% humidity, a 12/12 h light/dark cycle, and free feeding. Ethics approval came from Anhui University of Chinese Medicine (approval no. AHUCM-rabbits-2023195). Transcriptomics used three rabbit groups (n = 6 each): normal, model, and acupotomy. Mechanistic studies used five groups (n = 12 each): normal, model, acupotomy, TRPV4 inhibitor, and acupotomy + inhibitor. Rabbits were randomly allocated to experimental groups using a computer-generated random number sequence. The allocation was performed by an investigator not involved in the subsequent interventions or outcome assessment. Normal rabbits had no intervention across 63 days. Model rabbits received no treatment for 21 days after induction. Acupotomy occurred on days 49, 56, and 63 post-induction. For the inhibitor group, HC067047 (a TRPV4 channel inhibitor) was administered intragastrically once daily at 10 mg/kg from day 43 to day 63 after modeling (21 consecutive days). The dose of HC067047 (10 mg/kg) was selected based on previous studies demonstrating effective TRPV4 blockade in animal models, including paclitaxel-induced peripheral neuropathy31 . Acupotomy was performed on days 49, 56, and 63 after modeling. Thus, HC067047 treatment covered the entire acupotomy intervention period: the first, second, and third acupotomy sessions corresponded to treatment days 7, 14, and 21, respectively. In the combination group, HC067047 was administered using the same regimen, and acupotomy was performed on the same schedule (days 49, 56, and 63). Rabbits were selected because their knee joint size and cartilage thickness allow reliable acupotomy manipulation and reproducible histopathological/molecular sampling. Compared with rodents, the rabbit stifle joint provides a larger surgical field for standardized needle insertion and sweeping, and the modified Videman immobilization model in rabbits reproduces key KOA features including cartilage degeneration and chondrocyte apoptosis. These characteristics make the rabbit model particularly suitable for evaluating acupotomy and TRPV4-related mechanisms.
Details of all reagents, equipment, and software used in this study are provided in the Table of Materials.
Model establishment
KOA was induced by the modified Videman method29. Rabbits were fasted for 12 h before modeling and weighed. Anesthesia was induced by ear vein injection of 3% pentobarbital sodium at 30 mg/kg. The left hindlimb was shaved from the groin to the ankle. The procedure was performed by two experimenters: one restrained the rabbit by holding the forelimbs and the right hindlimb, while the other pulled the left hindlimb into full extension. Medical pressure-sensitive adhesive tape was first wrapped from the groin to the ankle to protect the fur. Double-sided foam tape was then applied over the adhesive tape from the groin to the middle of the lower leg. A polymer plaster bandage (softened in 65–85 °C water) was wrapped externally under pressure to ensure the knee was extended at 180° and the ankle dorsiflexed at 60°. The toes were left exposed for daily monitoring of blood supply. The gaps were filled with absorbent cotton, and the limb was further fixed with three layers of thick cardboard, polymer plaster bandage, and cloth strips to prevent biting. The bandage was inspected regularly for toe congestion or swelling. If swelling occurred, the cast was removed temporarily and re-applied after normal blood supply was restored. Immobilization lasted for 6 weeks. After the 6-week immobilization period, the casts and bandages were removed. Model success was verified by the modified Lequesne MG score30 (pain/gait/mobility: 0–3 each; swelling: 0–2; total 0–11). A score > 6 confirmed successful KOA induction.
During the immobilization period, the casts were checked daily, and the blood supply of the toes on the affected side was monitored. If mild toe swelling occurred due to an overly tight cast, the cast was adjusted appropriately. If swelling worsened, or ulceration or necrosis appeared, the cast was removed immediately and re-applied after swelling subsided. Because mobility was restricted, the contralateral limb was monitored to prevent it from becoming trapped in the metal cage; if trapped, it was gently released. Skin wounds were treated promptly, and the cast was examined for displacement or excessive tightness. If the cast became too loose to maintain knee extension, it was removed and replaced.
Acupotomy procedure
Acupotomy was performed on days 49, 56, and 63 post-induction. Rabbits were fixed on a treatment table. The following sites were selected as needle insertion points: the musculotendinous junction of the vastus medialis (superomedial to the patella), the musculotendinous junction of the vastus lateralis (superolateral to the patella), the musculotendinous junction of the rectus femoris (central superior patellar margin), the musculotendinous junction of the biceps femoris (below the lateral femoral condyle near the fibular head), the anserine bursa (anteroinferior to the medial femoral condyle, at the conjoined tendon of the sartorius, gracilis, and semitendinosus), and palpable cord-like nodules. Cord-like nodules were identified by palpation around the knee before each treatment; the two most prominent nodules were marked. All sites were localized with gentian violet, and the skin was shaved and disinfected routinely. A disposable sterile acupotomy needle (0.35 mm × 25 mm; vendor information provided in the Table of Materials) was used. The blade was oriented parallel to the longitudinal axis of the tendon or the selected muscle-tendon cord. The needle was inserted perpendicular to the skin surface until it reached the tendon-bone junction. Insertion depth was determined by the location of the cord-like structure or by a sense of resistance beneath the blade. Longitudinal sweeping was performed for 2–3 strokes, followed by transverse separation for 1–2 strokes. The needle was withdrawn when a sense of loosening was felt. After withdrawal, a sterile cotton ball was pressed briefly to stop bleeding. Acupotomy was performed once weekly for 3 weeks. The treatment endpoint was defined as a sense of loosening beneath the blade and completion of the scheduled three sessions.
Sample collection
On day 64, rabbits were euthanized by pentobarbital overdose (3%, ≥150 mg/kg, i.p.). Knee cartilage was collected and processed as follows: (1) fixed in 4% paraformaldehyde for histology; (2) fixed in 2.5% glutaraldehyde for electron microscopy; or (3) snap-frozen in liquid nitrogen and stored at –80 °C for molecular assays. The 18 rabbits used for transcriptomics were separate from the 60 rabbits used for mechanistic studies. For mechanistic studies, 60 rabbits were used (n = 12/group). Within each group, 6 rabbits were randomly selected for qRT-PCR, 3 for Western blotting, and 3 for histology/immunofluorescence/flow cytometry. All subset selections were performed randomly.
Histological scoring
Two pathologists who were blinded to the experimental group assignments independently evaluated the stained sections using the Mankin scoring system to assess the degree of cartilage degeneration27. Before evaluation, all sections were assigned random identification codes. The treatment allocations were not disclosed to the pathologists until all scoring was complete.
Transcriptomic analysis
Total RNA from pulverized cartilage was isolated using TRIzol (1 mL per 50–100 mg tissue), followed by chloroform extraction, isopropanol precipitation, 75% ethanol wash, and DEPC-water resuspension. Purity and integrity were checked on a spectrophotometer and bioanalyzer. For sequencing, mRNA was pulled down with oligo-dT beads, heat-fragmented (94 °C, 100–300 bp), and converted to cDNA. Libraries underwent end repair, A-tailing, adaptor ligation, and 12–15 cycles of PCR before paired-end sequencing (2 × 150 bp) on a high-throughput sequencing platform.
Trimming of raw reads used Trim Galore (Q < 10). Clean reads were mapped to the rabbit reference genome OryCun2.0 via HISAT2. FPKM values for gene expression were obtained with Cufflinks. To address the complexity of the rabbit genome—particularly expanded gene families and highly similar paralogous regions—multi-mapping reads (those aligning to more than one locus, NH:i > 1) were excluded and only uniquely mapped reads (NH:i = 1) were retained for gene-level quantification, which also prevents spurious allocation of paralogous reads to a single gene. Differential expression thresholds were p < 0.05 and |log₂FC| > 1 using cuffdiff. KEGG analysis was run with clusterProfiler (R, p < 0.05).
Hematoxylin and eosin (H&E) staining
Fixed cartilage (4% paraformaldehyde, 72 h) was decalcified in EDTA (4 weeks), dehydrated, embedded in paraffin, and sectioned (5 μm). Sections were stained with H&E and observed under a light microscope.
Transmission electron microscopy (TEM)
Specimens (~1 mm3) were fixed in 2.5% glutaraldehyde (6 h), post-fixed in 1% osmic acid (2 h), dehydrated in ethanol, stained with uranyl acetate (3 h or overnight), treated with propylene oxide (30 min), and embedded in epoxy resin. Ultrathin sections (70 nm) were stained with lead citrate and examined under a transmission electron microscope.
Measurement of intracellular calcium levels
Tissues were minced in PBS containing 1 mM EDTA, filtered (70 μm strainer), and centrifuged (300 × g, 10 min). After red blood cell lysis (10 min, room temperature), cells were incubated with 2 μM Fluo-4 AM (37 °C, 20 min, dark). After two washes with calcium-free PBS, fluorescence was measured by flow cytometry. For gating, chondrocytes were first gated on forward and side scatter, followed by single-cell gating. Unstained resting cells were used as a negative reference to set the calcium threshold. Data are expressed as the mean fluorescence intensity (MFI) ratio of Fluo-4 relative to the normal group. This threshold-based dynamic comparison reflects differences in calcium signaling among groups.
Immunofluorescence staining
After dewaxing and rehydration of paraffin sections, antigen retrieval was carried out. Sections were treated with 3% H₂O₂ for 10 min to inactivate endogenous peroxidase, followed by a 30-min block with goat serum for background reduction. Primary antibodies recognizing Calmodulin (1:200) and ACAN (1:300) were added and left overnight at 4 °C. Sections were then washed and reacted with HRP-conjugated secondary antibody (30 min), then with TSA-fluorescent dye for 10 min. DAPI was applied for 5 min to label nuclei. Fluorescence images were obtained using a standard fluorescence microscope. For quantification, three randomly selected fields per section were captured under identical exposure settings. The mean fluorescence intensity (MFI) of ACAN or Calmodulin was measured using ImageJ and normalized to the tissue area.
Western blotting
Protein extraction from cartilage (~0.1 g) was performed using RIPA buffer with 0.6 mM PMSF. Lysates were cleared by centrifugation (12,000 × g, 15 min). The resulting supernatants underwent SDS-PAGE separation (90 V, 90 min) and subsequent PVDF membrane transfer. Membranes were blocked with 5% non-fat milk (2 h) and then probed overnight at 4 °C with primary antibodies against Caspase-3 (1:5000), Caspase-9 (1:2000), Caspase-12 (1:2000), ACAN (1:1000), COL2 (1:2000), IFT88 (1:2000), and GAPDH (1:5000). After a 1.2-h incubation with HRP-conjugated secondary antibody (1:20,000), signals were detected by ECL and analyzed with image analysis software. The intensity of each target protein was normalized to GAPDH.
Quantitative real-time PCR (qRT-PCR)
Total RNA was reverse-transcribed using a reverse transcription kit. qPCR was performed using a SYBR Green-based qPCR master mix under the following conditions: 95 °C for 1 min; 40 cycles of 95 °C for 20 s and 60 °C for 1 min. Relative expression was calculated using the 2⁻ΔΔCt method, with β-actin as the reference gene. Primers are listed in Table 1.
Statistical analysis
Data are presented as mean ± SEM. Normality was assessed using the Shapiro–Wilk test and Q–Q plots, and homogeneity of variance was assessed using Levene’s test. Because some experiments had a small sample size (n = 3 per group), these assumption tests have limited statistical power; the parametric results should therefore be interpreted with this limitation in mind. For two-group comparisons, an unpaired two-tailed Student’s t-test was used. For multi-group comparisons, one-way ANOVA followed by Tukey’s post hoc test was used. A two-sided p < 0.05 was considered statistically significant. All analyses were performed using GraphPad Prism 10. No a priori sample-size calculation was performed; the sample size was based on previous studies using similar rabbit KOA models and on practical constraints.