Transcription and splicing regulation in human umbilical vein endothelial cells under hypoxic stress conditions by exon array
- Xingyi Hang†1,
- Peiyao Li†1,
- Zhifeng Li1,
- Wubin Qu1,
- Ying Yu1,
- Hualing Li1,
- Zhiyong Shen1,
- Hao Zheng2,
- Yan Gao1,
- Yonghong Wu1,
- Minghua Deng2,
- Zhixian Sun1 and
- Chenggang Zhang1Email author
© Hang et al; licensee BioMed Central Ltd. 2009
Received: 20 September 2008
Accepted: 25 March 2009
Published: 25 March 2009
The balance between endothelial cell survival and apoptosis during stress is an important cellular process for vessel integrity and vascular homeostasis, and it is also pivotal in angiogenesis during the development of many vascular diseases. However, the underlying molecular mechanisms remain largely unknown. Although both transcription and alternative splicing are important in regulating gene expression in endothelial cells under stress, the regulatory mechanisms underlying this state and their interactions have not yet been studied on a genome-wide basis.
Human umbilical vein endothelial cells (HUVECs) were treated with cobalt chloride (CoCl2) both to mimic hypoxia and to induce cell apoptosis and alternative splicing responses. Cell apoptosis rate analysis indicated that HUVECs exposed to 300 μM CoCl2 for 24 hrs were initially counterbalancing apoptosis with cell survival. We therefore used the Affymetrix exon array system to determine genome-wide transcript- and exon-level differential expression. Other than 1583 differentially expressed transcripts, 342 alternatively spliced exons were detected and classified by different splicing types. Sixteen alternatively spliced exons were validated by RT-PCR. Furthermore, direct evidence for the ongoing balance between HUVEC survival and apoptosis was provided by Gene Ontology (GO) and protein function, as well as protein domain and pathway enrichment analyses of the differentially expressed transcripts. Importantly, a novel molecular module, in which the heat shock protein (HSP) families play a significant role, was found to be activated under mimicked hypoxia conditions. In addition, 46% of the transcripts containing stress-modulated exons were differentially expressed, indicating the possibility of combinatorial regulation of transcription and splicing.
The exon array system effectively profiles gene expression and splicing on the genome-wide scale. Based on this approach, our data suggest that transcription and splicing not only regulate gene expression, but also carry out combinational regulation of the balance between survival and apoptosis of HUVECs under mimicked hypoxia conditions. Since cell survival following the apoptotic challenge is pivotal in angiogenesis during the development of many vascular diseases, our results may advance the knowledge of multilevel gene regulation in endothelial cells under physiological and pathological conditions.
The balance between endothelial cell (EC) survival and apoptosis is an important cellular process involved in preserving blood vessel integrity and vascular homeostasis [1–4]. Lining the surface of vascular structures, ECs should endure a variety of normal or abnormal stresses that are both chemical and physical in nature. Aberrant stresses may break the dynamic balance and contribute to irreversible endothelial dysfunctions due to EC apoptosis and vessel integrity defects [4–6]. Studies have demonstrated that modulating this balance is important in the initiation and development of many vascular diseases, e.g. stroke, diabetic retinopathies, thrombosis, and atherosclerosis [2, 7–9]. Therefore, identifying the regulatory mechanisms of the survival and apoptosis of ECs may provide opportunities to improve clinical therapies for the treatment of these vascular diseases.
Transcription has been well studied and has been shown to be of considerable importance in modulating EC apoptosis [10, 11]. Alternative splicing (AS), an important molecular mechanism increasing proteome diversity via the assembly of different exons, has been reported to regulate cellular processes in endothelial systems under stress. For example, a splicing isoform of platelet endothelial cell adhesion molecule-1 (PECAM-1, a suppressor of cell apoptosis) was proven to activate the EPH receptor B2 (EPHB2) in response to the early stages of shear stress . Splicing variants of vascular endothelial growth factor (VEGF) provide a balance of pro- and anti-angiogenic regulation, and they also act as determinants of tumor angiogenesis . Importantly, one study has reported that AS, like transcription, can enable rapid and specific changes in gene expression in response to stress . Thus, elucidating the transcriptional and splicing regulation that affects EC survival and apoptosis is critical for a better understanding of endothelial function under physiological and pathological stresses.
Although many studies have focused on transcriptional and proteome profiling of ECs under stress [15, 16], no study to date has addressed splicing and multilevel regulation from a genomic standpoint. Here, human umbilical vein endothelial cells (HUVECs) were treated with 300 μM CoCl2 for 24 hrs to mimic hypoxia [17–19] and to induce cell apoptosis and alternative splicing responses, as previously described [20, 21]. An Affymetrix Human Exon 1.0 ST array system containing over 1 million exon clusters and 5.5 million features was used to profile gene expression at both the transcriptional and splicing levels. After a comparative analysis of expression between treated and normal samples, Gene Ontology (GO) and protein annotation coupled with pathway analysis provided evidence illustrating the balance between cell survival and apoptosis. Furthermore, the classification of splicing patterns and the discovery of a group of genes affected by both transcription and splicing indicated multilevel regulations representing the response of HUVECs to stress. Our data may facilitate the development of new therapeutic approaches for vascular disease treatment.
Analysis of apoptosis in CoCl2-treated HUVECs
Differentially expressed genes and functional analysis
List of the top eight functional categories of DEGs, based on the protein annotations of the UniProt databases
Alternative splicing events and protein domain analysis
Using the "Splicing Index" algorithm described in the methods, 342 probe select regions (PSRs) labeled as "core" were identified as alternatively spliced exons (See Additional file 2), and these belonged to 293 alternatively spliced transcripts. Compared with normal HUVECs, 250 PSRs with higher expression were classified as "general exon inclusion" events, while the other 92 PSRs with lower expression were considered as "general exon skipping" events. Thirty-five percent (102/293) of the alternatively spliced transcripts are supported by experimental evidence based on the NCBI RefSeq database records. Since it is possible for multiple alternative splicing events to occur in the same transcript, 13% (37/293) of the transcripts were found to contain 25% (86/342) of the alternatively spliced exons. Therefore, there was an average of 2.3 (86/37) alternatively spliced exons per transcript under mimicked hypoxia conditions in our data. A typical example is ubiquitin-associated protein 2 (UBAP2), for which five alternatively spliced exons were detected, indicating a complicated pattern of splicing regulation of UBAP2 in HUVECs under hypoxic stress.
It is well known that splicing not only provides feedback to affect transcription, but also feeds forward to modulate protein function. An InterProScan search revealed 105 protein domains in 21% (71/342) of the coding regions of alternatively spliced exons, and 70 different domains were revealed after removing the redundant records (See Additional file 3).
RT-PCR and quantitative Real-time PCR validation
Relative quantitative comparison between the RT-qPCR technique and exon array system
Gene average CT
GAPDH average CT
ΔΔCT, ΔCT-ΔCT, 0 h
Fold change to 0 ha
based on exon array systemb
18.84 ± 0.47
16.67 ± 0.23
18.84 ± 0.47
0.00 ± 0.47
20.77 ± 0.37
15.33 ± 0.30
20.77 ± 0.37
1.93 ± 0.37
21.78 ± 0.33
16.67 ± 0.23
5.12 ± 0.41
0.00 ± 0.41
18.15 ± 0.09
15.33 ± 0.30
2.82 ± 0.31
-2.30 ± 0.31
Analysis of transcription and splicing pathways
Pathway enrichment of the up- and downregulated genes based on the KEGG database analysis
HSA04010: MAPK signaling pathway
HSA04612: Antigen processing and presentation
HSA04510: Focal adhesion
HSA04810: Regulation of actin cytoskeleton
HSA04110: Cell cycle
HSA00240: Pyrimidine metabolism
HSA04350: TGF-beta signaling pathway
Splicing patterns and complex regulation between transcription and splicing
Furthermore, we found that a large proportion of alternatively spliced genes overlapped with DEGs (46% of alternatively spliced genes, 134/293), and the overlapping genes differentially expressed at both the gene and exon levels. Afterwards, a functional analysis similar to Table 1 was performed on the alternatively spliced genes. Interestingly, parallel rankings of functional categories were found between alternatively spliced genes and DEGs because of the high proportion of overlaps between these two gene sets. The top three categories of alternatively spliced genes ("alternative splicing" with P = 2.5 × 10-13, "nuclear protein" with P = 5.0 × 10-7, and "phosphorylation" with P = 1.3 × 10-15, all Fisher's exact test) were very consistent with their counterparts at the gene level. We also found that alternatively spliced genes with multiple affected exons were largely included among the overlapping genes (54%, P = 2.9 × 10-53, Chi-square test). Finally, 17 alternatively spliced genes (See Additional file 7)were transcription factors supported by the publicly available TRANSFAC 7.0 database (Fisher's exact test, P-value = 2.6 × 10-12).
Functional and pathway analyses support the conflicting balance between HUVEC survival and apoptosis
Consistent with the results of cell apoptosis analysis, the gene function and pathway analyses of DEGs revealed the conflicting balance between HUVEC survival and apoptosis. IL-8, a gene known to directly promote endothelial cell survival and angiogenesis, has been demonstrated to be highly upregulated based on exon array system and RT-qPCR analysis . This suggests that chemokines may play an important role in resisting apoptosis in HUVECs. HUVEC survival may decrease with the expression of thrombospondin-1 (THBS1), which has been reported to induce endothelial cell apoptosis and inhibit angiogenesis . On the other hand, Yang et al. reported that loss of survivin (BIRC5) increased cellular sensitivity to apoptotic stimuli and caused spontaneous apoptosis . Our results indicated that downregulation of survivin in HUVECs is highly likely to result in apoptosis via this mechanism. It has also been reported that AIFM2 reduces cell survival signaling and contributes to the onset of apoptosis . The observed upregulation of AIFM2 suggests that this gene also plays a role in promoting cell apoptosis. These gene expression patterns indicated that HUVECs struggle to avoid apoptosis in order to survive under stress. In the results of the GO analysis, it is notable that the upregulated genes are significantly enriched in the "programmed cell death" functional annotation, demonstrating the ongoing apoptosis of HUVECs.
Since genes are usually functionally organized into pathways, it is necessary to explore the gene regulation in terms of the pathways involved. As shown in Table 3, the "Focal Adhesion" pathway is largely silenced, which is congruous with the fact that adhesion-dependent endothelial cell survival is regulated by focal adhesion kinase . This silenced pathway may result in the disorder of the cellular signaling that mediates the contact between endothelial cells and the extracellular matrix during apoptosis . In addition, Kulms et al. showed that disruption of the "Actin cytoskeleton" (in the downregulated pathway) is mediated via the activation of CD95 (Fas/APO-1) during the induction of apoptosis . With regard to the upregulated pathways, the "MAPK signaling" pathway was studied by inducing apoptosis in endothelial cells via phosphorylation . The upregulated "Antigen processing and presentation" pathway is supported by the expression of many antigens, especially platelet endothelial cell adhesion molecule-1 (PECAM-1/CD31 antigen), which provides survival signals to suppress apoptosis . However, the regulation of the "Proteasome" pathway is somewhat complex because proteasome inhibitors have dual functions, either facilitating or inhibiting apoptosis . In conclusion, the expression of genes in the examined pathways presents a comprehensive illustration of the state of homeostasis between cell survival and apoptosis.
Finally, a novel heat shock protein module composed of the Hsp27, Hsp70, Hsp105, and DnaJ subfamilies was discovered to underlie the functional modulation of biological networks under stress. These heat shock proteins have been individually demonstrated to resist apoptosis in response to a variety of stimuli including hypoxia [37–40]. Figure 5 shows that the 70 kDa heat shock protein 1A (HSPA1A) may function together with other heat shock proteins to form a protein complex that more effectively inhibits apoptosis. Notably, HSPA1A has been reported to confer resistance to apoptosis in conjunction with other heat shock proteins , which is consistent with its feature as a hub of the scale-free module. This novel module suggests that heat shock proteins and their collective regulation may be crucial to controlling HUVEC survival and apoptosis.
Complex regulation of transcription and splicing in stress induced HUVECs
Both dependent and independent regulations of transcription and splicing usually coexist under most physiological and pathological conditions. Based on the observation of a higher rate of overlapping between DEGs and alternatively spliced genes than that found in other studies [42, 43], we expect the possibility of combinatorial regulation between transcription and splicing in stress induced HUVECs. Although we also found that the general splicing patterns are highly correlated with gene expression levels, the exact molecular mechanism of the coupling regulation is still unknown. We hypothesize that splicing may modify the transcription activity or RNA stability [44, 45], while transcription may change the splicing efficiency . On one hand, alternative splicing of mRNA can change RNA stability, which in turn will probably affect the expression levels of the gene transcripts with different RNA stability . On the other hand, it is also possible that different expression levels of the upstream genes of splicing factors facilitate or inhibit splicing machinery by influencing spliceosome assembly or the cis-elements during the splicing process. These two aspects of regulations could both possibly result in the high degree of correlation between splicing patterns and transcriptional expression. Therefore, it is reasonable to speculate that HUVECs may utilize the combinatorial regulation of transcription and splicing to modulate the cellular response to stress finely and efficiently.
Transcription and splicing may be independent processes, but there are still possible correlations at specific spatiotemporal stages of the cellular response. In our results, 17 differentially expressed transcription factors (See Additional file 7) were detected as alternatively spliced genes(Fisher's exact test, P-value = 2.6 × 10-12). On the other hand, 15 splicing factors (See Additional file 7), including 6 SR proteins and 9 hnRNP proteins, were detected as DEGs (Fisher's exact test, P-value = 1.3 × 10-6). The existence of two possible regulatory mechanisms for these transcription factors and splicing factors can be conjectured: 1) the 17 alternatively spliced transcription factors are possible targets of splicing factors; 2) the 15 differentially expressed splicing factors are possible targets of transcription factors. If the differential expression of splicing factors directly influences the splicing efficiency and in turn triggers the alternative splicing of transcription factors, a loop of feedback regulation can then be established in response to stress. Since it is difficult to reveal the exact regulatory mechanisms underlying these trans-factors and their targets, further studies are needed to explain the regulatory model of the complex regulation under stress in the future.
Alternative splicing can influence biological networks through domain architectures
Since no significant enrichment of alternatively spliced genes was found in the KEGG pathways, splicing may follow a different set of regulatory rules than transcription in pathways. Alternative splicing can expand the protein repertoire and influence protein function by altering protein domains. Melissa et al. reported that 7,179 of 22,218 human genes in the Ensembl database encoded two or more different proteins. Of these, 2,229 genes encoded proteins with different PFAM domain architectures . The affected domains in the coding regions of alternatively spliced exons confirmed the existence of changes in the transcriptome and proteome resulting from alterations in the domain architecture of biological networks . We found that alternative splicing may influence transcription through the gain or loss of promoter binding domains. For example, the number of zinc finger domains (IPR007087) decreased in zinc finger protein 589 (ZNF589), whose transcription factor activity depends on the number of domain repeats. The same phenomenon was also found in the WD-40 repeat domain (IPR001680) of the SH3KBP1 and RRP9 genes. In our results, the DNA-binding domain HMG-I(Y) (IPR000637) was lost in the high mobility group AT-hook 2 (HMGA2). Previous studies have demonstrated that the domain HMG-I(Y) functions as part of a hypoxia-induced enhanceosome, promoting the transcription of COX-2 in HUVECs . Defects in the HMG-I(Y) DNA-binding domain will disorganize the transcriptional regulation under stress. The MAM domain (IPR000998) in neuropilin 1 (NRP1), representing adhesive function, may be altered to induce endothelial dysfunction in response to stress. These changes of domains were analyzed based on the coding regions of alternatively spliced exons (See Additional file 3).
In this study, HUVECs were incubated with 300 μM CoCl2 for 24 hrs to induce the balance between cell survival and apoptosis, followed by a genome-wide expression profiling of transcription and splicing by exon array system. Functional and pathway analyses of gene levels and exon levels demonstrated the importance of transcription and splicing regulation in cellular processes. Evidence from the splicing classifications and the overlap between the two levels suggested a combinatorial regulation. Because very few studies have investigated splicing regulation in endothelial cell survival and apoptosis, elucidating the underlying mechanisms associated with these phenomena is critical for a better understanding of vascular biology under normal and pathological conditions.
Cell culture and cell apoptosis analysis
HUVECs were purchased from Cascade Biologics (USA) and cultured in Medium 200 supplemented with Low Serum Growth Supplement (LSGS) (Cascade Biologics, USA) in a CO2 incubator (5% CO2) at 37°C. The cells were treated with different concentrations of CoCl2 (100 μM, 300 μM, 600 μM, 900 μM) (Sigma, USA) for 0, 12, 24, 36 and 48 hrs to mimic hypoxia. The cells were then incubated with fluorescein isothiocyanate-conjugated Annexin V (A-FITC) and propidium iodide (PI) using the Apoptest kit (Jiancheng, China) according to the manufacturer's instructions. Flow cytometry analysis was performed using the FACSCalibur system (Becton Dickinson, USA). The data were analyzed using CellQuest software to estimate the apoptosis rate at different time points.
Sample preparation and array hybridization
After being cultured under normoxia or mimicked hypoxia (300 μM CoCl2 for 24 hrs), total RNA was extracted from the HUVECs using the TRIzol reagent, according to the manufacturer's protocol (Invitrogen, USA). Total RNA was dissolved in an appropriate volume of DEPC-treated water following A260/A280 measurement, while the total RNA integrity was evaluated by electrophoresis in a denaturing gel. The RNA samples were further purified using DNase (TaKaRa, Japan). For each experimental condition, three independent replicate samples were obtained for exon array analysis. For each sample, 1 μg of RNA was processed using the Affymetrix GeneChip® Whole Transcript Sense Target Labeling Assay. The GeneChip® WT cDNA Synthesis Kit, the WT cDNA Amplification Kit, and the WT Terminal Labeling Kit (Affymetrix, Inc., Santa Clara, CA) were used for the sample preparation. 8 μg of cDNA were used for the second cycle cDNA reaction. Hybridization cocktails containing 3–4 μg of fragmented, end-labeled cDNA were applied to the GeneChip® Human Exon 1.0 ST arrays. Hybridization was performed for 16 hrs using the MES_EukGE-WS2v5_450-DEV fluidics wash and stain script. The arrays were scanned using the Affymetrix GCS 3000 7G and Gene-Chip Operating Software v1.3 to produce the intensity files.
RT-PCR and quantitative Real-time RT-PCR
1 μg of each RNA sample was used for first strand cDNA synthesis using SuperScript II reverse transcriptase (Invitrogen, USA) and a combination of random hexamer primers and oligo-dT in a total volume of 10 μl. PCR was carried out using 2 μl of cDNA, with specific primers flanking the constitutive exons, and ExTaq Polymerase (TaKaRa, Japan) in a volume of 25 μl. The conditions for PCR amplification were denaturation at 95°C for 5 min, 32 cycles of 95°C for 30 sec, 55°C for 30 sec, and 72°C for 45 sec, followed by a final elongation step at 72°C for 7 min. The PCR products were then separated on 1.5% agarose gels. The RT-PCR products were gel-purified using a PCR purification kit (Promega, USA) and subcloned into the pGEM-T Easy Vector (Promega, USA) for direct sequencing to validate the transcript variants.
1 μl of each cDNA product was used for quantitative real-time PCR amplification with SYBR Green PCR Master Mix (TianGen). The primers were designed and verified by the primer specificity-checking program MFEprimer http://biocompute.bmi.ac.cn/MFEprimer. PCR was carried out with an iCycler Real-time PCR detection system (Bio-Rad) under the following conditions (40 cycles): 95°C for 2 min, 95°C for 30 sec, 57°C for 30 sec, and 68°C for 30 sec. SYBR Green analyses were followed by dissociation curves in a temperature range of 60°C~90°C to assess the amplification specificity. Each sample was tested in triplicate and quantified according to the mean expression values obtained for both samples.
Low level analysis of the exon array
Low-level analysis of the optical intensity files of the exon array (".CEL" format) was performed by Affymetrix Power Tools (APT). Background noise was detected by the "Detection above Background (DABG)" algorithm. Normalization was performed using the "quantile normalization" algorithm for both the exon and gene levels. The "Probe Logarithmic Intensity Error Estimation" (PLIER) algorithm was used to estimate exon signals based on probe intensities. At the gene level, a variant algorithm called "Iter-PLIER" was used to summarize gene signals from probeset intensities. The "Iter-PLIER" algorithm can discard probesets with inconsistent signals to avoid low-weighted effects introduced by differentially expressed exons.
Hierarchical filtering was then performed to eliminate noise and outliers at both the gene and exon levels. At the exon level, only the probesets considered "Present" (DABG P < 0.05) in at least 50% of the samples in either group were reserved. At the gene level, only the "core" meta-probesets with high confidence were used to estimate gene signals. The differentially expressed genes were considered acceptable based on two principles. First, genes with more than 50% of the "core" exons designated as "Present" (DABG P < 0.05) should appear in more than 50% of the samples in both groups. Second, the gene signals needed to exceed 100. We subsequently removed the probesets labeled as potential cross-hybridization targets based on Affymetrix CSV annotation files.
Detection of differentially expressed genes and alternative splicing
The absolute value of SI represented the magnitude of difference of the exon inclusion rate between the two groups. To identify the significant alternatively spliced exons, a Student's t-test was used to compare TNS values between the two groups. Finally, the high proportion of true positives, with P-value < 0.015 and fold change magnitudes > 0.5, were retained as potential alternatively spliced exons.
The raw ".CEL" files and normalized data at both the gene and exon levels have been deposited in the Gene Expression Omnibus (GEO) of the National Center for Biotechnology Information http://www.ncbi.nlm.nih.gov/geo under GEO Series record GSE12546.
Visualization and classification of alternative splicingevents
Before validating the exon array data by various approaches, an expert investigation on gene structure and genomic context was carried out to assess the positions and surrounding mRNA/cDNA sequences of alternatively spliced exons. The Blat program  was used to map alternatively spliced exons in the UCSC Genome Browser http://genome.ucsc.edu/ referred to the mRNA/cDNA sequences (from the NCBI RefSeq and GenBank databases) or expressed sequence tags (ESTs). Alternatively spliced multi-exon genes were classified into six splicing patterns according to the relative positions of the affected probe selected regions (PSRs) in exons and genes based on the sequence mapping. These classifications were cassette exons, namely exon inclusion and exon skipping, alternative promoters, alternative polyadenylation, alternative donor sites, alternative acceptor sites and intron retention.
Function and pathway analysis
GO, protein function, and pathway enrichment analyses were carried out by the DAVID tool http://david.abcc.ncifcrf.gov/. DEGs and alternatively spliced genes were mapped to the KEGG database using GenMAPP software, in order to visualize their distributions in the pathways .
After detecting alternatively spliced exons, their sequences and gene annotations were obtained from the Affymetrix website http://www.affymetrix.com/support/technical/byproduct.affx?product=huexon-st. The protein sequences of the coding regions of alternatively spliced exons were extracted from the NCBI RefSeq database by a in-house developed Perl program . The InterProScan software was used to search protein domains via the interfaces of the PFAM, PROSITE, PRODOM, and SMART databases .
Literature mining for functional modules
The purpose of the analysis is to find functional modules from complex biological networks. The functional module was defined as a part of a biological network with specific functions and topological features . The nodes represent genes, and the links represent regulatory relationships between genes in the modules. A two-step literature mining strategy was performed on up- and downregulated genes to find activated functional modules in affected HUVECs. First, we used the cytoscape plugin "Agilent Literature Search" to construct the biological networks by a literature mining algorithm . Only direct regulatory relationships between genes were preserved in building the network. Second, some orphan nodes and fake links were manually removed by checking relevant sentences in the obtained literatures in the first step. During the manual module check, the nodes were annotated by description from NCBI Entrez Gene , and the links were classified by regulatory relationships stated in the sentences from the relevant literature.
The authors sincerely wish to thank the anonymous reviewers for their fruitful suggestions to improve this manuscript. We thank CapitalBio Inc. to provide microarray service. This work was supported by the National Basic Research Project of China (973 program) (2006CB504100, 2003CB715900); the General Program (30800196, 30771230, 30772293, 30600365) of the National Natural Science Foundation of China; the Major Program for Science and Technology Research of Beijing Municipal Bureau (7061004).
- Duval H, Harris M, Li J, Johnson N, Print C: New insights into the function and regulation of endothelial cell apoptosis. Angiogenesis. 2003, 6 (3): 171-183. 10.1023/B:AGEN.0000021390.09275.bc.View ArticlePubMedGoogle Scholar
- Winn RK, Harlan JM: The role of endothelial cell apoptosis in inflammatory and immune diseases. J Thromb Haemost. 2005, 3 (8): 1815-1824. 10.1111/j.1538-7836.2005.01378.x.View ArticlePubMedGoogle Scholar
- Santoro MM, Samuel T, Mitchell T, Reed JC, Stainier DYR: Birc2 (cIap1) regulates endothelial cell integrity and blood vessel homeostasis. Nat Genet. 2007, 39 (11): 1397-1402. 10.1038/ng.2007.8.View ArticlePubMedGoogle Scholar
- Dimmeler S, Zeiher AM: Endothelial cell apoptosis in angiogenesis and vessel regression. Circ Res. 2000, 87 (6): 434-439.View ArticlePubMedGoogle Scholar
- Kobayashi N, DeLano FA, Schmid-Schonbein GW: Oxidative stress promotes endothelial cell apoptosis and loss of microvessels in the spontaneously hypertensive rats. Arterioscler Thromb Vasc Biol. 2005, 25 (10): 2114-2121. 10.1161/01.ATV.0000178993.13222.f2.View ArticlePubMedGoogle Scholar
- Sudoh N, Toba K, Akishita M, Ako J, Hashimoto M, Iijima K, Kim S, Liang YQ, Ohike Y, Watanabe T, et al: Estrogen prevents oxidative stress-induced endothelial cell apoptosis in rats. Circulation. 2001, 103 (5): 724-729.View ArticlePubMedGoogle Scholar
- Ten VS, Pinsky DJ: Endothelial response to hypoxia: physiologic adaptation and pathologic dysfunction. Curr Opin Crit Care. 2002, 8 (3): 242-250. 10.1097/00075198-200206000-00008.View ArticlePubMedGoogle Scholar
- Carmeliet P: Angiogenesis in life, disease and medicine. Nature. 2005, 438 (7070): 932-936. 10.1038/nature04478.View ArticlePubMedGoogle Scholar
- Tricot O, Mallat Z, Heymes C, Belmin J, Leseche G, Tedgui A: Relation between endothelial cell apoptosis and blood flow direction in human atherosclerotic plaques. Circulation. 2000, 101 (21): 2450-2453.View ArticlePubMedGoogle Scholar
- Birdsey GM, Dryden NH, Amsellem V, Gebhardt F, Sahnan K, Haskard DO, Dejana E, Mason JC, Randi AM: Transcription factor Erg regulates angiogenesis and endothelial apoptosis through VE-cadherin. Blood. 2008, 111 (7): 3498-3506. 10.1182/blood-2007-08-105346.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang N, Stemerman MB: Ref-1 and transcriptional control of endothelial apoptosis. Circ Res. 2001, 88 (12): 1223-1225. 10.1161/hh1201.093162.View ArticlePubMedGoogle Scholar
- Sheibani N, Sorenson CM, Frazier WA: Differential modulation of cadherin-mediated cell-cell adhesion by platelet endothelial cell adhesion molecule-1 isoforms through activation of extracellular regulated kinases. Molecular Biology of the Cell. 2000, 11 (8): 2793-2802.PubMed CentralView ArticlePubMedGoogle Scholar
- Ladomery MR, Harper SJ, Bates DO: Alternative splicing in angiogenesis: The vascular endothelial growth factor paradigm. Cancer Letters. 2007, 249 (2): 133-142. 10.1016/j.canlet.2006.08.015.View ArticlePubMedGoogle Scholar
- Pleiss JA, Whitworth GB, Bergkessel M, Guthrie C: Rapid, transcript-specific changes in splicing in response to environmental stress. Mol Cell. 2007, 27 (6): 928-937. 10.1016/j.molcel.2007.07.018.PubMed CentralView ArticlePubMedGoogle Scholar
- Ning W, Chu TJ, Li CJ, Choi AMK, Peters DG: Genome-wide analysis of the endothelial transcriptome under short-term chronic hypoxia. Physiological Genomics. 2004, 18 (1): 70-78. 10.1152/physiolgenomics.00221.2003.View ArticlePubMedGoogle Scholar
- Scheurer SB, Rybak JN, Rosli C, Neri D, Elia G: Modulation of gene expression by hypoxia in human umbilical cord vein endothelial cells: A transcriptomic and proteomic study. Proteomics. 2004, 4 (6): 1737-1760. 10.1002/pmic.200300689.View ArticlePubMedGoogle Scholar
- Chachami G, Simos G, Hatziefthimiou A, Bonanou S, Molyvdas PA, Paraskeva E: Cobalt induces hypoxia-inducible factor-1alpha expression in airway smooth muscle cells by a reactive oxygen species- and PI3K-dependent mechanism. Am J Respir Cell Mol Biol. 2004, 31 (5): 544-551. 10.1165/rcmb.2003-0426OC.View ArticlePubMedGoogle Scholar
- Yuan Y, Hilliard G, Ferguson T, Millhorn DE: Cobalt inhibits the interaction between hypoxia-inducible factor-alpha and von Hippel-Lindau protein by direct binding to hypoxia-inducible factor-alpha. J Biol Chem. 2003, 278 (18): 15911-15916. 10.1074/jbc.M300463200.View ArticlePubMedGoogle Scholar
- Kim KS, Rajagopal V, Gonsalves C, Johnson C, Kalra VK: A novel role of hypoxia-inducible factor in cobalt chloride- and hypoxia-mediated expression of IL-8 chemokine in human endothelial cells. J Immunol. 2006, 177 (10): 7211-7224.View ArticlePubMedGoogle Scholar
- Malhotra R, Lin Z, Vincenz C, Brosius FC: Hypoxia induces apoptosis via two independent pathways in Jurkat cells: differential regulation by glucose. Am J Physiol Cell Physiol. 2001, 281 (5): C1596-1603.PubMedGoogle Scholar
- Munch C, Zhu BG, Mink A, Seefried U, Riepe MW, Ludolph AC, Meyer T: Chemical hypoxia facilitates alternative splicing of EAAT2 in presymptomatic APP23 transgenic mice. Neurochem Res. 2008, 33 (6): 1005-1010. 10.1007/s11064-007-9540-5.View ArticlePubMedGoogle Scholar
- Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang HZ, Lopez R, Magrane M, et al: UniProt: the Universal Protein knowledgebase. Nucleic Acids Research. 2004, 32: D115-D119. 10.1093/nar/gkh131.PubMed CentralView ArticlePubMedGoogle Scholar
- Blencowe BJ: Alternative splicing: New insights from global analyses. Cell. 2006, 126 (1): 37-47. 10.1016/j.cell.2006.06.023.View ArticlePubMedGoogle Scholar
- Minegishi N, Ohta J, Suwabe N, Nakauchi H, Ishihara H, Hayashi N, Yamamoto M: Alternative promoters regulate transcription of the mouse GATA-2 gene. J Biol Chem. 1998, 273 (6): 3625-3634. 10.1074/jbc.273.6.3625.View ArticlePubMedGoogle Scholar
- Xin D, Hu L, Kong X: Alternative promoters influence alternative splicing at the genomic level. PLoS ONE. 2008, 3 (6): e2377-10.1371/journal.pone.0002377.PubMed CentralView ArticlePubMedGoogle Scholar
- Bemmo A, Benovoy D, Kwan T, Gaffney DJ, Jensen RV, Majewski J: Gene expression and isoform variation analysis using Affymetrix Exon Arrays. BMC Genomics. 2008, 9 (529): 529-10.1186/1471-2164-9-529.PubMed CentralView ArticlePubMedGoogle Scholar
- Li AH, Dubey S, Varney ML, Dave BJ, Singh RK: IL-8 directly enhanced endothelial cell survival, proliferation, and matrix metalloproteinases production and regulated angiogenesis. J Immunol. 2003, 170 (6): 3369-76.View ArticlePubMedGoogle Scholar
- Nor JE, Mitra RS, Sutorik MM, Mooney DJ, Castle VP, Polverini PJ: Thrombospondin-1 induces endothelial cell apoptosis and inhibits angiogenesis by activating the caspase death pathway. Journal of Vascular Research. 2000, 37 (3): 209-218. 10.1159/000025733.View ArticlePubMedGoogle Scholar
- Yang D, Welm A, Bishop JM: Cell division and cell survival in the absence of survivin. Proceedings of the National Academy of Sciences of the United States of America. 2004, 101 (42): 15100-15105. 10.1073/pnas.0406665101.PubMed CentralView ArticlePubMedGoogle Scholar
- Gong M, Hay S, Marshall KR, Munro AW, Scrutton NS: DNA binding suppresses human AIF-M2 activity and provides a connection between redox chemistry, reactive oxygen species, and apoptosis. J Biol Chem. 2007, 282 (41): 30331-30340. 10.1074/jbc.M703713200.View ArticlePubMedGoogle Scholar
- Frisch SM, Vuori K, Ruoslahti E, ChanHui PY: Control of adhesion-dependent cell survival by focal adhesion kinase. Journal of Cell Biology. 1996, 134 (3): 793-799. 10.1083/jcb.134.3.793.View ArticlePubMedGoogle Scholar
- Cance WG, Golubovskaya VM: Focal adhesion kinase versus p53: apoptosis or survival?. Sci Signal. 2008, 1 (20): pe22-10.1126/stke.120pe22.PubMed CentralView ArticlePubMedGoogle Scholar
- Kulms D, Dussmann H, Poppelmann B, Stander S, Schwarz A, Schwarz T: Apoptosis induced by disruption of the actin cytoskeleton is mediated via activation of CD95 (Fas/APO-1). Cell Death and Differentiation. 2002, 9 (6): 598-608. 10.1038/sj.cdd.4401002.View ArticlePubMedGoogle Scholar
- Grethe S, Ares MP, Andersson T, Porn-Ares MI: p38 MAPK mediates TNF-induced apoptosis in endothelial cells via phosphorylation and downregulation of Bcl-x(L). Experimental Cell Research. 2004, 298 (2): 632-642. 10.1016/j.yexcr.2004.05.007.View ArticlePubMedGoogle Scholar
- Gao CJ, Sun WY, Christofidou-Solomidou M, Sawada M, Newman DK, Bergom C, Albelda SM, Matsuyama S, Newman PJ: PECAM-1 functions as a specific and potent inhibitor of mitochondrial-dependent apoptosis. Blood. 2003, 102 (1): 169-179. 10.1182/blood-2003-01-0003.View ArticlePubMedGoogle Scholar
- Wojcik C: Regulation of apoptosis by the ubiquitin and proteasome pathway. J Cell Mol Med. 2002, 6 (1): 25-48. 10.1111/j.1582-4934.2002.tb00309.x.View ArticlePubMedGoogle Scholar
- Hatayama T, Yamagishi N, Minobe E, Sakai K: Role of hsp105 in protection against stress-induced apoptosis in neuronal PC12 cells. Biochem Biophys Res Commun. 2001, 288 (3): 528-534. 10.1006/bbrc.2001.5802.View ArticlePubMedGoogle Scholar
- Mehlen P, SchulzeOsthoff K, Arrigo AP: Small stress proteins as novel regulators of apoptosis – Heat shock protein 27 blocks Fas/APO-1- and staurosporine-induced cell death. Journal of Biological Chemistry. 1996, 271 (28): 16510-16514. 10.1074/jbc.271.28.16510.View ArticlePubMedGoogle Scholar
- Terui K, Haga S, Enosawa S, Ohnuma N, Ozaki M: Hypoxia/re-oxygenation-induced, redox-dependent activation of STAT1 (signal transducer and activator of transcription 1) confers resistance to apoptotic cell death via hsp70 induction. Biochemical Journal. 2004, 380: 203-209. 10.1042/BJ20031891.PubMed CentralView ArticlePubMedGoogle Scholar
- Mailhos C, Howard MK, Latchman DS: Heat shock proteins hsp90 and hsp70 protect neuronal cells from thermal stress but not from programmed cell death. J Neurochem. 1994, 63 (5): 1787-1795.View ArticlePubMedGoogle Scholar
- Stankiewicz AR, Lachapelle G, Foo CPZ, Radicioni SM, Mosser DD: Hsp70 inhibits heat-induced apoptosis upstream of mitochondria by preventing Bax translocation. Journal of Biological Chemistry. 2005, 280 (46): 38729-38739. 10.1074/jbc.M509497200.View ArticlePubMedGoogle Scholar
- Gardina P, Clark T, Shimada B, Staples M, Yang Q, Veitch J, Schweitzer A, Awad T, Sugnet C, Dee S, et al: Alternative splicing and differential gene expression in colon cancer detected by a whole genome exon array. BMC Genomics. 2006, 7 (1): 325-10.1186/1471-2164-7-325.PubMed CentralView ArticlePubMedGoogle Scholar
- Ip JY, Tong A, Pan Q, Topp JD, Blencowe BJ, Lynch KW: Global analysis of alternative splicing during T-cell activation. RNA. 2007, 13 (4): 563-572. 10.1261/rna.457207.PubMed CentralView ArticlePubMedGoogle Scholar
- Nobuhiro F, Yoshitaka N, Popiel HA, Hiroki K, Masamitsu Y, Tatsushi T: Alternative splicing regulates the transcriptional activity of Drosophila heat shock transcription factor in response to heat/cold stress. FEBS letters. 2005, 579 (17): 3842-3848. 10.1016/j.febslet.2005.05.074.View ArticleGoogle Scholar
- Sellers RS, Luchin AI, Richard V, Brena RM, Lima D, Rosol TJ: Alternative splicing of parathyroid hormone-related protein mRNA: expression and stability. J Mol Endocrinol. 2004, 33 (1): 227-241. 10.1677/jme.0.0330227.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosonina E, Bakowski MA, McCracken S, Blencowe BJ: Transcriptional activators control splicing and 3'-end cleavage levels. J Biol Chem. 2003, 278 (44): 43034-43040. 10.1074/jbc.M307289200.View ArticlePubMedGoogle Scholar
- Hollams EM, Giles KM, Thomson AM, Leedman PJ: MRNA stability and the control of gene expression: implications for human disease. Neurochem Res. 2002, 27 (10): 957-980. 10.1023/A:1020992418511.View ArticlePubMedGoogle Scholar
- Melissa C, Mario A, Thomas L, Benno S: The consequences of alternative splicing on biological pathways. ISMB. 2006, [http://ismb2006.cbi.cnptia.embrapa.br/plos/13_Cline_10_July.pdf]Google Scholar
- Resch A, Xing Y, Modrek B, Gorlick M, Riley R, Lee C: Assessing the impact of alternative splicing on domain interactions in the human proteome. Journal of Proteome Research. 2004, 3 (1): 76-83. 10.1021/pr034064v.View ArticlePubMedGoogle Scholar
- Ji YS, Xu Q, Schmedtje JF: Hypoxia induces high-mobility-group protein I(Y) and transcription of the cyclooxygenase-2 gene in human vascular endothelium. Circulation Research. 1998, 83 (3): 295-304.View ArticlePubMedGoogle Scholar
- Qu W, Shen Z, Zhao D, Yang Y, Zhang C: MFEprimer: multiple factor evaluation of the specificity of PCR primers. Bioinformatics. 2009, 25 (2): 276-278. 10.1093/bioinformatics/btn614.View ArticlePubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proceedings of the National Academy of Sciences of the United States of America. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMed CentralView ArticlePubMedGoogle Scholar
- Yeo GW, Xu X, Liang TY, Muotri AR, Carson CT, Coufal NG, Gage FH: Alternative splicing events identified in human embryonic stem cells and neural progenitors. PLoS Comput Biol. 2007, 3 (10): 1951-1967. 10.1371/journal.pcbi.0030196.View ArticlePubMedGoogle Scholar
- Zheng H, Hang X, Zhu J, Qian M, Qu W, Zhang C, Deng M: REMAS: a new regression model to identify alternative splicing events from exon array data. BMC Bioinformatics. 2009, 10 (Suppl 1): S18-10.1186/1471-2105-10-S1-S18.PubMed CentralView ArticlePubMedGoogle Scholar
- Srinivasan K, Shiue L, Hayes JD, Centers R, Fitzwater S, Loewen R, Edmondson LR, Bryant J, Smith M, Rommelfanger C, et al: Detection and measurement of alternative splicing using splicing-sensitive microarrays. Methods. 2005, 37 (4): 345-359. 10.1016/j.ymeth.2005.09.007.View ArticlePubMedGoogle Scholar
- Clark TA, Sugnet CW, Ares M: Genomewide analysis of mRNA processing in yeast using splicing-specific microarrays. Science. 2002, 296 (5569): 907-910. 10.1126/science.1069415.View ArticlePubMedGoogle Scholar
- Kent WJ: BLAT – The BLAST-like alignment tool. Genome Research. 2002, 12 (4): 656-664.PubMed CentralView ArticlePubMedGoogle Scholar
- Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for annotation, visualization, and integrated discovery. Genome Biol. 2003, 4 (5): P3-10.1186/gb-2003-4-5-p3.View ArticlePubMedGoogle Scholar
- Salomonis N, Hanspers K, Zambon AC, Vranizan K, Lawlor SC, Dahlquist KD, Doniger SW, Stuart J, Conklin BR, Pico AR: GenMAPP 2: new features and resources for pathway analysis. BMC Bioinformatics. 2007, 8: 217-10.1186/1471-2105-8-217.PubMed CentralView ArticlePubMedGoogle Scholar
- Pruitt KD, Tatusova T, Maglott DR: NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Research. 2005, 33: D501-D504. 10.1093/nar/gki025.PubMed CentralView ArticlePubMedGoogle Scholar
- Zdobnov EM, Apweiler R: InterProScan – an integration platform for the signature-recognition methods in InterPro. Bioinformatics. 2001, 17 (9): 847-848. 10.1093/bioinformatics/17.9.847.View ArticlePubMedGoogle Scholar
- Sauro HM: Modularity defined. Mol Syst Biol. 2008, 4: 166-10.1038/msb.2008.3.PubMed CentralView ArticlePubMedGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T: Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13 (11): 2498-2504. 10.1101/gr.1239303.PubMed CentralView ArticlePubMedGoogle Scholar
- Sayers EW, Barrett T, Benson DA, Bryant SH, Canese K, Chetvernin V, Church DM, DiCuccio M, Edgar R, Federhen S, et al: Database resources of the National Center for Biotechnology Information. Nucl Acids Res. 2009, 37 (suppl_1): D5-15. 10.1093/nar/gkn741.PubMed CentralView ArticlePubMedGoogle 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.