PDF

Int J Med Sci 2022; 19(1):98-104. doi:10.7150/ijms.65343 This issue Cite

Research Paper

RNA isolation from micro-quantity of articular cartilage for quantitative gene expression by microarray analysis

Xiaowei Zhang, MD^1,2, Trevor J. McFarland, BS³, Kristina Vartanian, BS³, Yong Zhu, PhD⁴, Christina A. Harrington, PhD^3,5, Cong-Qiu Chu, MD, PhD^1,2

1. Division of Arthritis and Rheumatic Diseases, Oregon Health & Science University, Portland, Oregon 97239.
2. Section of Rheumatology, VA Portland Health Care System, Portland, Oregon 97239.
3. Gene Profiling Shared Resource, Oregon Health & Science University; Portland, Oregon 97239.
4. Vivoscript, Inc, P. O. Box 63025, Irvine, CA 92602.
5. Department of Molecular and Medical Genetics, Oregon Health & Science University; Portland, Oregon 97239.

✉ Corresponding author: chucedu

Abstract

Isolation of quality RNA from articular cartilage has been challenging due to low cellularity and the high abundance of extracellular matrix and proteoglycan proteins. Recently developed methods for isolation of high quality RNA from cartilage are more applicable to larger cartilage specimens typically weighing at least 25 mg. While these methods generate RNA suitable for analysis, they are less successful with smaller tissue inputs. For the study of small focal defect cartilage specimens an improved RNA extraction method is needed. Here we report a protocol for direct RNA isolation from less than 3 mg of wet weight rabbit articular cartilage for quantitative microarray gene profiling. This protocol is useful for identifying differentially expressed genes in chondrocytes following focal cartilage repair and can potentially be adopted for gene expression analysis of cartilage biopsy specimens from human joints.

Keywords: RNA isolation, articular cartilage, cartilage repair, microarray

Introduction

SRY-type high-mobility group box 9 (SOX9) is a master transcription factor of chondrogenesis. In previous studies we have shown that an engineered cell-permeable super-positively charged SOX9 (scSOX9) protein was able to improve the quality of repaired cartilage by microfracture [1, 2]. Microfracture is a common procedure for cartilage repair, but often generates only fibrocartilage, which is inferior to hyaline cartilage in function and durability. In a rabbit cartilage injury model, applied at the site of microfracture, scSOX9 successfully helped regenerate hyaline-like cartilage. We are interested in identifying gene expression profiles to elucidate the metabolism of cartilage regeneration and cartilage defect repair.

Isolation of intact RNA from adult articular cartilage is challenging. Cartilage is characterized by low cellularity and a high content of highly cross-linked extracellular matrix proteins with aggregating proteoglycans [3]. Both classes of molecules are extremely large and highly negatively charged. Direct isolation of RNA from articular cartilage is historically problematic with low yield, low purity, and poor integrity of the extracted RNA. Most gene expression studies are based on RNA extracted from cultured chondrocytes. In situ hybridization can represent an alternative technique to examine gene expressions in cartilage, but it is not quantitative [4]. Several developed protocols are capable of successfully isolating quality RNA from articular cartilage but require large amounts of tissue [5-8]. Zheng et al reported that good quality RNA can be isolated from 25 mg of tissue but this required pooling eight femoral heads of young mice (2-3 weeks of age) [8]. However, cartilage specimens harvested from individual repaired cartilage sites result in small quantities of tissue and hence these published protocols for RNA isolation are not suitable. Thus, our main goal was to develop a reliable protocol for the isolation of quality RNA for quantitative gene expression analysis from micro-quantities of tissue harvested from repaired cartilage at the site of the focal defect. With our modified method, we successfully isolated RNA from less than 3 mg of cartilage (wet weight) of sufficient quality for gene expression analysis.

Materials and Methods

Cartilage tissue

All animal studies were approved by Institutional Animal Care and Use Committee of the VAPORHCS, #3469-15). Mature female New Zealand white rabbits with body weights of 3.0 - 3.5 kg were used. A cylindrical osteochondral defect of 4 mm in diameter and 3 mm in depth is made at the patella groove of the femur. This is followed by microfracture [1] supplemented with collagen membrane or collagen membrane carrying scSOX9 or a loss-of-function mutant, scSOX9-A76E. Eight weeks following the procedure, full thickness cartilage plugs of 4 mm in diameter were harvested from the repaired cartilage in the defect area. All cartilage samples were weighed immediately upon harvest (weights ranged from 1.9 to 2.9 mg) and either snap frozen in liquid nitrogen or stored in 350 µL of TRIzol (Thermo Fisher Scientific). All samples were stored at -80°C until the RNA isolation step was performed.

RNA isolation

RNA isolation and microarray assays were performed at the Oregon Health & Science University Gene Profiling Shared Resource. Quality control of RNA isolation was performed according to the quality assessment procedures developed by the core laboratory [9]. To improve the RNA isolation process from this challenging tissue source two methods with various modifications were evaluated (see flow chart in Figure 1 and Table 1 for detailed methodology): 1. A modified TRIzol (Thermo Fisher Scientific) extraction followed by an RNeasy column purification (QIAGEN) [10] or 2. A modified RNAqueous (Ambion) extraction. Samples stored in TRIzol were thawed and disrupted using stainless steel beads and a TissueLyser II device (QIAGEN) set at 50Hz for 5 minutes. Samples that had been flash frozen were thawed after the addition of 500 µL of RNAqueous lysis buffer; these samples were also disrupted with the TissueLyser II as above. RNA was isolated per manufacturer's instructions with the following modifications: each protocol was tested with and without the addition of a Proteinase K (QIAGEN) digestion step using 10 µl of the enzyme (>600mAU, incubated at 55°C for 20 minutes) and with and without a viscosity-reducing homogenization step using a QIAshredder column (QIAGEN). Following the various tissue pre-processing steps, the lysates were passed through either an RNeasy (QIAGEN) or RNAqueous (Ambion) RNA binding column and washed with kit specific buffers. RNA was eluted with 30-50 µL of nuclease-free water (Figure 1 and Table 1). RNA concentration, yield, and purity ratios were measured by UV absorbance using a Nanodrop 1000 Spectrophotometer (ThermoFisher Scientific). RNA size distribution and integrity were determined using a 2100 Bioanalyzer (Agilent Technologies) and an RNA 6000 Pico chip. For samples assayed on the Bioanalyzer with the Total RNA Program, Agilent software uses an algorithm to assign an RNA integrity score based on the trace pattern. This score is a value between 1 and 10 and is referred to as an RNA integrity number (RIN) [11, 12]. Based on initial results obtained using the RNAqueous method, this procedure was tested on only a few samples due to the limited number of samples available. Ultimately, based on overall RNA quality and yield results (Table 2), Method 1a, the TRIzol method with Proteinase K digestion followed by QIAshredder column homogenization and RNeasy column purification (QIAGEN) (Figure 1 and Table 1), was selected for preparing samples for downstream expression profiling.

Microarray analysis

The input for microarray target preparation was 30 ng of total RNA. RNA was reverse-transcribed, amplified and labeled with biotin using the Ovation Pico WTA v.2 RNA amplification kit (NuGEN Technologies) with the Encore Biotin Labeling Module. Amplified target cDNA yield was determined by OD₂₆₀ absorbance and all samples passed vendor QC thresholds. Amplified and labeled cDNA targets were combined with hybridization solution and controls (ThermoFisher Scientific, formerly Affymetrix). Three micrograms of labeled target in 130 µl of hybridization solution was injected into a GeneChip Rabbit Gene 1.0 array cartridge (ThermoFisher Scientific, formerly Affymetrix) designed to interrogate over 20,000 rabbit genes. Arrays were incubated for 18 hours at 45 °C and processed as per manufacturer's instructions. Arrays were scanned using the GeneChip Scanner 3000 7G with autoloader. Image processing was performed using the Affymetrix GeneChip Command Console (AGCC) software. Each array file was then analyzed using Transcriptome Analysis Console v4.0.2 software (ThermoFisher Scientific) for differential gene expression.

 Figure 1 

Flow chart of protocols for RNA isolation from articular cartilage. Two methods with different workflow modifications were tested for isolation of quality RNA from articular cartilage for cDNA microarray assay for gene expression. Normal rabbit cartilage tissue was collected from distal femoral surface of knee joint and either snap-frozen or suspended in 350 µl and stored at -80^oC. The detailed protocols are described in Table 1.

Results

RNA quantity and quality

Yield and quality of RNA extracted from approximately 2-3 mg of rabbit cartilage using the different methods described is shown in Table 2. Based on RNA yield and Bioanalyzer RIN scores Method 1a performed best among all the methods tested and generated good quality RNA from normal cartilage. Therefore, subsequent RNA isolation for all the samples from repaired cartilage was carried out using Method 1a (Figure 1 and Table 1) which includes a Proteinase K treatment and QIAshredder homogenization.

Table 3 summarizes the quality and purity of RNA isolated from samples at 8 weeks post microfracture supplemented with scSOX9 or controls. Adequate yields of RNA for microarray assay were isolated from individual samples ranging in wet weight from 1.9 to 2.9 mg. Concentrations measured by UV absorbance were between 3.5 - 87.4 ng/µl and were generally consistent with Bioanalyzer estimated concentrations. The optical density (OD) 260/280 ratios ranged from 1.7- 2.0, indicating acceptable RNA purity for microarray assay.

The Bioanalyzer electropherogram traces of all samples displayed discrete ribosomal RNA peaks at 18S and 28S with several additional peaks between the 18S and 28S bands and just below the 18S band (Figure 2). RIN scores ranged from 6.5 - 8.3, indicating the extracted RNA was largely intact.

 Table 2 

Comparison of yield and quality of RNA isolated from rabbit cartilage using modified methodologies*

RNA Quality Comparison	Method 1a†	Method 1b	Method 1c	Method 2a	Method 2b	Method 2c
Method Details	Trizol/RNeasy Hybrid with Proteinase K and QIAshredder	Trizol/RNeasy Hybrid with Proteinase K	Trizol/RNeasy Hybrid without Proteinase K	Modified RNaqueous with Proteinase K and QIAshredder	Modified RNaqueous with Proteinase K	Modified RNaqueous without Proteinase K
Number of samples tested	24	9	1	1	2	2
UV absorbance 260/280 (average)	1.9	1.8	2.9	2.2	0.0‡	5.2‡
UV absorbance 260/230 (average)	1.2	0.9	1.8	1.0	0.0‡	0.0‡
Bioanalyzer Concentration (ng/μl) (average)	58	38	8.1	0.8	0.0	0.1
Bioanalyzer Total Yield (μg) (average)	1.7	1.1	0.3	0.03	0.0	0.0
Bioanalyzer RIN Score (average)	7.1	5.2	Not applicable	2.2	Not applicable	Not applicable
Representative Bioanalyzer Trace

*Normal rabbit articular cartilage (3 mg per sample) was used for RNA isolation in each method; ^†Method chosen for microarray analysis; ^‡UV data available for only one sample.

 Figure 2 

Quality assessment of RNA extracted from repaired cartilage. Bioanalyzer traces of RNA of selected samples from normal cartilage, collagen membrane only (Membrane), scSOX9 or scSOX9-A76E treated cartilage and electrophoresis graphs are shown. L: RNA ladder; normal cartilage (sample 5, 7, 8); Membrane only (sample 1, 11, 12); scSOX9 (sample 4, 9, 10); scSOX9-A76E (sample 2, 3, 6). M: RNA isolated from mouse brain is shown for comparison.

Quality assessment of microarray analysis

cDNA target yields prior to array hybridization were 9.4 ±1.552 µg (mean ± SD, n = 12) demonstrating that all RNA samples were capable of generating expected yields of cDNA for hybridization. Assessment of the sample hybridization performance metrics generated with Affymetrix GeneChip Command Console (AGCC) processing of array data indicated that RNA input quality was good and array hybridizations performed well according to GeneChip standard quality thresholds (9). Positive and negative Area Under the Curve (AUC) [13] measures how well the probe set signals separate the positive controls from the negative controls and is a useful indicator of overall array performance with each sample. AUC values of our samples were 0.89 ± 0.004 (mean ± SD, n = 12). AUC values approaching 0.9 and above are considered indicative of good sample performance [14].

Differential gene expression

The purpose of our study was to compare the gene expression profiles of cartilage regeneration induced by scSOX9 compared with mutant scSOX9-A76E and collagen membrane only and in reference to normal cartilage. The scSOX9 treated samples resulted in a total of 177 significantly differentially expressed genes (≤ -2 or ≥ 2 fold change, p < 0.05) compared with normal cartilage; 36 genes compared with membrane only and 37 genes compared with loss-of-function mutant scSOX9-A76E (Figure 3). Interestingly, the most prominently upregulated genes in scSOX9 vs. scSOX9-A96E were related to extracellular matrix homeostasis such as Beta-1,4-Galactosyltransferase 6 (B4GALT6) and C-type lectin domain family 4 member A (CLEC4A), bone differentiation such as gremlin 1, DAN family BMP antagonist (GREM1) and C-type lectin domain containing 5A (CLEC5A), and inflammation (interleukin-15 and CXCL10). More analyses of differential gene expression in specimens of repaired cartilage at 4, 8 and 12 weeks post microfracture are under way.

 Figure 3 

Changes of gene expression at 8 weeks post microfracture. Microarray assay detected differently expressed genes during cartilage repair induced by scSOX9 supplemented microfracture. (A) Total number and up-regulated and down-regulated genes with changes of ≤ - 2 or ≥ 2 folds (p < 0.05) comparing different treatment. (B) Volcano plots with highlight of changes of genes induced scSOX9 compared with collagen membrane only or with scSOX9-A76E.

Discussion

Direct isolation of sufficient amounts of high quality RNA from articular cartilage has been challenging largely due to the inherit nature of cartilage. Typically large amounts of cartilage tissue are required to achieve good RNA yields [15, 16]. Highly cross-linked proteoglycans present in cartilage can also interfere with RNA purity. This has been circumvented by extracting RNA from isolated chondrocytes derived from cartilage [17]. A recently published method successfully purified RNA with an average RIN value of 7.9 ± 03 from cartilage of a human osteoarthritis specimen. The quality of RNA isolated using this method met the requirement for RNA sequencing but required a 100 mg of cartilage [15]. Others have reported methods for isolation of quality RNA from cartilage with starting material as low as 25 mg in quantity from young cartilage [8]. Gehrsitz et al reported that 10 - 40 mg cartilage plugs obtained from human autopsy specimen was sufficient for isolation of quality RNA for RT-PCR [18]. Though high quality RNA is directly extracted from cartilage tissue, none of these protocols suits the requirement of our cartilage specimens due to the small size of the focal injury area. Our method took advantage of commercially available RNA isolation kits with specific modifications that yielded a sufficient amount of quality RNA for microarray analysis using micro quantities of cartilage tissue. In particular, the addition of Proteinase K digestion followed by QIAshredder homogenization yielded a high quantity of RNA with good quality. The RIN of RNA isolated from this method ranged from 6.5 - 8.3 and in 50% of the samples, an RIN score of at least 7 (as recommended for RNA sequencing) was achieved [19]. As evidenced by these results, it is our belief that complete digestion and homogenization is crucial for the release of RNA from this highly cross-linked protein-rich tissue source. This method is particularly suitable for RNA extraction from articular cartilage in osteoarthritis or cartilage injury and repair models with small animals where a large quantity of specimen is not available. This method may also be adopted for RNA isolation from cartilage biopsy specimens in human osteoarthritis studies.

SOX9 is the master transcription factor required for chondrogenesis and orchestrates the expression of many genes during this process [20]. However, much of our knowledge about SOX9 action during chondrogenesis is learnt from in vitro studies on cultured chondrocytes [20] and may not fully represent in vivo mechanisms. Using the recombinant scSOX9 in combination with microfracture, we have shown that scSOX9 is able to induce hyaline-like cartilage regeneration in a rabbit acute cartilage injury model [1]. This model offered a unique opportunity to understand the interaction of SOX9 with other molecules in vivo during cartilage repair by inducing bone marrow derived mesenchymal stem cells (MSCs). Since the target genes for scSOX9 during the process of cartilage repair are largely unknown, we employed microarray analysis for gene expression which allows simultaneous examination of over 20,000 genes [21]. The use of microarray technology to study thousands of genes involved in cartilage regeneration has not been reported, although RNA-sequencing has been recently used to study isolated chondrocytes [22, 23]. We, therefore, made an effort to modify methods to isolate quality RNA from a micro amount of cartilage specimen to satisfy the requirement for microarray assay. In our preliminary analysis, scSOX9 induced cartilage regeneration was associated with genes involved in cell proliferation and differentiation, inflammation, energy metabolism, cell survival, and extracellular matrix protein homeostasis.

Limitations of our study must be noted. First, good quality RNA as judged by RIN was isolated from micro quantity of cartilage by our modified method. Although it is widely used for assessment of total RNA integrity, RIN does not necessarily reflect the integrity of mRNA [12]. For example, inconsistencies between RIN and corresponding RNA electropherogram profiles and lack of correlation of RIN with RNA concentration were observed in postmortem human brain tissue were observed [24]. Second, bulk RNA isolation followed by microarray analysis cannot precisely assign the differentially expressed genes to specific cell types. Single cell RNA-sequencing may be helpful in achieving that end and is being considered for future studies. Third, the cartilage samples obtained from repaired cartilage tissue may not represent a pure chondrocyte cell type but may also contain MSCs which are in transition status. Forth, using collagen membrane to carry scSOX9 in our experimental model may not be the ultimate method to apply for human cartilage repair; and collagen membrane might have impact on scSOX9 induced gene expression.

Acknowledgements

The project was supported by VA Merit Review grant, IBX002858 and partially by a grant from Vivoscript, Inc. RNA isolation and microarray assays were performed in the OHSU Gene Profiling Shared Resource.

Author Contributions

All authors critically revised the manuscript for intellectual content and provided final approval of the version submitted for publication.

Study conception and design: Yong Zhu, Christina Harrington, Cong-Qiu Chu; Acquisition of data: Xiaowei Zhang, Trevor J. McFarland, Kristina Vartanian; Analysis and interpretation of data: Yong Zhu, Kristina Vartanian, Christina Harrington, Cong-Qiu Chu.

Competing Interests

YZ is employee of Vivoscript, Inc. which produced scSOX9 used in the previous and ongoing studies. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1. Zhang X, Wu S, Naccarato T, Prakash-Damani M, Chou Y, Chu CQ. et al. Regeneration of hyaline-like cartilage in situ with SOX9 stimulation of bone marrow-derived mesenchymal stem cells. PLoS One. 2017;12:e0180138

2. Zhang X, Wu SL, Zhu Y, Chu CQ. Long term durable repaired cartilage Induced by SOX9 in situ with bone marrow-derived mesenchymal stem cells. Internation Journal of Medical Research. 2021;18:1399-405

3. Carney SL, Muir H. The structure and function of cartilage proteoglycans. Physiol Rev. 1988;68:858-910

4. Aigner T, Vornehm SI, Zeiler G, Dudhia J, von der Mark K, Bayliss MT. Suppression of cartilage matrix gene expression in upper zone chondrocytes of osteoarthritic cartilage. Arthritis Rheum. 1997;40:562-9

5. Cs-Szabo G, Melching LI, Roughley PJ, Glant TT. Changes in messenger RNA and protein levels of proteoglycans and link protein in human osteoarthritic cartilage samples. Arthritis Rheum. 1997;40:1037-45

6. Melching LI, Cs-Szabo G, Roughley PJ. Analysis of proteoglycan messages in human articular cartilage by a competitive PCR technique. Matrix Biol. 1997;16:1-11

7. Wilson T, Kaur N, Davis J, Ali SA. Tissue Collection and RNA Extraction from the Human Osteoarthritic Knee Joint. J Vis Exp. 2021

8. Zheng H, McAlinden A. RNA Isolation from Articular Cartilage Tissue. Methods Mol Biol. 2021;2245:121-33

9. Harrington CA, Winther M, Garred MM. Use of bioanalyzer electropherograms for quality control and target evaluation in microarray expression profiling studies of ocular tissues. J Ocul Biol Dis Infor. 2009;2:243-9

10. Vangsness CT Jr, Farr J 2nd, Boyd J, Dellaero DT, Mills CR, LeRoux-Williams M. Adult human mesenchymal stem cells delivered via intra-articular injection to the knee following partial medial meniscectomy: a randomized, double-blind, controlled study. J Bone Joint Surg Am. 2014;96:90-8

11. Mueller O, Lightfoot S, Schroeder A. RNA integrity number (RIN) - standardization of RNA quality control 2016. p. 1-8.

12. Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M. et al. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006;7:3

13. Burdick JA, Mauck RL, Gerecht S. To Serve and Protect: Hydrogels to Improve Stem Cell-Based Therapies. Cell Stem Cell. 2016;18:13-5

14. AffymetrixGeneChip®. Quality Assessment of Exon and Gene Arrays Affymetrix GeneChip® Gene and Exon Array Whitepaper Collection. 2007; http://tools.thermofisher.com/content/sfs/brochures/exon_gene_arrays_qa_whitepaper.pdf.

15. Le Bleu HK, Kamal FA, Kelly M, Ketz JP, Zuscik MJ, Elbarbary RA. Extraction of high-quality RNA from human articular cartilage. Anal Biochem. 2017;518:134-8

16. McKenna LA, Gehrsitz A, Soder S, Eger W, Kirchner T, Aigner T. Effective isolation of high-quality total RNA from human adult articular cartilage. Anal Biochem. 2000;286:80-5

17. Ali SA, Alman B. RNA extraction from human articular cartilage by chondrocyte isolation. Anal Biochem. 2012;429:39-41

18. Gehrsitz A, McKenna LA, Soder S, Kirchner T, Aigner T. Isolation of RNA from small human articular cartilage specimens allows quantification of mRNA expression levels in local articular cartilage defects. J Orthop Res. 2001;19:478-81

19. Li S, Tighe SW, Nicolet CM, Grove D, Levy S, Farmerie W. et al. Multi-platform assessment of transcriptome profiling using RNA-seq in the ABRF next-generation sequencing study. Nat Biotechnol. 2014;32:915-25

20. Liu CF, Lefebvre V. The transcription factors SOX9 and SOX5/SOX6 cooperate genome-wide through super-enhancers to drive chondrogenesis. Nucleic Acids Res. 2015;43:8183-203

21. Minnier J, Pennock ND, Guo Q, Schedin P, Harrington CA. RNA-Seq and Expression Arrays: Selection Guidelines for Genome-Wide Expression Profiling. Methods Mol Biol. 2018;1783:7-33

22. Duan X, Cai L, Schmidt EJ, Shen J, Tycksen ED, O'Keefe RJ. et al. RNA-seq analysis of chondrocyte transcriptome reveals genetic heterogeneity in LG/J and SM/J murine strains. Osteoarthritis Cartilage. 2020;28:516-27

23. Kelly NH, Huynh NPT, Guilak F. Single cell RNA-sequencing reveals cellular heterogeneity and trajectories of lineage specification during murine embryonic limb development. Matrix Biol. 2020;89:1-10

24. Sonntag KC, Tejada G, Subburaju S, Berretta S, Benes FM, Woo TU. Limited predictability of postmortem human brain tissue quality by RNA integrity numbers. J Neurochem. 2016;138:53-9

Author contact

Corresponding author: chucedu

Received 2021-7-25
Accepted 2021-10-19
Published 2022-1-1

Citation styles

©2024 Ivyspring International Publisher. Terms of use