Next Article in Journal
Deep-Learning Algorithm and Concomitant Biomarker Identification for NSCLC Prediction Using Multi-Omics Data Integration
Next Article in Special Issue
An Efficient Bayesian Method for Estimating the Degree of the Skewness of X Chromosome Inactivation Based on the Mixture of General Pedigrees and Unrelated Females
Previous Article in Journal
Proteome Profile Changes Induced by Heterologous Overexpression of Mycobacterium tuberculosis-Derived Antigens PstS-1 (Rv0934) and Ag85B (Rv1886c) in Mycobacterium microti
Previous Article in Special Issue
Application of the Dynamical Network Biomarker Theory to Raman Spectra
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 12
Issue 12
10.3390/biom12121837
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systems Biology Analysis of Temporal Dynamics That Govern Endothelial Response to Cyclic Stretch

by
Michael W. Lai
1,
Nathan Chow
1,
Antonio Checco
1,
Balvir Kunar
2,
David Redmond
2,
Shahin Rafii
2 and
Sina Y. Rabbany
1,2,*
1
Bioengineering Program, DeMatteis School of Engineering and Applied Science, Hofstra University, New York, NY 11549, USA
2
Division of Regenerative Medicine, Ansary Stem Cell Institute, Weill Cornell Medicine (WCM), New York, NY 10065, USA
*
Author to whom correspondence should be addressed.
Biomolecules 2022, 12(12), 1837; https://0-doi-org.brum.beds.ac.uk/10.3390/biom12121837
Submission received: 21 October 2022 / Revised: 23 November 2022 / Accepted: 28 November 2022 / Published: 8 December 2022
(This article belongs to the Special Issue Systems Biology and Omics Approaches for Complex Human Disease)
Browse Figures
Versions Notes

Abstract

:
Endothelial cells in vivo are subjected to a wide array of mechanical stimuli, such as cyclic stretch. Notably, a 10% stretch is associated with an atheroprotective endothelial phenotype, while a 20% stretch is associated with an atheroprone endothelial phenotype. Here, a systems biology-based approach is used to present a comprehensive overview of the functional responses and molecular regulatory networks that characterize the transition from an atheroprotective to an atheroprone phenotype in response to cyclic stretch. Using primary human umbilical vein endothelial cells (HUVECs), we determined the role of the equibiaxial cyclic stretch in vitro, with changes to the radius of the magnitudes of 10% and 20%, which are representative of physiological and pathological strain, respectively. Following the transcriptome analysis of next-generation sequencing data, we identified four key endothelial responses to pathological cyclic stretch: cell cycle regulation, inflammatory response, fatty acid metabolism, and mTOR signaling, driven by a regulatory network of eight transcription factors. Our study highlights the dynamic regulation of several key stretch-sensitive endothelial functions relevant to the induction of an atheroprone versus an atheroprotective phenotype and lays the foundation for further investigation into the mechanisms governing vascular pathology. This study has significant implications for the development of treatment modalities for vascular disease.

1. Introduction

Vascular endothelial cells (ECs) line the luminal surfaces of blood vessels and experience dynamic mechanical forces from their microenvironment. These stimuli result from the hemodynamic parameters exerted by blood flow and pressure and play a role in regulating the gene expression and physiological responses of ECs—known as mechanotransduction. Two important concurrent hemodynamic forces are shear stress and cyclic stretch. Shear stress is a result of the friction that is generated between blood flow and the endothelium, while cyclic stretch arises due to the pulsatile nature of blood flow [1]. Subsequent morphological changes and atheroprone or atheroprotective EC phenotypes result in a response to these stimuli.
The differentiation of atheroprone versus atheroprotective markers and morphologies is necessary for the investigation of the fine balance that ECs exhibit within the body, where subtle hemodynamic changes over the course of years can precipitate adverse health outcomes if not kept within ideal parameters. A wide range of sequelae from aberrant EC stress remain significant health concerns, such as atherosclerosis, thrombus formation, aneurysm, aortic dissection, and nephropathy. In the context of thrombi, for example, decreased EC vitality leads to diminished barrier function and the exposure of sub-endothelial collagen to platelets, leading to aggregation and vessel occlusion [2].However, the most prominent of these concerns is atherosclerosis, due to its driving role in vascular disease, which despite advances remains the number one cause of death in the US [3].
Both laminar shear stress and physiological cyclic stretch have been demonstrated to be essential in the maintenance of vascular homeostasis and are thus associated with an atheroprotective phenotype [4]. Meanwhile, disturbed blood flow imparts oscillatory shear stress and pathological cyclic stretch to the endothelium, disrupting appropriate EC mechanotransduction, and thus is associated with an atheroprone phenotype [5,6,7]. The mechanisms by which laminar and oscillatory shear stress induce an atheroprotective and atheroprone phenotype, respectively, have been well characterized [8,9,10]. For cyclic stretch, several alterations to the cytoskeleton and downstream cellular signaling have been reported [11]. These alterations occur due to signal transduction by the mechanosensitive receptors associated with ECs, including ion channels known as the transient receptor potential (TRP), which have a demonstrated role in stretch-mediated calcium influx, integrins, and the platelet endothelial cell adhesion molecule (PECAM-1) [12,13,14,15]. Transduction by these mechanoreceptors induces several hallmarks of the endothelial response to cyclic stretch, including actin stress fiber formation regulated by the small GTPase Rho and stretch-mediated angiogenesis [16,17]. The endothelial atheroprone phenotype brought about by pathological cyclic stretch disrupts vascular homeostasis through the increased production of reactive oxygen species (ROS) and increased inflammation [18,19]. However, the molecular events governing the progression from an atheroprotective to an atheroprone phenotype in response to cyclic stretch require further investigation. To better understand how the magnitude of the cyclic stretch affects cellular responses, we examined the level of mechanical stimuli in a quantitative manner to discern the relative contribution of different regulatory genes.
Identifying markers signifying EC vitality is a prerequisite for the analysis of the factors that influence it and lead to vascular disease. EC vitality can be monitored by EC proliferation, stiffness, inflammatory response, mTOR expression, and fatty acid metabolism. Proliferation ensures viable endothelium and barrier function. Stiffness refers to the EC’s ability to absorb hemodynamic stress, with excessive stiffness being a function of poor resilience. The EC modulation of inflammatory responses protects or promotes immune deposition-forming plaques, while fatty acid metabolism influences lipid deposition. Molecular Target of Rapamycin (mTOR) is linked with vasodilation, inflammation, cell survival, and cell proliferation and provides a comprehensive view of EC function. Therefore, each of these variables is useful in the study of EC vitality or the lack thereof [20].
The physiological and pathological responses of ECs to cyclic stretch have been primarily elucidated using experimental models that impose uniaxial cyclic stretch in vitro at a set magnitude [21]. For more physiological relevance, it is necessary for equibiaxial cyclic stretch to be applied to ensure greater strain homogeneity [22,23]. In this study, primary human umbilical vein endothelial cells (HUVECs) were subject to equibiaxial cyclic stretch at magnitudes in the range of 5% ∆r to 25% ∆r. Stretch magnitudes of 10% ∆r and 20% ∆r were utilized for the final RNA-sequencing (RNA- seq) analyses, which are often used in vitro to represent physiological and pathological conditions [24]. HUVECs were chosen for our experimental system given that they have been used previously to model cardiac pathologies due to their in vitro durability [8]. Smooth muscle cells are prone to differentiation when cultured in vitro, in the absence of natural stimuli such as high shear and pressure [25]. Once the HUVECs are cultured, these cells assume the molecular and cellular features of “generic endothelial cells”. We then performed a series of post-sequencing analyses to determine the molecular networks that drive the key endothelial responses to cyclic stretch (cell cycle regulation, inflammatory response, fatty acid metabolism, and mTOR signaling). Here, we utilize a systems biology-based approach to present a detailed mapping of the molecular networks which govern the dynamic transition from an atheroprotective to an atheroprone phenotype in response to cyclic stretch.

2. Materials and Methods

2.1. Endothelial Cell Culture

Primary HUVECs freshly isolated from two different vascular beds were cultured in Medium 199 (Corning, Corning, NY, USA), supplemented with l-glutamine (Corning, Corning, NY, USA), antibiotic antimycotic (Corning, Corning NY, USA), HEPES buffer (Corning, Corning, NY, USA), heparin (Sigma Aldrich, St. Louis, MO, USA), and 20% fetal bovine serum (Hyclone, Logan, UT, USA). The HUVECs were transfected with the E4ORF1 of the Ad E4 gene complex, which has previously been demonstrated to support the long-term survival of ECs through sustained Akt phosphorylation [26]. The cells were maintained within a trigas humidified Heracell 150i incubator (Thermo Scientific, Waltham, MA, USA) at 37 °C and 5% CO2/O2. For experiments to inhibit the mTOR’s kinase activity, the cells were pre-treated with rapamycin (Cell Signaling Technology, Danvers, MA, USA) at a concentration of 100 nM for 24 h prior to stretch.

2.2. Cyclic Stretch

A CellScale MCB1 (Cell Scale Biomaterials Testing, Ontario, Canada) device was adapted to deliver the equibiaxial cyclic stretch to the HUVECs seeded onto custom-built circular polydimethylsiloxane (PDMS) molds. To allow for the application of the equibiaxial stretch, the HUVECs were cultured on circular molds which were mounted within the stretch chamber of the device. The molds were cast using a Sylgard 184 Silicone Elastomer Kit (Dow Chemical, Midland, MI, USA) at a standard 10:1 ratio of base to curing agent. The PDMS molds were then cured for 24 h and sterilized under ultraviolet light. To assist the adhesion of the HUVECs to the PDMS, the molds were plasma treated for one minute using a PE-25 plasma cleaner (Plasma Etch, Carson City, NV, USA) and coated with fibronectin (Thermo Fisher Scientific, Waltham, MA, USA) at a concentration of 20 µg/mL. After 12 h of culture, the equibiaxial cyclic stretch was applied to the cells at physiologically relevant (10% change of radius (∆r)) and pathologically relevant (20%∆r) magnitudes, with a frequency of 0.5 Hz for a duration of 24 h, to recapitulate the mechanical stretches experienced by the cells in vivo. At the end of each experiment, RNA was collected for downstream RNA-Seq and qPCR applications or utilized for imaging by atomic force microscopy. To verify that the experimental system delivered sufficient cyclic stretch, the cells were stained for actin (conjugated phalloidin fluorescein isothiocyanate, Sigma-Aldrich, St. Louis, MO, USA) and DAPI nuclear counterstain following the 24 h stretch period.

2.3. RNA Sequencing

Total RNA was extracted from the HUVECs which had undergone both 10% and 20% equibiaxial stretch using an RNEasy Mini Kit (Qiagen, Hilden, Germany). For RNA-seq analysis, the RNA extracts were submitted to the Weill Cornell Medicine Genomic Core Facility. The RNA quality was assessed through RNA integrity numbers (RINs) using an Agilent Technologies 2100 Bioanalyzer. All the RNA samples used in this study had a RIN of 10 out of 10. The RNA library preps were generated and multiplexed using Illumina’s TruSeq RNA Library Preparation Kit v2 (non-stranded and poly-A selection). Ten nanometers of cDNA was used as input for high-throughput sequencing via Illumina’s HiSeq 4000, producing 51 bp paired-end reads. The sequencing reads were de-multiplexed (bcl2fastq v2.17), checked for quality (FastQC v0.11.5), and trimmed/filtered when appropriate (Trimmomatic v0.36). The resultant high-quality reads were mapped (TopHat2 v2.1.0; Bowtie2 v2.2.6) to the transcriptome sequence reference of the UCSC mm10 genome build. The gene counts were quantified using the Python package HTSeq (v0.11.1). The transcript abundance measures (FPKM values) were quantified using Cufflinks (v2.2.1). The gene counts were then analyzed for differential expression using edgeR [27]. Differentially expressed genes were selected across biological replicates with an adjusted p-value of <0.05. Counts per million (CPM) values were calculated and utilized to filter out lowly expressed genes. The genes with CPM < 1 were excluded. The selection of differentially expressed genes was performed through pairwise comparisons of 20% stretch vs. control (non-stretch), 20% stretch vs. 10% stretch, and 10% stretch vs. control.

2.4. Pathway Analysis

Gene set enrichment analysis (GSEA) was utilized to perform pathway enrichment analysis on the genes that were found to be differentially expressed in each of the pairwise comparisons [28]. Enrichment analysis was performed against the curated gene sets within the Molecular Signatures Database (Broad Institute). Cytoscape version 3.7 was then utilized to map the interaction networks within the enriched gene sets. To create network maps, the Cytoscape plugin GeneMANIA (University of Toronto) was utilized. GeneMANIA employs the guilt-by-association approach to map interactions among related genes [29]. Functional pathway diagrams were then constructed through a combination of the Cytoscape connectivity values and manual literature searches.

2.5. Immunofluorescence and Western Blotting

To functionally validate the molecular changes reported by the RNA-seq and pathway analysis, immunofluorescence and Western blotting were performed. The cells were stained for actin (conjugated phalloidin fluorescein isothiocyanate, Sigma-Aldrich, St. Louis, MO, USA), VE-cadherin (goat IgG mouse VE-cadherin, R&D Systems, Minneapolis, MN, USA), and DAPI nuclear counterstain and visualized using a Zeiss Axio Observer Z1. Western blotting was performed with primary Rabbit anti-VEGFR2 (Cell Signaling Technology #2479), Rabbit anti-Akt (Cell Signaling Technology #9272), and Rabbit anti-GAPDH (Cell Signaling Technology #2118) and secondary Goat anti-Rabbit IgG (Li-Cor, Lincoln, NE, USA).

2.6. Quantitative PCR

The expression analysis of target genes by quantitative PCR (qPCR) was performed as previously described [30]. The primer sequences are listed in Supplementary Table S1.
Relative gene expression was calculated according to the 2−ΔΔCt method and was normalized to the expression of the housekeeping gene GAPDH. Pairwise comparisons of 20% stretch vs. control and 10% stretch vs. control were performed to determine standard deviation and mean gene expression, in which the stretch samples were normalized to the control. A one-way analysis of variance (ANOVA), followed by a post-hoc Fisher’s LSD multiple comparison test, was utilized to determine statistical significance. A p-value < 0.05 was deemed statistically significant.

2.7. Atomic Force Microscopy

An MFP-3D-Bio atomic force microscope (Asylum Research, Santa Barbara, CA, USA) was used for all the atomic force microscopy (AFM) measurements. All the AFM measurements were carried out on live cells with the probe immersed in a cell growth medium kept at a temperature of 37 °C using an AFM Petri dish heater (Asylum Research, Santa Barbara, CA, USA). Young’s moduli of live cell monolayers were obtained using a silicon nitride cantilever (Novascan, Ames, IA) with a spherical borosilicate glass probe of 5 μm in diameter. The deflection sensitivity of the cantilever was measured using glass as a stiff surface, while the cantilever spring constant was determined through analysis of the power spectral density (PSD) of the cantilever thermal fluctuations. The AFM software (Asylum Research, Santa Barbara, CA, USA) was utilized to record the PSD and fit the theoretical PSD of a simple harmonic oscillator to the experimental data. The fit described the data accurately, yielding a cantilever stiffness of k = 0.015 N/m, which was consistent with the nominal stiffness (0.02 N/m) reported by the manufacturer. The cell monolayers were positioned under the AFM tip, and the force versus displacement curves were captured on a 32 × 32 grid spanning 80 × 80 μm. Each force curve was recorded at a rate of 4 μm/s and indented at a maximum distance of 200 nm. Low-quality indentation curves (i.e., curves with undefined contact point or non-monotonous slope) were rejected.
The indentation was sufficiently small compared to the height of the cell, allowing the Hertz–Sneddon (HS) model to be used to fit the indentation curves using the AFM software. According to the HS model, the force (denoted as F) versus indentation (denoted as δ) experienced by a paraboloidal indenter with curvature radius R is F = 4/3[E/(1 − ν2)]R1/2 δ 3/2, where E and ν denote the cell’s Young’s modulus and Poisson’s ratio, respectively. The data were then fitted from the point of contact between the spherical probe and the sample (chosen to be the point at which the derivative of the force distance curve became nonzero) to the maximum indentation value. The best fit yielded the contact point height and the Young’s modulus of the cell, assuming a cell’s Poisson ratio = 0.45, and the nominal indenter curvature radius R = 2.5 μm. The fitted Young’s modulus was utilized to construct the stiffness maps, whereas the contact point height was utilized to map the sample topography. The topography data were used to locate the regions of the stiffness maps where the sample thickness was <1 μm. These regions were occupied by the substrate and the cell periphery and were excluded from the statistical analysis carried out to assess the average stiffness of the cell’s body. The statistical distribution of the AFM measurements was displayed using box plots, in which the horizontal line is representative of the median, the box is representative of the data from the first to third quartile, and the whiskers are representative of the data minimum and maximum. Statistical analyses were conducted using paired t-tests, in which p < 0.05 was indicative of statistical significance.

3. Results

3.1. Cyclic Stretch Validation

To validate that the equibiaxial stretch was delivered to the cellular system, a calibration of the cyclic stretch apparatus was performed by sectioning a PDMS membrane into five evenly distributed locations and by determining the mean displacement induced by stretch in the x and y directions across the five technical replicates (Figure 1A). In each of the five locations, there was a mean displacement of 0.0635 cm in both the x and y directions after stretch, indicating that equibiaxial cyclic stretch was imposed on the system (Figure 1B). To verify that the ECs experienced cyclic stretch within the cellular system, F-actin staining was performed on the ECs prior to and after stretch. In the static ECs, the staining is primarily concentrated on the cortical actin, with a limited appearance of actin stress fibers (Figure 1C). However, Figure 1D demonstrates the clear appearance of actin stress fibers and an elongated morphology in the stretched ECs, which is a hallmark of the actin remodeling which occurs in response to cyclic stretch [31].

3.2. Cell Cycle Regulation

It has previously been reported that there is increased endothelial proliferation under higher magnitudes of cyclic stretch [32]. The GSEA analysis of the 10% sequencing data demonstrated that there is a significant upregulation of the genes associated with cell cycle progression in the 20% stretch relative to the 10% stretch, specifically those downstream of the E2F family of transcription factors (Figure 2A). To validate the proliferative response, Western blot analysis of VEGFR2 and Akt were performed (Figure 2B). The protein expressions of both were upregulated in response to the 20% stretch relative to the 10% stretch and the static conditions.
The pathway analysis was performed to delineate the central regulatory network driving the cell cycle progression under the 20% stretch. To do this, the relevant connections between the target genes and the central network nodes were mapped (Figure 3A). Of the E2F target genes, TUBG1, POLD2, CSE1L, DEK, HMGA1, HMGB2, MCM6/7, PA2G4, RAN, BIRC5, CDC45, CDC7, CKS2, DBF4, EZH2, MCM3/4, and POLD1 maintained connections to CDK2, which was the central node with the greatest network connectivity (Table 1). MCM2 had the second greatest network connectivity, maintaining connections to MCM3/4, MCM6/7, CDC45, CDC7, CKS2, and DBF4. MCM2 is a minichromosome maintenance protein that, along with MCM3-7, makes up the protein complex which forms the core of the DNA helicase during DNA replication [33]. The upregulation of the MCM proteins with greater magnitudes of cyclic stretch may provide a molecular basis for the increased rate of DNA synthesis that underlies the cellular response to pathological stretch [34]. The activation of RPA1, a replication protein, and PCNA, a cofactor for DNA polymerase-∂, in response to cyclic stretch may also play a role in this process, as both genes were connected to target genes associated with the enzymes involved in DNA replication (MCM3/4, MCM6/7, POLD1, POLD2, and LIG1). Although RPA1 had the third greatest number of connections in the cell cycle connectivity network as a whole (Supplementary Figure S1), CDK1 had the third most connections to the E2F target genes, emphasizing the importance of CDKs in E2F pathway regulation. CDK1 was connected to HMGA1, BIRC5, CKS2, EZH2, and MCM4.

3.3. Inflammatory Response

Continuous cyclic stretch at high magnitudes may induce the inflammatory secretion of cytokines [35]. Interferons are antiviral factors that are responsible for immune activation through cytokine release [36]. Of the relevant genes associated with interferon response identified by the GSEA analysis, STAT1, HIF1A, BTG1, and CCL2 were differentially expressed at 20% stretch relative to 10% stretch (Figure 2). The lesser-known CCL2 is a cytokine and likely modulates inflammatory response, while BTG1 is an antiproliferative factor which should be considered a result of the inflammatory response. Each of these genes has been implicated in the interferon-mediated inflammatory response [37,38]. Pathway analysis was conducted to identify the relevant connections between the interferon-associated genes and the overall regulatory network. STAT1, HIF1A, BTG1, and CCL2 each maintained a connection to the transcription factor JUN, while HIF1A also maintained a connection to the transcription factor MYC, and CCL2 maintained a connection to the transcription factor JUNB (Figure 3B). The genes identified as being involved in the cyclic stretch-mediated cholesterol homeostasis in Figure 2, namely SREBF2, a sterol regulatory element binding transcription factor, and ATF3, an activating transcription factor 3, have also been implicated in the endothelial inflammatory response [39,40]. SREBF2 was connected to CREB1, while ATF3 was connected to JUNB, JUN, ATF2, and CEBPG (Table 2).

3.4. Fatty Acid Metabolism

A pathological vascular phenotype has been associated with an increase in glucose metabolism and a decrease in the metabolic utilization of fatty acids [33,41]. The GSEA analysis was utilized to determine the differentially expressed genes associated with the metabolic changes in response to the pathological magnitudes of the cyclic stretch. Of those identified in the gene set, DLD and HSP90AA1 were upregulated, while ACOT8, IDH3G, ALDOA, and S100A10 were downregulated (Figure 2). ATF2 and MYC were the central regulators of these metabolic genes; ATF2 maintained connections to DLD, ALDOA, and S100A10, while MYC maintained connections to ACOT8, IDH3G, and HSP90AA1 (Table 2).

3.5. mTOR Signaling

Vascular mTOR signaling has been implicated in both physiological and pathological processes [42]. Using the GSEA analysis, we delineated the mTOR targets which were specifically involved in the pathological mTOR signaling. Of these, CXCR4, SSR1, NUP205, CCT6A, HSPD1, and PLK1 were upregulated in the 20% stretch, while SERPINH1, NFYC, TPI1, ACLY, and PLOD2 were downregulated (Figure 2). MYC was found to be the primary central regulator of the mTOR target genes, maintaining connections to CXCR4, SSR1, NUP205, CCT6A, HSPD1, PLK1, SERPINH1, and NFYC. ATF2 and TERF1 each had three connections; ATF2 was connected to HSPD1, SERPINH1, and TPI1, while TERF1 was connected to TPI1, ACLY, and PLOD2. JUN and NR3C1 each had one connection, to CXCR4 and HSPD1, respectively (Table 2).

3.6. Regulation of Cell Stiffness by Cyclic Stretch

To determine the relationship between cyclic stretch and cell stiffness, we utilized AFM to measure the stiffness of the ECs subjected to 10% and 20% stretch. We found that the cyclic stretch decreased the EC stiffness in a dose-dependent manner, with a median stiffness decreasing from 606 Pa to 405 Pa in response to the 10% stretch and to 306 Pa in response to 20% stretch (Figure 4A). We then treated ECs with rapamycin, a known inhibitor of the mTOR kinase activity, which has previously been demonstrated to play a role in cytoskeletal reorganization, and found a reversal of the previously observed trend, with median stiffness increasing to 461 Pa in response to the 10% stretch and to 421 Pa in response to the 20% stretch (Figure 4B) [43]. The initial decrease in stiffness from 10% to 20% stretch correlated with an increased modulation of the cytoskeletal component genes in response to the 20% stretch (Figure 4C). Of these genes, we mapped the expression of those downstream of mTOR and found a greater number of mTOR target genes in response to the 20% stretch relative to the 10% stretch (Figure 4D). We then constructed a connectivity network consisting of mTOR and the cytoskeletal component genes (Figure 5A). AURKA and CDK1 were identified as central regulators, both of which were connected to mTOR.

3.7. Cytoskeletal Remodeling

The ECs showed increased localization of VE-cadherin (red) to cell junctions when stretched (Figure 5C) compared to the unstretched controls (Figure 5B). In the unstretched controls, VE-cadherin was preferentially located in the cytoplasm. Differential cell morphology, as highlighted by the actin stain (green), was observed compared with the controls, and it also varied as a function of the distance from the center of the substrate, with greater cell polarity observed further away from center.

3.8. Shear Stress Comparison

Shear stress and cyclic stretch are complementary hemodynamic forces that act on the endothelium. To determine whether these forces impart similar downstream expressional changes, we compared the expression of the cyclic stretch effector genes identified in Figure 2 in response to the 20% cyclic stretch, the 10% cyclic stretch, the oscillatory shear stress, and the pulsatile shear stress (Figure 6A). Previously published oscillatory and pulsatile shear stress data (GEO Accession Number GSE103672) were reanalyzed for the purposes of this study [8]. These data were obtained from HUVECs exposed to a shear stress of 12 ± 5 dyn/cm2 (pulsatile shear stress) or 0.5 ± 5 dyn/cm2 (oscillatory shear stress). Among the mapped effector genes, the 10% cyclic stretch expression clustered with the pulsatile shear expression, while the 20% cyclic stretch expression clustered with the oscillatory shear expression, suggesting that the 10% cyclic stretch imposes downstream expressional changes similar to those of the pulsatile shear stress, while the 20% cyclic stretch imposes downstream expressional changes to those of the oscillatory shear stress. We then conducted pairwise comparisons to determine whether there are shared transcriptional regulators between each clustered pair (Figure 6B). Seventeen transcription factors were commonly upregulated, and one transcription factor was commonly downregulated between the 20% cyclic stretch and the oscillatory shear stress, while eleven transcription factors were commonly upregulated and downregulated between the 10% cyclic stretch and the pulsatile shear stress.

3.9. qPCR Expression Validation

The expression of four regulatory target genes, E2F1, STAT1, JUNB, and ATF2 was evaluated by quantitative PCR to verify the RNA-seq dataset. For E2F1, JUNB, and ATF2, there was a significant increase in expression in response to the 20% stretch (Figure 7A–C). Meanwhile, there was a significant decrease in STAT1 expression in response to the 20% stretch (6D), thus verifying the differential expression observed in the RNA-seq data. The non-significant changes in expression in response to the 10% stretch serve as a further validation, as these target genes were only differentially expressed at the magnitude of the 20% stretch.

4. Discussion

This study utilized a systems biology-based approach to map the molecular networks and functional responses that govern the endothelial phenotype induced by the physiological and pathological magnitudes of cyclic stretch. Previous studies have investigated the effects of cyclic stretch by using experimental systems that primarily deliver uniaxial stretch at set magnitudes [21]. However, these studies provide limited insight into the differential endothelial responses to the increasing magnitudes of equibiaxial cyclic stretch. This study utilized RNA-seq data collected from ECs subjected to 10% and 20% equibiaxial cyclic stretch to conduct functional pathway analyses and to visualize the gene regulatory networks, providing comprehensive insights into the temporal relationships between endothelial gene expression and the relevant phenotypic responses.
Low-magnitude cyclic stretch (5–10% ∆r) has largely been considered a physiological stretch, associated with an atheroprotective phenotype, while high-magnitude cyclic stretch (>20% ∆r) has largely been considered a pathological stretch, associated with an atheroprone phenotype [24]. The pathological stretch coupled with oscillatory shear stress, imparted by disturbed blood flow, has been implicated in the pathogenesis of vascular disease [44,45,46]. Thus, the aim of this study was twofold. First, we aimed to determine the molecular signatures which constitute the endothelial state during the pathological stretch. Second, we aimed to compare these molecular signatures to the molecular signatures that govern the endothelial state during oscillatory shear stress.

4.1. Endothelial Response to Cyclic Stretch

Following overexpression analysis, we identified four key endothelial responses to cyclic stretch: cell cycle regulation, interferon response, fatty acid metabolism, and mTOR signaling. Cell cycle hyperactivity has been implicated in the pathological progression of vascular disease [47]. In this study, we identified a cohort of cell cycle-related genes, all downstream targets of the E2F family of the transcription factors that are involved in driving the progression of the cell cycle under the pathological cyclic stretch.
Outside of their role in the cell cycle, the E2F transcription factors have also been implicated in the hydrogen peroxide (H2O2)-induced downregulation of the tetrahydrobiopterin salvage enzyme DHFR (dihydrofolate reductase), causing eNOS (endothelial NO synthase) uncoupling and a subsequent increase in blood pressure [48]. Because of this dual role, the targeting of the E2F transcription factors may be a viable therapeutic option to mitigate the progression of vascular disease. The role of hyperactive cell cycle activity and proliferation in the pathological progression of vascular disease was confirmed by protein expression analysis as the expression of both VEGFR2 and Akt, which are known markers of endothelial proliferation modulated by E2F activity, increased in response to the pathological magnitudes of cyclic stretch [49,50].
The genes identified in the interferon response (STAT1, HIF1A, BTG1, CCL2) and cholesterol homeostasis (ATF3, SREBF2) have previously been implicated in the inflammatory response, which occurs during the onset of vascular disease [37,38,39,40,51]. ATF3 and SREBF2 have both been demonstrated to mediate the effects of high-density lipoprotein (HDL) and low-density lipoprotein (LDL), respectively [40,52]. The downregulation of ATF3 in response to pathological stretch may suggest an impairment of the ATF3-dependent atheroprotective effects of HDL during the onset of vascular disease [53]. The development of primary aldosteronism is another factor contributing to vascular disease progression and is the most common cause of reversible vascular disease [54]. Overactive mTOR signaling has been demonstrated to contribute to primary aldosteronism; mTORC1 complex inhibition has been associated with decreased levels of plasma aldosterone and decreased blood pressure [55,56]. Thus, the upregulation of mTOR pathway components in response to pathological stretch may represent another underlying molecular contributor driving vascular disease.
Vascular disease is also characterized by the preferential utilization of glucose for metabolism rather than fatty acids [33,41]. We found that this metabolic profile may be due in part to the upregulation of DLD and the downregulation of ACOT8 in response to pathological stretch. DLD forms a subunit in the pyruvate dehydrogenase complex, a highly regulated enzyme preceding the tricarboxylic acid cycle during glucose metabolism [57]. Meanwhile, ACOT8 plays a role in ß-oxidation during the metabolism of fatty acids [58].
Each of the identified functional responses was driven by a central regulatory network. This regulatory network consisted of a cohort of transcription factors (JUN, MYC, ATF2, TERF1, CREB1, CEBPG, NR3C1, and JUNB) that were specifically upregulated in response to the pathological cyclic stretch relative to the physiological cyclic stretch. Of these functional responses, we found that the inflammatory response was mediated by JUN, MYC, and JUNB; cholesterol homeostasis was mediated by CREB1, JUNB, JUN, ATF2, and CEBPG; fatty acid metabolism was mediated by ATF2 and MYC; and mTOR signaling was mediated by MYC, ATF2, TERF1, JUN, and NR3C1.
The combined activity of these transcription factors may play an Important role in coordinating the molecular transition that occurs as ECs shift from an atheroprotective to an atheroprone phenotype in response to the increasingly pathological magnitudes of cyclic stretch, culminating in vascular disease.

4.2. Effect of Cyclic Stretch on the Endothelial Cytoskeleton

The endothelial cytoskeleton has been shown to be involved in vascular disease [59]. We aimed to evaluate whether cyclic stretch impacted the EC cytoskeletal reorganization by determining cell stiffness in response to physiological and pathological cyclic stretch. We found that the cytoskeletal response is dose-dependent and that the ECs become softer at higher magnitudes of cyclic stretch. Thus, a decrease in EC stiffness may contribute to the development of an atheroprone phenotype. This decrease in stiffness was reversed by the treatment with rapamycin, a known inhibitor of mTOR [60]. We found that there was a greater number of differentially expressed cytoskeletal component genes in response to the 20% stretch relative than to the 10% stretch, signifying that the 20% stretch leads to an increase in the regulation of microtubule, microfilament, and intermediate filament genes. This provides a molecular basis for the larger stiffness decrease that was observed, given that alternations in cell stiffness depend on the corresponding alterations in the cytoskeleton [61]. There were also a greater number of mTOR target cytoskeletal genes that were differentially expressed in the 20% stretch, potentially explaining why the rapamycin treatment had a more significant impact on the stiffness in response to the 20% stretch rather than the 10% stretch. Furthermore, the VE-cadherin mobilization in response to stretch suggests a highly involved cellular response, with a notable example in the fact that VEGFR2 is known to associate with VE-cadherin near cell membranes in a process that modulates receptor functionality [62].
To establish a mechanism for stiffness regulation by cyclic stretch, we constructed a connectivity network consisting of differentially expressed cytoskeletal component genes and mTOR. Two central regulators were identified, AURKA and CDK1, which were both connected to mTOR. AURKA has been shown to drive mTOR activation, while CDK1 has previously been demonstrated to help modulate downstream mTOR effectors [63,64]. Both of these proteins are regulators of microtubule behavior [45,65]. Thus, we hypothesize a potential regulatory mechanism in which AURKA is stimulated by an upstream mechanosensitive pathway to cyclic stretch and leads to downstream mTOR activation. CDK1, which was connected to an additional 21 downstream cytoskeletal component genes (CDC27, PRC1, BUB1, FBXO5, INCENP, BIRC5, KIF11, DLGAP5, TOP2A, CCNB2, ESPL1, MKI67, KIF20B, NDE1, KRT18, LMNB1, VIM, AMPH, CALD1, CCP110, and ACTR3) may then regulate the downstream effects of mTOR kinase activity on the cytoskeleton. The activation of this mechanosensitive pathway is likely to be magnitude-dependent, with a greater magnitude of stretch leading to greater pathway activation, which would explain the increased number of differentially expressed mTOR target genes found in response to the 20% stretch.

4.3. Comparison of Endothelial Response to Cyclic Stretch and Shear Stress

Because of the complementary nature of cyclic stretch and shear stress as hemodynamic forces, we wanted to determine whether similar molecular networks regulate the endothelial response to each. We compared the previously published data that mapped the molecular response of ECs to pulsatile and oscillatory shear stress, as being representative of atheroprotective and atheroprone conditions, respectively [8]. After mapping the expression of the identified effector genes, we found that the expressional pattern of ECs in response to the 10% stretch was most closely related to that of the ECs in response to pulsatile shear stress, while the expressional pattern in response to the 20% stretch was most closely related to that of the ECs in response to oscillatory shear stress. This suggests that the physiological and pathological magnitudes of cyclic stretch and shear stress impart similar molecular and, in turn, functional responses to the endothelium to coordinate an atheroprotective and atheroprone phenotype, respectively.
To determine whether common transcriptional regulators exist in the responses to cyclic stretch and shear stress, we identified transcription factors that were differentially expressed in response to both the 10% stretch and pulsatile shear stress and the 20% stretch and oscillatory shear stress, respectively. We found 22 transcription factors that were commonly expressed between the 10% stretch and the pulsatile shear stress, 11 of which were upregulated (ZNF692, E2F1, CENPA, ATF3, MYBL2, FOSL1, DNMT1, ZNF678, HMGA1, TFAP4, and IRF1) and 11 of which were downregulated (LRP5, ZNF189, ZNF449, SOX17, ZNF41, AEBP1, NR2F2, XBP1, EPAS1, ZBTB41, and ZNF467). In addition, we found 17 transcription factors commonly upregulated between the 20% stretch and the oscillatory shear stress (SNAI1, MYBL2, JUNB, BHLHE40, FOSL1, ZBTB24, BCL6B, E2F8, CSRNP1, ZNF697, HIC1, NFATC2, E2F1, DPF3, KLF3, HMGA1, and RARA) and 1 transcription factor that was commonly downregulated (TSC22D1). The common upregulation of SNAI1 in response to both the 20% stretch and the oscillatory shear stress is particularly of interest, given its central role as a transcriptional regulator of the endothelial to mesenchymal transition (EndMT) and the increased speculation surrounding the role of EndMT in vascular disease progression [66]. Although our analysis is preliminary in scope, the transcription factors provided here provide a basis for future investigations into the shared regulatory networks that drive the molecular response to cyclic stretch and shear stress. Specifically, future genomic analyses evaluating the potential synergistic effects of pathological cyclic stretch and pathological shear stress should be conducted using in vitro systems capable of delivering concurrent cyclic stretch and shear stress, as previously established [67].

5. Conclusions

In this study, we present a thorough overview of the functional and molecular endothelial responses to pathological and physiological cyclic stretch. Through detailed network analysis, we delineated the regulatory transcriptional networks driving the progression from an atheroprotective to an atheroprone endothelial phenotype induced by equibiaxial cyclic stretch. The genes identified here provide a novel insight as an invaluable predicative framework for the dissection of future clinical data on cohorts of patients with vascular disease; a few of these (S100A10, FOS) have already been identified in clinical genomic analyses of vascular disease [68]. Thus, this study verifies and expands upon the regulatory mechanisms which are currently understood and provides gene targets that lay the foundation for further investigation into the pathogenesis of vascular disease.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biom12121837/s1.

Author Contributions

M.W.L.: Methodology, Software, Validation, Formal analysis, Writing—Original Draft, Writing—Review & Editing. N.C.: Conceptualization, Methodology, Formal analysis, Investigation, Writing—Review & Editing, A.C.: Atomic Force Microscopy experiments, Validation. B.K.: Methodology, Investigation. D.R.: Methodology pertaining to RNA Seq Analysis. S.R.: Conceptualization, Resources, Writing—Review & Editing, Funding acquisition. S.Y.R.: Conceptualization, Methodology, Validation, Investigation, Resources, Writing—Original Draft, Supervision, Project administration, Funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

We are grateful for funding support from Hofstra University and the National Institute of Health (NIH) Director’s Transformative award R01 HL128158.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

S.R. is the founder and a non-paid consultant to Angiocrine Bioscience.

Conflicts of Interest

The authors declare no potential conflict of interest.

References

  1. Charbonier, F.; Zamani, M.; Huang, N. Endothelial Cell Mechanotransduction in the Dynamic Vascular Environment. Adv. Biosyst. 2018, 3, 1800252. [Google Scholar] [CrossRef] [PubMed]
  2. Alberio, L.; Dale, G.L. Platelet-collagen interactions: Membrane receptors and intracellular signalling pathways. Eur. J. Clin. Investig. 1999, 29, 1066–1076. [Google Scholar] [CrossRef] [PubMed]
  3. Ma, J.; Ward, E.; Siegel, R.; Jemal, A. Temporal Trends in Mortality in the United States, 1969–2013. JAMA 2015, 314, 1731. [Google Scholar] [CrossRef] [Green Version]
  4. Baeyens, N.; Bandyopadhyay, C.; Coon, B.; Yun, S.; Schwartz, M. Endothelial fluid shear stress sensing in vascular health and disease. J. Clin. Investig. 2016, 126, 821–828. [Google Scholar] [CrossRef] [PubMed]
  5. Chistiakov, D.; Orekhov, A.; Bobryshev, Y. Effects of shear stress on endothelial cells: Go with the flow. Acta Physiol. 2016, 219, 382–408. [Google Scholar] [CrossRef] [PubMed]
  6. Chiu, J.; Chien, S. Effects of Disturbed Flow on Vascular Endothelium: Pathophysiological Basis and Clinical Perspectives. Physiol. Rev. 2011, 91, 327–387. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Gimbrone, M.; García-Cardeña, G. Endothelial Cell Dysfunction and the Pathobiology of Atherosclerosis. Circ. Res. 2016, 118, 620–636. [Google Scholar] [CrossRef] [Green Version]
  8. Ajami, N.; Gupta, S.; Maurya, M.; Nguyen, P.; Li, J.; Shyy, J.; Chen, Z.; Chien, S.; Subramaniam, S. Systems biology analysis of longitudinal functional response of endothelial cells to shear stress. Proc. Natl. Acad. Sci. USA 2017, 114, 10990–10995. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, J.; Friedman, M. Adaptive response of vascular endothelial cells to an acute increase in shear stress magnitude. Am. J. Physiol.-Heart Circ. Physiol. 2012, 302, H983–H991. [Google Scholar] [CrossRef]
  10. Zhou, J.; Li, Y.; Chien, S. Shear Stress-Initiated Signaling and Its Regulation of Endothelial Function. Arterioscler. Thromb. Vasc. Biol. 2014, 34, 2191–2198. [Google Scholar] [CrossRef]
  11. Figueroa, D.; Kemeny, S.; Clyne, A. Glycated Collagen Decreased Endothelial Cell Fibronectin Alignment in Response to Cyclic Stretch Via Interruption of Actin Alignment. J. Biomech. Eng. 2014, 136, 101010. [Google Scholar] [CrossRef] [PubMed]
  12. Adapala, R.; Talasila, P.; Bratz, I.; Zhang, D.; Suzuki, M.; Meszaros, J.; Thodeti, C. PKCα mediates acetylcholine-induced activation of TRPV4-dependent calcium influx in endothelial cells. Am. J. Physiol.-Heart Circ. Physiol. 2011, 301, H757–H765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Katsumi, A.; Naoe, T.; Matsushita, T.; Kaibuchi, K.; Schwartz, M. Integrin Activation and Matrix Binding Mediate Cellular Responses to Mechanical Stretch. J. Biol. Chem. 2005, 280, 16546–16549. [Google Scholar] [CrossRef] [Green Version]
  14. Osawa, M.; Masuda, M.; Kusano, K.; Fujiwara, K. Evidence for a role of platelet endothelial cell adhesion molecule-1 in endothelial cell mechanosignal transduction. J. Cell Biol. 2002, 158, 773–785. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Thodeti, C.; Matthews, B.; Ravi, A.; Mammoto, A.; Ghosh, K.; Bracha, A.; Ingber, D. TRPV4 Channels Mediate Cyclic Strain–Induced Endothelial Cell Reorientation Through Integrin-to-Integrin Signaling. Circ. Res. 2009, 104, 1123–1130. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kaunas, R.; Nguyen, P.; Usami, S.; Chien, S. From the Cover: Cooperative effects of Rho and mechanical stretch on stress fiber organization. Proc. Natl. Acad. Sci. USA 2005, 102, 15895–15900. [Google Scholar] [CrossRef] [Green Version]
  17. Wilkins, J.; Pike, D.; Gibson, C.; Kubota, A.; Shiu, Y. Differential effects of cyclic stretch on bFGF- and VEGF-induced sprouting angiogenesis. Biotechnol. Prog. 2014, 30, 879–888. [Google Scholar] [CrossRef] [PubMed]
  18. Ali, M.; Pearlstein, D.; Mathieu, C.; Schumacker, P. Mitochondrial requirement for endothelial responses to cyclic strain: Implications for mechanotransduction. Am. J. Physiol.-Lung Cell. Mol. Physiol. 2004, 287, L486–L496. [Google Scholar] [CrossRef] [Green Version]
  19. Xiong, Y.; Hu, Z.; Han, X.; Jiang, B.; Zhang, R.; Zhang, X.; Lu, Y.; Geng, C.; Li, W.; He, Y.; et al. Hypertensive stretch regulates endothelial exocytosis of Weibel-Palade bodies through VEGF receptor 2 signaling pathways. Cell Res. 2013, 23, 820–834. [Google Scholar] [CrossRef] [Green Version]
  20. Contreras, A.; Dormond, O.; Edelbauer, M.; Calzadilla, K.; Hoerning, A.; Pal, S.; Briscoe, D. mTOR—Understanding the Clinical Effects. Transplant. Proc. 2008, 40, S9–S12. [Google Scholar] [CrossRef]
  21. Wang, J.; Goldschmidt-Clermont, P.; Wille, J.; Yin, F. Specificity of endothelial cell reorientation in response to cyclic mechanical stretching. J. Biomech. 2001, 34, 1563–1572. [Google Scholar] [CrossRef] [PubMed]
  22. Huang, A.; Balestrini, J.; Udelsman, B.; Zhou, K.; Zhao, L.; Ferruzzi, J.; Starcher, B.C.; Levene, M.J.; Humphrey, J.D.; Niklason, L.E. Biaxial Stretch Improves Elastic Fiber Maturation, Collagen Arrangement and Mechanical Properties in Engineered Arteries. Tissue Eng. Part C Methods 2016, 22, 524–533. [Google Scholar] [CrossRef]
  23. Ursekar, C.; Teo, S.; Hirata, H.; Harada, I.; Chiam, K.; Sawada, Y. Design and Construction of an Equibiaxial Cell Stretching System That Is Improved for Biochemical Analysis. PLoS ONE 2014, 9, e90665. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Anwar, M.; Shalhoub, J.; Lim, C.; Gohel, M.; Davies, A. The Effect of Pressure- Induced Mechanical Stretch on Vascular Wall Differential Gene Expression. J. Vasc. Res. 2012, 49, 463–478. [Google Scholar] [CrossRef] [PubMed]
  25. Armulik, A.; Genové, G.; Betsholtz, C. Pericytes: Developmental, Physiological, and Pathological Perspectives, Problems and Promises. Dev. Cell 2011, 21, 193–215. [Google Scholar] [CrossRef] [Green Version]
  26. Seandel, M.; Butler, J.; Kobayashi, H.; Hooper, A.; White, I.; Zhang, F.; Vertes, E.L.; Kobayashi, M.; Zhang, Y.; Shmelkov, S.V.; et al. Generation of a functional and durable vascular niche by the adenoviral E4ORF1 gene. Proc. Natl. Acad. Sci. USA 2008, 105, 19288–19293. [Google Scholar] [CrossRef] [Green Version]
  27. Robinson, M.; McCarthy, D.; Smyth, G. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  28. Subramanian, A.; Tamayo, P.; Mootha, V.; Mukherjee, S.; Ebert, B.; Gillette, M.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Montojo, J.; Zuberi, K.; Rodriguez, H.; Kazi, F.; Wright, G.; Donaldson, S.; Morris, Q.; Bader, G.D. GeneMANIA Cytoscape plugin: Fast gene function predictions on the desktop. Bioinformatics 2010, 26, 2927–2928. [Google Scholar] [CrossRef] [Green Version]
  30. Guo, P.; Poulos, M.; Palikuqi, B.; Badwe, C.; Lis, R.; Kunar, B.; Ding, B.-S.; Rabbany, S.Y.; Shido, K.; Butler, J.M.; et al. Endothelial jagged-2 sustains hematopoietic stem and progenitor reconstitution after myelosuppression. J. Clin. Investig. 2017, 127, 4242–4256. [Google Scholar] [CrossRef]
  31. Tondon, A.; Kaunas, R. The Direction of Stretch-Induced Cell and Stress Fiber Orientation Depends on Collagen Matrix Stress. PLoS ONE 2014, 9, e89592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, L.; Bao, H.; Wang, K.; Zhang, P.; Yao, Q.; Chen, X.; Huang, K.; Qi, Y.; Jiang, Z. Secreted miR- 27a Induced by Cyclic Stretch Modulates the Proliferation of Endothelial Cells in Hypertension via GRK6. Sci. Rep. 2017, 7, srep41058. [Google Scholar] [CrossRef]
  33. Polak-Iwaniuk, A.; Harasim-Symbor, E.; Gołaszewska, K.; Chabowski, A. How Hypertension Affects Heart Metabolism. Front. Physiol. 2019, 10, 435. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Hamada, M.; Nishio, I.; Baba, A.; Fukuda, K.; Takeda, J.; Ura, M.; Hano, T.; Kuchii, M.; Masuyama, Y. Enhanced DNA synthesis of cultured vascular smooth muscle cells from spontaneously hypertensive rats Difference of response to growth factor, intracellular free calcium concentration and DNA synthesizing cell cycle. Atherosclerosis 1990, 81, 191–198. [Google Scholar] [CrossRef]
  35. Kobayashi, S.; Nagino, M.; Komatsu, S.; Naruse, K.; Nimura, Y.; Nakanishi, M.; Sokabe, M. Stretch-induced IL-6 secretion from endothelial cells requires NF-κB activation. Biochem. Biophys. Res. Commun. 2003, 308, 306–312. [Google Scholar] [CrossRef] [PubMed]
  36. Kopitar-Jerala, N. The Role of Interferons in Inflammation and Inflammasome Activation. Front. Immunol. 2017, 8, 873. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Palazon, A.; Goldrath, A.; Nizet, V.; Johnson, R. HIF Transcription Factors, Inflammation and Immunity. Immunity 2014, 41, 518–528. [Google Scholar] [CrossRef] [Green Version]
  38. Rauch, I.; Müller, M.; Decker, T. The regulation of inflammation by interferons and their STATs. JAK-STAT 2013, 2, e23820. [Google Scholar] [CrossRef] [Green Version]
  39. Xiao, H.; Lu, M.; Lin, T.; Chen, Z.; Chen, G.; Wang, W.; Marin, T.; Shentu, T.-P.; Wen, L.; Gongol, B.; et al. Sterol Regulatory Element Binding Protein 2 Activation of NLRP3 Inflammasome in Endothelium Mediates Hemodynamic-Induced Atherosclerosis Susceptibility. Circulation 2013, 128, 632–642. [Google Scholar] [CrossRef] [Green Version]
  40. Zhang, H.; Reilly, M. Anti-Inflammatory Effects of High-Density Lipoprotein Through Activating Transcription Factor 3. Arterioscler. Thromb. Vasc. Biol. 2014, 34, e11–e12. [Google Scholar] [CrossRef]
  41. Talati, M.; Hemnes, A. Fatty Acid Metabolism in Pulmonary Arterial Hypertension: Role in Right Ventricular Dysfunction and Hypertrophy. Pulm. Circ. 2015, 5, 269–278. [Google Scholar] [CrossRef] [Green Version]
  42. Sciarretta, S.; Forte, M.; Frati, G.; Sadoshima, J. New Insights into the Role of mTOR Signaling in the Cardiovascular System. Circ. Res. 2018, 122, 489–505. [Google Scholar] [CrossRef]
  43. Liu, L.; Luo, Y.; Chen, L.; Shen, T.; Xu, B.; Chen, W.; Zhou, H.; Han, X.; Huang, S. Rapamycin Inhibits Cytoskeleton Reorganization and Cell Motility by Suppressing RhoA Expression and Activity. J. Biol. Chem. 2010, 285, 38362–38373. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Cheng, C.; Tempel, D.; van Haperen, R.; van der Baan, A.; Grosveld, F.; Daemen, M.; Krams, R.; De Crom, R. Atherosclerotic Lesion Size and Vulnerability Are Determined by Patterns of Fluid Shear Stress. Circulation 2006, 113, 2744–2753. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Katayama, H.; Sasai, K.; Kloc, M.; Brinkley, B.; Sen, S. Aurora kinase-A regulates kinetochore/chromatin associated microtubule assembly in human cells. Cell Cycle 2008, 7, 2691–2704. [Google Scholar] [CrossRef] [Green Version]
  46. Thubrikar, M.; Robicsek, F. Pressure-induced arterial wall stress and atherosclerosis. Ann. Thorac. Surg. 1995, 59, 1594–1603. [Google Scholar] [CrossRef] [PubMed]
  47. Tanner, F.; Greutert, H.; Barandier, C.; Frischknecht, K.; Lüscher, T.F. Different Cell Cycle Regulation of Vascular Smooth Muscle in Genetic Hypertension. Hypertension 2003, 42, 184–188. [Google Scholar] [CrossRef] [Green Version]
  48. Li, H.; Li, Q.; Zhang, Y.; Liu, W.; Gu, B.; Narumi, T.; Siu, K.L.; Youn, J.Y.; Liu, P.; Yang, X.; et al. Novel Treatment of Hypertension by Specifically Targeting E2F for Restoration of Endothelial Dihydrofolate Reductase and eNOS Function under Oxidative Stress. Hypertension 2019, 73, 179–189. [Google Scholar] [CrossRef]
  49. Bakker, W.; Weijts, B.; Westendorp, B.; de Bruin, A. HIF proteins connect the RB- E2F factors to angiogenesis. Transcription 2013, 4, 62–66. [Google Scholar] [CrossRef] [Green Version]
  50. Chaussepied, M.; Ginsberg, D. Transcriptional Regulation of AKT Activation by E2F. Mol. Cell 2004, 16, 831–837. [Google Scholar] [CrossRef]
  51. De Miguel, C.; Rudemiller, N.; Abais, J.; Mattson, D. Inflammation and Hypertension: New Understandings and Potential Therapeutic Targets. Curr. Hypertens. Rep. 2014, 17, 507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Miserez, A.; Muller, P.; Barella, L.; Barella, S.; Staehelin, H.; Leitersdorf, E.; Kark, J.D.; Friedlander, Y. Sterol-regulatory element-binding protein (SREBP)-2 contributes to polygenic hypercholesterolaemia. Atherosclerosis 2002, 164, 15–26. [Google Scholar] [CrossRef]
  53. De Nardo, D.; Labzin, L.; Kono, H.; Seki, R.; Schmidt, S.; Beyer, M.; Xu, D.; Zimmer, S.; Lahrmann, C.; Schildberg, F.A.; et al. High- density lipoprotein mediates anti-inflammatory reprogramming of macrophages via the transcriptional regulator ATF3. Nat. Immunol. 2013, 15, 152–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Aronova, A.; Fahey, T., III; Zarnegar, R. Management of hypertension in primary aldosteronism. World J. Cardiol. 2014, 6, 227. [Google Scholar] [CrossRef]
  55. Su, H.; Gu, Y.; Li, F.; Wang, Q.; Huang, B.; Jin, X.; Ning, G.; Sun, F. The PI3K/AKT/mTOR Signaling Pathway Is Overactivated in Primary Aldosteronism. PLoS ONE 2013, 8, e62399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Trinh, B.; Hepprich, M.; Betz, M.; Burkard, T.; Cavelti-Weder, C.; Seelig, E.; Meienberg, F.; Kratschmar, D.V.; Beuschlein, F.; Reincke, M.; et al. Treatment of Primary Aldosteronism with mTORC1 Inhibitors. J. Clin. Endocrinol. Metab. 2019, 104, 4703–4714. [Google Scholar] [CrossRef] [PubMed]
  57. Lacroix, M.; Rodier, G.; Kirsh, O.; Houles, T.; Delpech, H.; Seyran, B.; Gayte, L.; Casas, F.; Pessemesse, L.; Heuillet, M.; et al. E4F1 controls a transcriptional program essential for pyruvate dehydrogenase activity. Proc. Natl. Acad. Sci. USA 2016, 113, 10998–11003. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Westin, M.; Hunt, M.; Alexson, S. The Identification of a Succinyl-CoA Thioesterase Suggests a Novel Pathway for Succinate Production in Peroxisomes. J. Biol. Chem. 2005, 280, 38125–38132. [Google Scholar] [CrossRef] [Green Version]
  59. Atwal, G.; Sun, Y.; Chen, J.; Komarova, Y. Role of Microtubule Cytoskeleton as a Novel Target Against Hypoxia Induced Pulmonary Hypertension. J. Am. Coll. Cardiol. 2017, 69, 1882. [Google Scholar] [CrossRef]
  60. Barilli, A.; Visigalli, R.; Sala, R.; Gazzola, G.; Parolari, A.; Tremoli, E.; Bonomini, S.; Simon, A.; Closs, E.; Dall’Asta, V.; et al. In human endothelial cells rapamycin causes mTORC2 inhibition and impairs cell viability and function. Cardiovasc. Res. 2008, 78, 563–571. [Google Scholar] [CrossRef]
  61. Gavara, N.; Chadwick, R. Relationship between cell stiffness and stress fiber amount, assessed by simultaneous atomic force microscopy and live-cell fluorescence imaging. Biomech. Model. Mechanobiol. 2015, 15, 511–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Tian, Y.; Gawlak, G.; O’Donnell, J.; Birukova, A.; Birukov, K. Activation of Vascular Endothelial Growth Factor (VEGF) Receptor 2 Mediates Endothelial Permeability Caused by Cyclic Stretch. J. Biol. Chem. 2016, 291, 10032–10045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Pal, S.; He, M.; Tong, T.; Wu, H.; Liu, X.; Lau, C.; Wang, J.-H.; Warden, C.; Wu, X.; Signoretti, S.; et al. RNA-seq Reveals Aurora Kinase–Driven mTOR Pathway Activation in Patients with Sarcomatoid Metastatic Renal Cell Carcinoma. Mol. Cancer Res. 2014, 13, 130–137. [Google Scholar] [CrossRef] [Green Version]
  64. Ramirez-Valle, F.; Badura, M.; Braunstein, S.; Narasimhan, M.; Schneider, R. Mitotic Raptor Promotes mTORC1 Activity, G2/M Cell Cycle Progression, and Internal Ribosome Entry Site-Mediated mRNA Translation. Mol. Cell. Biol. 2010, 30, 3151–3164. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Fourest-Lieuvin, A.; Peris, L.; Gache, V.; Garcia-Saez, I.; Juillan-Binard, C.; Lantez, V.; Job, D. Microtubule Regulation in Mitosis: Tubulin Phosphorylation by the Cyclin-dependent Kinase Cdk1. Mol. Biol. Cell 2006, 17, 1041–1050. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Kovacic, J.; Dimmeler, S.; Harvey, R.; Finkel, T.; Aikawa, E.; Krenning, G.; Baker, A. Endothelial to Mesenchymal Transition in Cardiovascular Disease. J. Am. Coll. Cardiol. 2019, 73, 190–209. [Google Scholar] [CrossRef]
  67. Meza, D.; Musmacker, B.; Steadman, E.; Stransky, T.; Rubenstein, D.; Yin, W. Endothelial Cell Biomechanical Responses are Dependent on Both Fluid Shear Stress and Tensile Strain. Cell. Mol. Bioeng. 2019, 12, 311–325. [Google Scholar] [CrossRef]
  68. Huan, T.; Esko, T.; Peters, M.; Pilling, L.; Schramm, K.; Schurmann, C.; Chen, B.H.; Liu, C.; Joehanes, R.; Johnson, A.D.; et al. A Meta-analysis of Gene Expression Signatures of Blood Pressure and Hypertension. PLoS Genet. 2015, 11, e1005035. [Google Scholar] [CrossRef]
Biomolecules 12 01837 g001 550
Figure 1. Calibration of the biaxial stretch system. (A) Schematic of the sectioned PDMS membrane utilized for calibration. Each box (A–E) is 0.25 cm × 0.25 cm. (B) Membrane displacement after stretch at each of the five locations indicated in (A). (C) Fluorescent image demonstrating staining of actin (green) and DAPI counterstain (blue) after 24 h of static culture on a PDMS membrane at 43× magnification. (D) Fluorescent image demonstrating actin staining (green) and DAPI counterstain (blue) after 24 h of biaxial stretch at 43× magnification. Scale bar represents 20 µm.
Figure 1. Calibration of the biaxial stretch system. (A) Schematic of the sectioned PDMS membrane utilized for calibration. Each box (A–E) is 0.25 cm × 0.25 cm. (B) Membrane displacement after stretch at each of the five locations indicated in (A). (C) Fluorescent image demonstrating staining of actin (green) and DAPI counterstain (blue) after 24 h of static culture on a PDMS membrane at 43× magnification. (D) Fluorescent image demonstrating actin staining (green) and DAPI counterstain (blue) after 24 h of biaxial stretch at 43× magnification. Scale bar represents 20 µm.
Biomolecules 12 01837 g001
Biomolecules 12 01837 g002 550
Figure 2. Endothelial response to pathological cyclic stretch. (A) Schematic summarizing the endothelial response to pathological magnitudes of cyclic stretch (20%). Green nodes list sets of upstream regulatory genes that correspond to the downstream cellular responses subsequently listed in red nodes. Downstream cellular responses were identified by GSEA analysis. The false discovery rate (FDR) of each identified gene set is listed. An FDR < 0.25 was utilized as the threshold for further investigation. Relevant target genes within each response are listed below the respective node. Target genes colored red are upregulated at 20% stretch while genes colored blue are downregulated at 20% stretch, relative to the static (no stretch) control. (B) Western blot analysis of protein modulators involved in the E2F response to pathological cyclic stretch. S represents expression after 24 h of static culture, 10% represents expression after 24 h of 10% cyclic stretch, and 20% represents expression after 24 h of 20% cyclic stretch.
Figure 2. Endothelial response to pathological cyclic stretch. (A) Schematic summarizing the endothelial response to pathological magnitudes of cyclic stretch (20%). Green nodes list sets of upstream regulatory genes that correspond to the downstream cellular responses subsequently listed in red nodes. Downstream cellular responses were identified by GSEA analysis. The false discovery rate (FDR) of each identified gene set is listed. An FDR < 0.25 was utilized as the threshold for further investigation. Relevant target genes within each response are listed below the respective node. Target genes colored red are upregulated at 20% stretch while genes colored blue are downregulated at 20% stretch, relative to the static (no stretch) control. (B) Western blot analysis of protein modulators involved in the E2F response to pathological cyclic stretch. S represents expression after 24 h of static culture, 10% represents expression after 24 h of 10% cyclic stretch, and 20% represents expression after 24 h of 20% cyclic stretch.
Biomolecules 12 01837 g002
Biomolecules 12 01837 g003 550
Figure 3. Regulatory gene networks governing the endothelial response to cyclic stretch. Networks represent (A) endothelial cell cycle regulation and (B) endothelial transcriptional regulation in response to pathological magnitudes of cyclic stretch (20%). Green nodes represent regulatory genes that were identified due to their high number of network connections, which in turn are signified by red nodes representing target genes. Node size of regulatory genes correlates with the number of direct network connections for the respective node. Network was constructed using the GeneMANIA Cytoscape plugin, which utilizes a guilt-by-association approach to map physical interactions between input genes. The input gene list was curated from gene sets identified by GSEA analysis. Only direct first neighbor connections to the regulatory genes are shown here. The full network maps can be seen in Supplementary Figures S1 and S2.
Figure 3. Regulatory gene networks governing the endothelial response to cyclic stretch. Networks represent (A) endothelial cell cycle regulation and (B) endothelial transcriptional regulation in response to pathological magnitudes of cyclic stretch (20%). Green nodes represent regulatory genes that were identified due to their high number of network connections, which in turn are signified by red nodes representing target genes. Node size of regulatory genes correlates with the number of direct network connections for the respective node. Network was constructed using the GeneMANIA Cytoscape plugin, which utilizes a guilt-by-association approach to map physical interactions between input genes. The input gene list was curated from gene sets identified by GSEA analysis. Only direct first neighbor connections to the regulatory genes are shown here. The full network maps can be seen in Supplementary Figures S1 and S2.
Biomolecules 12 01837 g003
Biomolecules 12 01837 g004a 550Biomolecules 12 01837 g004b 550
Figure 4. Regulation of the cytoskeleton by cyclic stretch. Cell stiffness of (A) HUVECs and (B) rapamycin-treated HUVECs as determined through AFM. * indicates p < 0.05 between respective groups. (C) Heatmap demonstrating the expression of cytoskeletal component genes grouped into their respective classifications in response to 10% and 20% stretch. (D) Heatmap demonstrating the expression of cytoskeletal mTOR-target genes grouped into their respective classifications in response to 10% and 20% stretch. For both (C,D), expression was normalized to the static (no stretch) control. Cytoskeletal genes were identified through the Molecular Signatures Database.
Figure 4. Regulation of the cytoskeleton by cyclic stretch. Cell stiffness of (A) HUVECs and (B) rapamycin-treated HUVECs as determined through AFM. * indicates p < 0.05 between respective groups. (C) Heatmap demonstrating the expression of cytoskeletal component genes grouped into their respective classifications in response to 10% and 20% stretch. (D) Heatmap demonstrating the expression of cytoskeletal mTOR-target genes grouped into their respective classifications in response to 10% and 20% stretch. For both (C,D), expression was normalized to the static (no stretch) control. Cytoskeletal genes were identified through the Molecular Signatures Database.
Biomolecules 12 01837 g004aBiomolecules 12 01837 g004b
Biomolecules 12 01837 g005a 550Biomolecules 12 01837 g005b 550
Figure 5. Cytoskeletal activation by cyclic stretch. (A) Regulatory network highlighting the connections between cytoskeletal component genes in response to 20% stretch. Green nodes represent central regulators, while red and yellow nodes represent first and second degree connections downstream of CDK1, respectively. (B) Fluorescent image demonstrating staining of actin (green), VE-cadherin staining (red), and DAPI counterstain (blue) after 24 h of static culture on a PDMS membrane at 43× magnification. Scale bar represents 20 µm. (C) Fluorescent image demonstrating actin staining (green), VE-cadherin staining (red), and DAPI counterstain (blue) after 24 h of biaxial stretch at 43× magnification.
Figure 5. Cytoskeletal activation by cyclic stretch. (A) Regulatory network highlighting the connections between cytoskeletal component genes in response to 20% stretch. Green nodes represent central regulators, while red and yellow nodes represent first and second degree connections downstream of CDK1, respectively. (B) Fluorescent image demonstrating staining of actin (green), VE-cadherin staining (red), and DAPI counterstain (blue) after 24 h of static culture on a PDMS membrane at 43× magnification. Scale bar represents 20 µm. (C) Fluorescent image demonstrating actin staining (green), VE-cadherin staining (red), and DAPI counterstain (blue) after 24 h of biaxial stretch at 43× magnification.
Biomolecules 12 01837 g005aBiomolecules 12 01837 g005b
Biomolecules 12 01837 g006 550
Figure 6. Comparison of gene expression changes in response to shear stress and cyclic stretch. (A) Heatmap profiling expression changes of network target genes identified in Figure 3 in response to oscillatory shear, pulsatile shear (PS), 20% stretch, and 10% stretch. OS and PS expression data were compiled from previously published data by Ajami et al., GEO Accession Number GSE103672. For each condition, expression was normalized to the static (no shear/stretch) control and scaled by row. Hierarchical clustering analysis was performed on each column and row. (B) Identification of transcription factors that demonstrated changes in expression between both pulsatile shear and 10% stretch and between both oscillatory shear and 20% stretch. Transcription factors colored in red are upregulated in response to the respective conditions, while transcription factors colored in blue are downregulated.
Figure 6. Comparison of gene expression changes in response to shear stress and cyclic stretch. (A) Heatmap profiling expression changes of network target genes identified in Figure 3 in response to oscillatory shear, pulsatile shear (PS), 20% stretch, and 10% stretch. OS and PS expression data were compiled from previously published data by Ajami et al., GEO Accession Number GSE103672. For each condition, expression was normalized to the static (no shear/stretch) control and scaled by row. Hierarchical clustering analysis was performed on each column and row. (B) Identification of transcription factors that demonstrated changes in expression between both pulsatile shear and 10% stretch and between both oscillatory shear and 20% stretch. Transcription factors colored in red are upregulated in response to the respective conditions, while transcription factors colored in blue are downregulated.
Biomolecules 12 01837 g006
Biomolecules 12 01837 g007 550
Figure 7. qPCR validation of RNA-seq. Gene expression of (A) E2F1, (B) JUNB, (C) ATF2, and (D) STAT1 was determined in response to 24 h of 10% and 20% cyclic stretch. Fold change was calculated by normalizing expression to static controls (fold change = 1). * indicates p < 0.05 relative to the static control using one-way ANOVA followed by a post-hoc Fisher’s LSD multiple comparison test, while # indicates p < 0.05 between stretch groups.
Figure 7. qPCR validation of RNA-seq. Gene expression of (A) E2F1, (B) JUNB, (C) ATF2, and (D) STAT1 was determined in response to 24 h of 10% and 20% cyclic stretch. Fold change was calculated by normalizing expression to static controls (fold change = 1). * indicates p < 0.05 relative to the static control using one-way ANOVA followed by a post-hoc Fisher’s LSD multiple comparison test, while # indicates p < 0.05 between stretch groups.
Biomolecules 12 01837 g007
Table
Table 1. Cell cycle regulatory network interaction list. Connectivity value is included in parentheses next to each regulatory gene and represents the number of node connections associated with each regulatory gene.
Table 1. Cell cycle regulatory network interaction list. Connectivity value is included in parentheses next to each regulatory gene and represents the number of node connections associated with each regulatory gene.
CDK2 (49)PCNA (40)RPA1 (39)CDK1 (35)BRCA1 (32)MCM2 (32)
PA2G4USP1TIPINTIPINORC2RPA3
PAICSTIMELESSEXOSC8EXOSC8MSH6ORC2
MYBL2POLD3SMC4SMC4CCNB1PLK1
HMGA1POLD2TMPOTMPOH2AFXCDK6
CCNE2POLD1CHTF18CHTF18HMMRORC5
HMGB2RPA1UNGUNGCDK4L3MBTL1
CDKN3RPA2NOP56NOP56TUBBMCM7
CKS2RPA3RFC2RFC2EZH2ORC4
CCNA1RFC2RFC3RFC3MYCMCM5
CDC25ARFC3SNRPBSNRPBUBE2TGINS3
GMNNCHTF18RPA3RPA3MLH1CDC6
CKS1BHUS1PAN2PAN2DBF4CDK4
CDC20MLH1RPA2RPA2MSH3ORC6
CDC7DSCC1POLD1POLD1CHEK2MCM3
DBF4POLEPCNAPCNARBBP7DBF4
MCM4UNGTOP2ATOP2ADDX39ACCNB2
MCM6PMS2MSH3MSH3ORC3CDC45
CDC45TOP2AHUS1HUS1BARD1MCM4
ORC2MSH3RAD9ARAD9ACDKN2ACKS2
MCM7MSH6MSH6MSH6TOP2ACDKN2A
ORC1RAD9AMSH2MSH2RPA1RPA1
CDC6DNMT1ATRATRTP53ASF1B
CDT1MSH2BRCA1BRCA1TUBG1MCM10
CDKN1ACHEK1TK1TK1MSH2RPA2
CCNB1BARD1CDK1CDK1ATRATR
CDK1CDK6CCNB1CCNB1CHEK1ORC1
TUBG1LIG1CCNA1CCNA1CDK1CCNA1
LMNB1CDKN2AORC6ORC6CCNA1CDC7
CCNB2TP53MCM4MCM4CCND1MMS22L
CCND1CDK4MCM6MCM6AURKACCND1
MCM2CCND3MCM7MCM7CDK2MCM6
MCM3CCND1MCM2MCM2RFC2CDK2
POLD1 POLD2CDK2 CDC6MCM3 ORC2MCM3 ORC2
RPA1CDKN1ACDK2CDK2
PCNACCNB2CDC5LCDC5L
BRCA1CDK1TP53TP53
TUBB CDC5LCCNB1 CDT1TUBB RRM2TUBB RRM2
BIRC5L3MBTL1
CCND3
CDK6
MSH2 SMC4
Table
Table 2. Transcriptional regulatory network interaction list. Connectivity value is included in parentheses next to each regulatory gene and represents the number of node connections associated with each regulatory gene.
Table 2. Transcriptional regulatory network interaction list. Connectivity value is included in parentheses next to each regulatory gene and represents the number of node connections associated with each regulatory gene.
JUN (39)MYC (37)TERF1 (20)ATF2 (29)NR3C1 (15)JUNB (15)CREB1 (16)CEBPG (15)
JUNWARSHSP90AA1NR4A1SAT1SREBF2CEBPE
DDIT3HIF1APGMIFOSL1MVPDDIT3ZNF451DDIT3
TRAF2
CCL2		HSP90AA1FOXJ3CEBPACLUATF3JUNATF2
FOSL1MCM2BATF3SMARCC1CCL2MYCATF7
ATF3
TGIF1
JUNB		CEBPAHMOX1JUNCREB1SMAD4NR3C1BATF
CTNNB1TBPL1ATF7POU2F1JUNATF7BATF3
GSK3BDDX39AFOSCEBPABATFDR1CEBPA
DBP
BATF
MAFG		RUNX1ENO2CTNNB1FOSMAFGDLDATF3
SERPINH1ACAT2CEBPGJUNBATF3ETS1ATF5
TCF12MVKATF3MAFFMAFKFOSMAFK
BATF3
MAFK
HIF1A		KDM5BZNF281ETS1ETS1ATF2NFIL3DBP
HSPD1LDHADDIT3KMT2AATF7POU2F1FOS
NMIGAPDHBACH1RARAFOSL1MTF2FOSL1
HDAC9
CEBPE
NFATC2		CREB1TPI1BATFNR3C2FOSRFX3NFIL3
CDK6HSPH1JUNBHSPD1ETS1ZNF436JUN
SMAD4
FOS
ETS1		CXCR4NR4A1SERPINH1 ZNF92
TOP2APLOD2SMAD4
NFYCCASP7
ATF7
ATF2
FOSL1		SMARCC1ACLYCANX
HSPH1ALDOACD44
PIM1 NFYA
ZNF281
CCT6A HSPH1
HSP90AA1
CEBPA
CEBPG		ETV3 DLD
IDH3B HSPD1
CTNNB1
GSK3B
RUNX1		GCDH ALDOA
ACOT8 TPI1
IKBKB
NR3C1
STAT1		ZNF121 LDHA
NUP205 GAPDH
POU2F1
CXCR4
TOP2A		PTPN11 S100A10
IDH3G ACADVL
SSR1
LGALS1
NFYA
MYC
CREB1		ERCC3
PLK1
GCLM
TCF4
BTG1		ADNP
RBPJ
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lai, M.W.; Chow, N.; Checco, A.; Kunar, B.; Redmond, D.; Rafii, S.; Rabbany, S.Y. Systems Biology Analysis of Temporal Dynamics That Govern Endothelial Response to Cyclic Stretch. Biomolecules 2022, 12, 1837. https://0-doi-org.brum.beds.ac.uk/10.3390/biom12121837

AMA Style

Lai MW, Chow N, Checco A, Kunar B, Redmond D, Rafii S, Rabbany SY. Systems Biology Analysis of Temporal Dynamics That Govern Endothelial Response to Cyclic Stretch. Biomolecules. 2022; 12(12):1837. https://0-doi-org.brum.beds.ac.uk/10.3390/biom12121837

Chicago/Turabian Style

Lai, Michael W., Nathan Chow, Antonio Checco, Balvir Kunar, David Redmond, Shahin Rafii, and Sina Y. Rabbany. 2022. "Systems Biology Analysis of Temporal Dynamics That Govern Endothelial Response to Cyclic Stretch" Biomolecules 12, no. 12: 1837. https://0-doi-org.brum.beds.ac.uk/10.3390/biom12121837

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop