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

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.

A gene expression signature of CD34+ cells to predict major cytogenetic response in chronic-phase chronic myeloid leukemia patients treated with imatinib.

In chronic-phase chronic myeloid leukemia (CML) patients, the lack of a major cytogenetic response (< 36% Ph(+) metaphases) to imatinib within 12 months indicates failure and mandates a change of therapy. To identify biomarkers predictive of imatinib failure, we performed gene expression array profiling of CD34(+) cells from 2 independent cohorts of imatinib-naive chronic-phase CML patients. The learning set consisted of retrospectively selected patients with a complete cytogenetic response or more than 65% Ph(+) metaphases within 12 months of imatinib therapy. Based on analysis of variance P less than .1 and fold difference 1.5 or more, we identified 885 probe sets with differential expression between responders and nonresponders, from which we extracted a 75-probe set minimal signature (classifier) that separated the 2 groups. On application to a prospectively accrued validation set, the classifier correctly predicted 88% of responders and 83% of nonresponders. Bioinformatics analysis and comparison with published studies revealed overlap of classifier genes with CML progression signatures and implicated beta-catenin in their regulation, suggesting that chronic-phase CML patients destined to fail imatinib have more advanced disease than evident by morphologic criteria. Our classifier may allow directing more aggressive therapy upfront to the patients most likely to benefit while sparing good-risk patients from unnecessary toxicity.

Gene expression profiling of whole blood: comparison of target preparation methods for accurate and reproducible microarray analysis.

BACKGROUND: Peripheral blood is an accessible and informative source of transcriptomal information for many human disease and pharmacogenomic studies. While there can be significant advantages to analyzing RNA isolated from whole blood, particularly in clinical studies, the preparation of samples for microarray analysis is complicated by the need to minimize artifacts associated with highly abundant globin RNA transcripts. The impact of globin RNA transcripts on expression profiling data can potentially be reduced by using RNA preparation and labeling methods that remove or block globin RNA during the microarray assay. We compared four different methods for preparing microarray hybridization targets from human whole blood collected in PAXGene tubes. Three of the methods utilized the Affymetrix one-cycle cDNA synthesis/in vitro transcription protocol but varied treatment of input RNA as follows: i. no treatment; ii. treatment with GLOBINclear; or iii. treatment with globin PNA oligos. In the fourth method cDNA targets were prepared with the Ovation amplification and labeling system. RESULTS: We find that microarray targets generated with labeling methods that reduce globin mRNA levels or minimize the impact of globin transcripts during hybridization detect more transcripts in the microarray assay compared with the standard Affymetrix method. Comparison of microarray results with quantitative PCR analysis of a panel of genes from the NF-kappa B pathway shows good correlation of transcript measurements produced with all four target preparation methods, although method-specific differences in overall correlation were observed. The impact of freezing blood collected in PAXGene tubes on data reproducibility was also examined. Expression profiles show little or no difference when RNA is extracted from either fresh or frozen blood samples. CONCLUSION: RNA preparation and labeling methods designed to reduce the impact of globin mRNA transcripts can significantly improve the sensitivity of the DNA microarray expression profiling assay for whole blood samples. While blockage of globin transcripts during first strand cDNA synthesis with globin PNAs resulted in the best overall performance in this study, we conclude that selection of a protocol for expression profiling studies in blood should depend on several factors, including implementation requirements of the method and study design. RNA isolated from either freshly collected or frozen blood samples stored in PAXGene tubes can be used without altering gene expression profiles.

Mapping quantitative trait loci that influence femoral cross-sectional area in mice.

Size and shape are critical determinants of the mechanical properties of skeletal elements and can be anticipated to be highly heritable. Moreover, the genes responsible may be independent of those that regulate bone mineral density (BMD). To begin to identify the heritable determinants of skeletal geometry, we have examined femoral cross-sectional area (FCSA) in male and female mice from two inbred strains of mice with divergent FCSA (C57BL/6 [B6] and DBA/2 [D2]), a large genetically heterogeneous population (n = 964) of B6D2F2 mice and 18 BXD recombinant inbred (RI) strains derived from their F2 cross. Femora were harvested from 16-week-old mice and FCSA (bone and marrow space enclosed within the periosteum) was measured at the midshaft by digital image analysis. In all mouse populations examined, FCSA was positively correlated with body weight and weight-corrected FCSA (WC-FCSA) values were normally distributed in the BXD-RI and F2 populations, suggesting polygenic control of this trait. Genome-wide quantitative trait locus (QTL) analysis of the B6D2F2 population revealed regions on four different chromosomes that were very strongly linked to WC-FCSA (chromosomes 6, 8, 10, and X) in both genders. Evidence of gender-specific genetic influences on femoral geometry was also identified at three other chromosomal sites (chromosomes 2, 7, and 12). Supporting evidence for the WC-FCSA QTLs on chromosomes 2, 7, 8, 10, and 12 also was present in the RI strains. Interestingly, none of these WC-FCSA QTLs were identified in our previous QTL analysis of whole body BMD in the same B6D2F2 population. Thus, the genetic determinants of bone size appear to be largely, if not entirely, distinct from those that regulate BMD attainment. The identification of the genes responsible for geometric differences in bone development should reveal fundamentally important processes in the control of skeletal integrity.

Computational detection of alternative exon usage.

BACKGROUND: With the advent of the GeneChip Exon Arrays, it is now possible to extract "exon-level" expression estimates, allowing for detection of alternative splicing events, one of the primary mechanisms of transcript diversity. In the context of (1) a complex trait use case and (2) a human cerebellum vs. heart comparison on previously validated data, we present a transcript-based statistical model and validation framework to allow detection of alternative exon usage (AEU) between different groups. To illustrate the approach, we detect and confirm differences in exon usage in the two of the most widely studied mouse genetic models (the C57BL/6J and DBA/2J inbred strains) and in a human dataset. RESULTS: We developed a computational framework that consists of probe level annotation mapping and statistical modeling to detect putative AEU events, as well as visualization and alignment with known splice events. We show a dramatic improvement (∼25 fold) in the ability to detect these events using the appropriate annotation and statistical model which is actually specified at the transcript level, as compared with the transcript cluster/gene-level annotation used on the array. An additional component of this workflow is a probe index that allows ranking AEU candidates for validation and can aid in identification of false positives due to single nucleotide polymorphisms. DISCUSSION: Our work highlights the importance of concordance between the functional unit interrogated (e.g., gene, transcripts) and the entity (e.g., exon, probeset) within the statistical model. The framework we present is broadly applicable to other platforms (including RNAseq).

Randomization in laboratory procedure is key to obtaining reproducible microarray results.

The quality of gene expression microarray data has improved dramatically since the first arrays were introduced in the late 1990s. However, the reproducibility of data generated at multiple laboratory sites remains a matter of concern, especially for scientists who are attempting to combine and analyze data from public repositories. We have carried out a study in which a common set of RNA samples was assayed five times in four different laboratories using Affymetrix GeneChip arrays. We observed dramatic differences in the results across laboratories and identified batch effects in array processing as one of the primary causes for these differences. When batch processing of samples is confounded with experimental factors of interest it is not possible to separate their effects, and lists of differentially expressed genes may include many artifacts. This study demonstrates the substantial impact of sample processing on microarray analysis results and underscores the need for randomization in the laboratory as a means to avoid confounding of biological factors with procedural effects.