Variable promoter methylation contributes to differential expression of key genes in human placenta-derived venous and arterial endothelial cells
BMC Genomics volume 14, Article number: 475 (2013)
The endothelial compartment, comprising arterial, venous and lymphatic cell types, is established prenatally in association with rapid phenotypic and functional changes. The molecular mechanisms underpinning this process in utero have yet to be fully elucidated. The aim of this study was to investigate the potential for DNA methylation to act as a driver of the specific gene expression profiles of arterial and venous endothelial cells.
Placenta-derived venous and arterial endothelial cells were collected at birth prior to culturing. DNA methylation was measured at >450,000 CpG sites in parallel with expression measurements taken from 25,000 annotated genes. A consistent set of genomic loci was found to show coordinate differential methylation between the arterial and venous cell types. This included many loci previously not investigated in relation to endothelial function. An inverse relationship was observed between gene expression and promoter methylation levels for a limited subset of genes implicated in endothelial function, including NOS3, encoding endothelial Nitric Oxide Synthase.
Endothelial cells derived from the placental vasculature at birth contain widespread methylation of key regulatory genes. These are candidates involved in the specification of different endothelial cell types and represent potential target genes for environmentally mediated epigenetic disruption in utero in association with cardiovascular disease risk later in life.
Mounting evidence linking environmental exposures in early life to later risk of cardiovascular disease has led to intense interest in the process of vasculature development in utero[1–3]. Primarily, arteries and veins are defined by physiologic factors such as the direction and pressure of blood flow and by functional and anatomical differences such as the arrangement of smooth muscle cells around the vessels. In general, arteries carry oxygenated blood, and have tighter endothelial junctions, whilst veins carry deoxygenated blood and have looser endothelial junctions . The feto-placental vascular system differs from that of most other human organs, because the arteries carry the deoxygenated blood coming from the fetus whilst the veins carry the oxygen enriched blood.
Not surprisingly, given their different physiologic functions, the identity of arterial and venous endothelial cells is established before the onset of circulation [5–7] in association with distinct gene expression signatures that define endothelial cell identity [8–10]. This includes endothelial nitric oxide synthase (eNOS) encoded by NOS3, von Willebrand factor (vWF), vascular endothelial cadherin (VE-cadherin, CDH5), intracellular adhesion molecular-2 (ICAM-2), endothelial growth factor receptor tyrosine kinases VEGF-R1 (Flt-1) and VEGF-R2 (Flk-1, KDR), angiopoietin receptors Tie-1 and Tie2 (TEK) and NOTCH4 [Reviewed in [11, 12]]. In addition, a number of transcription factors have been shown to be preferentially expressed in endothelial progenitor cells and mature endothelial cells and have been argued to orchestrate the expression of such genes . Despite this, the associated epigenetic regulators determining this specific gene expression profile are largely unclear.
Numerous lines of evidence suggest that the endothelial compartment represents a potential target tissue for transmission of environmentally mediated risk of complex disease, potentially via a process of epigenetic disruption in early development [2, 13]. Epigenetic mechanisms, including DNA methylation, are now widely accepted to underpin developmentally regulated changes in cell morphology and function , however little is known about the role of such modifications, and their relative plasticity, during endothelial cell development. It is clear that one of the key regulators of vascular function, endothelial nitric oxide synthase (eNOS) is under epigenetic control by several mechanisms, including DNA methylation . Indeed the chromatin structure of the eNOS promoter is transcriptionally permissive only in cells of the endothelial compartment, such as those isolated from human umbilical vein, mouse aortic, and pooled human dermal microvascular endothelial cells [15, 16]. Altered chromatin has also been reported for other endothelial-specific genes including vWF, NOTCH4. Furthermore, studies have shown changes in epigenetic marks induced by oxidative stress (i.e. hypoxia) in VEGF and NOS3, supporting the potential epigenetic roles in tissue-specific regulation of those genes. Recent genome-scale comparison of DNA methylation in dermal derived lymphatic and blood-derived endothelial cells has highlighted the role of differential methylation in the specification of endothelial function , but beyond these limited data, little is known about the relationship between epigenetic modifications and gene expression underpinning endothelial phenotype or developmental stage at the genome-wide level. As a first step towards addressing these questions, we have examined the relationship between genome-scale DNA methylation and gene expression in purified human placental arterial and venous endothelial cells (HPAEC and HPVEC).
Results and Discussion
Placenta-derived venous and arterial endothelial cells play especially important functional roles (e.g. nutrition, cholesterol delivery) as part of the fetal-maternal supply line [22, 23]. Furthermore, several physiological differences exist between these two cell types  and they differ in their degree of maturity. Whereas placental arterial endothelial cells have a mature arterial phenotype (fully differentiated), venous derived cells show a juvenile (less differentiated) phenotype, potentially representing a pool of endothelial progenitor cells .
We compared DNA methylation in 9 purified cell populations each of Human Placental Arterial Endothelial Cells (HPAEC) and venous equivalents (HPVEC). After removing probes that did not pass a quality control cutoff (i.e. probes with >0.05 p-value detection) and probes on sex chromosomes, a total of 351,952 probes remained that were common to all datasets for inclusion in subsequent analyses. β-values (DNA methylation values between 0 and 1, approximating 0-100% methylation) were calculated for each probe in HPAEC and HPVEC samples.
General hypomethylation in placental venous relative to arterial endothelial cells
As a first comparison of methylation within the two cell types, we calculated an average β-value across the entire dataset as a proxy for global methylation levels in venous and arterial cell samples as previously described . Interestingly, we found clear evidence for a general hypomethylation in HPVEC (average β = 0.434) relative to HPAEC (average β = 0.479), with the latter comparable to peripheral blood average methylation, whereas HPVEC levels were more similar to that of buccal cells taken at birth. Both HPAEC and HPVEC were generally more methylated on average than placental tissue (Figure 1).
Widespread epigenetic differences in placental arterial and venous endothelial cell compartments
Importantly, where tested, HM450 DNA methylation levels showed strong correlations with those derived using locus-specific SEQUENOM MassARRAY EpiTYPER platform (Additional file 1: Figure S3), validating the HM450 platform, in accordance with other recent studies [26–29]. Unsupervised hierarchical clustering of variable probes in this dataset (SD/mean > 0.4; 151453 probes) clearly discriminated arterial from venous samples (Additional file 1: Figure S1) highlighting the genome-wide DNA methylation differences between these endothelial compartments. HPAEC were more variable as a group than HPVEC in terms of DNA methylation profile. This discriminatory capacity was maintained following unsupervised clustering of the 1000 most highly variable probes (Figure 2). Further, a similar comparison with recently reported methylation datasets generated from cultured adult dermal blood and lymphatic endothelial cell populations  with our dataset, along with unpublished Human Umbilical Vein Endothelial cell (HUVEC) data, highlights the distinct DNA methylation profile of endothelial source from different vascular compartments. HPVEC show the most unique methylation profile, clustering separately from lymphatic and blood derived endothelial cells, HUVECs and HPAEC cells (Additional file 1: Figure S2).
To identify sets of differentially methylated probes (DMPs) between HPAEC and HPVEC, average β-values were calculated for each probe in HPAEC and HPVEC groups. Linear regression analysis revealed a large set of significantly differentially methylated probes (95,266 in total; Benjamini-Hochberg adjusted p-value ≤ 0.05). A total of 15,115 probes, associated with 5,142 genes, showed elevated methylation in HPVEC relative to HPAEC (average β-value difference ≥0.10; Additional file 2: Table S2) whereas 53,725 probes associated with 12,659 genes showed the reciprocal pattern (Additional file 3: Table S3). These data highlight the coordinated nature of methylation change, with blocks of probes (differentially methylated regions – DMRs), showing a common shift in the direction and magnitude of methylation change between cell types.
In order to examine the potential for altered DNA methyltransferase activity to explain the observed global DNA methylation differences described above, we extracted methylation and gene expression data for found evidence for differential methylation of both the DNMT1 and DNMT3A (but not DNMT3B) genes, with elevated average methylation in both promoter regions in HPVEC relative to HPAEC (Additional file 4: Table S1). Further, this was associated with a reciprocal pattern of gene expression specifically for DNMT3A consistent with methylation induced down regulation of DNMT activity in the placental venous compartment (Additional file 4: Table S1).
Of particular interest to endothelial cell biology, many DMP/DMRs were identified in genes previously implicated in endothelial functioning, including nitric oxide synthase 2 (iNOS) [12, 30] and NOS3 (encoding eNOS) [15, 31], von Willebrand Factor (VWF) , NOTCH4, VEGFA, VEGFC, SELE, FLT, and KDR genes, although the number of DMPs and direction of methylation difference between cell types was variable (Additional file 2: Tables S2 and Additional file 3: Table S3).
Interestingly, many genes showed a reciprocal effect on methylation profile according to genomic location relative to the transcription start site. For example, NOS3, showed two regions of highly significant differential methylation. The first, located in the major promoter previously described to drive endothelial expression, encompassed 4 probes with higher mean methylation in HPVEC relative to HPAEC cells, while a reciprocal mean methylation pattern was observed for 5 probes in a downstream exon-associated region in a previously described, placenta-specific variant of this gene (Figure 3). Locus-specific analysis of methylation confirmed highly differential methylation across each of these regions (Figure 3D), however the regulatory function (if any) of this second region remains to be determined. Such complex patterns of methylation within individual genes have recently emerged as a common feature in gene regulation in humans [36–38].
NOS3 is constitutively expressed in the vascular endothelium and plays a critical role in cardiovascular physiology as evidenced by systemic and pulmonary hypertension, abnormal vasculature, defective angiogenesis, poor healing in response to injury and impaired mobilization of stem cells, in eNOS-null mice [reviewed in ]. Disruption of methylation in either of these regions early in development could be a modifiable risk factor for disease in later life.
Complex interplay between endothelial DNA methylation and gene expression
Extensive epigenetic studies have highlighted the complex relationship between DNA methylation and gene expression levels. In addition to the well documented role of elevated methylation at promoter-associated CpG islands in down-regulation of gene expression, it is also apparent that some methylation marks at these (and other) sites may be independent of gene expression status, or even predict active expression, as observed in several DMRs of imprinted genes or in the case of gene-body methylation [36–38]. Thus, in order to examine the relationship between DNA methylation and gene expression, we integrated the methylation and gene expression data sets. For each HM450 probe, differences in β-values between the two groups (Δβ) were calculated and these were compared with corresponding gene expression differences (Additional file 1: Figure S5). By setting cutoffs in average methylation difference higher than 10% (Δβ > 0.1) and log2 fold change in gene expression of ~1.4 (0.5 logFC), we identified a subset of genes that showed a reciprocal relationship between methylation and expression level in HPVEC and HPAEC (Figure 4; Additional file 1: red dots in Figure S5), consistent with a role of increasing methylation in down-regulation of gene expression. This includes 866 HM450 probes showing lower methylation in association with elevated expression of 306 genes in arterial cells (Additional file 5: Table S6), and 2388 probes covering 513 genes showing the same relationship in venous cells (Additional file 5: Table S7).
In order to test for any pattern regarding genomic location and likelihood of regulating gene expression, we examined the proportion of each probe type showing coordinate changes in differential methylation and gene expression, according to annotation (promoter-associated, location within 1500 bp of transcription start site, gene-body, and enhancer-associated). This revealed that the majority of probes were located in enhancer or upstream regulatory (TSS1500, promoter) regions (Table 1, Additional file 1: Figure S6).
Interestingly, endothelial genes likely to be directly regulated by DNA methylation in our dataset included NOS3, vWF, APOLD1, ANGPTL2 genes. Many endothelial genes show both differential methylation at upstream regulatory regions and gene expression differences, but not always in the anticipated reciprocal direction (Table 2A, Additional file 5: Table S8; Additional file 1: Figure S5).
The FLT1 and KDR genes code for two different subtypes of VEGF receptors. These genes have previously been shown to be expressed in most endothelial cells including aortic, vein, microvessel and in many tumors, albeit variably across tumors and individuals [40–42]. FLT1 and KDR are postulated to play a significant role in angiogenesis and/or tumor progression by influencing the activity of VEGF[43, 44]. These genes are shown to be hypermethylated in some cancer cell lines (e.g. colon, stomach) . In our study, FLT1 showed differential promoter methylation and gene expression in HPAEC and HPVEC, however expression and methylation were positively correlated rather than anti-correlated (Table 2). In contrast, the KDR gene was differentially expressed between two cell types but this was independent of any DNA methylation differences in the vicinity of the transcription start site. These findings are interesting given that previous studies have shown down regulation of each of these genes in association with promoter methylation in various cell lines, while cell lines expressing KDR (including HUVECs) appear universally unmethylated [41, 45].
The TLR2 gene has previously been shown to be expressed in endothelial cells in a cell-type specific manner . We found no significant gene expression difference in TLR2 in HPVEC relative to HPAEC (average log FC < 0.073) despite significant differences in promoter methylation levels, with HPVEC having higher regional methylation (average Δβ = 0.213), spanning seven HM450 probes, located proximal to the transcription start site. The activity of this gene in HUVECs has previously been shown to be regulated by inflammation in the absence of altered DNA methylation in HUVECs .
In combination these data highlight the differential and complex role of epigenetics in the regulation of gene expression in endothelial cells of different origin, with several genes likely to be regulated by alternative epigenetic mechanisms (such as histone modification and/or non-coding RNAs) in a cell context-dependent manner. Differential utilization of genomic regulatory elements or expression of upstream transcription factors may also contribute to the observed pattern of expression.
Epigenetic regulation of miR expression in endothelial cells?
Multiple examples of miRNA gene differential methylation were apparent in HPAEC and HPVEC cells (Additional file 5: Table S9), including miRs previously implicated in regulation of endothelial function and/or angiogenesis such as miR-126 and miR-130A . Interestingly miR-200 family members (miR200-A, –B and –C), involved in regulation of VEGF signaling in endothelial cells , were found to be more highly methylated in HPAEC than HPVEC. Furthermore, miR-10A-associated probes, previously linked with arthero-function, were specifically hypermethylated in HPAEC relative to HPVEC . In contrast, the gene encoding miR-125B, a translational suppressor of VE-cadherin , were hypermethylated in HPVEC, relative to HPAEC. Despite the clear evidence of differential methylation, an examination of miR expression profiles failed to reveal a corresponding difference in gene expression between the two cell types (data not shown). Thus, the biological significance of the methylation differences remains.
Ontology of genes associated with DMRs
In order to more fully understand any biological/cellular functions subject to coordinated regulation in endothelial cells, we performed gene ontology and Ingenuity Pathway Analyses (http://www.ingenuity.com) on differentially methylated genes. Prior to the analysis, the list of DMPs was selectively culled to focus on those probes predicted to play a role in gene regulation according to genomic location (associated with gene promoters, enhancers, 5’ UTR, or TSS associated regions as specified by the HM450 manifest annotation file HM450_V1.2). Ingenuity Pathway Analysis revealed an enrichment of genes involved in biological functions such as “Tissue Development”, “Cardiovascular System Development and Function”, and “Tissue Morphology” as differentially methylated in HPAEC and HPVEC (Additional file 5: Table S4). Many of these genes warrant further investigation in the development of arterial or venous endothelial cell phenotype. Gene Ontology Analysis performed by DAVID (http://david.abcc.ncifcrf.gov/, ) highlighted related biological processes associated with the “Extracellular matrix”, previously shown to be critical for proper vascular development , and “Plasma membrane” and “Vasculature Development” (Additional file 5: Table S5). Importantly, coordinated regulation of genes involved in specific cellular pathways such as “cardiovascular pathway genes” was apparent (Additional file 1: Figure S4).
Placenta-derived arterial and venous endothelial cells differ in both their functional characteristics  and differentiation state . Whereas arterial cells are mature and fully differentiated, their venous counterparts have been regarded as immature, representing a juvenile phenotype with a high degree of plasticity. We speculate that differential phenotype of these cells is largely driven by distinct gene expression changes, many of which are mediated by promoter methylation differences in key genes identified here. At present it is not possible to distinguish which epigenetic changes are associated with the degree of maturity of the cells, and which drive more general distinct differences in venous vs arterial characteristics. This will require further investigation.
The two cell types examined show a distinct difference in global DNA methylation level, with the HPVEC hypomethylated relative to HPAEC, and other endothelial cells from a variety of tissues. The low average methylation seen in the venous compartment is reminiscent of the hypomethylation seen in placental tissue and may be one of the underlying mechanisms by which environmental cues modulate their phenotype to adapt to the microenvironment. Indeed, the HPVECs are more sensitive to changes in their local environment in vitro relative to the HPAEC counterparts .
To the best of our knowledge this is the first study comparing global gene expression and DNA methylation of primary arterial and venous endothelial cells isolated from the same organ. In general our data have revealed a highly coordinated series of DNA methylation events, many of which are directly implicated in regulating underlying gene expression levels, while the functional significance of others is less apparent. Included in our dataset are numerous genes previously not studied in relation to endothelial function, that warrant such an investigation in future studies. The combination of DNA methylation and gene expression profiling of early life endothelial cells represents a powerful approach to identify candidate loci potentially subject to environmentally mediated epigenetic disruption, in association with modified risk of later complex diseases involving cardiovascular dysfunction.
HPAEC and HPVEC: isolation and culturing
Primary HPAEC and HPVEC were isolated from arteries and veins, respectively dissected from placentas after uncomplicated vaginal delivery as described previously . Each of the 9 venous and arterial endothelial cell pairs used here was isolated from the same vascular loop of 9 full term human placentas. Cells were cultured on 1% (v/v) gelatin-coated plates using Endothelial Basal Medium (EBM, Cambrex, Clonetics™, Walkersville, MD) supplemented with the EGM™-MV BulletKit (Clonetics™). They were characterized by internalization of acetylated low-density-lipoprotein and immunohistochemical staining for the endothelial cell marker von Willebrand factor and negative staining for fibroblast-specific antigen and smooth muscle actin . Only HPAEC and HPVEC pairs, i.e. isolated from the same placentas, were used in order to minimize variance.
Nucleic acid purification and QC
Total RNA from HPAEC and HPVEC was isolated with RNeasy mini kit (QIAGEN, Dusseldorf, Germany) and scrutinized for quality on the BioAnalyzer BA2100 (Agilent, Santa Clara, CA, USA) with the RNA 6000 Nano Chip Kit (Agilent, Cat No 5067–1511). The RIN (RNA Integrity Number) of the samples ranged between 8.7 and 10. Genomic DNA from HPAEC and HPVEC was isolated using phenol/chloroform density gradient centrifugation method as described previously .
Genome-scale DNA methylation analysis: data acquisition and processing
A total of 1 μg of DNA isolated from 9 HPAEC and 9 HPVEC cell populations was bisulphite converted using the MethylEasy™ bisulphite modification kit (Human Genetic Signatures, Sydney, Australia), according to the manufacturer’s instructions. Unpublished HUVEC data used for unsupervised hierarchical clustering were obtained from purified cells isolated as previously described . Unpublished HM450 data from buccal cells, peripheral blood and placental tissue, used for a comparative analysis of average β values across over 300,000 genomic loci, was kindly provided by Drs Jeff Craig and Boris Novakovic, Murdoch Childrens Research Institute. Conversion efficiency was assessed by bisulphite-specific PCR. Hybridization of bisulphite-treated samples to Illumina Infinium Human Methylation450 (HM450) Beadchips was performed at the Australian Genome Research Facility (AGRF). Raw data files were exported from Genome Studio (Illumina, San Diego, CA) into the R statistical environment (http://cran.r-project.org/index.html). Infinium HM450 data was normalised using the SWAN method from the minfi package available from Bioconductor [56, 57]. This has been specifically designed for such platform where a bias from the two types of probes is apparent. M-values were calculated after removing probes on the sex chromosomes to eliminate any potential gender bias and any poor performing probes, defined as those with a detection p-value cut-off > 0.05 in any sample. β-values were derived from intensities as defined by the ratio of methylated to unmethylated probes given by β =M / (U+M) and were used as a measure of effect size.
Gene expression array analysis: data acquisition and processing
Total RNA was labeled using Affymetrix GeneChip® Whole Transcript (WT) Sense Target Labeling Kit (Affymetrix, Santa Clara, CA, USA; Cat No. 900652) and then prepared for hybridization. For expression analysis RNA was hybridized against GeneChip® Human 1.0 ST arrays (Affymetrix, Cat No. 901087) according to the manufacturer’s instructions. Labeling and hybridization controls were evaluated with Expression Console EC 1.1. Hybridizations and analysis were carried out at the Division Core Facility for Molecular Biology at the Centre of Medical Research at the Medical University of Graz. Microarray data were analysed with Partek Genomic Suite v6.4 software (Partek Inc, St Louis, MO, USA). The import process of the CEL files contained RMA normalization (robust multi-chip average) including background correction, quantile normalization across all arrays and median polished summarization based on log transformed expression values. Significant different genes were extracted with FDR5% / p<0.005 / p<0.05 using ANOVA. Annotations were obtained from NetAffx (Affymetrix).
Statistical analysis and bioinformatics
The Benjamini–Hotchberg method was used for controlling the false discovery rate and correct for multiple testing when comparing HPAEC and HPVEC methylation . The HM450 data (M-value) underwent unsupervised hierarchical clustering analysis using the lumi package . Linear regression analysis was performed using the limma package . For combined gene expression and DNA methylation analysis, Δβ values (average β-values of HPVEC subtracted from average values of HPAEC) were plotted against the average Log2 fold gene expression change of HM450 gene-associated probes. To determine the association between genomic location and differential methylation status, HM450 probes were exclusively assigned to one of the following groups based on the HM450 manifest annotation file version 1.2: promoter associated; enhancer associated; within 1500 bp upstream of a transcription start site (TSS); within 1500 bp upstream of a gene, 5’UTR and first exon; and probes located in gene body and 3’UTR only. Gene Ontology and pathway analysis was performed using the DAVID bioinformatics tool (david.abcc.ncifcrf.gov/)  and Ingenuity Pathway Analysis (http://www.ingenuity.com) under the default settings.
Gene specific DNA methylation
NOS3 locus-specific methylation was performed using the Sequenom EpiTYPER MassARRAY platform (Sequenom, San Diego, USA) as previously described [62, 63]. Amplicons were designed using EpiDesigner software (http://www.epidesigner.com/) and amplification conditions were as follows: 95°C for 10 min; 95°C for 10 s, 56°C for 30 s, and 72°C for 1 min 30 s for 40 cycles; 72°C for 7 min. Primer and target sequences, along with amplification cleavage product patterns and analyzable CpG units are provided in Additional file 1: Figure S7.
Human placental tissue was obtained at term of gestation from uncomplicated pregnancies. All women were lean (BMI 20–24.9), non-smokers and with blood pressure in the normal range. All underwent an oral Glucose Tolerance Test between weeks 24 and 28 and their blood pressure was measured at each visit. Values were in the normal range (ie. did not exceed threshold levels to classify the women as being diabetic/gestational diabetic or having pregnancy-induced hypertension/pre-eclampsia). Informed consent was obtained and ethical approval was granted by the ethics committee of the Medical University of Graz. This study meets the principles of the Declaration of Helsinki.
- β (values):
Beta (values), relative DNA methylation values
Delta Beta, difference in Beta values
Differentially methylated region
Differentially methylated probe
Human placental arterial endothelial cells
Human placental venous endothelial cells
Human umbilical vein endothelial cells
RNA integrity number
Transcription start site.
Barker DJ, Gluckman PD, Godfrey KM, Harding JE, Owens JA, Robinson JS: Fetal nutrition and cardiovascular disease in adult life. Lancet. 1993, 341 (8850): 938-941. 10.1016/0140-6736(93)91224-A.
Harding JE: The nutritional basis of the fetal origins of adult disease. Int J Epidemiol. 2001, 30 (1): 15-23. 10.1093/ije/30.1.15.
Hanson M, Gluckman P: Endothelial dysfunction and cardiovascular disease: the role of predictive adaptive responses. Heart. 2005, 91 (7): 864-866. 10.1136/hrt.2004.047381.
Aitsebaomo J, Portbury AL, Schisler JC, Patterson C: Brothers and sisters: molecular insights into arterial-venous heterogeneity. Circ Res. 2008, 103 (9): 929-939. 10.1161/CIRCRESAHA.108.184937.
Lawson ND, Scheer N, Pham VN, Kim CH, Chitnis AB, Campos-Ortega JA, Weinstein BM: Notch signaling is required for arterial-venous differentiation during embryonic vascular development. Development. 2001, 128 (19): 3675-3683.
Torres-Vazquez J, Kamei M, Weinstein BM: Molecular distinction between arteries and veins. Cell Tissue Res. 2003, 314 (1): 43-59. 10.1007/s00441-003-0771-8.
Swift MR, Weinstein BM: Arterial-venous specification during development. Circ Res. 2009, 104 (5): 576-588. 10.1161/CIRCRESAHA.108.188805.
Kume T: Specification of arterial, venous, and lymphatic endothelial cells during embryonic development. Histol Histopathol. 2010, 25 (5): 637-646.
Lang I, Schweizer A, Hiden U, Ghaffari-Tabrizi N, Hagendorfer G, Bilban M, Pabst MA, Korgun ET, Dohr G, Desoye G: Human fetal placental endothelial cells have a mature arterial and a juvenile venous phenotype with adipogenic and osteogenic differentiation potential. Differentiation. 2008, 76 (10): 1031-1043. 10.1111/j.1432-0436.2008.00302.x.
Ho M, Yang E, Matcuk G, Deng D, Sampas N, Tsalenko A, Tabibiazar R, Zhang Y, Chen M, Talbi S, et al: Identification of endothelial cell genes by combined database mining and microarray analysis. Physiol Genomics. 2003, 13 (3): 249-262.
Fish JE, Marsden PA: Endothelial nitric oxide synthase: insight into cell-specific gene regulation in the vascular endothelium. Cell Mol Life Sci. 2006, 63 (2): 144-162. 10.1007/s00018-005-5421-8.
Matouk CC, Marsden PA: Epigenetic regulation of vascular endothelial gene expression. Circ Res. 2008, 102 (8): 873-887. 10.1161/CIRCRESAHA.107.171025.
Napoli C, Hayashi T, Cacciatore F, Casamassimi A, Casini C, Al-Omran M, Ignarro LJ: Endothelial progenitor cells as therapeutic agents in the microcirculation: an update. Atherosclerosis. 2011, 215 (1): 9-22. 10.1016/j.atherosclerosis.2010.10.039.
Bennett-Baker PE, Wilkowski J, Burke DT: Age-associated activation of epigenetically repressed genes in the mouse. Genetics. 2003, 165 (4): 2055-2062.
Chan Y, Fish JE, D'Abreo C, Lin S, Robb GB, Teichert AM, Karantzoulis-Fegaras F, Keightley A, Steer BM, Marsden PA: The cell-specific expression of endothelial nitric-oxide synthase: a role for DNA methylation. J Biol Chem. 2004, 279 (33): 35087-35100. 10.1074/jbc.M405063200.
Yan MS, Matouk CC, Marsden PA: Epigenetics of the vascular endothelium. J Appl Physiol. 2010, 109 (3): 916-926. 10.1152/japplphysiol.00131.2010.
Peng Y, Jahroudi N: The NFY transcription factor inhibits von Willebrand factor promoter activation in non-endothelial cells through recruitment of histone deacetylases. J Biol Chem. 2003, 278 (10): 8385-8394. 10.1074/jbc.M213156200.
Wu J, Iwata F, Grass JA, Osborne CS, Elnitski L, Fraser P, Ohneda O, Yamamoto M, Bresnick EH: Molecular determinants of NOTCH4 transcription in vascular endothelium. Mol Cell Biol. 2005, 25 (4): 1458-1474. 10.1128/MCB.25.4.1458-1474.2005.
Johnson AB, Barton MC: Hypoxia-induced and stress-specific changes in chromatin structure and function. Mutat Res. 2007, 618 (1–2): 149-162.
Fish JE, Yan MS, Matouk CC, St Bernard R, Ho JJ, Gavryushova A, Srivastava D, Marsden PA: Hypoxic repression of endothelial nitric-oxide synthase transcription is coupled with eviction of promoter histones. J Biol Chem. 2010, 285 (2): 810-826. 10.1074/jbc.M109.067868.
Bronneke S, Bruckner B, Peters N, Bosch TC, Stab F, Wenck H, Hagemann S, Winnefeld M: DNA methylation regulates lineage-specifying genes in primary lymphatic and blood endothelial cells. Angiogenesis. 2012, 15 (2): 317-329. 10.1007/s10456-012-9264-2.
Palinski W: Maternal-fetal cholesterol transport in the placenta: good, bad, and target for modulation. Circ Res. 2009, 104 (5): 569-571. 10.1161/CIRCRESAHA.109.194191.
Stefulj J, Panzenboeck U, Becker T, Hirschmugl B, Schweinzer C, Lang I, Marsche G, Sadjak A, Lang U, Desoye G, et al: Human endothelial cells of the placental barrier efficiently deliver cholesterol to the fetal circulation via ABCA1 and ABCG1. Circ Res. 2009, 104 (5): 600-608. 10.1161/CIRCRESAHA.108.185066.
Scholler M, Wadsack C, Metso J, Chirackal Manavalan AP, Sreckovic I, Schweinzer C, Hiden U, Jauhiainen M, Desoye G, Panzenboeck U: Phospholipid transfer protein is differentially expressed in human arterial and venous placental endothelial cells and enhances cholesterol efflux to fetal HDL. J Clin Endocrinol Metab. 2012, 97 (7): 2466-2474. 10.1210/jc.2011-2969.
Novakovic B, Yuen RK, Gordon L, Penaherrera MS, Sharkey A, Moffett A, Craig JM, Robinson WP, Saffery R: Evidence for widespread changes in promoter methylation profile in human placenta in response to increasing gestational age and environmental/stochastic factors. BMC Genomics. 2011, 12: 529-10.1186/1471-2164-12-529.
Martino D, Maksimovic J, Joo JH, Prescott SL, Saffery R: Genome-scale profiling reveals a subset of genes regulated by DNA methylation that program somatic T-cell phenotypes in humans. Genes Immun. 2012, 13 (5): 388-398. 10.1038/gene.2012.7.
Roessler J, Ammerpohl O, Gutwein J, Hasemeier B, Anwar SL, Kreipe HH, Lehmann U: Quantitative cross-validation and content analysis of the 450k DNA methylation array from Illumina. Inc. BMC Res Notes. 2012, 5 (1): 210-10.1186/1756-0500-5-210.
Bibikova M, Barnes B, Tsan C, Ho V, Klotzle B, Le JM, Delano D, Zhang L, Schroth GP, Gunderson KL, et al: High density DNA methylation array with single CpG site resolution. Genomics. 2011, 98 (4): 288-295. 10.1016/j.ygeno.2011.07.007.
Dedeurwaerder S, Defrance M, Calonne E, Denis H, Sotiriou C, Fuks F: Evaluation of the Infinium Methylation 450K technology. Epigenomics. 2011, 3 (6): 771-784. 10.2217/epi.11.105.
Wilcox JN, Subramanian RR, Sundell CL, Tracey WR, Pollock JS, Harrison DG, Marsden PA: Expression of multiple isoforms of nitric oxide synthase in normal and atherosclerotic vessels. Arterioscler Thromb Vasc Biol. 1997, 17 (11): 2479-2488. 10.1161/01.ATV.17.11.2479.
Fish JE, Matouk CC, Rachlis A, Lin S, Tai SC, D'Abreo C, Marsden PA: The expression of endothelial nitric-oxide synthase is controlled by a cell-specific histone code. J Biol Chem. 2005, 280 (26): 24824-24838. 10.1074/jbc.M502115200.
Gale NW, Yancopoulos GD: Growth factors acting via endothelial cell-specific receptor tyrosine kinases: VEGFs, angiopoietins, and ephrins in vascular development. Genes Dev. 1999, 13 (9): 1055-1066. 10.1101/gad.13.9.1055.
Edelstein LC, Pan A, Collins T: Chromatin modification and the endothelial-specific activation of the E-selectin gene. J Biol Chem. 2005, 280 (12): 11192-11202. 10.1074/jbc.M412997200.
Morishita K, Johnson DE, Williams LT: A novel promoter for vascular endothelial growth factor receptor (flt-1) that confers endothelial-specific gene expression. J Biol Chem. 1995, 270 (46): 27948-27953. 10.1074/jbc.270.46.27948.
Eichmann A, Yuan L, Moyon D, Lenoble F, Pardanaud L, Breant C: Vascular development: from precursor cells to branched arterial and venous networks. Int J Dev Biol. 2005, 49 (2–3): 259-267.
Jjingo D, Conley AB, Yi SV, Lunyak VV, Jordan IK: On the presence and role of human gene-body DNA methylation. Oncotarget. 2012, 3 (4): 462-474.
Aran D, Toperoff G, Rosenberg M, Hellman A: Replication timing-related and gene body-specific methylation of active human genes. Hum Mol Genet. 2011, 20 (4): 670-680. 10.1093/hmg/ddq513.
Ball MP, Li JB, Gao Y, Lee JH, LeProust EM, Park IH, Xie B, Daley GQ, Church GM: Targeted and genome-scale strategies reveal gene-body methylation signatures in human cells. Nat Biotechnol. 2009, 27 (4): 361-368. 10.1038/nbt.1533.
Tsutsui M, Shimokawa H, Otsuji Y, Yanagihara N: Pathophysiological relevance of NO signaling in the cardiovascular system: novel insight from mice lacking all NO synthases. Pharmacol Ther. 2010, 128 (3): 499-508. 10.1016/j.pharmthera.2010.08.010.
Masood R, Cai J, Zheng T, Smith DL, Hinton DR, Gill PS: Vascular endothelial growth factor (VEGF) is an autocrine growth factor for VEGF receptor-positive human tumors. Blood. 2001, 98 (6): 1904-1913. 10.1182/blood.V98.6.1904.
Kim JY, Hwang JH, Zhou W, Shin J, Noh SM, Song IS, Kim JY, Lee SH, Kim J: The expression of VEGF receptor genes is concurrently influenced by epigenetic gene silencing of the genes and VEGF activation. Epigenetics. 2009, 4 (5): 313-321. 10.4161/epi.4.5.9160.
Shibuya M: Structure and dual function of vascular endothelial growth factor receptor-1 (Flt-1). Int J Biochem Cell Biol. 2001, 33 (4): 409-420. 10.1016/S1357-2725(01)00026-7.
Ferrer FA, Miller LJ, Lindquist R, Kowalczyk P, Laudone VP, Albertsen PC, Kreutzer DL: Expression of vascular endothelial growth factor receptors in human prostate cancer. Urology. 1999, 54 (3): 567-572. 10.1016/S0090-4295(99)00156-9.
Veikkola T, Karkkainen M, Claesson-Welsh L, Alitalo K: Regulation of angiogenesis via vascular endothelial growth factor receptors. Cancer Res. 2000, 60 (2): 203-212.
Quentmeier H, Eberth S, Romani J, Weich HA, Zaborski M, Drexler HG: DNA methylation regulates expression of VEGF-R2 (KDR) and VEGF-R3 (FLT4). BMC Cancer. 2012, 12: 19-10.1186/1471-2407-12-19.
Fitzner N, Clauberg S, Essmann F, Liebmann J, Kolb-Bachofen V: Human skin endothelial cells can express all 10 TLR genes and respond to respective ligands. Clinical and vaccine immunology: CVI. 2008, 15 (1): 138-146. 10.1128/CVI.00257-07.
Diesel B, Ripoche N, Risch RT, Tierling S, Walter J, Kiemer AK: Inflammation-induced up-regulation of TLR2 expression in human endothelial cells is independent of differential methylation in the TLR2 promoter CpG island. Innate Immun. 2012, 18 (1): 112-123. 10.1177/1753425910394888.
Wu F, Yang Z, Li G: Role of specific microRNAs for endothelial function and angiogenesis. Biochem Biophys Res Commun. 2009, 386 (4): 549-553. 10.1016/j.bbrc.2009.06.075.
Choi YC, Yoon S, Jeong Y, Yoon J, Baek K: Regulation of vascular endothelial growth factor signaling by miR-200b. Mol Cells. 2011, 32 (1): 77-82. 10.1007/s10059-011-1042-2.
Fang Y, Shi C, Manduchi E, Civelek M, Davies PF: MicroRNA-10a regulation of proinflammatory phenotype in athero-susceptible endothelium in vivo and in vitro. Proc Natl Acad Sci USA. 2010, 107 (30): 13450-13455. 10.1073/pnas.1002120107.
Muramatsu F, Kidoya H, Naito H, Sakimoto S, Takakura N: MicroRNA-125b inhibits tube formation of blood vessels through translational suppression of VE-cadherin. Oncogene. 2012, 32 (4): 414-421.
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.
Lassance L, Miedl H, Konya V, Heinemann A, Ebner B, Hackl H, Desoye G, Hiden U: Differential response of arterial and venous endothelial cells to extracellular matrix is modulated by oxygen. Histochem Cell Biol. 2012, Epub
Novakovic B, Wong NC, Sibson M, Ng HK, Morley R, Manuelpillai U, Down T, Rakyan VK, Beck S, Hiendleder S, et al: DNA methylation-mediated down-regulation of DNA methyltransferase-1 (DNMT1) is coincident with, but not essential for, global hypomethylation in human placenta. J Biol Chem. 2010, 285 (13): 9583-9593. 10.1074/jbc.M109.064956.
Gordon L, Joo JE, Powell JE, Ollikainen M, Novakovic B, Li X, Andronikos R, Cruickshank MN, Conneely KN, Smith AK, et al: Neonatal DNA methylation profile in human twins is specified by a complex interplay between intrauterine environmental and genetic factors, subject to tissue-specific influence. Genome Res. 2012, 22 (8): 1395-1406. 10.1101/gr.136598.111.
Hansen KD, Aryee M: Minfi: Analyze Illumina's 450k methylation arrays. R package version 100.
Maksimovic J, Gordon L, Oshlack A: SWAN: Subset quantile Within-Array Normalization for Illumina Infinium Human Methylation450 BeadChips. Genome Biol. 2012, 13 (6): R44-10.1186/gb-2012-13-6-r44.
Reiner A, Yekutieli D, Benjamini Y: Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics. 2003, 19 (3): 368-375. 10.1093/bioinformatics/btf877.
Du P, Kibbe WA, Lin SM: lumi: a pipeline for processing Illumina microarray. Bioinformatics (Oxford, England). 2008, 24 (13): 1547-1548. 10.1093/bioinformatics/btn224.
Smyth GK: Limma: linear models for microarray data. In: Bioinformatics and Computational Biology Solutions using R and Bioconducto. Edited by: Gentleman R, Carey V, Dudoit S, Huber W, Irizarry R. 2005, New York: Springer, 397-420.
da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.
Ollikainen M, Smith KR, Joo EJ, Ng HK, Andronikos R, Novakovic B, Abdul Aziz NK, Carlin JB, Morley R, Saffery R, et al: DNA methylation analysis of multiple tissues from newborn twins reveals both genetic and intrauterine components to variation in the human neonatal epigenome. Hum Mol Genet. 2010, 19 (21): 4176-4188. 10.1093/hmg/ddq336.
Wong NC, Novakovic B, Weinrich B, Dewi C, Andronikos R, Sibson M, Macrae F, Morley R, Pertile MD, Craig JM, et al: Methylation of the adenomatous polyposis coli (APC) gene in human placenta and hypermethylation in choriocarcinoma cells. Cancer Lett. 2008, 268 (1): 56-62. 10.1016/j.canlet.2008.03.033.
The authors wish to thank Dr Alicia Oshlack for advice on data analysis. This research received funding from the European Union’s Seventh Framework Programme (FP7/2007-2013), project Early Nutrition under grant agreement n °289346 and was funded by a National Health and Medical Research Council (NHMRC): EU application (APP1038018) to RS and the Victorian Government's Operational Infrastructure Support Program. RS is also supported by a NHMRC Senior Research Fellowship (APP1045161).
The authors declare that they have no competing interests.
RS and GD conceived the idea and designed the study. JEJ performed the DNA methylation experiments and the analysis. UH, LL and GD performed the isolation of endothelial cells and the gene expression experiments. LG and DM contributed to the bioinformatics analysis. JEJ and RS wrote the majority of the manuscript. All authors read and contributed to the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1: Unsupervised hierarchical clustering of probes showing variable DNA methylation levels (coefficient of variation >0.4) in HPAEC and HPVEC samples. Figure S2. Unsupervised hierarchical clustering of probes showing variable DNA methylation levels (coefficient of variation >0.4) in endothelial cells derived from different tissue compartments. Figure S3. Correlation between the HM450 and Sequenom EpiTYPER. Infinium HumanMethylation450 methylation accurately reflects DNA methylation levels in HPAEC and HPVEC. Figure S4. Coordinated gene expression and DNA methylation in “Cardiovascular System Development and Function, Connective Tissue Development and Function, Skeletal and Muscular System Development and Function” pathway genes in HPAEC and HPVEC. Figure S5. Scatterplot showing relationship between DNA methylation and gene expression in venous and arterial cells. Figure S6. Proportion of probes associated with specific gene expression change by genomic location. Figure S7. NOS3 Sequenom Assays used to measure regional methylation in HPVEC and HPAEC for Assay 1 (A), and Assay 2 (B). (PPT 939 KB)
Additional file 5: Table S4: IPA (Ingenuity Pathway Analysis) in differentially methylated probes by 10% (DNA methylation) and differentially expressed by 0.5 Log FC (located near transcription start sites). Table S5. Gene Ontology analysis in differentially methylated probes by 10% (DNA methylation) and differentially expressed by 0.5 Log FC (located near transcription start sites). Table S6. List of probes showing coordinated methylation (less methylated) and gene expression (higher expression) in HPAEC. Table S7. List of probes showing coordinated methylation (less methylated) and gene expression (higher expression) in HPVEC. Table S8. List of selected probes associated in genes shown in Table 1. Table S9. List of probes on selected miRs. (XLS 568 KB)
Authors’ original submitted files for images
About this article
Cite this article
Joo, J.E., Hiden, U., Lassance, L. et al. Variable promoter methylation contributes to differential expression of key genes in human placenta-derived venous and arterial endothelial cells. BMC Genomics 14, 475 (2013). https://doi.org/10.1186/1471-2164-14-475