Jump to

Abstract

Objective—Endothelial cells are central to the initiation of atherosclerosis, yet there has been limited success in studying their gene expression in the mouse aorta. To address this, we developed a method for determining the global transcriptional changes that occur in the mouse endothelium in response to atherogenic conditions and applied it to investigate inflammatory stimuli.

Approach and Results—We characterized a method for the isolation of endothelial cell RNA with high purity directly from mouse aortas and adapted this method to allow for the treatment of aortas ex vivo before RNA collection. Expression array analysis was performed on endothelial cell RNA isolated from control and hyperlipidemic prelesion mouse aortas, and 797 differentially expressed genes were identified. We also examined the effect of additional atherogenic conditions on endothelial gene expression, including ex vivo treatment with inflammatory stimuli, acute hyperlipidemia, and age. Of the 14 most highly differentially expressed genes in endothelium from prelesion aortas, 8 were also perturbed significantly by ≥1 atherogenic conditions: 2610019E17Rik, Abca1, H2-Ab1, H2-D1, Pf4, Ppbp, Pvrl2, and Tnnt2.

Conclusions—We demonstrated that RNA can be isolated from mouse aortic endothelial cells after in vivo and ex vivo treatments of the murine vessel wall. We applied these methods to identify a group of genes, many of which have not been described previously as having a direct role in atherosclerosis, that were highly regulated by atherogenic stimuli and may play a role in early atherogenesis.

Introduction

The cells that compose a healthy artery, particularly the vascular endothelial cells (ECs) that line the lumen, are critical to the development of atherosclerotic plaques.1,2 The activation of ECs leads to the recruitment of leukocytes, such as monocytes, and the increased permeability of the endothelial layer permits the entry of low-density lipoproteins (LDL) to the intima, where they are taken up by macrophages to form foam cells.3 As atherosclerosis progresses, smooth muscle cells (SMCs) that reside in the media of the artery migrate to the developing lesion, where they proliferate and join the accumulating foam cells and cellular debris in the growing atheroma. Although the general events that occur in the vessel wall during the transition of a healthy artery to an atherogenic vessel have been established, many of the specific molecular changes that take place during disease progression remain unclear. Determining the precise genetic perturbations that occur in the endothelium at the initiation of atherogenesis, before lipid accumulation and leukocyte infiltration, would be vital in understanding the events that underlie atherosclerosis initiation.

Microarray analysis of ECs from the carotid arteries of human atherosclerosis samples provides valuable insight into the gene expression differences present in early and advanced lesions.4 However, these ECs are dissected from vessels that have already developed atherosclerosis, making it difficult to separate the transcriptional profile of the initial response to atherogenic conditions from downstream perturbations. Genome-wide transcriptional studies of in vivo ECs have been performed in healthy swine, comparing atherosusceptible disturbed flow regions of the vessel wall with atheroprotective undisturbed laminar flow regions to identify genetic pathways that differ between these conditions.5,6 These studies have shown that the atherosusceptible endothelium exhibits elevated levels of proinflammatory and endoplasmic reticulum stress genes, whereas the atheroprotective endothelium has elevated levels of antioxidant genes. Similar studies on ECs isolated from swine maintained on a brief hypercholesterolemic diet supported these findings and reported an upregulation of ABCA1 in response to diet across all regions tested.7 Although swine is a useful model of human atherosclerosis and its size allows for successful in vivo endothelial expression studies, there are significant limitations for their use in genetic studies. Although the small size of the mouse aorta does not allow for the spatial resolution achieved in swine, inbred mouse strains allow for genetic manipulation, biological replicates, and more high-throughput approaches.8,9 Furthermore, the atherosclerotic disease pathology of an inbred mouse strain has been studied in detail,10,11 and the healthy and prelesional disease states can be clearly identified, predicted, and studied.

The mouse is the most commonly used animal model of human atherosclerosis,9,12 yet no studies have been published that examine the global transcriptional changes in the endothelium during atherosclerosis. This may be partially because of the difficulty in obtaining pure, differentiated cell cultures of mouse aortic ECs (MAECs), as cell culture conditions require growing cells away from their native environment and may lead to rapid dedifferentiation.13 The study of cells isolated directly from the aorta as opposed to cell culture offers significant advantages: the data gathered would closely represent the in vivo physical state of these cells because the cell–cell interactions remain intact up until the moment of cell isolation, and it would be possible to assess the transcriptional response of the cells to in vivo conditions, such as hyperlipidemia. However, commonly used cell-isolation techniques, such as fluorescence-activated cell sorting and laser capture microdissection, cannot be used to successfully isolate MAECs in a high-throughput manner because of the tendency of ECs to adhere to one another, the irregular shape of this cell type, and the curvature of the mouse aorta.14 Previous studies on preparations described as mouse aortic intimal cells have explored differential gene expression in disturbed compared with laminar flow regions.15,16 Studies have also examined the effect of hyperlipidemia on the expression of previously identified candidate genes in MAECs,17–20 but a whole-genome hypothesis-generating study has not yet been reported for this cell type.

We characterized cells obtained by a method to isolate the intimal layer from disease-free mouse aortas15 as ECs, and adapted this method to allow for the treatment of these cells with proinflammatory agents ex vivo. These methods were then used to compare the gene expression profiles of MAECs isolated from aortas predisposed to atherosclerosis as compared with those from healthy vessels to determine the major transcriptionally perturbed genes. We used the inbred strain C57BL/6, or BL6, fed a chow diet for these studies; as BL6 is both the most commonly used inbred strain and one of the most atherosusceptible,9,12 it is valuable to analyze the genetic perturbations that occur in the endothelium during atherosclerosis initiation of this widely studied and useful strain. Wild-type BL6 mice fed a normal chow diet do not develop atherosclerosis, whereas BL6 mice with an apolipoprotein-deficient (ApoE−/−) mutation on a chow diet are hyperlipidemic and do develop atherosclerosis. By comparing the gene expression profiles between the healthy wild-type and atherosclerotic ApoE−/− mice, we were able to assess the transcriptional differences that occur during the earliest stages of atherosclerosis initiation. We further characterized the top differentially expressed genes in the endothelium of the prelesion aorta by observing their transcriptional activity in the following atherogenic conditions: proinflammatory stimuli, age, and acute hyperlipidemia. In addition, as the roles of inflammatory mediators in atherosclerosis initiation remain unclear, we applied our methods to provide a global picture of the genetic perturbations caused by 3 inflammatory agents: 1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphocholine (oxPAPC), oxidized LDL (oxLDL), and the more general inflammatory molecule lipopolysaccharide (LPS). We present a comprehensive study of genes that are differentially regulated between healthy and atherogenic mice in discrete MAEC cell preparations and provide valuable insight into the genetics of atherosclerosis initiation.

Materials and Methods

Materials and Methods are available in the online-only Supplement.

Results

Transcript levels of the EC-specific markers von Willebrand factor (vWF), platelet/endothelial cell adhesion molecule 1 (CD31), cadherin 5 (VE Cad), and endothelial nitric oxide synthase (eNOS) were enriched significantly in healthy aortic intimal cell preparations from wild-type BL6 mice compared with RNA extracted from whole mouse aortas (Figure 1A; Table I in the online-only Data Supplement), indicating that the intimal preparations are predominantly ECs. Expression levels of macrophage, leukocyte, and adipose markers were below detection thresholds, indicating the lack of these cell types in our preparations. Previous studies have reported ≈215 CD68+ and 8 CD3+ cells in the entire ascending aorta of the BL6 mouse strain at the ages of 3 to 6 months and an EC density of ≈2500 ECs per millimeter squared in the murine aorta.15 As the aortic arch has a surface area ≥9 mm,2 leukocytes would represent <1% of >22 500 cells in this region; therefore, we conclude that our preparations are >99% pure ECs. Studies in swine report 0.72% presence of leukocytes in preparations taken from the porcine aortic arch endothelium using a scraping method for cell isolation, further indicating that the presence of leukocytes adherent to the EC layer is <1%.6

Cells isolated from the intima of healthy mouse aortas are endothelial cells (ECs). A, Graph represents relative expression of RNA normalized to the housekeeping gene B2M. RNA was isolated from mouse aortic endothelial cell (MAEC) preparations and whole aortas, from 3 animals each, and expression values are shown for individual mice. Normalized gene expression for the 2 EC markers von Willebrand factor (vWF) and platelet/endothelial cell adhesion molecule 1 (CD31) show significant enrichment in the intimal preparations and are nearly undetectable in the whole aorta samples, whereas the smooth muscle cell (SMC) protein marker ACTA2 is expressed to a much higher level in the whole aorta samples. Error bars show SD for 2 technical replicates. B, Relative transcription levels of 8 different smooth muscle cell protein markers and the EC marker CD31 in MAECs and human aortic endothelial cell (HAEC) cultures as determined by reverse transcription quantitative polymerase chain reaction, normalized to the housekeeping gene B2M. mRNA transcript levels for all SMC protein markers are present in the HAEC cultures, which have been confirmed as SMC free using protein markers and serve as a positive control for our purposes. SMC protein transcript levels for the mouse homologs of these genes show similar expression levels in our intimal preparations, indicating a confidence level of EC purity equal to that for HAEC cultures. Error bars show SE for 6 biological replicates for each sample set. t test: **P <0.01; and ***P <0.005.

RNA expression levels of the SMC protein marker, Acta2, were significantly lower in the intimal preparations compared with whole aortas (Figure 1A). Because the amount of cells obtained does not yield sufficient amounts of protein for Western blots, we were unable to assess directly protein expression of the cell preparations. However, we were able to determine indirectly the protein purity levels through comparison with human aortic EC (HAEC) cultures (Figure 1B). These HAEC cultures have been characterized previously as at ≥99% pure as judged by endothelial protein markers and do not show expression of smooth muscle proteins.21 We compared the mRNA expression levels of the EC marker CD31 and 8 different SMC protein markers in our MAEC preparations with those present in HAECs. We found that CD31 levels were equivalent between mouse and human samples, and all SMC protein markers tested are expressed as mRNA transcripts in HAEC cultures as well. Thus, ECs express several transcripts characteristic of SMCs, but these are not translated into proteins, a finding that has been reported in previous studies on ECs.22 These data indicate that cells isolated from the intimal layer of healthy mouse aortas are highly enriched for ECs and can be used as such with the same level of confidence as HAEC cultures.

The RNA obtained consistently represents cells from both high and low areas of shear stress in the aortic arch. We chose to include the entire aortic arch region primarily for ease of reproducibility across the mice used for these studies and to collect enough RNA of high quality. Although separating different spatial regions on a larger animal such as swine is feasible, this is more challenging in a smaller animal such as the mouse, and day-to-day variability could result in identification of false-positive differential expression based on aortic region and not the environmental effects studied.

Amplification of RNA Isolated From MAECs Does Not Affect Relative Abundance of Transcripts

As the MAEC isolation process yielded ≈3 ng of RNA per mouse, we wanted to identify a RNA amplification procedure that would not alter the relative expression levels of transcript in the samples. Although the manufacturer and additional independent research laboratories have demonstrated that the single primer isothermal amplification method used by NuGEN (San Carlos, CA) yields expression data that are concordant with native RNA expression,23,24 we confirmed this in our specific data set using individual genes relevant to our studies. We tested 4 samples: 2 that had been treated with oxPAPC before isolation and 2 control samples incubated in media only. After isolation, half the RNA was converted to cDNA and the other half was amplified using the NuGEN RNA amplification kit. Both the cDNA and the amplified RNA product were assayed for transcript levels of vWF, heme oxygenase 1 (Hmox1), vascular cell adhesion molecule 1 (Vcam1). Both the raw and the normalized reverse transcription quantitative polymerase chain reaction (PCR) values of the native cDNA and amplified RNA products correspond closely, indicating that the expression data acquired from amplified RNA will provide relative expression levels consistent with the native RNA sample (Figure II in the online-only Data Supplement).

MAECs From Prelesion Aortas Show a Differential Expression Pattern Compared With MAECs From Healthy Aortas

To understand the changes in EC gene expression in the initiation and later stages of atherosclerosis, we used 4-week-old BL6 wild-type mice as on a chow diet as healthy controls and 4-week-old hyperlipidemic BL6 ApoE−/− mice on a chow diet as prelesioned mice. BL6 ApoE−/− mice at 24 weeks on a chow diet develop lesions, whereas BL6 wild-type mice at 24 weeks on a chow diet do not develop lesions (Figure III in the online-only Data Supplement).25 At 4 weeks, BL6 ApoE−/− mice on a chow diet have not yet developed lipid deposits observable by oil-red-O staining (Figure IIIA in the online-only Data Supplement), which is why we term these ApoE−/− mice prelesioned aortas as opposed to the healthy BL6 wild-type mice that will not develop lesions under these experimental conditions. High-resolution, freeze-etch electron microscopy indicates that monocytes have not yet adhered to the endothelium in a 4-week-old BL6 ApoE−/− chow-fed mouse because these cells are not observable until 5 weeks.26 We did not detect macrophage transcripts in unamplified RNA collected from wild-type mice (Table I in the online-only Data Supplement) and do not see a difference in the transcription levels of the macrophage markers CD68 and Msr1, which are barely detectable even in amplified RNA, in BL6 wild-type compared with BL6 ApoE−/− 4-week-old mice (Figure IV in the online-only Data Supplement). We wanted to study the intima at a time point before lipid accumulation and monocyte adhesion and determined the 4-week time point to be ideal based on our observations and previous reports.

MAEC RNA was isolated, amplified, and subjected to gene expression microarray analysis. A total of 797 genes were differentially expressed between ECs isolated from the 4-week-old healthy and prelesion aortas at a false-discovery rate of 10% (Figure 2; Table 1), and 32 genes were differentially expressed at the stringent standards of 2-fold change at 5% false-discovery rate (Table 1). Functional annotation analysis of the 797 differentially expressed genes (Table II in the online-only Data Supplement) identified the plasma membrane and immune response–related categories as the most enriched.

To study the response of MAECs to proatherogenic stimuli in the context of the vessel wall, we applied a protocol used for treating HAEC cultures to whole mouse aortas27 (Figure I in the online-only Data Supplement). Neither incubation of the aortas in media before EC isolation affects significantly the expression levels of EC markers (Figure 3A) nor did treatment with the oxidized phospholipid oxPAPC, oxLDL, or LPS (Figure 3B–3D). We next determined whether endothelium isolated from treated vessels displayed gene regulation previously observed in HAEC cell cultures. We found that expression of Hmox-1, a gene induced during oxidative stress and the inflammatory response in HAECs, is upregulated ≈10-fold in response to oxPAPC (Figure 3B), a response also seen in HAEC cultures.28 A similar response was seen after treatment with oxLDL (Figure 3C). LPS treatment (Figure 3D) showed increased expression of proinflammatory marker Vcam1 but no induction of Hmox-1, also consistent with prior observations in HAECs.29 Because gene regulation was consistent with past studies, we proceeded with the expression array analysis.

Differential gene expression of mouse aortic endothelial cells (MAECs) isolated from prelesion compared with healthy aortas and the overlap with treated MAECs. Heatmap displaying differential expression of gene expression transcripts in MAECs. The top bar displays differential expression of MAECs isolated from healthy to prelesion aortas with 10% false-discovery rate. The red bars indicate upregulation of expression from healthy to prediseased, whereas green bars indicate downregulation. The next 3 rows show differential expression data of the same 797 differentially expressed prelesion transcripts in MAECs treated with 1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphocholine (oxPAPC), oxidized low-density lipoprotein (oxLDL), or lipopolysaccharide (LPS), respectively. Gene expression cut-off for treatment is 1.2-fold and is maximized at an absolute value of 4 for purposes of graph visualization. Numbers below treatments indicate the number of transcripts that are regulated in the same direction as prelesion samples.

Patterns of Differential Expression Elicited by oxPAPC Exhibits the Highest Similarity With ECs From Prelesion Aortas

Because many of the early EC changes in atherogenesis are thought to be affected by activation of the inflammatory response, we compared the endothelial gene expression profile after treatment with proinflammatory agents to the gene expression of MAECs in the early stages of atherosclerosis as determined in the previous section. An additional set of MAEC expression data was generated by treating aortas for 4 hours with medium alone or media containing oxPAPC, oxLDL, or LPS before cell isolation and then amplifying the RNA and performing microarray analysis. Differentially expressed genes were determined using MAECs collected from aortas incubated in control media as the baseline. Of the 797 differentially expressed transcripts in ECs from healthy to prelesion aortas, nearly half are regulated in the same direction after treatment with oxPAPC, followed closely with treatment of LPS and oxLDL (Figure 2). Although there is overlap among the 3 different proinflammatory treatments, there are also distinct sets of transcripts that differ. There are also smaller sets of transcripts that are differentially expressed in opposite direction in the prelesion samples compared with the treated samples and a subset of transcripts in ECs from prelesion aortas whose regulation is not seen in the treatments tested, represented by the black areas on the heat map presented in Figure 2.

Differential Expression Analysis Identifies New Candidates for a Role in Atherosclerosis

To identify genes that may play a critical role in EC-mediated atherogenesis, we used a more stringent threshold to identify the top differentially expressed gene transcripts from healthy to MAECs from prelesion aortas as determined by fold change in the microarray data. The expression levels of the genes represented by these microarray probe sets were tested by reverse transcription quantitative PCR. Genes that were determined to have significant differential expression in both microarray and reverse transcription quantitative PCR data have a correlation coefficient of 0.91 and are shown in Table 2. The ApoE−/− mutation commonly used to obtain hyperlipidemic mice could affect gene expression separate from the perturbations that result from the proatherogenic environment. To gain additional evidence that the most regulated genes in the endothelium of prelesion aortas resulted from an atherogenic stimulus, we determined their transcriptional response to additional atherogenic stimuli. These experimental comparisons were all performed using BL6 wild-type mice and did not include the ApoE−/− mutation as a variable. In addition to the oxPAPC, oxLDL, and LPS treatments previously described, we examined transcriptional differences of these highly regulated genes in MAECs isolated from aged mice and mice with acute hyperlipidemia. We compared MAEC expression from BL6 wild-type mice at 4 to 24 weeks to represent aging because BL6 wild-type mice on a chow diet >1 year begin to show lipid accumulation in the aortic root, but not in the aortic arch (Figure IIID and IIIE in the online-only Data Supplement). Although these mice are not hyperlipidemic, as atherosclerosis increases with age in humans,30 they represent the aspects of atherosclerosis that may be regulated by advanced age. We also collected MAEC RNA from mice that had been injected intravenously with human LDL twice during a 5-day time period to represent acute hyperlipidemia (Table III in the online-only Data Supplement). Cell preparations from these treatments were compared with MAECs of prelesion aortas from BL6 ApoE−/− mice at 4 weeks. Eight of the top 14 differentially expressed genes in the prelesioned condition were also perturbed by these separate atherogenic treatments: 2610019E17Rik,H2-Ab1, H2-D1, Pf4, Ppbp, Pvrl2, and Tnnt2 (Figure 4; Table 2). None of the top differentially expressed genes were responsive to LPS treatment, the most nonspecific of the 3 proinflammatory substances tested.

Discussion

We characterized fully and adapted a method15 to isolate ECs from mouse aortas that allowed us to profile gene expression in ECs during the early stages of atherosclerosis. This method provides a quick and efficient way to collect MAEC RNA, generally yielding RNA from a mouse in <1 hour after harvesting, and has several advantages instead of cell culture methods. Not only do we avoid the passage and dedifferentiation issues that occur during cell culture procedures but we are able to directly assay the transcriptional state of the cells that closely represents their in vivo status. In these experiments, we were able to obtain the transcriptional state of ECs in prelesion aortas, and through the use of 1 RNA amplification, we were able to perform microarray analysis on MAECs taken from a single mouse. Moreover, we demonstrate that the whole aorta, with its endothelium intact, can be treated in cell culture medium similar to conventional cell culture with the added advantage of not disrupting the natural environment and cell–cell communication. This method does not allow for the isolation of protein for downstream applications such as Western blots, and because the cells die after the hematoxylin staining step, this procedure would not be appropriate for downstream cell-sorting techniques. We also found that the small size of the mouse aorta made it difficult to consistently gather cells from discrete flow regions. A limitation of the resulting data is that cells from both atherosusceptible and atheroresistant laminar flow regions have been collected, resulting in a heterogeneous RNA sample representing both types of cells. As some of the aortic regions from which we have collected RNA will never develop atherosclerosis31 and are, therefore, not representative of prelesion cells, the resulting differential expression profile may contain false-positives. Further experiments that compare ECs that are isolated exclusively from atherosusceptible aortic regions of the mouse would refine these results. However, we envision that our characterization of this method and its potential applications for various biological treatments will be of use to vascular biologists who require MAEC preparations that are highly enriched and provide a reliable representation of in vivo transcriptional activity.

This method was used to assay the transcriptional profile of MAECs from prelesion aortas for the first time. We found that ≈800 genes are perturbed in the in vivo endothelium before lipid accumulation, a significant alteration to the transcriptional profile of healthy ECs. The most highly regulated genes were enriched for immune response–related categories (Table II in the online-only Data Supplement), a finding that is consistent with prior studies in swine,5 and plasma membrane/cell surface categories, perhaps indicative of the endothelial activation and dysfunction that occur during atherosclerosis initiation.32 We also assessed the role of biologically relevant proinflammatory substances in atherogenic transcriptional regulation through ex vivo treatment. We found that the oxidized phospholipid oxPAPC has the most overlap with ECs from prelesion aortas, and that none of the 14 top differentially expressed genes in these ECs are regulated by LPS treatment, whereas some are regulated by oxPAPC and oxLDL. We also observed that there are sets of genes that are both up- and downregulated in endothelium of prelesion aortas that are unaffected by inflammatory mediators. These findings indicate that although oxidized lipids may be more biologically relevant than LPS in the inflammatory response leading to the initiation of atherosclerosis, these inflammatory mediators are not likely to be solely responsible for the transcriptional changes seen in the early hyperlipidemic endothelium.

We identified 14 highly differentially expressed genes in the prelesion endothelium. We were able to further characterize the transcriptional activity of these novel genes in response to relevant proinflammatory agents using our ex vivo treatment method. Eight of the top differentially expressed genes in prelesion MAECs were also differentially expressed in response to additional disease-relevant stimuli and represent top candidates for further study (Figure 4; Table 2). Among these genes are some that have been described previously as playing a role in the atherosclerotic endothelium, particularly Abca1. This gene has been identified as upregulated in the swine endothelium in response to a brief hypercholesterolemic diet,7 a finding that is consistent with its regulation in our ApoE−/− hyperlipidemic and acute hyperlipidemic mice. The protein encoded by Abca1 has been shown to modulate the transport of phospholipids and cholesterol to apo AI through the endothelial layer and is atheroprotective in vascular ECs because of this role in cholesterol efflux.33,34 The upregulation of this gene that we observe in the hyperlipidemic prelesion intima and the acute hyperlipidemic mice may reflect the response to the influx of LDL in these mice. Pvrl2, a gene expressed by ECs and involved with the immune response,35 shows significant differential expression as well. Its closely related cell surface ligand, Pvr, has been shown to regulate the movement of leukocytes across the endothelium36,37 and interact with vascular endothelial growth factor 2 (VEGFR2),38 suggesting that Pvrl2 may also be involved with endothelial transmigration and VEGF-induced angiogenesis. This gene may also affect EC permeability because it has been shown to control cell–cell adhesion.39 The strong upregulation we observe in both the hyperlipidemic prelesion and the acute hyperlipidemic endothelium, a similar pattern to Abca1, suggests that this gene is strongly responsive to LDL levels.

Our study also identified several genes that have not been linked previously to the endothelial role in atherogenesis, such as the sarcomere protein gene Tnnt2.40 The Tnnt2 transcript, which is homologous to human TNNT2, encodes the protein cardiac muscle troponin T, also known as cTnT. Circulating cardiac troponin T is the most widely used biomarker for the detection of myocardial injury,41 and the extent of coronary atherosclerosis is associated with increasing circulating levels of high-sensitive cardiac troponin T.42,43 We observe elevation of this transcript in mice susceptible for atherosclerosis before developing the disease, as they advance in age, and after treatment with oxLDL. It is possible that this biomarker or its mRNA transcript could indicate endothelial activation during atherogenesis at a much earlier stage of coronary artery disease than myocardial injury. A pair of related transcripts, the platelet genes Pf4 and Ppbp,44 is downregulated in both the prelesion and the aged mouse endothelium. Pf4 has been shown to inhibit EC proliferation and angiogenesis in vitro and in vivo.45 Because atherosclerosis has high levels of angiogenesis,46 downregulation of an inhibitory gene, such as Pf4, could be a key step in allowing the progression of atherosclerosis to proceed once initiation has occurred. Ppbp may also be involved in angiogenesis47 although less is known about its specific role in the process.

We also observe the perturbation of 4 histocompatibility genes: H2-Aa, H2-Ab1, H2-D1, and H2-Ea. These genes are all located in the same cluster in both mouse and human, and a Genome Wide Association Study (GWAS) meta-analysis identified a single-nucleotide polymorphism within this region that is associated with coronary artery disease48 (Table 3). Transcripts for the mouse homolog of HLA-B, H2-D1, were downregulated in response to treatment with the oxidized phospholipid oxPAPC and in the prelesion endothelium. When transfected into ECs, a pulmonary hypertension risk allele of HLA-B influences the production of endothelin -1 (Edn1) and endothelial nitric oxide synthase (eNOS) and induces both ER stress and the unfolded protein response,49 all processes that are central to endothelial activation during atherosclerosis.50 The human homologs to mouse H2-Aa, H2-Ab1, and H2-Ea, HLA-DR and HLA-DQ, are major histocompatibility complex class II members that regulate T-cell–dependent immune responses. Vascular SMCs in atherosclerotic lesions express HLA-DR,51 and as studies have shown that hypoxia induces both HLA-DR expression and secretion in EC cultures,52 it is likely that MAEC expresses this transcript as part of the inflammatory response in atherogenesis as we have seen here. Although the HLA-DQ transcripts are also expressed in ECs,53 little is known about the HLA-DQ subtypes in the context of atherosclerosis and why they would be downregulated in atherosusceptible conditions as we show here. The previously unstudied gene 2610019E17RIK is downregulated in the prelesion endothelium, by oxLDL treatment, and in older mice. It is a noncoding RNA of 494 base pairs located on mouse chromosome 17 and may merit further characterization as noncoding RNAs have been shown to play central roles in disease.54

In summary, our study provides several valuable leads for the study of the genes mediating atherosclerosis in the vessel wall. Our characterization of a vascular endothelium isolation method and demonstration that it is applicable to treatment of cells in the aorta provide new tools to researchers performing individual gene follow-up studies in the mouse vasculature.55 Future studies could apply these methods to observe differences in discrete flow regions and compare them from atheroresistant and atherosusceptible mouse strains to refine these data and identify additional genetic factors. The application of this method identifies several genes that are differentially expressed in the endothelium during the early stages of atherogenesis. Although most of these genes have been described previously as having roles in inflammation, ECs, or heart disease, this is the first time the majority of these transcripts have been correlated directly to vascular ECs in atherosclerosis and provide exciting leads for further study.

Acknowledgments

We thank Jan Danciger for her advice on endothelial cells, Ladan Vakili for the oxidized low-density lipoprotein, and Xuping Wang and Judy Wu for sectioning and staining the aortic arch samples.

. Troponin measurements during drug development–considerations for monitoring and management of potential cardiotoxicity: an educational collaboration among the Cardiac Safety Research Consortium, the Duke Clinical Research Institute, and the US Food and Drug Administration.Am Heart J. 2011;162:64–73.

Significance

Vascular endothelial cells are central to the initiation and progression of atherosclerosis, but there have been few studies that directly assess the genetic changes that occur in vivo during these events. As the mouse is the most commonly used animal model of human atherosclerosis, isolating and studying the endothelial cells from the mouse aorta would provide valuable information on their role in disease. This study characterizes cells isolated from the intima of the mouse aorta as endothelial cells and applies this method to identify a set of genes that are highly differentially expressed in endothelial cells of mice susceptible to atherosclerosis compared with mice that are resistant. The expression of these genes is also studied in response to several atherogenic stimuli. On the basis of these data, we identify 8 candidate genes for a role in the endothelium during atherosclerosis initiation, many of which have not been described previously in this context.