A model of anti-angiogenesis: differential transcriptosome profiling of microvascular endothelial cells from diffuse systemic sclerosis patients
- Betti Giusti1,
- Gabriella Fibbi2,
- Francesca Margheri2,
- Simona Serratì2,
- Luciana Rossi1,
- Filippo Poggi1,
- Ilaria Lapini1,
- Alberto Magi1,
- Angela Del Rosso3,
- Marina Cinelli3,
- Serena Guiducci3,
- Bashar Kahaleh4,
- Laura Bazzichi5,
- Stefano Bombardieri5,
- Marco Matucci-Cerinic3,
- Gian Franco Gensini1, 6,
- Mario Del Rosso2Email author and
- Rosanna Abbate1Email author
© Giusti et al.; licensee BioMed Central Ltd. 2006
Received: 17 October 2005
Accepted: 30 June 2006
Published: 19 July 2006
The objective of this work was to identify genes involved in impaired angiogenesis by comparing the transcriptosomes of microvascular endothelial cells from normal subjects and patients affected by systemic sclerosis (SSc), as a unique human model disease characterized by insufficient angiogenesis. Total RNAs, prepared from skin endothelial cells of clinically healthy subjects and SSc patients affected by the diffuse form of the disease, were pooled, labeled with fluorochromes, and hybridized to 14,000 70 mer oligonucleotide microarrays. Genes were analyzed based on gene expression levels and categorized into different functional groups based on the description of the Gene Ontology (GO) consortium to identify statistically significant terms. Quantitative PCR was used to validate the array results. After data processing and application of the filtering criteria, the analyzable features numbered 6,724. About 3% of analyzable transcripts (199) were differentially expressed, 141 more abundantly and 58 less abundantly in SSc endothelial cells. Surprisingly, SSc endothelial cells over-express pro-angiogenic transcripts, but also show up-regulation of genes exerting a powerful negative control, and down-regulation of genes critical to cell migration and extracellular matrix-cytoskeleton coupling, all alterations that provide an impediment to correct angiogenesis. We also identified transcripts controlling haemostasis, inflammation, stimulus transduction, transcription, protein synthesis, and genome organization. An up-regulation of transcripts related to protein degradation and ubiquitination was observed in SSc endothelial cells. We have validated data on the main anti-angiogenesis-related genes by RT-PCR, western blotting, in vitro angiogenesis and immunohistochemistry. These observations indicate that microvascular endothelial cells of patients with SSc show abnormalities in a variety of genes that are able to account for defective angiogenesis.
Systemic sclerosis (SSc) affects the connective tissue of the skin and internal organs, such as gastrointestinal tract, lungs, heart and kidneys. Disease progression involves the immune system, extracellular matrix (ECM) deposition and the microvasculature . In the later stages of the disease, the vessel walls are thickened and hyalinized and their lumen is narrowed, leading to devascularization and tissue ischemia, which is not counterbalanced by active neo-angiogenesis. Angiogenesis, the process of new blood vessel generation from capillary or post-capillary venules, requires gross changes in endothelial cell function. In this process, an endothelial cell modifies the interaction with its basement membrane, remodels and migrates through ECM, proliferates, and differentiates. The final effect is the formation of endothelial tubules with a lumen, which are capable of transporting blood . Newly expressed molecules or hyper-expression of pre-existing ones are coordinately required in this series of events, including proteolytic enzymes that are believed to be critical to ECM remodeling , growth factor activation  and release of ECM-trapped regulatory molecules [5, 6].
While gene-expression profiling using microarray technologies is available for skin biopsies  and cultured fibroblasts from individuals with a diagnosis of SSc [8, 9], a global portrait of gene expression of microvascular endothelial cells (MVECs) has not been reported in the literature. In order to better understand whether dysregulated genes may contribute to the pathogenesis of defective angiogenesis, we have undertaken studies of gene expression in MVECs isolated from the lesional skin of patients affected by the diffuse form of SSc and matched healthy controls, using a 14,000 oligonucleotide (70 mer) microarray. After the identification of differentially expressed genes by a Bayesian empirical model [10, 11], genes were annotated on the basis of biological process ontology and statistically significant gene ontology terms were evaluated.
The results show that of the several thousands genes that passed filtering criteria, 199 genes are differentially expressed, 141 being up-regulated and 58 down-regulated in SSc endothelial cells. We observed that SSc endothelial cells overexpress pro-angiogenic and anti-angiogenic transcripts, and down-regulation of genes critical to cell migration and proliferation (including tissue kallikreins (KLKs)) , adhesion and capillary differentiation. We have validated the data on the main anti-angiogenesis-related molecules by RT-PCR and have focused functional experiments on differentially expressed molecules that have recently been shown to be relevant to endothelial cell physiology, such as plexin B1, pent(r)axin 3 and desmoglein (DSG) 2. Plexin B1, which we found to be down-regulated in MVECs of SSc patients, has been reported to bind and mediate the pro-angiogenic signal of semaphorin 4D . Pent(r)axin 3, which we found to be up-regulated in MVECs of SSc patients, inhibits the pro-angiogenic effect of Fibroblast Growth Factor-2 (FGF2), including that produced by autocriny of endothelial cells . Desmoglein 2 is a calcium-binding trans-membrane protein of the cadherin cell adhesion molecule superfamily that mediates homophilic cell adhesion, and has been identified as a structural component of endothelial cell intercellular junctions . Here we show that DSG2 down-regulation in MVECs of SSc subjects associates with an anti-angiogenic phenotype.
We also identified differential expression of transcripts controlling haemostasis, inflammation, stimulus transduction, transcription, protein synthesis and genome organization. Other up-regulated genes are markers of cellular stress, such as those of the ubiquitin-proteasoma family. Taken together, our results show that over-expression of some genes in SSc MVECs indicates a response to a powerful pro-angiogenic environment, while over-expression of others may render them unable to respond to angiogenic stimuli by over-expression of anti-angiogenic and hypo-expression of pro-angiogenic molecules.
Materials and methods
Patients, controls, tissue biopsies and endothelial cell preparation
Patients with diffuse SSc (submitted to skin biopsies for diagnostic purposes at the Departments of Medicine, Division of Rheumatology, Florence and Pisa Universities) and healthy controls were used as sources of MVECs. All patients (20 females and 10 males, with a mean disease duration of 9 years (range 3.1 to 12)) fulfilled the American College of Rheumatology criteria for the classification of SSc . Only patients classified as having the diffuse cutaneous SSc were admitted to the study (sclerosis of both distal and proximal extremities, with or without truncal involvement). Patients with overlap symptoms to other connective diseases were excluded from the study, as well as patients affected by other diseases involving the vascular system. Biopsies were performed on the dorsal involved skin of the hands. Fifteen healthy patients undergoing surgery for traumatic events involving the hands were subjected to the MVEC isolation procedure, after punch biopsies of the dorsal skin of the hands, which were processed as skin biopsies of SSc patients. The study was approved by the local Ethical Committee and patient consent was obtained from each subject enrolled. Ethics approval and patient consent were granted for this manuscript.
The patients were not on steroids, cyclophosphamide, D-penicillamine, relaxin or other disease-modifying drugs. Calcium channel blockers were stopped ten days before the biopsy. Only proton pump inhibitors and cisapride were allowed. Briefly, skin biopsies have been mechanically cleaned of epidermis and adipose tissue in order to obtain a pure specimen of vascularized dermis, and treated as described elsewhere [17, 18]. In some cases, clusters of round-shaped cells were squeezed from microvessels and formed colonies composed by polygonal elements. Such colonies were detached with EDTA, and CD31-positive cells were subjected to immuno-magnetic isolation with Dynabeads CD31 (Dynal ASA, Oslo, Norway) . Isolated cells were further identified as MVECs by labeling with anti-factor VIII-related antigen and by re-probing with anti-CD31 antibodies. Cells were maintained in complete MCDB medium, supplemented with 30% FCS, 20 μg/ml endothelial cell growth supplement (ECGS), 10 μg/ml hydrocortisone, 15 UI/ml heparin, and antibiotics (100 UI/ml penicillin, 100 μg/ml streptomycin, 50 μg/ml amphotericin). MVECs from normal individuals and from SSc patients are referred to as N-MVECs and SSc-MVECs, respectively, and were used between the third and seventh passage in culture. A biopsy specimen from each subject was formalin fixed and paraffin embedded for immunohistochemistry assays. Each case has been stained with hematoxilin and eosin to assess the original diagnosis.
Since the success rate for isolation of SSc-MVECs is lower than 20%, compared to a success rate of more than 50% for N-MVEC, controls were matched by age and sex to the SSc cases that yielded MVECs. Therefore, total RNA was prepared from six N-MVEC and six SSc-MVEC pellets using the RNeasy Minikit (Qiagen, Hilden, Germany) according to the manufacturer's protocol. Equal quantities of total RNA from each of the six N-MVEC and six SSc-MVEC pellets were pooled to give a N-MVEC pool and a SSc-MVEC pool.
Microarray based gene expression analysis
The setting and the subsequent hybridization of microarrays were performed as described in a previous paper . Briefly, we used poly-L-lysine (Sigma Chemical Company, St. Louis, MO, USA) coated arrays representing 14,000 genes (70 mer oligonucleotides; Human Array-Ready Oligo set version 1.1, Operon Technologies, Inc., Alameda, CA, USA). We reverse transcribed and labeled 20 μg of the N-MVEC pool and 20 μg of the SSc-MVEC pool with NHS-cyanine dyes (Cy3 and Cy5; Amersham Biosciences, Amersham Place, England). The two labeled probes were hybridized with the array for 16 h at 65°C. Arrays were scanned by using a 4000B Scanner (Axon Instruments, Union City, CA, USA). Due to difficulty in growing SSc-MVECs, we performed two replicates of the microarray experiment: one with N-MVECs Cy3 labeled and SSc-MVECs Cy5 labeled and one with N-MVECs Cy5 labeled and SSc-MVECs Cy3 labeled (dye swap).
Image processing and statistical analysis
Each hybridization produced a pair of 16-bit images, which were processed using the GenePix Pro 4.1 software (Axon Instruments). Poorly spotted genes, expressing weak or distorted signals, were automatically discarded by the GenePix Pro 4.1 software and manually by visual inspection. In order to reduce the identification of false positive differentially expressed genes, spots exhibiting at the same time low Cy3 and Cy5 fluorescent signal intensities (<100) were discarded from consideration by pre-processing the raw data in Axon .gpr format using Microsoft Excel. For each microarray, we performed a local intensity-dependent normalization using an lowess scatter plot smoother to remove dye and spatial (print-tip) effects . After single-slide normalization we applied a dye-swap normalization as proposed by Kerr and colleagues [11, 20] to correct the different properties of the dyes on a gene by gene basis . The data obtained from normalization was analyzed by Newton algorithm . More details on image processing and statistical analysis have been previously reported . The full data set is available at ArrayExpress .
Gene ontology data analysis
In our study we used one of the three ontologies produced by the Gene Ontology (GO) consortium, the biological process ontology. The term 'biological process' should be interpreted as a biological function to which the gene product contributes. The actual mappings of genes to GO terms are provided by the Gene Ontology Annotation Database [23, 24]. The mappings were downloaded from .
In brief, given a set of genes and one ontology, we first found the set of all unique GO terms within the ontology that are associated with one or more of the genes of interest. Next, we determined how many of the selected 199 differentially expressed genes are annotated at each term and how many of the genes that were assayed (all the genes represented on the microarray) are annotated at the term. The test evaluated if there are more interesting genes at the term than one might expect by chance. Due to the small number of genes in some categories, significant terms were inferred by two-sided Fisher's exact test . The statistical analyses were implemented in the R environment using Bioconductor packages .
Criteria based only on GO terms were not sufficient to classify a gene as positively or negatively involved in the regulation of angiogenesis. Therefore, we included the biological processes obtained by GO into the following families: angiogenesis, apoptosis, haemostasis, inflammation and immunity, stress and ubiquitination, transductions, DNA/RNA organization, transcription, protein synthesis, and mitochondrial functions. In particular, in order to be classified as pro-angiogenic, a gene must play a significant role in endothelial cell adhesion, invasion, proliferation, and differentiation.
Real Time RT-PCR based gene expression analysis
Time-fold up- or down-expression of genes analyzed by real time PCR in SSc-MVECs versus N-MVECs
Time-fold up- or down-expression in SSc-MVECs relative to N-MVECs
For immunohistochemistry, tissue sections were 3 to 5 μm thick and placed on pretreated glass slides, dewaxed and treated to block endogenous peroxidase activity. The following primary antibodies were employed: rabbit anti-human KLK9 (catalog n K005-12, raised against a synthetic peptide corresponding to amino acids 239 to 250 of the human KLK9 protein) and anti-human KLK12 (catalog n K005-15, raised against a synthetic peptide corresponding to amino acids 236–248 of the human KLK12 protein), and mouse anti-human KLK11 (catalog n K005-14, raised against human recombinant KLK11), all from US Biological (Swampscott, MA, USA); anti-DSG mouse monoclonal antibody (Chemicon, Temecula, CA, USA); anti-pentraxin 3 rat monoclonal antibody MNB4 (Alexis Biochemical, Lausen, Switzerland); anti-plexin B1 and anti-CTGF rabbit polyclonal antibodies, both from Santa Cruz Biotechnology (Santa Cruz, CA, USA); murine monoclonal antibody 5B4 (mAb5B4), which recognizes the kringle domain of the A chain of PLAU, a kind gift of Dr ML Nolli (Lepetit Research Center, Varese, Italy); and anti-IL8 rabbit polyclonal antibodies (Chemicon, Temecula, CA, USA). All the primary antibodies were diluted 1:40 and incubated overnight with tissue sections in a moist chamber at 4°C. A standard streptavidin-biotin detection system (Vector, Burlingame, CA, USA) was carried out. Isotype Ig controls were used in parallel with primary antibodies to assess the specificity of the staining. Primary antibody bound to antigen was visualized by diaminobenzidine staining and a nuclear counterstaining with hematoxilin was performed. Immunohistochemistry was performed on the skin biopsies of the six normal subjects and six SSc patients whose MVECs were used for RNA preparation. Immunohistochemistry quantification was performed by image analysis using the ScnImage program .
N-MVECs and SSc-MVECs were grown to 70% confluence and serum-starved overnight in MCDB supplemented with 2% FCS. Cells were then suspended in lysis buffer (10 mM Tris-HCl, pH 7.4, containing 150 mM NaCl, 1% Triton X-100, 15% glycerol, 1 mM sodium orthovanadate, 1 mM NaF, 1 mM EDTA, 1 mM phenylmethylsulfonyl fluoride, and 10 μg/ml aprotinin). We electrophoresed 60 μg of the cell extract proteins on 12% SDS-PAGE under reducing conditions and then blotted to a polyvinylidene difluoride membrane (Hybond-C Extra; Amersham Biosciences, Piscataway, NJ, USA) for 3 h at 35 V. The membrane was incubated with 5% skim milk in 20 mM Tris buffer, pH 7.4, for 1 h at room temperature to block non-specific binding and then probed with antibodies directed against pentraxin 3, or plexin B1 or DSG2 overnight at 4°C. After incubation with horseradish peroxidase-conjugated donkey anti-mouse IgG (1:5,000) for 1 h (Amersham Biosciences), immune complexes were detected with the enhanced chemiluminescence detection system (Amersham Biosciences). The membranes were exposed to autoradiographic films (Hyperfilm MP; Amersham Biosciences) for 1 to 30 minutes.
A Boyden chamber was used to evaluate spontaneous and stimulated invasion (chemoinvasion) through Matrigel-coated porous filters, as described . For spontaneous invasion, 50 μl of cell suspension (6.25 × 103 cells) were placed in the upper compartment of the chamber and migration was allowed to occur for 6 h. To inhibit the activity of the relevant molecules, specific antibodies (anti-pentraxin 3 and anti-plexin B1, each at 3 μg/ml final concentration) were added to both the upper and lower compartment of the migration chamber. Irrelevant mouse IgGs were used to verify the specificity of the effect. The number of cells moving across the filter measured mobilization. Experiments were performed in triplicate. Migration was expressed as mean ± standard deviation (SD) of the number of total cells counted per filter or as the percentage of basal response.
Preparation of SSc-MVEC conditioned medium
Confluent cultures of SSc-MVECs were washed twice with phosphate-buffered saline and incubated overnight in the presence of MCDB medium supplemented with 2% FCS. The culture supernatant was centrifuged at 1,500 rpm for 10 minutes and either used immediately or stored at -20°C.
In vitrocapillary morphogenesis assay
Matrigel (0.5 ml; 10 to 12 mg/ml) was pipetted into 13 mm diameter tissue culture wells and polymerized for 30 minutes to 1 h at 37°C, as described . N-MVECs were plated (60 × 103/ml) in complete MCDB medium supplemented with 30% FCS and 20 μg/ml endothelial cell growth supplement. Capillary morphogenesis was also performed in the presence of 3 μg/ml of anti-pentraxin 3 or anti-plexin B1 antibody. Irrelevant mouse IgGs were used as negative control. Plates were photographed at 24 h. Results were quantified by measuring the percentage of the photographic field occupied by endothelial cells by image analysis. Six to nine photographic fields from three plates were scanned for each point.
Results are expressed as means ± SD for (n) experiments. Multiple comparisons were performed by the Student-Newman-Keuls test, after demonstration of significant differences among medians by nonparametric variance analysis according to Kruskal-Wallis.
Microarray, gene ontology analysis, and class distribution of differentially expressed genes
Of the 14,000 transcripts represented on our arrays, after data processing and application of the filtering criteria, the analyzable features numbered 6,724. The full list of the 150 most expressed genes, independent of the cellular source (N-MVEC and SSc-MVEC), is available as Additional file 1a. We used a Newton algorithm after single slide and dye-swap normalization to assess the 6,724 genes expressed by MVECs for differential expression between SSc-MVECs and N-MVECs.
Genes found differentially expressed between SSc-MVECs and N-MVECs numbered 199 (3% of the total transcripts analyzed; Additional files 2 to 7). Of these, 141 transcripts were expressed more abundantly and 58 less abundantly in the SSc-MVECs.
List of all the Gene Ontology significant terms with more than two annotated genes on the array (N > 2)
GO category number
EIF4G1(↑), ITGB4BP(↑), LAMR1(↑), RPL7(↑), RPL9(↑), RPL10(↑), RPL12(↑), RPL23A(↑), RPS5(↑), RPS10(↑), RPS20(↑), RPL14(↑), NOLA2(↑), RPL10A(↑), RPL38(↑)
Induction of positive chemotaxis
Copper ion homeostasis
Proteolysis and peptidolysis
ANPEP(↑), AZU1(↓), CAPN2(↑), NEDD8(↑), PLAT(↑), PLAU, PRSS1(↓), ADAM15(↑), NPEPPS(↑), KLK11(↓), SUPT16H(↑), CASP14(↓), KLK12(↓)
Rho protein signal transduction
RHOA(↑), ARHGDIB(↑), CFL1(↑)
Actin cytoskeleton organization and biogenesis
RHOA(↑), ARHGDIB(↑), CFL1(↑), PFN1(↑)
ACTN4(↑), CTGF(↑), IL8(↑), SPOCK(↑), ACTR3(↑), ARPC2(↑)
Positive regulation of apoptosis
Regulation of translation
MKNK2(↑), LAMR1(↑), PPP2CA(↑)
Negative regulation of apoptosis
Regulation of cell adhesion
Response to oxidative stress
ATOX1(↑), DUSP1(↑), GPX2(↓)
GP9(↓), PLAT(↑), PLAU(↑), THBD(↑)
APOA1(↓), HMGA1(↑), LBP(↓)
Response to wounding
Regulation of transcription, DNA-dependent
BTF3L3(↑), RUNX2(↓), ENO1(↓), GATA6(↓), HMGA1(↓), FOXA2(↓), NFKB2(↓), SSRP1(↑), UBE2V1(↑), NFAT5(↓), SIX4(↓)
Positive regulation of transcription
Protein kinase cascade
MKNK2(↑), STK17A(↓), MAP4K1(↓)
Perception of sound
GJA1(↑), PMP22(↑), WDR1(↑), TIMM8B(↑)
The GO biological processes include many of those known to be required to fulfill an angiogenic program. Of particular interest are the genes involved in proteolysis and peptidolysis, cell migration and cell motility, Rho protein signal transduction, regulation of cell adhesion, blood coagulation, and mitosis.
However, many of the differentially expressed genes have multiple functions, each one often required for angiogenesis, and some recognized pro- or anti-angiogenic properties of several genes are not yet available in the GO biological processes.
Class distribution of differentially expressed genes with a log odds ratio >1
↑ Metallothionein 1A
↑ Plasminogen activator, tissue
↑ Angiopoietin-like 4
↑ Connective tissue growth factor
↑ Enolase 1
↑ Pentaxin-related gene, rapidly induced by IL-1 beta
↑ Laminin receptor 1
↑ Profilin 1
↑ Metallothionein 2A
↑ Sparc/osteonectin (testican)
↑ Metallothionein 1E
↑ Cyclin D1
↑ Cofilin 1
↑ Calpain 2, large subunit
↑ Macrophage migration inhibitory factor
↑ WD repeat domain 1
↑ Calsyntenin 1
↑ A disintegrin and metalloproteinase domain 15
↑ EGF-containing fibulin-like extracellular matrix protein 1
↓ Kallikrein 9
↓ Plexin B1
↓ Kallikrein 12
↓ Kallikrein 11
↓ Desmoglein 2
↑ Tumor necrosis factor receptor superfamily, member 6b, decoy
↑ Proteoglycan 1, secretory granule
↑ Cofilin 1
↑ Calpain 2, (m/II) large subunit
↑ Serum/glucocorticoid regulated kinase
↑ Macrophage migration inhibitory factor
↑ Phosphoprotein enriched in astrocytes 15
↑ Mitochondrial carrier homolog 1
Haemostasis, inflammation and immunity
↑ Plasminogen activator, tissue
↑ Enolase 1
↑ Pentaxin-related gene
↑ Interferon induced transmembrane protein 2
↑ Interleukin enhancer binding factor 2
↓ Lipopolysaccharide binding protein
↓ Prostaglandin D2 synthase 21 kDa
↓ Azurocidin 1
↑ ATX1 antioxidant protein 1 homolog
↑ Ubiquitin carboxy-terminal esterase L1
↑ Anaphase promoting complex subunit 11 homolog
↑ Ferritin, light polypeptide
↑ Peripheral myelin protein 22
↑ Solute carrier family 38, member 2
↑ Proteasome (prosome, macropain) 26S subunit, ATPase, 1
↑ FK506 binding protein 4
↑ Proteasome (prosome, macropain) 26S subunit, non-ATPase, 13
↑ Peroxiredoxin 1.
↓ Anaphase promoting complex subunit 10
↑ Gap junction protein, alpha 1, 43 kDa (connexin 43)
↑ Hippocalcin-like 1
↑ Cofilin 1
↑ Ras-related C3 botulinum toxin substrate 2
↑ Serum/glucocorticoid regulated kinase
↑ Interleukin 6 signal transducer (gp130, oncostatin M receptor)
↑ Rho GDP dissociation inhibitor beta
↑ Guanine nucleotide binding protein, gamma 10
↑ Peroxiredoxin 4
↑ Ras homolog gene family, member A
↓ Plexin B1
↓ Mitogen-activated protein kinase/kinase/kinase/kinase 1
↑ Shwachman-Bodian-Diamond syndrome
↑ Enolase 1
↑ H3 histone, family 3B
↓ TruB pseudouridine synthase homolog 2
↓ Arginine vasopressin-induced 1
↓ Chromosome 1 open reading frame 34
↑ Basic transcription factor 3, like 3
↑ Interleukin enhancer binding factor 2
↓ GATA binding protein 6
↓ SRY (sex determining region Y)-box 9
↓ Runt-related transcription factor 2
↓ Nuclear factor of activated T-cells 5
↓ SRY (sex determining region Y)-box 5
↓ Forkhead box A2
Table 3, which reports differentially expressed genes with a log odds ratio (LOR) >1.0 (see also Additional file 2), where gene function was classified according to GO and to information available in the NCBI web site and related links [29, 30], indicates that many genes that mediate endothelial cell migration/invasion, proliferation, cytoskeletal remodeling and capillary differentiation (angiogenesis section) are up-regulated in SSc-MVECs. However, the angiogenesis inhibitor pent(r)axin-related gene (PTX3) is also up-regulated, while other genes critical for the angiogenic process, such as plexin B1 (PLXNB1, semaphorin receptor), tissue kallikreins KLK9, KLK11, and KLK12) , and DSG2 (a cadherin that mediates homophilic cell adhesion), undergo down-regulation in SSc-MVECs. Apoptosis-related genes (Table 3 and Additional file 3) were variously altered, including down-regulation of BCL2 in SSc-MVECs, which also exhibited a general pro-fibrinolytic pattern, coupled with over-expression of PTX3, which increases tissue factor expression and stimulates generation of inflammation mediators, . SSc-MVECs also show up-regulation of genes related to a response to oxidative, osmotic, and shear stress, and of genes linked with protein ubiquitination and proteasoma activation (Table 3 and Additional file 4). Table 3 and Additional file 5 indicate an overall perturbation of signal transductions mediated by small GTPase proteins, and down-regulation of MAP4K in SSc-MVECs, and of genes involved in nucleosome and chromatin remodeling and in regulation of transcription, including down-regulation of GATA6, which controls transcription of von Willebrand factor, and of RUNX2, which controls endothelial cell migration and invasion. Table 3 and additional file 6 show up-regulation in SSc-MVECs of a large number of structural constituents of ribosomes and of genes engaged in oxidative phosphorylation and related ATP production, indicating an intense protein synthesis and energy production in SSc-MVECs. A series of 36 differentially expressed genes whose functions are unknown or cannot be included within a class is available as Additional file 7.
Validated expression of selected genes by RT-PCR
To validate the results of the cDNA microarray analysis, the mRNA expression of nine selected genes was independently examined with real time RT-PCR. We selected nine differentially expressed transcripts, including many of those that are functional to the main hypothesis of the present work (KLK9, AF135026; KLK11, AB012917; KLK12, AF135025; IL8, M17017; PLAU, X02419; PTX3, M31166; PLXNB1, AJ011414; DSG2, NM_001943; CTGF, NM_001901); among these were genes exhibiting a significant decrease (KLK9, KLK11, KLK12, PLXNB1, and DSG2) or a significant increase (IL8, PLAU, PTX3, and CTGF) in expression in SSc-MVECs in comparison to N-MVECs. We evaluated these in the same total RNAs used for comparative microarray experiments. Real time RT-PCR analysis confirmed the data obtained by microarray technology (Table 1). To reinforce the data on the genes reported in Table 1, we also performed RT-PCR determinations on single RNA preparations (from six SSc patients and six healthy subjects), as previously described for single KLKs . The values obtained from these determinations were similar to those obtained from the RNA pools: a mean fold increase for PTX3 (1.72, range 1.18 to 4.73, p = 0.041), IL8 (3.3, range 1.51 to 7.9, p = 0.039), PLAU (2.76, range 1.51 to 5.2, p = 0.02), and CTGF (1.79, range 1.19 to 3.09, p = 0.026); and a mean fold decrease for PLXNB1 (1.96, range 1.35 to 5.91, p = 0.042), DSG2 (29.91, range 10.62 to 68.9, p = 0.02), KLK9 (20.69, range 3.82 to 75.0, p = 0.022), KLK11 (34.48, range 5.83 to 150.0, p = 0.021), and KLK12 (24.26, range 2.64 to 118.0, p = 0.020).
Functional studies on the angiogenic effects of pentraxin 3, plexin B1 and DSG2
To date, this is the first study that compares the differential transcriptosome of MVECs isolated from the skin of normal subjects and from the lesional skin of SSc patients affected by the diffuse form of the disease in the avascular phase.
The number of 6,724 genes that passed the filtering criteria is in agreement with those previously obtained in other human cell lines by microarray analysis  and by massively parallel signature sequencing . We have observed that the majority of genes are expressed at similar levels in N-MVECs and SSc-MVECs and that, interestingly, 3.2% of the total filtered transcripts (199 genes) were differentially expressed. Considering the correctness of comparing data limited to transcripts, we show that, in SSc-MVECs, the dysregulation involves only a small number of genes controlling a large number of processes that are critical to the biology of endothelial cells.
Surprisingly, we have found that SSc-MVECs exhibit a pro-angiogenic gene expression pattern (Table 3). In SSc, the angiogenic process is severely impaired in the late phases of the disease , independent of the increased levels of circulating Vascular Endothelial Growth Factor (VEGF) and FGF2 [34, 35]. Therefore, it is likely that some critical checkpoints in the control of angiogenesis are altered at the MVEC level. Our data indicate down-regulation of PLXNB1, a receptor for semaphorin that tracks the pathway to migrating endothelial cells by activation of the MET oncogene and by stimulating Rho-initiated pathways [13, 36, 37], of tissue KLK9, KLK11 and KLK12, shown to be required for MVEC migration and proliferation , and of DSG2, which is positively involved in homophilic cell-cell interaction [15, 38]. At the same time, PTX3, an angiogenesis inhibitor that acts by binding FGF2 , is up-regulated (Table 3). We propose to interpret these data as a stabilization of a pro-angiogenic pattern dictated by angiogenesis factors such as VEGF and FGF2 that is blocked or rendered ineffective by the strong down-regulation of the critical adhesion/invasion/proliferation systems and by up-regulation of the angiogenesis inhibitor PTX3 (Table 3). Due to a common technical problem during the spotting procedure, microarray data were not available for matrix metallo-protease-12 (MMP12), which in a previous study on the same SSc-MVECs was found to be up-regulated and responsible for urokinase-type plasminogen activator receptor (uPAR) truncation and subsequent angiogenesis impairment . In this previous work we suggested the possibility that the alterations we observed in SSc-MVECs may initially be stimulated by environment factors and then become the product of the hypoxia-induced selection of an endothelial cell population more suitable to survive in the hypoxic micro-environment typical of the disease. Still, unexpectedly, in SSc-MVECs we have observed an overall up-regulation of many components of several transduction systems. The large majority of such alterations (Table 3, Additional file 5) deal with transduction by small GTPase proteins, which couple signals from ECM molecules to the cell cytoskeleton, inducing alternating states of cell contraction/relaxation . However, up-regulation of Rho GDP dissociation inhibitor beta (ARHGDIB) is likely to provide a critical 'bottle-neck' to GTPase protein transductions by impairing the substitution of GDP with GTP, thus preventing the receptor to enter its activation state . Also, PLXNB1 down-regulation inhibits semaphorin-directed MVEC migration by blunting Rho-initiated pathways [36, 37]. Further, an overall impairment to cell proliferation in SSc-MVECs must also be related to down-regulation of MAP4K1, a serine/threonine kinase involved in a variety of cellular signaling cascades .
Over-expression of genes of proteasome and ubiquitin pathways suggests the possibility that the observed gene alterations are the effect of a cell adaptation to a particularly hostile environment. It is interesting to underline the induction of 3 metallothionein genes in SSc-MVECs. Metallothioneins belong to a family of stress-induced proteins that regulate Zn and Cu availability and also modulate the amount and activity of NF-KB, a transcription factor for genes involved in apoptosis, immune response and inflammation .
Since RNA harvested from MVECs was obtained from cells between the fourth and seventh passage, this could raise concerns about MVEC stability and/or possible selection of cells more likely to survive in culture. However, we previously showed that RT-PCR, performed with RNA isolates of cells from single patients at early and late culture passages to validate tissue KLK expression, gave similar results, in agreement with microarray data obtained with pooled RNA . Moreover, the demonstration of down-regulation of KLK9, KLK11, KLK12 and DSG2 in skin biopsies of SSc patients by immunocytochemistry suggests that at least some of the reported alterations pre-exist to the isolation and culture propagation techniques. Nonetheless, a possible selection bias, responsible for some of the reported differences, cannot be ruled out. There is evidence that inference from most genes is not adversely affected by pooling, such that pooling is recommended when fewer than three arrays are used in each condition . Additionally, due to the small number of replicates, we decided to apply the dye swap design to minimize the gene-specific dye bias , which is the major source of experimental variability between replicates. In the present study we have applied markedly stringent criteria for feature extraction and data normalization, which could blunt the identification of differentially expressed genes. Nevertheless, we believe that our findings, related to single genes, or classes of genes, may provide hints for future research and are worthy of a more in-depth study to identify possible ways to correct some critical molecular defects and to recover, at least partially, the angiogenic attitude of SSc-MVECs.
In this study we analyzed gene expression profiles of skin MVECs isolated from healthy subjects and diffuse SSc patients, using a microarray approach and validating data on a series of selected genes by quantitative RT-PCR. Based on gene ontology and other classification criteria, we found that SSc-MVECs, while being unable to perform angiogenesis, over-express a large variety of transcripts positively involved in angiogenesis. However, such a pro-angiogenic pattern of SSc-MVECs is counterbalanced by up-regulation of pent(r)axin 3, a powerful angiogenesis inhibitor, by down-regulation of a few critical pro-angiogenic transcripts (three tissue KLKs, plexin-B1, and DSG), and by alteration of transcripts involved in signal transduction pathways. This different expression profile and identification of a few molecules able to account for impaired angiogenesis in SSc provide the focus for future attempts to correct anti-angiogenesis in SSc.
= Rho GDP dissociation inhibitor beta
= connective tissue growth factor
= extracellular matrix
= fetal calf serum
= glyceraldehyde-3-phosphate dehydrogenase
= gene ontology
= log odds ratio
= microvascular endothelial cell
= urokinase type plasminogen activator
= reverse transcription PCR
= standard deviation
= systemic sclerosis.
This work was supported by grants from Scleroderma Foundation (USA), Italian MIUR (Progetti di Ricerca di Interesse Nazionale, PRIN), and Ente Cassa di Risparmio di Firenze (Italy). Dr. Francesca Margheri was the recipient of a fellowship from the italian FIRC (Fondazione Italiana per la Ricerca sul Cancro).
- Haustein UF: Systemic sclerosis-scleroderma. Dermatol Online J. 2002, 8: 3-PubMedGoogle Scholar
- Carmeliet P, Collen D: Transgenic mouse models in angiogenesis and cardiovascular disease. J Pathol. 2000, 190: 387-405. 10.1002/(SICI)1096-9896(200002)190:3<387::AID-PATH595>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Mignatti P, Rifkin DB: Biology and biochemistry of proteinases in tumor invasion. Physiol Rev. 1993, 73: 161-195.PubMedGoogle Scholar
- Rifkin DB, Mazzieri R, Munger JS, Noguera I, Sung J: Proteolytic control of growth factor availability. APMIS. 1999, 107: 80-85.View ArticlePubMedGoogle Scholar
- Plouet J, Moro F, Bertagnolli S, Coldeboeuf N, Mazarguil H, Clamens S, Bayard F: Extracellular cleavage of the vascular endothelial growth factor 189-amino acid form by urokinase is required for its mitogenic effect. J Biol Chem. 1997, 272: 13390-13396. 10.1074/jbc.272.20.13390.View ArticlePubMedGoogle Scholar
- Bergers G, Brekken R, McMahon G, Vu TH, Itoh T, Tamaki K, Tanzawa K, Thorpe P, Itohara S, Werb S, Hanahan D: Matrix metalloproteinase-9 triggers the angiogenic switch during carcinogenesis. Nat Cell Biol. 2000, 2: 737-744. 10.1038/35036374.PubMed CentralView ArticlePubMedGoogle Scholar
- Whitfield ML, Finlay DR, Murray JI, Troyanskaya OG, Chi JT, Pergamenschikov A, McCalmont TH, Brown PO, Botstein D, Connolly MK: Systemic and cell type-specific gene expression patterns in scleroderma skin. Proc Natl Acad Sci USA. 2003, 100: 12319-12324. 10.1073/pnas.1635114100.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou X, Tan FK, Xiong M, Milewicz DM, Feghali CA, Fritzler MJ, Reveille JD, Arnett FC: Systemic sclerosis (scleroderma): specific autoantigen genes are selectively overexpressed in scleroderma fibroblasts. J Immunol. 2001, 167: 7126-7133.View ArticlePubMedGoogle Scholar
- Tan FK, Hildebrand BA, Lester MS, Stivers DN, Pounds S, Zhou X, Wallis DD, Milewicz DM, Reveille JD, Mayes MD, et al: Classification analysis of the transcriptosome of nonlesional cultured dermal fibroblasts from systemic sclerosis patients with early disease. Arthritis Rheum. 2005, 52: 865-876. 10.1002/art.20871.View ArticlePubMedGoogle Scholar
- Newton MA, Kendziorski CM, Richmond CS, Blattner FR, Tsui KW: On differential variability of expression ratios: improving statistical inference about gene expression changes from microarray data. J Comput Biol. 2001, 8: 37-52. 10.1089/106652701300099074.View ArticlePubMedGoogle Scholar
- Kerr K, Martin M, Churchill G: Analysis of variance for gene expression microarray data. J Comput Biol. 2000, 7: 819-837. 10.1089/10665270050514954.View ArticlePubMedGoogle Scholar
- Giusti B, Serratì S, Margheri F, Papucci L, Rossi L, Poggi F, Magi A, Del Rosso A, Cinelli M, Guiducci S, et al: The anti-angiogenic tissue kallikrein pattern of endothelial cell in systemic sclerosis. Arthritis Rheum. 2005, 52: 3618-3628. 10.1002/art.21383.View ArticlePubMedGoogle Scholar
- Conrotto P, Valdembri D, Corso S, Serini G, Tamagnone L, Comoglio PM, Bussolino F, Giordano S: Sema4 induces angiogenesis through recruitment by Plexin B1. Blood. 2005, 105: 4321-4229. 10.1182/blood-2004-07-2885.View ArticlePubMedGoogle Scholar
- Rusnati M, Camozzi M, Moroni E, Bottazzi B, Peri G, Indraccolo S, Amadori A, Mantovani A, Presta M: Selective recognition of fibroblast growth factor-2 by the long pentraxin PTX3 inhibits angiogenesis. Blood. 2004, 104: 92-99. 10.1182/blood-2003-10-3433.View ArticlePubMedGoogle Scholar
- Suda K, Rothen-Rutishauser B, Gunthert M, Wunderli-Allenspach H: Phenotypic characterization of human umbilical vein endothelial (ECV304) and urinary carcinoma (T24) cells: endothelial versus epithelial features. In vitro Cell Dev Biol Anim. 2001, 37: 505-514. 10.1290/1071-2690(2001)037<0505:PCOHUV>2.0.CO;2.View ArticlePubMedGoogle Scholar
- Subcommittee for Scleroderma Criteria of the American Rheumatism Association Diagnostic and Therapeutic Criteria Committee: Preliminary criteria for the classification of systemic sclerosis (scleroderma). Arthritis Rheum. 1980, 23: 581-590.View ArticleGoogle Scholar
- Manconi F, Markham R, Fraser IS: Culturing endothelial cells of microvascular origin. Methods Cell Sci. 2000, 22: 89-99. 10.1023/A:1009895723488.View ArticlePubMedGoogle Scholar
- D'Alessio S, Fibbi G, Cinelli M, Guiducci S, Del Rosso A, Margheri F, Serratì S, Pucci P, Kahaleh B, Fan P, et al: Matrix metalloproteinase 12-dependent cleavage of urokinase receptor in systemic sclerosis microvascular endothelial cells results in impaired angiogenesis. Arthritis Rheum. 2004, 50: 3275-3285. 10.1002/art.20562.View ArticlePubMedGoogle Scholar
- Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30: e15-10.1093/nar/30.4.e15.PubMed CentralView ArticlePubMedGoogle Scholar
- Kerr K, Churchill G: Experimental design for gene expression microarrays. Biostatistics. 2001, 2: 183-201. 10.1093/biostatistics/2.2.183.View ArticlePubMedGoogle Scholar
- System Biology and Bioinformatics website of theUniversity of Rostock. [http://www.sbi.uni-rostock.de]
- ArrayExpress. [http://www.ebi.ac.uk/arrayexpress/]
- Camon E, Magrane M, Barrell D, Lee V, Dimmer E, Maslen J, Binns D, Harte N, Lopez R, Apweiler R: The Gene Ontology Annotation (GOA) Database: sharing knowledge in Uniprot with Gene Ontology. Nucleic Acids Res. 2004, 32 (Database issue): D262-266. 10.1093/nar/gkh021.PubMed CentralView ArticlePubMedGoogle Scholar
- Camon E, Barrell D, Lee V, Dimmer E, Apweiler R: The Gene Ontology Annotation (GOA) Database – an integrated resource of GO annotations to the UniProt Knowledgebase. In Silico Biol. 2004, 4: 5-6.PubMedGoogle Scholar
- FTP site of the National Center for BiotechnologyInformation (NCBI). [ftp://ftp.ncbi.nlm.nih.gov/refseq/LocusLink/ARCHIVE/Loc2GO]
- Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S, et al: GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol. 2003, 4: R28-10.1186/gb-2003-4-4-r28.PubMed CentralView ArticlePubMedGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: Open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMed CentralView ArticlePubMedGoogle Scholar
- Scion Image for Windows. [http://www.scioncorp.com/pages/scion_image_windows.htm]
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralView ArticlePubMedGoogle Scholar
- Pruitt KD, Maglott DR: RefSeq and LocusLink: NCBI gene-cemtered resources. Nucleic Acid Res. 2001, 29: 137-140. 10.1093/nar/29.1.137.PubMed CentralView ArticlePubMedGoogle Scholar
- Napoleone E, Di Santo A, Bastone A, Peri G, Mantovani A, de Gaetano G, Donati MB, Lorenzet R: Long pentraxin PTX3 upregulates tissue factor expression in human endothelial cells: a novel link between vascular inflammation and clotting activation. Arterioscler Thromb Vasc Biol. 2002, 22: 782-787. 10.1161/01.ATV.0000012282.39306.64.View ArticlePubMedGoogle Scholar
- Jongeneel CV, Iseli C, Stevenson BJ, Riggins GJ, Lal A, Mackay A, Harris RA, O'Hare MJ, Neville AM, Simpson AJ, Strausberg RL: Comprehensive sampling of gene expression in human cell lines with massively parallel signature sequencing. Proc Natl Acad Sci USA. 2003, 100: 4702-4705. 10.1073/pnas.0831040100.PubMed CentralView ArticlePubMedGoogle Scholar
- Cutolo M, Grassi W, Matucci Cerinic M: Raynaud's phenomenon and the role of capillaroscopy. Arthritis Rheum. 2003, 48: 3023-3030. 10.1002/art.11310.View ArticlePubMedGoogle Scholar
- Distler O, Distler JH, Scheid A, Acker T, Hirth A, Rethage J, Michel BA, Gay RE, Muller-Ladner U, Matucci-Cerinic M, et al: Uncontrolled expression of vascular endothelial growth factor and its receptors leads to insufficient skin angiogenesis in patients with systemic sclerosis. Circ Res. 2004, 95: 109-116. 10.1161/01.RES.0000134644.89917.96.View ArticlePubMedGoogle Scholar
- Kuwana M, Okazaki Y, Yasuoka H, Kawakami Y, Ikeda Y: Defective vasculogenesis in systemic sclerosis. Lancet. 2004, 364: 603-610. 10.1016/S0140-6736(04)16853-0.View ArticlePubMedGoogle Scholar
- Basile JR, Barac A, Zhu T, Guan KL, Gutkind JS: Class IV semaphorins promote angiogenesis by stimulating Rho-initiated pathways through plexin-B. Cancer Res. 2004, 64: 5212-5224. 10.1158/0008-5472.CAN-04-0126.View ArticlePubMedGoogle Scholar
- Giordano S, Corso S, Conrotto P, Artigiani S, Gilestro G, Barberis D, Tamagnone L, Comoglio PM: The semaphorin 4D receptor controls invasive growth by coupling with Met. Nat Cell Biol. 2002, 4: 720-724. 10.1038/ncb843.View ArticlePubMedGoogle Scholar
- Schafer S, Koch PJ, Franke WW: Identification of the ubiquitous human desmoglein, Dsg2, and the expression catalogue of the desmoglein subfamily of desmosomal cadherins. Exp Cell Res. 1994, 211: 391-399. 10.1006/excr.1994.1103.View ArticlePubMedGoogle Scholar
- Bhattacharya M, Babwah AV, Ferguson SS: Small GTP-binding protein-coupled receptors. Biochem Soc Trans. 2004, 32: 1040-1044. 10.1042/BST0321040.View ArticlePubMedGoogle Scholar
- Guillemot JC, Kruskal BA, Adra CN, Zhu S, Ko JL, Burch P, Nocka K, Seetoo K, Simons E, Lim B: Targeted disruption of guanosine diphosphate-dissociation inhibitor for Rho-related proteins, GDID4: normal hematopoietic differentiation but subtle defect in superoxide production by macrophages derived from in vitro embryonal stem cell differentiation. Blood. 1996, 88: 2722-2731.PubMedGoogle Scholar
- Champion A, Picaud A, Henry Y: Reassessing the MAP3K and MAP4K relationships. Trends Plant Sci. 2004, 9: 123-129. 10.1016/j.tplants.2004.01.005.View ArticlePubMedGoogle Scholar
- Butcher HL, Kennette WA, Collins O, Zalups RK, Koropatnick J: Metallothionein mediates the level and activity of nuclear factor kB in murine fibroblasts. J Pharmacol Exp Ther. 2004, 310: 589-598. 10.1124/jpet.104.066126.View ArticlePubMedGoogle Scholar
- Kendziorski C, Irizzarry RA, Chen KS, Haag JD, Gould MN: On the utility of pooling biological samples in microarray experiments. Proc Natl Acad Sci USA. 2005, 102: 4252-4257. 10.1073/pnas.0500607102.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin-Magnette ML, Aubert J, Cabannes E, Daudin JJ: Evaluation of the gene-specific dye bias in cDNA microarray experiments. Bioinformatics. 2005, 21: 1995-2000. 10.1093/bioinformatics/bti302.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.