1 Human Immunology, Hanson Institute, Adelaide, Australia
2 Bionomics Limited, Thebarton, Australia
3 School of Mathematics, Adelaide, Australia
4 Department of Medicine, The University of Adelaide, Adelaide, Australia
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
quantification; gene expression; mitogen-activated protein kinase pathway
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Approaches to identify genes that are differentially regulated in angiogenesis include colony arrays, cDNA filter arrays, microarrays, gene-calling, SAGE, differential display, and suppression subtractive hybridization (2, 21, 22, 27, 48). These approaches have been successful in identifying genes involved in endothelial cell proliferation, migration, or differentiation, but have been limited by the fact that the majority of studies have utilized cells cultured in a restricted two-dimensional (2-D) milieu.
We utilized our in vitro three-dimensional (3-D) collagen gel model of capillary tube formation, in which the time course of morphological stages is very well characterized (20) and which recapitulates many of the critical individual processes in angiogenesis (36). Importantly, because these stages are highly synchronized, whereby the majority of cells at any one time are undergoing the same morphological rearrangements, crucial transcriptional events in the angiogenic process are demonstrable. We produced an "angiogenic" microarray chip from libraries generated using suppression subtractive hybridization. This allowed enrichment for regulated cDNAs relevant to the angiogenic process and amplification, via normalization, of a number of rare and lowly regulated mRNAs.
To date, the majority of microarray-based analyses of gene expression in angiogenic models have performed single or duplicate hybridizations per experiment and then applied a typical twofold or more cutoff (1, 2, 26, 56). However, this results in the identification of only highly regulated genes. By increasing the number of hybridizations, performing time course analyses, and applying rigorous statistical analyses, we identified lowly differentially expressed genes and entire signaling pathways, which have proven biologically relevant. Furthermore, we demonstrated the coordinated regulation of successive members of signaling pathways, in particular the mitogen-activated protein kinase (MAPK) and leukemia inhibitory factor (LIF) pathways, as an additional strength of this approach. In addition, we also showed that "virtual Northern" (VN) blot analysis is an alternative approach to quantification, being up to 2,500 times more saving on RNA than Northern blots and having benefits over quantitative RT-PCR (Q-RT-PCR) without jeopardizing sensitivity.
![]() |
MATERIALS AND METHODS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
The capillary tube formation assay was performed as previously described (20). Briefly, 6.4 x 104 cells/160 µl HUVEC medium were plated onto 100 µl gelled bovine type I collagen (Celtrix Laboratories; Palo Alto, CA) in 96-well flat-bottomed microtiter plates. Capillary tube formation was stimulated by the addition of 20 ng/ml phorbol myristate acetate (PMA) and an antibody against 2ß1-integrin (RMAC11), which promotes the formation of complex multicellular tubes.
Generation of subtracted cDNA libraries.
HUVECs from three different umbilical cords were plated onto collagen gels to induce capillary tube formation (15 wells/cord). At times of 0, 0.5, 3, 6, and 24 h after plating, total RNA was isolated using the TRIzol procedure (Invitrogen; Carlsbad, CA). These times were chosen because they approximate time points at which distinct angiogenic morphological events are known to occur (20, 36). TRIzol reagent was added directly to each well and gently agitated to aid in dissolution. RNA was extracted according to the manufacturer's protocol except that two TRIzol extractions and a further two phenol-chloroform-isoamyl alcohol (25:24:1) extractions were required to completely remove the collagen. The RNA samples, derived from tubes from each of the three HUVEC lines, corresponding to the same time points, were pooled to minimize between-subject variation. An additional precipitation was performed in which 3 vol of 4 M sodium acetate were added to the RNA, and it was pelleted at 13,000 g for 30 min. This removed other nonproteinaceous contaminants from the RNA. Yields of 713 µg of total HUVEC RNA were obtained. To achieve enough mRNA for subtraction experiments, the "switching mechanism at the 5'-end of RNA transcript" (SMART) PCR cDNA synthesis procedure (Clontech; Palo Alto) was used to convert the mRNA to cDNA and to preferentially amplify full-length cDNAs. The PCR Select cDNA Subtraction Kit (Clontech) was then used to perform suppression subtractive hybridization between each adjacent time point in both the forward and reverse orientation. Four forward and four reverse libraries were then generated by digesting the subtracted cDNAs with EagI and cloning them into the NotI site in pBluescriptIISK+ plasmids (Stratagene; La Jolla, CA). These libraries were used to transform XL-10 Gold ultracompetent bacteria (Stratagene).
Generation and probing of microarrays.
Generation of microarrays was performed by the Australian Genome Research Facility (Melbourne, Australia). Colonies from the forward subtracted libraries (00.5 h, 3,200 colonies; 0.53 h, 3,000 colonies; 36 h, 2,800 colonies; and 624 h, 1,000 colonies) were robotically picked, and glycerols stocks were generated. A total of 10,000 clones was picked. Each clone was then PCR amplified using primers to the T7 and T3 promoters in the vector that flank the cloning site. The PCR products were microarrayed onto glass slides in duplicate along with controls known to be regulated in this angiogenic model [e.g., prostaglandin endoperoxidase synthase 2 (PTGS2), matrix metalloproteinase 10 (MMP10), and connective tissue growth factor (CTGF)] together with genes predicted not to be regulated [thymosin-ß4 (TMSB4X) and heat shock protein 150 kDa (HSPH1)]. A dilution series of a pool of the 10,000 clones was also incorporated, along with cDNAs from other species as negative controls. We refer to these slides as "angiogenic" microarray chips.
Because of the small amounts of RNA available, the Message Amp aRNA Kit (Ambion; Austin, TX) was used to convert 0.5 µg of the total RNA from each time point to cDNA and then linearly amplify the mRNA 400-fold. The amplified RNA (aRNA) was generated incorporating 5-(3-aminoallyl)-UTP and subsequently coupled to Cy3 and Cy5 dyes (Amersham Biosciences; Buckinghamshire, UK).
aRNA (3.5 µg) was desiccated and resuspended in 0.1 M sodium bicarbonate (pH 9.0) for 15 min. Solubilized aRNA was coupled to fluorescent CyDyes and incubated at room temperature for 1 h. Uncoupled CyDyes were removed using a PCR Purification Kit (Qiagen). aRNA was eluted in RNase-free water (60 µl). Blocking solutions [60 µl human Cot-1 DNA (1 mg/ml, Invitrogen); 0.96 µl yeast tRNA (25 mg/ml, Sigma); 1.5 µl poly dA (8 mg/ml, Sigma); 1.5 µl of 50x Denhardt's solution; and 4 µl salmon sperm DNA (10 mg/ml, Invitrogen)] were added to the aRNA and desiccated to almost complete dryness. aRNA was resuspended in hybridization buffer (16 µl of 100% deionized formamide and 16 µl of 12.5x SSC) and denatured at 100°C for 3 min. Denatured aRNA was chilled on ice, and 0.4 µl of 10% SDS were added. The hybridization mix was then added to slide arrays that had prehybridized at 42°C in 25% formamide, 5x SSC, 0.1% SDS, and 10 mg/ml BSA in a Coplin jar for 45 min. Hybridizations were performed in a hybridization chamber (Corning) at 42°C for 16 h. The slides were then washed in 0.5x SSC and 0.01% SDS for 1 min, 0.5x SSC for 3 min, and 0.06x SSC for 3 min and rinsed once in filtered MQ water. All slides were dried in a bench centrifuge at 750 rpm for 5 min. A probing strategy in which 3.5 µg of labeled aRNA at each time point was compared with the zero time point was used. All comparisons were performed in quadruplicate including duplicate dye swaps and two independent probe amplifications.
Analysis of microarray data and gene identification.
The microarray slides were scanned using a GenePix 4000B Scanner driven by GenePix Pro 4.0 (Axon Instruments; Foster City, CA). All analyses were performed using freely available statistical programming and graphics environment R (http://cran.r-project.org). Image processing was done using the "SPOT" software package (7) with a seed-growing algorithm for spot image segmentation and a morphological opening algorithm for background subtraction. The extracted information was normalized using the "Statistics in Microarray Analysis" package (http://www.stat.berkeley.edu/terry/zarray/software/smacode.html). Within-slide normalization was based on intensity-dependent Lowess, performed for each print tip group. The data from four individual arrays used for each time point were scaled to make them directly comparable with each other. Bayesian analysis was used on the data from quadruplicate slides and duplicate spots to identify and rank clones showing consistent differential expression both within and between arrays (33). For confirmation, "moderated t-statistics" were calculated (46). In addition, profiles of fold induction versus time were generated for all pairs of spots. With the use of Bayesian analysis, moderated t-statistics, and profile analysis, the clones most likely to be regulated were selected, and the inserts were sequenced by the Australian Genome Research Facility.
Analysis of the gene sequences was performed using two programs from the freeware Staden Package (http://www.mrc-lmb.cam.ac.uk/pubseq/staden_home.html), PreGAP4 and GAP4. With the use of PreGAP4, the gene sequences were first converted from the ABI file format to SCF format. The quality of the gene sequences were analyzed (Phred), and the sequences were stripped of vector and adapter sequences. This information was then stored in a file output. GAP4 was then used to analyze and assemble the various sequences into contigs. The analysis from GAP4 produces five potential different categories including contigs, singlets (if not assembled), vector (no inserts), Escherichia coli sequences, and failed read (bad sequences). Each analysis was stored in a separate file. All the files were then merged into a single FASTA file, and BlastAll [http://www.ncbi.nlm.nih.gov/; National Center for Biotechnology Information (NCBI)] was used to blast all the sequences against the RefSeq database (NCBI). For sequences unable to be verified in this way, a BLAT search (http://genome.ucsc.edu/cgi-bin/hgBlat; University of California-Santa Cruz Genome Bioinformatics) was used to position the sequences onto chromosomes and, where possible, to link these sequences to a neighbouring gene by overlapping expressed sequence tags. NCBI UniGene nomenclature was used for naming the genes. The MIAME compliant data is available at the NCBI Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/) [platform: GPL555; series: GSE779 (all experiments), GSE767 (0.5 h), GSE771 (3 h), GSE777 (6 h), GSE778 (24 h); sample: GSM1195911962, 1203412035, 1219112194, and 1242812435].
VN blot analysis.
VN blots were performed using the SMART PCR cDNA synthesis procedure (Clontech). In brief, full-length double-stranded cDNAs were generated from 1 µg of total RNA by the SMART procedure according to the manufacturer's protocol, ensuring optimization of the number of PCR cycles (typically 1821) to remain within the linear range of amplification. Approximately 1 µg of cDNA from each time point was run on a 1% agarose gel. The cDNA was transferred by capillary transfer to nylon membranes (Hybond-N, Amersham Biosciences) and cross-linked (Stratalinker, Stratagene) to the membrane. For probe generation, cloned cDNA fragments were PCR amplified using T7 and T3 primers, and the products were labeled with [-32P]dATP using either a MegaPrime Kit (Amersham) or StripEZ Kit (Ambion). Hybridization was performed in ExpressHyb solution (Clontech) for 2 h, and blots were washed in 2x SSC and 0.1% SDS at 25°C (twice) and 0.1x SSC and 0.1% SDS at 60°C for 20 min. Blots were visualized and quantified using a Typhoon 9410 PhosporImager and Imagequant 3.3 software (Molecular Dynamics and Amersham Biosciences). All blots were exposed for the appropriate time to ensure that the signal was within the linear range of the machine for quantification.
Q-RT-PCR.
Total RNA (3 µg) was converted to first-strand cDNA using an oligo(dT) primer (T18VN) and SuperscriptII reverse transcriptase (Invitrogen) according to the manufacturer's protocol. Residual contaminating DNA was removed using a DNA-free Kit (Ambion). Q-RT-PCR was performed on cDNA equivalent to 40 ng of starting total RNA with gene-specific primers and 0.1x SYBR green I nucleic acid gel stain (Molecular Probes) using AmpliTaq Gold DNA polymerase (Applied Biosystems) according to the Applied Biosystems' protocol. The amplification parameters were optimized for each primer pair so that appropriate-size single gene products were generated. After each real-time PCR amplification, melt curves were examined to confirm the generation of single product amplicons. Amplification parameters were typically 95°C for 10 min and cycling at 95°C for 20 s, 6065°C for 20 s, and 72°C for 30 s for a maximum of 40 cycles using a RoboCycler thermal cycler (Corbett Research; Mortlake, New South Wales, Australia). Data were analyzed using Rotor-Gene Analysis Software version 4.6. Primer pairs were designed within the 3'-untranslated region using PrimerExpress 1.5 software (Applied Biosystems) or Primer3 freeware and synthesized by GeneWorks (Adelaide, Australia). The sequences are given in Table 1.
|
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
|
VN blot analysis as a powerful and informative alternative to Q-RT-PCR for confirmation of microarray data. As with all microarray data, it is important and routine to verify the regulation using alternative more quantifiable methods such as Northern blotting, RNase protection, or Q-RT-PCR. This is particularly important in light of the data above demonstrating a wide range of magnitudes of induction even for the same gene. In the in vitro HUVEC 3-D collagen model of angiogenesis, cell numbers are limited, and this model is difficult to scale up while maintaining the morphological synchrony of the cells. This left Northern blot analysis out of the question, and Q-RT-PCR and VN analysis as the only options. Here, we investigated the feasibility of using VN blots as a method for confirming microarray results, particularly because no prior knowledge of the gene sequence is necessary. VN blots have been used previously to confirm differential expression (15, 32, 43, 44, 58), although, in the majority of cases, no numerical quantification of the levels of regulation were determined.
We chose five genes to analyze in detail as they were expressed to varying degrees in our libraries and had different regulation profiles. These were MMP10, hairy/enhancer-of- split related with YRPW motif 1 (HEY1), early growth response 1 (EGR1), PTGS2, and CTGF, all of which are well known to be involved in angiogenesis (5, 8, 16, 23, 45, 52). A comparison between microarray profiles (Fig. 3A), Q-RT-PCR (Fig. 3B), and VN blot analysis (Fig. 3, C and D) is shown. Microarray profiles of clones representing each gene have been overlaid. The maximal fold of regulation for each gene using all three methods of analysis is compiled in Fig. 3E. An excellent correlation was seen between the VN profiles and those derived by microarray and Q-RT-PCR in terms of profile shape. Over 100 VN blots were performed on genes identified as differentially expressed on microarrays, and their profiles closely resembled those obtained by microarray analysis (data not shown).
|
VN blots show a number of advantages over other quantification procedures. First, they enable visualization of different mRNA transcripts. One example is testis enhanced gene transcript (TEGT; Fig. 3G), which shows two differentially expressed transcripts due to utilization of alternative polyadenylation sites (54). Second, VN blots provide superior sensitivity and resolution compared with Northern blots (data not shown). Finally, they can be reused multiple times (up to 8 in our hands).
The comparisons between methods also underscore interesting results. The MMP10 cDNA fragment that gave a 1.2-fold maximal induction by microarray (Figs. 2B and 3A) showed 80- to 120-fold maximal induction when measured by VN analysis and Q-RT-PCR, respectively (Fig. 3E). Thus lowly regulated genes, as determined by microarray, can in fact be highly regulated when measured by these alternative methods.
Gene expression profiles during capillary tube formation.
The sequenced regulated genes were placed, where possible, into ontological groups (Fig. 4). A number of interesting and previously unreported observations was found.
|
The generation of the lumen requires extensive vesicle/vacuole formation, fusion, and directional targeting to the apical surface (13, 19, 36). In keeping with this, several RabGTPases and an ArfGTPase (ARF6), known for roles in vesicle trafficking and membrane fusion in the cell, were upregulated. Expression of these genes peaked at times where vacuole expansion and fusion is leading to lumen formation. RAB5A promotes circular ruffle formation that is required for 3-D migration rather than lamellipodia, which is associated with cell motility/spreading on 2-D matrixes (50). RAB6C, which is phosphorylated by PKC (17) and plays a role in vesicle trafficking and secretion, was induced from 6 to 24 h. RAB11A and RAB21 are involved in endosomal and vesicular trafficking to the apical surfaces of epithelial cells, respectively (38), the membrane surface to which fused vacuoles travel in the forming lumen (41). Also upregulated were TBC1D4 (a RabGAP) and DDEF1 (an ArfGAP for ARF1, -5, and -6) (6), GAPs linked to vesicle trafficking, the focal adhesion protein paxillin, and fibronectin-dependent migration, respectively (39). In all, 19 of 26 genes from the cytoskeleton/vesicular trafficking ontological group peaked in expression at 36 h, consistent with the majority of cell migration and vesicular expansion and trafficking occurring during this time.
PTGS1 and PTGS2 (cyclooxygenase 1 and 2) are involved in lipid metabolism, both catalyzing similar reactions in the conversion of arachidonic acid to prostaglandins. PTGS2 is well known to be differentially expressed during angiogenesis, mitogenesis, and inflammation, whereas PTGS1 is generally considered to be a constitutively expressing gene (52, 55). Whereas PTGS2 expression peaked at 36 h and then returned to basal levels by 24 h, PTGS1 expression remained constant for 6 h before dropping (>4-fold) by 24 h, indicating that overall cyclooxygenase activity rises dramatically early in capillary tube formation and then drops markedly upon maturation of the endothelial cells. In one study (52) of colon cancer, PTSG2 modulated the production of growth factors in the cancer cells, whereas PTGS1 regulated angiogenesis in endothelial cells. The finding that PTGS1 and PTGS2 were differentially expressed in endothelial cells is an intriguing one, suggesting that both may play a role in capillary tube formation and are coordinately regulated.
Other ontological groups of genes that were regulated include the nuclear transport genes RAN, RANBP2/4/9, KPNA4, and KPNB1. A highly regulated group of genes, with respect to magnitude, were the proteases/protease inhibitors, indicating the importance of extracellular matrix modulation in the generation of capillary tubes. Other groups such as the transcription factors, growth factors/hormones and receptors, and transporters/channels also had numerous highly regulated genes. Not surprisingly, the majority of these highly regulated genes have been reported previously in other in vitro and in vivo angiogenic model systems. However, we also identified a host of more lowly and subtly differentially expressed genes, which may shed new light on important processes or pathways in angiogenesis.
Transcriptional control of MAPK signal transduction pathway genes.
Classifying the regulated genes according to their functional role highlighted that the genes of several signal transduction pathways contained members that were significantly regulated. Figure 5A shows the microarray profiles of four sequential members of the MAPK signal transduction pathway. In this model system of angiogenesis, the levels of MAP3K1 (MEKK1) peak at 3 h and then return toward uninduced levels. Q-RT-PCR demonstrates that this level of induction is 13-fold (Fig. 5B). MAP2K1 (MEK1) and MAPK1 (ERK2) levels peak at 36 h, both rising only approximately twofold, but remaining elevated even out to 24 h. In marked contrast, the MAPK1/2 phosphatase (MKP1) dual specificity phosphatase 1 (DUSP1), which is a negative regulator of the MAPK pathway, peaked very early and rapidly declined to well below basal levels by 3 h. (Note that all other genes identified that peaked at 0.5 h are transcription factors; Fig. 4.) Within this time period, there was an 18-fold repression of the DUSP1 gene expression from its peak to lowest level. This highly dynamic regulation of the MAPK pathway members appears to be specific for the tube forming milieu because HUVECs, plated on a flat 2-D matrix of collagen and stimulated in the same way with PMA and anti-2ß1-integrin antibody but where cells do not form tubes, did not show significant regulation of MAPK pathway members (Fig. 5C).
|
A number of downstream target genes of the MAPK pathway were also identified, including the immediate-early genes Fos, JUNB, and EGR1 as well as phospholipase A2 (PLA2), signal transducer and activator of transcription 3 (STAT3), and the c-Myc binding partner MAX. Here, we demonstrated that the level of gene expression at the mRNA level of consecutive members of the MAPK signaling pathway is coordinately regulated in addition to its well-documented phosphorylation status.
Transcriptional control of the LIF signaling pathway.
LIF is a cytokine involved in the maintenance of embryonic stem cell pluripotency and inhibits angiogenesis and endothelial cell proliferation and migration (30). Microarray analysis has revealed the upregulation of members of this pathway during capillary tube formation (Fig. 5F). LIF, together with the LIF, oncostatin M, and IL-6 common receptor chain (IL6ST or gp130), peaked at 36 h. The levels of the downstream effectors JAK1 and STAT3 transcription show a sustained increase from 3 h. Targets of STAT3 that have been identified in this screen include Fos and IL-8, thus showing differential gene expression in the whole LIF signal transduction cascade from the extracellular ligand to the nucleus including targets in the nucleus.
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
We were able to identify on a global scale genes not only highly regulated but also many whose transcript levels vary by only 2- to 4-fold, representing a change of only 1.2- to 2-fold in microarray terms. In the vast majority of microarray studies reported to date on angiogenesis, a standard cutoff of greater than twofold has been used (1, 2, 26, 56). From our studies, this precludes the identification of a host of potentially important genes. For instance, at this cutoff, we would have listed only 7 rather than the 28 genes in the cytoskeleton/vesicular trafficking ontological group. Our recent report (49) on the role of ARHGAP24 in endothelial cell function in vitro and angiogenesis in vivo highlights the significance of identifying lowly regulated but biologically important genes. A meager twofold difference in gene expression is well recognized for numerous genes including VEGF and apolipoprotein B, where heterozygous knockout mice expressing half the wild-type levels of gene transcript are either nonviable or significantly affected (9, 24).
Categorizing regulated genes into ontological groupings is a useful way of visualizing the types of genes involved in a particular biological situation. The MAPK pathway is known to be crucial in the angiogenic process, and inhibitors of it prevent capillary tube formation and proliferation of endothelial cells (3). Activation of each of the kinases in this pathway is via rapid and cascading phosphorylation, whereas DUSP1/MKP-1 transcription and protein stability is upregulated by MAPK1 (ERK2) (and MAPK2/ERK1) (4, 25). DUSP1 then acts back on MAPK1/2 to dephosphorylate and inactivate the pathway as well as playing a possible role in the anchoring of MAPK1/2 within the nucleus (40). In most systems and experimental models, the levels of total proteins for these kinases remain quite consistent even over periods of up to 24 h. Here, we showed, for the first time, that the level of gene expression of multiple consecutive members of this pathway is regulated in a 3-D setting. The timing of this regulation takes hours, in contrast to the rapid signal transduction of this pathway by phosphorylation, which occurs in minutes. Interestingly, the three kinases of the MAPK signal transduction cascade peak at a similar time, with MAP2K1 and MAPK1 remaining elevated throughout the whole 24-h differentiation process. At least for MAP3K1 (MEKK1), this peak in expression occurs at a time when endotheial cells are actively migrating into the collagen gel. In fibroblasts, MAP3K1 orchestrates rear-end detachment of migrating cells through activation of calpain (12), and, in a human breast cancer cell line (T47D), protein levels increase threefold during mitosis and microtubule disruption, consistent with its role as a regulator of cell responses with respect to changes in microtubule integrity (57).
In contrast to the other members of the MAPK pathway, the expression profile of DUSP1, which inactivates MAPK1/2, is markedly downregulated after an initial spike of transcription. This may indicate the importance of ensuring that the MAPK pathway remains activated or activatable throughout the angiogenic process. A previous study (42) in rat glomerular mesangial cells demonstrated that transient or weak activation of MAPK1 leads to proliferation, whereas chronic activation, as might be expected here, results in differentiation. In other studies (28, 29), DUSP1 is differentially regulated at the transcript and protein levels in HUVECs, peaking at 1 h after adhesion to fibronectin or treatment with the cardiovascular hormone atrial natriuretic peptide. This results in a shortened activation of MAPK1/2 and p38 MAPK, respectively.
Long-term changes in the expression of ERK1 have been observed during neuronal regeneration (31), whereas upregulation of the murine Erk2 promoter occurs during differentiation of pC19 endothelial cells with retinoic acid (51). Upregulation of both the murine Erk1 and Erk2 genes was also reported during osteoblast differentiation from mesenchymal progenitors with BMP2 (34). Therefore, our results for upregulation of ERK2 transcript and phosphorylation levels during endothelial differentiation, as occurs in our 3-D model but not in a 2-D system, extend the observations of long-term ERK1/2 changes occurring during differentiation processes. Consistent with this are increases in markers of differentiation or maturation in this angiogenic model such as FLT1, TFPI2, ITGB1, ANGPT2, CXCR4, STC1, and HEY1. Each of these genes have previously been identified as important in angiogenesis, and their differential regulation at the mRNA level has been shown by others (2, 27).
In a similar fashion to the MAPK pathway, we found members of the LIF signaling pathway that were coordinately regulated. In this instance, differential expression is seen in members of the entire pathway from ligand through receptor and into the nucleus. This pathway is implicated in the inhibition of angiogenesis and endothelial cell proliferation and migration (30). Individual members of each of these pathways have been reported previously to be differentially expressed in other systems including the oncostatin M, LIF, and IL-6 common receptor (IL6ST). Consistent with our results, JAK1 was also regulated in a similar 3-D collagen HUVEC model (2). Interestingly, both LIF and oncostatin M, but not IL-6, promote angiogenesis in vivo (53).
A significant technical finding from this study was the large variation in fold induction seen for different spots representing not only the same gene but also the same cDNA fragment. This variation may be due to target fragment length, GC content, concentration, distance from the polyA+ tail, or ability to be PCR amplified before being spotted onto glass slides. It may also be due to abundance of transcript, alternatively processed variants, probe amplification, or signal-to-noise ratio. Preliminary studies showed a trend where target fragments most distant to the polyA+ tail resulted in low magnitudes, although this was not true for all genes. We are currently pursuing which of these factors plays a role in such variation, as this may facilitate future design of target sequences on microarrays. What is clear is that, until the rules for generating optimal cDNA target sequences are better understood, it remains crucial in microarray design to provide multiple targets and target concentrations for each gene transcript. As demonstrated here, a subtractive hybridization approach is one way to achieve this while providing the advantage of identifying as-yet-unrecognized novel cDNA sequences or splice variants.
Q-RT-PCR has become the most commonly used method to confirm differential expression of genes identified using microarray. In our system, sample size is small with total RNA available being <10 µg. Ten micrograms of total RNA is sufficient for 1 Northern blot, 510 RNase protection assays, 500 Q-RT-PCRs, or 2,500 VN blots. In addition, VN blots are 1050 times more sensitive than Northern blots and, because they contain cDNA rather than labile RNA, can be probed multiple times without a significant loss of signal. We find that both VN analysis and Q-RT-PCR display similarly sensitivity (data not shown) and give similar magnitudes of differential expression. Whereas Q-RT-PCR requires knowledge of the gene sequence, the design and synthesis of primer pairs, and may require optimization of conditions, VN analysis needs no knowledge of gene sequence, utilizes simple probe preparation by excision of inserts or their PCR amplification using common vector-derived primers, and requires no expensive equipment. Hence, VN blot analysis provides a reproducible, sensitive, rapid, scalable, and cost-effective approach to quantify gene transcript levels.
The combination of precision and utilization of multiple statistical analysis in this study has uncovered biologically relevant but lowly regulated genes and entire pathways that are coordinately regulated at the gene transcript level. Thus, in a situation where a whole system is being altered, as takes place during angiogenesis, analysis of entire pathways may be crucial to our ability to both understand and ultimately manipulate such processes.
![]() |
GRANTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
ACKNOWLEDGMENTS |
---|
Present address of C. J. Drogemuller: Applied Genetic Technology Corporation, 12085 Research Dr., Alachua, FL 32615.
Present address of S. R. Waterman: Dept of Microbiology, Univ. of Texas Southwestern Medical Center, 5323 Harry Hines Blvd., Dallas, TX 73590-9048.
![]() |
FOOTNOTES |
---|
Address for reprint requests and other correspondence: C. N. Hahn, Vascular Biology Laboratory, Hanson Institute, Frome Rd., Adelaide SA 5000, Australia (E-mail: chris.hahn{at}imvs.sa.gov.au).
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Visit Other APS Journals Online |