|
|
||||||||
Technical Developments |
1 From the Lucas MR Spectroscopy and Imaging Research Center (S.G., Y.S.Y., G.S., D.Y.L., M.D.B.) and Department of Radiology (S.G., Y.S.Y., G.S., D.Y.L., K.C.P.L., M.D.B.), Stanford University School of Medicine, 1201 Welch Rd, P260, Stanford, CA 94305-5488. Received July 31, 2002; revision requested September 7; revision received October 11; accepted December 10. Supported in part by the Lucas Foundation, the Phil Allen Trust, and National Institutes of Health grant P41 RR09784. Y.S.Y. and S.G. supported by a research fellowship of the National Cancer Institute at the Lucas Center, Stanford University. Address correspondence to M.D.B. (e-mail: mark@s-word.stanford.edu).
| ABSTRACT |
|---|
|
|
|---|
© RSNA, 2003
Index terms: Genes and genetics Magnetic resonance (MR), experimental studies, **.12143
| INTRODUCTION |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
The animals were sacrificed once the tumors reached the stage of progression suitable for further analysis on the basis of the MR contrast-enhancement pattern (see the Image Analysis section for details) or once the tumor began to affect the mobility of the animals. Regions at the periphery and center of the tumor were marked with dissecting pins on the basis of the identifiable anatomic landmarks on the MR images. The tumors were removed, and these regions were sampled at biopsy. Then, the specimens were frozen rapidly in liquid nitrogen. The volume of these samples from the periphery and center of the tumors ranged from approximately 50 to 200 mm3 depending on the size of the tumor at different stages of tumor progression. These samples were processed for gene expression analysis, reverse transcriptasepolymerase chain reaction, immunohistochemistry, in situ hybridization, Western blot analysis, and two-dimensional gel electrophoresis.
The animals were chosen for each technique on the basis of the following criteria: Tumors with similar T1- and T2-weighted MR image patterns were grouped together (see the Image Analysis section for details). Tissue samples from different tumors in the same group were used for these techniques. The sample sizes for each of the protocols were determined on the basis of previous literature precedents (1822).
MR Imaging
A clinical 3.0-T MR imager (Signa Horizon LX8.2; GE Medical Systems, Milwaukee, Wis) was used with a custom-designed quadrature high-pass birdcage coil tuned to 127 MHz for signal reception. Fifty-four animals were imaged. The following MR imaging protocols were performed in the transverse plane with a field of view of 8 cm, 256 x 192-pixel matrix, section thickness of 1.5 mm, and one signal acquired: (a) precontrast T1-weighted spin echo (repetition time msec/echo time msec of 300/13), (b) precontrast T2-weighted fast spin echo (4,000/85, echo train length of eight), and (c) postcontrast T1-weighted spin echo (300/13) at 2 minutes after injection of gadopentetate dimeglumine (200 µL of 0.5 mol/L solution administered via a tail vein at a rate of 40 µL/sec with no flushing solution). Postprocessing was performed with the software package in the MR imager.
Image Analysis
Images were reviewed by two clinicians (K.C.P.L., D.Y.L.). The reviewers were blinded to information about the stage of tumor progression. The animals were then categorized into three groups on the basis of the MR images: group 1 (18 mice) had homogeneously enhanced postcontrast T1-weighted MR images; group 2 (18 mice) had heterogeneously enhanced postcontrast T1-weighted MR images; and group 3 (18 mice) had heterogeneously enhanced postcontrast T1-weighted MR images, with elevated signal intensity in the T2-weighted MR images that corresponded to the contrast-enhanced areas in the T1-weighted images. The relative tumor pixel intensities of the center and peripheral areas on the T2-weighted images were analyzed (MRVISION, version 1.5.4b; MRVision, Menlo Park, Calif). For groups 1 and 2, the mean relative pixel intensity was 1.0 ± 0.1 (SD). However, the center areas with elevated signal intensity in the more advanced tumors on heterogeneous T2-weighted MR images had higher pixel intensities relative to those in the peripheral region of the tumor (1.6 ± 0.1) (group 3). The former tumors were excluded from the current study.
Hematoxylin-Eosin Staining
Two mice from each group were used for hematoxylin-eosin staining. Tissue samples were preserved in 10% formalin solution for 2448 hours. Then they were embedded in paraffin, sectioned, stained with hematoxylin-eosin, and mounted on glass slides. Morphologic features used to identify viable nuclei included the absence of overall shrinkage and homogenous dark basophilia. The slides were reviewed by a veterinary pathologist at Stanford University Medical School.
Oligonucleotide Microarray Analysis
Four mice were used from each group for microarray analysis. Total RNA was isolated (GStract Total RNA Isolation Kit; Maxim Biotech, South San Francisco, Calif). The procedure involves (a) homogenizing tissue or cells, (b) adding salt solution to precipitate the DNA and RNA, (c) extracting the RNA from the precipitate by using phenol chloroform, and (d) washing the RNA pellet to eliminate all contaminants. The polyadenylated messenger RNA (mRNA) was then purified from the total RNA (Oligotex mRNA Kit; Qiagen, Valencia, Calif). The purified mRNA was handled for analysis by using the GeneChip Expression Analysis protocol (Affymetrix, Santa Clara, Calif): (a) The preparation was targeted by synthesizing double-stranded complementary DNA (cDNA) from purified mRNA by using reverse transcriptase. (b) The purified cDNA was used to produce biotin-labeled complementary RNA (cRNA) by means of in vitro transcription reaction. (c) The target biotin-labeled cRNA was fragmented and hybridized with the control cRNA to the oligonucleotide probes on the probe array for 16 hours at 45°C. (d) The probe array was washed and stained with streptavidin phycoerythrin conjugate immediately after hybridization. (e) The probe array was scanned with a confocal microscope scanner (Hewlett-Packard, Palo Alto, Calif) at the excitation wavelength of 488 nm. The amount of light emitted at 570 nm was proportional to that of the bound target at each location on the probe array. Algorithms (GeneChip Expression Analysis; Affymetrix) were used to analyze data generated from the oligonucleotide microarrays (Mu11subA and Mu11subB; Affymetrix) (total of 11,000 genes) to determine the intensity of expressed genes. Data were analyzed (Microarray Suite, version 4.0, Affymetrix; Gene Spring, version 4.0, Silicon Genetics, Redwood City, Calif).
Quantitative Reverse TranscriptasePolymerase Chain Reaction Analysis
With three mice from each group, quantitative reverse transcriptasepolymerase chain reaction was performed according to standard procedures. RNA was extracted from tumor tissues with the following steps: Tissue was homogenized in 5 mL of lysis buffer (6 mol/L urea, 3 mol/L lithium chloride, 50 mmol/L sodium acetate, 200 µg/mL heparin, and 0.1% sodium dodecyl sulfate). The homogenized tissue was centrifuged at a rate of 13,100 rpm (16,000 g) for 20 minutes, extracted twice with an equal volume of phenol and chloroform, and precipitated with ethanol. The RNA pellet was air dried and dissolved in water treated with diethylpyrocarbonate. Reverse transcriptasepolymerase chain reaction was performed (DNA Thermal Cycler 480; Perkin-Elmer; Oak Brook, Ill). The chain reaction involved incubation of the samples at (a) 42°C for 15 minutes, 99°C for 5 minutes, and 5°C for 5 minutes for one cycle; (b) 95°C for 45 seconds and 60°C for 45 seconds for 35 cycles; and (c) 72°C for 7 minutes for one cycle. The products were checked in 2% agarose gel, along with the 100base pair ladder (Promega, Madison, Wis).
Statistical Analysis
Statistical analysis was performed by means of the Student t test. Differences with a P value of less than .05 were considered statistically significant. This analysis was performed by comparing the relative gene expression levels (fold change [Table 1]) of genes in the homogeneously enhanced tumors with those same genes in the heterogeneously enhanced tumors.
|
Two-dimensional Electrophoresis Argorase Gel and Sequencing
Three mice from each group were used for two-dimensional gel electrophoresis and sequencing. Gel electrophoresis was performed with standard protocol (23): Cytosolic fractions of the periphery (contrast-enhanced) and center (nonenhanced) regions of the tumor lysate (250 µg) were subjected to isoelectric focusing on impedance plethysmographic strips (isoelectric point at pH of 47); they were separated by 4%20% at polyacrylamide gel electrophoresis with sodium dodecyl sulfate. The gels were then stained with silver. The band of interest (marked by arrow on the right side in Fig 1, B) was digested with trypsin. The resulting fragments were analyzed by means of matrix-assisted laser desorption ionization time-of-flight mass spectrometry to determine their molecular weights. These fragments were identified as MGP (Profound Peptide Mapping software, version 4.10.5; Rockefeller University, New York, NY).
|
Western Blot Analysis
Two mice from each group were used for Western blot analysis, which was performed according to standard protocols (24). Membrane fractions of the peripheral and center regions of the tumor lysate (50 µg) were mixed with sample-loading buffer, heated for 5 minutes at 95°C in 50 mmol/L dithiothreitol, separated at polyacrylamide gel electrophoresis with sodium dodecyl sulfate, and transferred to nitrocellulose. The nitrocellulose filters were blocked with 5% reconstituted dry milk and incubated with an antiPDGF-receptor
chain antibody. Bound antibodies were visualized by using a chemiluminescence detection kit (Amersham Biosciences, England). The quantity of the expression level was determined by means of densitometric analysis.
| Results |
|---|
|
|
|---|
|
Hematoxylin-Eosin Staining
With hematoxylin-eosin staining, we observed no distinguishable difference between the periphery and center regions of the homogeneously enhanced tumors (Fig 3, A, B) or between the contrast-enhanced and nonenhanced regions of the heterogeneously enhanced tumors (Fig 3, C, D). At this stage of tumor growth, neither region showed any necrosis or substantial inflammatory infiltrates (Fig 2, C), as was reflected in the T2-weighted MR image (Fig 2, E). In contrast, hematoxylin-eosin staining for the nonenhanced region (Fig 3, F) in heterogeneously enhanced tumor showed a different pattern compared with those in the enhanced area (Fig 3, E). Significant generalized necrosis is evident with the presence of inflammatory infiltrates and karyorrhexis.
|
cysteine protease (28), extracellular matrix-associated protein Sc1 (29), insulin-like growth factor binding protein-3 (30), and apolipoprotein E (31). Three additional genes (
-globin, ß-globin complex, and ß-1-globin) are associated with globin synthesis and hypoxia (Table 1, entries 79) (32,33), and one gene is associated with vasculogenesis (PDGF receptor: Table 1, entry 10) (34). We observed one elevated mRNA level in the nonenhanced region that corresponded to mast cell protease-12 (Table 1, entry 11) (35). Interestingly, many genes such as fibronectin were highly up-regulated in the tumor but were also spatially homogeneous throughout the mass (Table 1, entry 12, compare mean intensity, contrast enhanced [3,822] and nonenhanced [3,639] with fold change [x1.08]) (36).
To ensure that differences in gene expression profiles are related to the changes we observed at MR imaging, we also analyzed four tumors that were homogeneously enhanced in the postcontrast MR images. In this case, tissue samples were taken from the periphery and center of the tumor, similar to the regions sampled at biopsy in the heterogeneously enhanced tumors (Fig 2, C, thin and thick arrows); the samples were analyzed in the same way as were the heterogeneous tumors. The only gene that was significantly up-regulated (more than twofold) in the homogeneously enhanced tumors was
-actin (x3.2 ± 0.4).
Statistical Analysis
A Student t test was performed to compare the mean fold change in gene expression profile between the homogeneously and heterogeneously enhanced tumor sets (Table 2). The results showed a statistically significant difference (P < .05) for MGP, PDGF receptor, and apolipoprotein E. No statistical significance was observed between the mean fold change of fibronectin in the homogeneously and heterogeneously enhanced tumors (P = .483).
|
Target Validation
To further confirm the heterogeneous spatial distribution of these genes, we performed in situ hybridization studies, immunohistochemical analysis, two-dimensional gel electrophoresis, protein sequencing, and Western blot analysis on MGP and PDGF receptor.
MGP: In Situ Hybridization, Two-dimensional Gel Electrophoresis, and Protein Sequence
We performed in situ hybridization for MGP, which showed the highest difference in mRNA expression in both microarray and quantitative reverse transcriptasepolymerase chain reaction analysis. As demonstrated in Figure 1, A, the periphery of the tumor showed intense chromogen staining compared with that in the center of the tumor, which reflects a higher degree of probe hybridization in the periphery. This pattern of mRNA up-regulation in the tumor periphery was consistent with findings on the corresponding contrastenhanced MR image (Fig 2, F). Next, protein purification and sequence identification was performed for MGP by means of mass spectroscopy. This approach was used as opposed to a Western blot analysis because MGP may be distinctly altered by
-carboxylation, which makes antibody-based detection less reliable. These forms of MGP were resolved at two-dimensional gel electrophoresis (Fig 1, B). A series of bands at 11,000 Da were observed from three distinct tissue samples obtained from the contrast-enhanced region of the tumors. These bands represented similar molecular weights but different isoelectric points consistent with the expected carboxylation state of MGP (Fig 1, B, right). The bands and patterns were distinctly different for the tissue obtained from the nonenhanced region in the same tumor (Fig 1, B, left). The most intense band present in the enhanced region of the tumor (arrow on Fig 1, B, left) was isolated and digested with trypsin. The resulting fragments were analyzed by means of matrix-assisted laser desorption ionization time-of-flight mass spectrometry and were identified as mouse MGP.
PDGF Receptor: Immunohisto-chemical and Western Blot Analysis
For PDGF receptor, we performed immunohistochemical analysis by means of a biotinylated monoclonal antibody (BD Bioscience) on tissue sections. In the heterogeneously enhanced tumors, the results showed a distinctly darker stain in the periphery of the tumor compared with that in the center, which represented chromogen deposition related to the presence of PDGF-receptor
(Fig 4, A). The staining pattern corresponded precisely with areas of contrast enhancement seen at MR imaging (Fig 2, F). To further confirm that the differences in mRNA expression in the contrast-enhanced regions of the tumor are related to protein expression, Western blot analysis was performed to determine the level of protein expression. Membrane fractions of the peripheral and central regions of the tumor were analyzed by means of immunoblot analysis with the antiPDGF-receptor
chain antibody. Results of densitometric analysis showed that PDGF receptor was up-regulated in the peripheral region by x2.7 ± 0.2 (n = 3) (Fig 4, B).
|
| Discussion |
|---|
|
|
|---|
MGP is an extracellular protein that is dependent on vitamin-K
-carboxylation. Findings in previous studies identified MGP in metastatic breast cancer cell lines (26); therefore, it may be involved in the process of metastasis in cancer (22). The overproduction of PDGF is associated with autocrine and paracrine growth stimulation of human tumors, and activation of their corresponding receptors may be critically important in tumor progression (34). Both PDGF and its receptors are implicated in the process of angiogenesis and neovascular maturation (34). Fibronectin, an extracellular matrix protein, is a ligand for the integrin family of receptors and functions in the mediation of cell adhesion (36). Evidence from both genomic analysis and immunohistologic examination show correlations between fibronectin, tumorigenesis, and metastasis (22). Apolipoprotein E is believed to be synthesized and secreted by astrocytes and has been of great interest as a participant in the pathogenesis of Alzheimer disease. This gene is expressed in a variety of tumors, but its role in tumor biology has not been defined (31).
Quantitative reverse transcriptase polymerase chain reaction was performed for these four genes, and the results were consistent with the results at oligonucleotide microarray analysis. MGP, apolipoprotein E, and PDGF receptor all had a higher level of expression in the contrast-enhanced area compared with the nonenhanced area, while fibronectin showed a similar expression level in both regions. Thus, we conclude that the findings at quantitative reverse transcriptasepolymerase chain reaction supported the microarray data: The levels of MGP, PDGF receptor, and apolipoprotein E are associated with the differential MR contrast-enhancement pattern.
In addition, the up-regulation of MGP in the contrast-enhanced regions of the tumor was confirmed with in situ hybridization, two-dimensional electrophoresis, and protein sequencing. On the basis of these results, we conclude that MGP was expressed differentially and was present mostly in the contrast-enhanced region of the tumor. Immunohistochemical examination and Western blot analysis were performed to demonstrate the up-regulation of PDGF receptor in the contrast-enhanced area of the tumor. We conclude that the level of PDGF receptor
was more than twofold the protein level in the contrast-enhanced region relative to the nonenhanced region in the heterogeneously enhanced tumors.
Thus, we hypothesize that strongly contrast-enhanced regions in tumors are associated with changes in the extracellular matrix (the "milieu") within the environment of the tumor. This matrix seems to contain factors that increase permeability to small molecules, such as gadopentetate dimeglumine, or has properties to retain these molecules long enough for them to be observed with standard clinical imaging sequences. This milieu is probably a combination of proteins from both the tumor cells themselves and the cells associated with the extracellular matrix. In the present study, we showed that the mRNA of proteins associated with metastatic cell lines and tumors, such as MGP and PDGF receptor, are up-regulated in the contrast-enhanced areas. Thus, the environment to which tumor cells are exposed may be as important as the cell taxonomy when tumors are staged.
New imaging techniques can also provide discrimination of tissue samples for genomic analysis. With dynamic contrast-enhanced MR imaging, the kinetics of contrast enhancement can be studied to provide more parameters for characterizing tissues in vivo (13). Similarly, MR imaging can be performed with tissue-specific contrast agents and molecular imaging techniques that provide even more information about the tissue of interest (37,38). In addition, clinical imaging modalitiessuch as computed tomography (CT), positron emission tomography, and single photon emission CTcan provide new ways of correlating imaging findings with genomic analysis findings. In addition, optical imaging methods combined with laser-capture microdissection can provide discrimination to the cellular level of samples obtained at image-guided biopsy (39).
In conclusion, we demonstrated that contrast-enhanced MR imaging can serve as a noninvasive tool for characterizing different regions of tumors and can provide correlates for genomic analysis with high spatial and temporal resolution. The differential expression of genes and proteins related to the heterogeneity within solid tumors as depicted at MR imaging has important consequences. First, an in vivo evaluation with noninvasive imaging of the tumor in its unperturbed state can guide tissue sampling to areas of the tumor with distinct pathophysiology. Second, by comparing the differences in gene and protein expression between spatially distinct regions, important targets can be identified that may otherwise be overlooked. Finally, since the clinical imaging protocols used in this study are implemented routinely in cancer patients, they can be easily added as a parameter to annotate gene expression analysis. MR imaging can be used as a noninvasive tool to characterize the distinct pathophysiology within cancers for genomic and proteomic analysis.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Abbreviations: MGP = matrix Gla protein, mRNA = messenger RNA, PDGF = platelet-derived growth factor
Author contributions: Guarantors of integrity of entire study, K.C.P.L., M.D.B.; study concepts, K.C.P.L., M.D.B.; study design, K.C.P.L., M.D.B., S.G., Y.S.Y.; literature research, S.G., Y.S.Y.; experimental studies, S.G., Y.S.Y., G.S.; data acquisition and analysis/interpretation, S.G., Y.S.Y., G.S.; statistical analysis, S.G., Y.S.Y.; manuscript preparation, S.G., Y.S.Y., M.D.B.; manuscript definition of intellectual content, K.C.P.L., M.D.B., S.G., Y.S.Y., D.Y.L.; manuscript editing, all authors; manuscript revision/review, K.C.P.L., M.D.B., S.G., Y.S.Y., D.Y.L., G.S.; manuscript final version approval, S.G., Y.S.Y., M.D.B.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R. G. Pautler Mouse MRI: Concepts and Applications in Physiology Physiology, August 1, 2004; 19(4): 168 - 175. [Abstract] [Full Text] [PDF] |
||||
![]() |
K C P Li Angiogenesis imaging in the post-genomic era Br. J. Radiol., December 1, 2003; 76(suppl_1): S1 - S2. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |