Multiomics analysis reveals the mechanical stress-dependent changes in trabecular meshwork cytoskeletal-extracellular matrix interactions

Trabecular meshwork (TM) tissue is subjected to constant mechanical stress due to the ocular pulse created by the cardiac cycle. This brings about alterations in the membrane lipids and associated cell–cell adhesion and cell–extracellular matrix (ECM) interactions, triggering intracellular signaling responses to counter mechanical insults. A loss of such response can lead to elevated intraocular pressure (IOP), a major risk factor for primary open-angle glaucoma. This study is aimed to understand the changes in signaling responses by TM subjected to mechanical stretch. We utilized multiomics to perform an unbiased mRNA sequencing to identify changes in transcripts, mass spectrometry- (MS-) based quantitative proteomics for protein changes, and multiple reaction monitoring (MRM) profiling-based MS and high-performance liquid chromatography (HPLC-) based MS to characterize the lipid changes. We performed pathway analysis to obtain an integrated map of TM response to mechanical stretch. The human TM cells subjected to mechanical stretch demonstrated an upregulation of protein quality control, oxidative damage response, pro-autophagic signal, induction of anti-apoptotic, and survival signaling. We propose that mechanical stretch-induced lipid signaling via increased ceramide and sphingomyelin potentially contributes to increased TM stiffness through actin-cytoskeleton reorganization and profibrotic response. Interestingly, increased phospholipids and diacylglycerol due to mechanical stretch potentially enable cell membrane remodeling and changes in signaling pathways to alter cellular contractility. Overall, we propose the mechanistic interplay of macromolecules to bring about a concerted cellular response in TM cells to achieve mechanotransduction and IOP regulation when TM cells undergo mechanical stretch.

Introduction

Glaucoma is an aging disease that causes irreversible blindness worldwide, affecting ∼80 million people across the globe (Tham et al., 2014; Susanna et al., 2015). Primary open-angle glaucoma (POAG) is the most common form of glaucoma, affecting over 57.5 million people worldwide (Allison et al., 2020). Elevated intraocular pressure (IOP) is the principal and sole treatable risk factor in POAG (Weinreb et al., 2016), which is regulated by aqueous humor (AH) production and drainage.

The conventional aqueous humor (AH) outflow pathway comprises trabecular meshwork (TM), juxtacanalicular connective tissue (JCT), and Schlemm’s canal (SC), which ends in the collector channels and aqueous veins (Johnstone et al., 2021). This drainage structure aids in the removal of 60%–80% of AH in humans (Hann and Fautsch, 2015). The TM tissue comprises beams of endothelial-type cells fenestrated by extracellular matrix (ECM) and is highly contractile, exhibiting a smooth muscle-like property (Pattabiraman and Toris, 2016). The TM is constantly stretched due to the ocular pulse created by the cardiac cycle, blinking, and eye movement inducing physiological changes in IOP (Johnstone, 2004a; Xin et al., 2017). The IOP is regulated due to the passage of AH across the TM-SC via the transcellular and/or paracellular flow (Raviola and Raviola, 1981; Johnson and Erickson, 2000). Though ECM imparts the resistance to fluid flow via the TM outflow pathway, it is not completely believed to be the only mechanistic impedance to AH drainage (Acott and Kelley, 2008; Pattabiraman and Toris, 2016; Acott et al., 2021; Johnstone et al., 2021). The distention and contraction phenomenon of TM-SC tissues is regulated by the elastic properties of the tissue (Stumpff and Wiederholt, 2000; Wiederholt et al., 2000). For the efficient function of TM, the elastic nature of the TM must be maintained. Additionally, in metabolically active tissues such as TM and SC (Reina-Torres et al., 2020), acute stress should induce protective measures to maintain tissue homeostasis. Failure to achieve homeostasis can lead to decreased AH flow and increased AH retention leading to elevated IOP. In fact, mechanistic evidence points out that experimental elevation of IOP can drive the cells toward an increased contractile phenotype by activating Rho GTPases, which is known to increase resistance to AH drainage (Pattabiraman et al., 2015a; Pattabiraman et al., 2015b). In glaucoma, the TM has a greater stiffness, and the elastic nature of the TM-SC is compromised (Last et al., 2011; Xin et al., 2016; Wang et al., 2017; Vahabikashi et al., 2019). The cardiac cycle sets up the pulsatile motion of the TM (Johnstone, 2004b). A loss of mechanosensation or transduction of such pulsatile mechanical stress can aid in TM stiffness, which in turn affects the pulsatile motion of the TM. A recent study using phase-sensitive optical coherence tomography (PhS-OCT) showed that the eyes with POAG had a significantly lower maximum velocity of the pulsatile motion and the motion amplitude of the TM compared to the healthy eye (Gao et al., 2020). This can be attributed to the increased stiffness and less elasticity of TM tissue. Thus, indicating a direct relationship between the mechanical stress on TM, TM elasticity, TM stiffness, and regulation of AH drainage.

A well-accepted model to study the mechanical stress on TM is the usage of cyclic mechanical stretch (CMS) (Matsuo et al., 1996; Liton et al., 2005; Vittal et al., 2005; Luna et al., 2009; Ramos et al., 2009; Youngblood et al., 2020). Studies have tested the effects of CMS on the changes in specific genes, proteins, and metabolites, including prostaglandins (Matsuo et al., 1996), matrix metalloproteases (Okada et al., 1998; Bradley et al., 2001; Keller et al., 2007), oculomedin (Sato et al., 1999), transforming growth factor beta-1 (TGF-β1) (Liton et al., 2005), and lysophosphatidic acid (Ho et al., 2018), as well as unbiased gene expression (Luna et al., 2009) and RNA sequencing studies (Youngblood et al., 2020). These studies show effective changes of TM under stretch in the processes, including calcium influx (Lakk et al., 2021), induction of autophagy (Shim et al., 2020; Shim et al., 2021), changes in actin networks (Duffy and O'Reilly, 2018), and upregulation of tyrosine phosphorylation of ECM protein in the TM (Maddala et al., 2017). Besides these studies toward understanding the response by TM subjected to CMS, studies involving global analysis of protein and lipid changes have not been performed, and an integrated analysis to investigate the underlying response mechanisms in TM subjected to mechanical stress is unavailable. Toward addressing the above lacunae in the field, we utilized a multiomics approach to better understand the physiological and pathological relevance of changes in the levels of RNA, proteins, and lipids exerted by CMS by combining mRNA sequencing, proteomics, and lipidomic analysis.

Material and methods

Materials

Reagents, antibodies, and materials

The reagents, antibodies, and materials used in this study were purchased as follows: Life Technologies (Austin, TX, United States)—TRIzol^TM LS Reagent (10296028); Zymo Research (Irvine, CA, United States)—Direct-zolTM RNA MiniPrep (R2050S); primary antibodies: Novus Biologicals LLC (Centennial, CO, United States)—anti-MKL2 (NBP1-46209) and anti-SREBP2 (NBP1-54446); Santa Cruz Biotechnology, (Dallas, TX, USA)—anti-HMGCS1 (sc-166763), anti-PAI-1 (sc-5297), anti-COL1A2 (sc-393573), and anti-GAPDH (sc-25778); rabbit anti-laminin HB-4 (anti-laminin-1 sera recognizing the total laminin including α1, β1, and γ1) (LeMosy et al., 1996) (gift from Dr. Harold Erickson, Duke University). Secondary antibodies conjugated with horseradish peroxidase (Jackson Immuno Research, West Grove, USA)—Donkey anti-mouse IgG (715-035-150) and anti-rabbit IgG (715-035-144). Cell Signaling Technology (Danvers, MA, United States) RIPA Buffer (9806S). GE Healthcare (Chicago, IL, United States)—nitrocellulose membrane (10600003).

Primary trabecular meshwork cell culture

Primary human TM (HTM) cells were cultured from TM tissue isolated from the leftover donor corneal rings after they had been used for corneal transplantation at the Indiana University Clinical Service, Indianapolis, as described previously (Keller et al., 2018). HIPPA compliance guidelines have been adhered to for the use of human tissues. The usage of donor tissues was exempt from the DHHS regulation and the IRB protocol (1911117637), approved by the Indiana University School of Medicine IRB review board. The age, race, and sex of the donors were obtained from eye banks, which provided the corneas (Supplementary Figure S1). Briefly, TM tissue extracted from the corneal ring was chopped into fine pieces and placed in a 2% gelatin-coated 6-well tissue culture plate sandwiched by a coverslip. The tissues were grown in OptiMEM (Gibco, #31985-070), containing 20% FBS and penicillin-streptomycin-glutamine solution (Gibco, #10378-016). The expanded population of HTM cells was sub-cultured after 1–2 weeks in DMEM, containing 10% FBS and characterized by the detection of dexamethasone-induced myocilin (Supplementary Figure S1) using rabbit anti-myocilin antibody (generously provided by Dr. Daniel Stamer, Duke University).

Cyclic mechanical stretch

Up to four biological replicates of HTM cells were seeded (∼40,000 cells/well) in collagen-coated plates (BioFlex plates, BF-3001C Flexcell International, Burlington, NC, United States) after counting using a hemocytometer. Once cells attain 90% confluency, cells were serum-starved for 3 h. Plates were placed on a Flexcell FX-6000 Tension System (Flexcell International, NC, USA). The cells were subjected to CMS at 0.69 Hz frequency with 15% stretching as published earlier (Luna et al., 2009; Shim et al., 2020; Youngblood et al., 2020) for 8 h for mRNA sequencing (HTM lines 1–4) and 24 h for proteomics and lipidomic analyses (HTM lines 1, 2, and 4). Control cells were cultured in the same conditions but were not subjected to CMS.

mRNA sequencing

Sample preparation for mRNA sequencing

Post CMS, the cells were washed with 1XPBS and collected and homogenized in TRI reagent. The total RNA was extracted and purified using Direct-zol^TM RNA MiniPrep kit following the manufacturer’s protocol.

KAPA mRNA HyperPrep method for mRNA sequencing

Purified total RNA was first evaluated for its quantity and quality using Agilent Bioanalyzer 2100. For RNA quality, a RIN number of 7 or higher was utilized. One hundred nanograms of total RNA were used for preparing a cDNA library that includes mRNA purification/enrichment, RNA fragmentation, cDNA synthesis, ligation of index adaptors, and RNA amplification by following the KAPA mRNA Hyper Prep Kit Technical Data Sheet, KR1352—v4.17 (Roche Diagnostics, Basel, Switzerland). Each resulting indexed library was quantified, its quality was assessed by Qubit and Agilent Bioanalyzer, and multiple libraries were pooled in equal molarity. The pooled libraries were then denatured and neutralized before loading to NovaSeq 6000 sequencer at 300 pM final concentration for 100 b paired-end sequencing (Illumina, Inc.). Approximately 30–40 M reads per library were generated. A Phred quality score (Q score) was used to measure the sequencing quality. More than 90% of the sequencing reads reached Q30 (99.9% base call accuracy).

Mapping QC and data analysis

The sequencing data were first assessed using FastQC (Babraham Bioinformatics, Cambridge, United Kingdom) for quality control. Then, all sequenced libraries were mapped to the human genome (hg38) using STAR RNA-seq aligner (v.2.5) (Dobin et al., 2013) with the following parameter: “--outSAMmapqUnique 60.” The reads distribution across the genome was assessed using bamutils (from NGSUtils v.0.5.9) (Breese and Liu, 2013). Uniquely mapped sequencing reads were assigned to hg38 refGene genes using feature Counts (from subread v.1.5.1) (Liao et al., 2014) with the following parameters: “-s 2 –p–Q 10.” Differential expression analysis was performed using edgeR (Robinson et al., 2010), (McCarthy et al., 2012). Counts were normalized to counts per million reads (CPM) for each sample. Data were examined by Multidimensional Scaling in the edgeR package (Robinson et al., 2010) to detect outliers. The data were normalized using the TMM (trimmed mean of M values) method. False discovery rates (FDR) were calculated using the Benjamini & Hochberg method (Benjamini and Hochberg, 1995) within edgeR. All raw and processed data are available via Gene expression omnibus (GEO).

The GEO accession: GSE195756

Proteomics

Sample preparation for proteomics

Cell pellets were lysed in 8 M urea (Bio-Rad Laboratories Inc., CA, United States, 161-0731), 100 mM Tris-HCl, pH 8.5 (Sigma-Aldrich, St. Louis, MO, United States, 10812846001) by sonication in 1.5 ml Micro Tubes (TPX Plastic for Sonication from Diagende Inc.) using a Bioruptor® sonication system (Diagenode Inc., New Jersey, United States, B01020001) with 30 s/30 s on/off cycles for 15 min in a water bath at 4°C. After subsequent centrifugation at 14,000 rcf for 20 min, protein concentrations were determined by Bradford protein assay (BioRad, Hercules, 5000006). A 20 µg equivalent of protein from each sample was reduced with 5 mM tris(2-carboxyethyl) phosphine hydrochloride (TCEP, Sigma-Aldrich, C4706) for 30 min at room temperature, and the resulting free cysteine thiols were alkylated with 10 mM chloroacetamide (CAA, Sigma-Aldrich, C0267) for 30 min at room temperature in the dark. Samples were diluted with 50 mM Tris HCl, pH 8.5, to a final urea concentration of 2 M for Trypsin/Lys-C based overnight protein digestion at 37°C (1:100 protease: substrate ratio, mass spectrometry grade, Promega Corporation, Madison, WI, V5072).

Peptide purification and labeling

Digestions were acidified with trifluoracetic acid (TFA, 0.5% v/v) and desalted on Sep-Pak® Vac cartridges (Waters^TM, Milford, Massachusetts, United States, WAT054955) with a wash of 1 ml 0.1% TFA followed by elution in 70% acetonitrile 0.1% formic acid (FA). Peptides were dried by speed vacuum and resuspended in 29 µL of 50 mM triethylammonium bicarbonate. Peptide concentrations were checked by Pierce Quantitative colorimetric assay (Thermo Fisher Scientific, Waltham, Massachusetts, United States, 23275). The same amount of peptide from each sample was then labeled for 2 hours at room temperature with 0.2 mg of Tandem Mass Tag reagent (Thermo Fisher Scientific, TMT™ Isobaric Label Reagent Set, 90309, lot no. VI307195B). Labeling reactions were quenched by adding 0.3% hydroxylamine (v/v) to the reaction mixtures at room temperature for 15 min. Labeled peptides were then mixed and dried by speed vacuum.

High pH basic fractionation

For high pH basic fractionation, peptides were reconstituted in 0.1% TFA and fractionated using methodology and reagents from Pierce™ High pH reversed-phase peptide fractionation kit (8 fractions, ThermoFisher, A32993).

Nano-LC-MS/MS analysis

Nano-LC-MS/MS analyses were performed on an EASY-nLC HPLC system (SCR:014993, Thermo Fisher Scientific) coupled to Orbitrap Exploris 480™ mass spectrometer (Thermo Fisher Scientific) with a FAIMS pro interface. About 1/8 of each global peptide fraction and 1/4 of each phosphopeptide fraction were loaded onto a reversed-phase EasySpray™ C18 column (2 μm, 100 Å, 75 μm × 25 cm, Thermo Scientific, ES902A) at 400 nl/min. Peptides were eluted from 4%–30% with mobile phase B (mobile phases A: 0.1% FA, water; B: 0.1% FA, 80% acetonitrile (Thermo Fisher Scientific, LS122500) over 160 min, 30%–80% B over 10 min, and dropping from 80%–10% B over the final 10 min. The mass spectrometer was operated in positive ion mode with a 4 s cycle time data-dependent acquisition method with advanced peak determination. The FAIMS CV was maintained at −50 V. Precursor scans (m/z 375-1600) were done with an orbitrap resolution of 60,000, RF lens% 40, maximum inject time 50 ms, normalized AGC target 300%, and MS2 intensity threshold of 5e4, including charges of 2–6 for fragmentation with 60 s dynamic exclusion. MS2 scans were performed with a quadrupole isolation window of 0.7 m/z, 35% HCD CE, 45,000 resolution, 200% normalized AGC target, auto maximum IT, and fixed first mass of 110 m/z.

Data analysis

The resulting RAW files were analyzed in Proteome Discover™ 2.4 (Thermo Fisher Scientific) with a Homo sapiens UniProt FASTA (last modified 021517) plus common contaminants. SEQUEST HT searches were conducted with a maximum number of two missed cleavages, precursor mass tolerance of 10 ppm, and a fragment mass tolerance of 0.02 Da. Static modifications used for the search were as follows: 1) carbamidomethylation on cysteine (C) residues, 2) TMT6plex label on N-peptide N-termini, and 3) TMT6plex label on lysine (K) residues. Dynamic modifications used for the search were oxidation of methionines and acetylation, Met-loss, or Met-loss plus acetylation of protein N-termini. The percolator false discovery rate was set to a strict setting of 0.01 and a relaxed setting of 0.050. In the consensus workflow, peptides were normalized by total peptide amount with no scaling. Co-isolation of 50% and average reporter ion S/N 10 were used as thresholds for quantification. Resulting normalized abundance values for each sample type, abundance ratio and log₂(abundance ratio) values (Supplementary Figure S2), and respective p-values (ANOVA) from Proteome Discover™ were exported to Microsoft Excel. All raw and processed data are available via ProteomeXchange.

Project accession: PXD031347

Project DOI: 10.6019/PXD031347

Lipidomics

Cellular lipid extraction and multiple reaction monitoring profiling-based MS for lipidomics

Samples represented by pellets containing around 50,000 cells were processed using the Bligh and Dyer protocol (Bligh and Dyer, 1959). Immediately after stretch, the cell lysis was achieved by adding 200 μl of ultrapure water and repeated pipetting for 1 minute. After that, 550 μl of methanol and 250 μl of chloroform were added, and the samples were incubated for 15 min at 4°C. Then, 250 μl of ultrapure water and 250 μl of chloroform were added, causing the solution to become biphasic. Phase separation was improved by centrifugation at 10,000 rpm for 5 min. The bottom phase (lipid extract) was then transferred to a clean microtube, and the extracts were dried using a SpeedVac centrifuge. Dried lipid extracts were diluted in 5 μl of chloroform and 45 μl of injection solvent (acetonitrile/methanol/ammonium acetate 300 mM 3:6.65:0.35 [v/v]) to obtain a stock solution. The stock solution was diluted ten times into injection solvent spiked with 0.1 ng/μl of EquiSPLASH Lipidomics (Avanti Polar Lipids, AL, United States #330731) for sample injection. The MRM-profiling methods and instrumentation used have been recently described in previous reports (de Lima et al., 2018) (Dipali et al., 2019) (Suárez-Trujillo et al., 2020) (Clyde-Brockway et al., 2021). In summary, data acquisition was performed using flow injection (no chromatographic separation) from 8 μl of the diluted lipid extract stock solution delivered using a micro-autosampler (G1377A) to the ESI source of an Agilent 6410 triple quadrupole mass spectrometer (Agilent Technologies, Santa Clara, CA, United States). A capillary pump was connected to the autosampler and operated at a flow rate of 7 μl/min and a pressure of 100 bar. The capillary voltage on the instrument was 5 kV, and the gas flow was 5.1 L/min at 300°C. The MS data obtained from these methods were processed using an in-house script to obtain a list of MRM transitions with their respective sum of absolute ion intensities over the acquisition time.

Cell culture media lipid extraction and high-performance liquid chromatography-mass spectrometry

For extracellular lipid analysis, control and stretched HTM cell culture media were collected (around 6 ml), and the speed vacuum was concentrated (SpeedVac DNA 130, Thermo Fisher Scientific) to generate the pellets. The lipids were extracted using the Bligh& Dyer protocol as described above (Bligh and Dyer, 1959). Lipid samples were eluted using reverse phase liquid chromatography with a Hypersil Gold C18 column (particle size 1.9 uM, 150 × 2.1 mm ID, Thermo Fisher Scientific). The column was used with the Vanquish Horizon UHPLC System (Thermo Fisher Scientific) equipped with an autosampler and binary pump combination. The HPLC system was coupled to a Q-Exactive mass spectrometer (Thermo Fisher Scientific) for high-resolution mass spectrometry analysis. The column temperature was 55°C, and the injection volume was 8 μl. A gradient of 30 min with a flow rate of 0.260 ml/min was run from 10% to 100% solvent B. Solvent A was a 50:50 ratio of acetonitrile:water +0.1% FA + 5 mM ammonium formate, and Solvent B was 88:10:2 ratio of isopropanol:acetonitrile:water +0.1% formic acid + 5 mM ammonium formate. A heated electrospray ionization (HESI) source was used as an ionization method and coupled to the Q-Exactive instrument. The conditions for the HESI were as follows: spray voltage was 4.00 kV for positive mode and 2.50 kV for negative mode. The heated capillary temperature was set to 325°C and the auxiliary heater to 275°C. The sheath gas flow rate was set to 35 units and the aux gas flow rate to 15 units. The S-lens RF level was 70. The resolution for the full scan was set to 70,000 at a mass range of 250–1,200 m/z. The automatic gain control target was 1 × 10^6, and the maximum injection time (IT) was 240 ms. The data-dependent acquisition had a resolution of 17,500 and an AGC target of 1 × 10^5 with a maximum IT of 50 ms. The isolation window was 1.0 m/z, and the collision energies were set to 20, 30, and 40 eV.

Raw files from the mass spectrometer were analyzed for lipid identification using LipidSearch software version 4.2. The following parameters were used for identification: parent and product search tolerance of 5 ppm. Filters were set at the top rank, main isomer peak, and FA priority, with an m-score threshold of 5.0 and a c-score threshold of 2.0. Quantification was set with an m/z tolerance of 5 ppm and retention time tolerance of 0.5 min. The following adducts were allowed in positive mode: [M + H]+, [M + NH4]+, [M + 2H]2+, and [M − H]−, [M − 2H]2−, [M − CH3] − in negative mode. All lipid classes were selected for search except for glycoglycerolipids and derivatized lipids.

All potential lipid matches for each sample were aligned together for the positive and negative modes. All technical replicates were aligned together. The retention time set for the alignment was 0.1 min, with top-ranked filtered and main isomer peak selected. All peaks with the same annotated lipid species were merged in the result.

Lipidomics data analysis

Statistical analysis was performed utilizing MetaboAnalyst 5.0 (http://www.metaboanalyst.ca/). Data on the relative amounts were auto-scaled to obtain a normal distribution and evaluated by principal component analysis (PCA) to find the directions that best explain the variance in the dataset, cluster analysis/heatmap to visualize the concentration values in the data table, and paired t-test to compare different lipid classes content in control and stretched TM cells/media. Informative lipids were analyzed according to class, fatty acyl residue chain length, and unsaturation level.

Data integration analysis

The RNA-seq, proteomics, and lipidomics data were integrated based on the differentially expressed genes and proteins and regulated lipid classes. The common up- and downregulated genes were matched with the proteins using gene names. For pathway analysis, comparison between differentially expressed genes and proteins based on pathways involving cytoskeleton organization (actin- and tubulin-related), ECM organization, apoptosis, oxidative stress (redox), lipid metabolism, autophagy, and proteostasis (chaperones and protein quality control) were filtered. These were fed on STRING search (https://string-db.org/) to obtain predicted interaction maps for each pathway. Lipid classes known to be involved in the regulation of the pathway mentioned above were compared with differentially expressed genes and proteins.

Western blotting

The cell lysates containing total protein were prepared using 1X RIPA buffer composed of 50 mM Tris-HCl (pH 7.2), 150 mM NaCl, 1% NP-40, 0.1% SDS, 1 mM EDTA, and 1 mM PMSF with protease and phosphatase inhibitors and then sonicated. The protein concentration was determined using Bradford Assay Reagent. A total of 25–40 μg of the protein sample was mixed with 4X Laemmli buffer and separated on 8%–10% SDS polyacrylamide gel. Following the gel run, the proteins were transferred to a 0.45 µM pore size nitrocellulose membrane. Ponceau S staining of the membrane was performed to document protein loading after transfer. Membranes were blocked in either 5% non-fat dry milk or 5% BSA in Tris-buffered saline with 0.1% Tween for 2 h, followed by respective primary antibodies overnight at 4°C (∼16 h) and then horseradish peroxidase-conjugated secondary antibodies (Jackson Immuno Research). The blots were washed with 1X TBST, and the immunoreactivity was detected using Western Lightning Plus Enhanced Chemiluminescence (ECL) Substrate (Perkin Elmer, Shelton, United States) and imaged in a ChemiDoc MP imaging system (Bio-Rad). Blots were stripped using mild stripping buffer if required to reprobe for the loading control and multiple proteins within the same molecular weight range. Data were normalized to GAPDH. Semi-quantitative analyses and fold changes were calculated from the band intensities measured using ImageJ software.

Statistical analysis

All data are presented as the mean ± standard error of the mean (SEM) of biological replicates of four independent observations for RNA-seq and three independent observations for proteomics and lipidomics mass spectrometry analysis. For immunoblotting, four independent biological replicates were used. GraphPad Prism 8 was used to generate graphs. Quantitative data were analyzed by the student’s paired t-test. p-value ≤ 0.050 was considered statistically significant.

Results

Differential expression of genes in TM subjected to CMS

Primary HTM cells from four biological replicates were subjected to CMS for 8 h to evaluate the effects of mechanical stress on the gene expression profile in TM. We focused on the 8 h window because the estimated median mRNA half-life in human cells is approximately 10 h (Yang et al., 2003). The mRNA sequencing analysis identified 14,942 transcripts. Based on the statistical significance (p ≤ 0.050) and false discovery rate (FDR) < 6%, 211 genes were identified to be differentially expressed. A systematic evaluation of the RNA-seq data is shown in Figure 1. The volcano plot (Figure 1A) with all genes significantly modulated due to CMS—red dots represent downregulated genes, blue dots represent upregulated genes, and grey dots represent the unchanged genes—compared to the controls. Among them, based on absolute log₂(fold change), 148 genes were upregulated (≥0.3 log₂FC), and 63 genes were downregulated (≤0.3 log₂FC) (Supplementary Table S1) in HTM cells subjected to CMS. Figure 1B represents the heatmap of the top 20 upregulated and downregulated genes. The PCA correlation matrix (Figure 1C, left), scree plot representing eigenvalues (proportion of variance) (Figure 1C middle), and the elbow nature of the graph indicating PCA worked well on the data and 2D PCA loading plots for the top two components are shown in Figure 1C (right). The top 20 upregulated genes are provided in Table 1.

FIGURE 1

FIGURE 1. Transcriptomics analysis of RNA-seq data. (A) Volcano plot of all genes significantly regulated by cyclic mechanical stress. (B) Heatmaps of the top 20 upregulated (left) and top 20 downregulated (right) transcripts abundances based on Log₂ FPKM (fragments per kilobase of transcript per million mapped reads). (C) PCA on the correlation matrix (left), scree plot representing the eigenvalues and the proportion of variance accounted for by the principal components (middle), and 2D PCA loading plots for the top 2 (PC1 and PC2) components (right).

TABLE 1

TABLE 1. Top 20 upregulated genes.

Interestingly, the most upregulated gene under CMS was mutS homolog 4 (MSH4) (2.54-log₂FC, 5.83-FC), which is involved in meiotic recombination (Kneitz et al., 2000) and protein ubiquitination pathway (Xu and Her, 2013). Several stress-related genes were differentially expressed: a) heat shock protein 70 family (Radons, 2016), HSPA5 (0.60-log₂FC; 1.51-FC), HSPA1B (0.92-log₂FC; 1.90-FC); b) DnaJ/heat shock protein 40 family, DNAJA1 (0.46-log₂FC; 1.37-FC), DNAJB1 (0.45-log₂FC; 1.36-FC), DNAJB9 (0.76-log₂FC; 1.70-FC); c) heat shock protein family H (Hsp110) member 1 (HSPH1) (0.73-log₂FC; 1.66-FC); d) heat shock protein family B (small) member 8 (HSPB8) (0.43-log₂FC; 1.35-FC); and e) co-chaperone for HSP90 cysteine and histidine-rich domain containing 1 (CHORDC1) (0.66-log₂FC; 1.583-FC) were upregulated significantly. Under CMS, there was significant upregulation of actin-cytoskeleton associated genes including SLIT-ROBO rho GTPase-activating protein 1 (SRGAP1) (0.39-log₂FC; 1.31-FC), spectrin-β chain, non-erythrocytic 5 (SPTBN5) (1.40-log₂FC; 2.64-FC), formin-like 2 (FMNL2) (0.38-log₂FC; 1.30-FC), and RAS like proto-oncogene A (RalA) (0.43-log₂FC; 1.348-FC). Additionally, both frizzled receptor class, FZD4 (0.65-log₂FC; 1.57-FC) and FZD8 (1.04-log₂FC; 2.06-FC) genes, were upregulated, which are involved in Wnt/β-catenin canonical signaling pathway (Gurney et al., 2012). On the contrary, Dickkopf WNT signaling pathway inhibitor 1 (DKK1) gene, a Wnt signaling inhibitor (Semënov et al., 2008), was also increased (0.54-log₂FC; 1.45-FC), indicating a homeostatic transcriptional control of Wnt signaling when TM is under mechanical stress. Further, looking at genes involved in ECM and fibrosis regulation, sphingosine kinase (SphK1) (0.83-log₂FC; 1.78-FC), a rate-limiting enzyme in the synthesis of profibrotic sphingosine-1-phosphate (S1P) (Fyrst and Saba, 2010) was augmented, and countering fibrosis increased in membrane-type metalloprotease 16 (MMP16) (0.36-log₂FC; 1.29-FC) involved in ECM degradation and Mothers against decapentaplegic homolog 7 (Smad7) (0.5-log₂FC; 1.42-FC), a negative regulator of the transforming growth factor β and fibrosis (Zhu et al., 2011).

Intriguingly, we found that genes involved in transmembrane transport were upregulated: a) calcium ion transport (solute carrier family 8 member A2) (SLC8A2) (2.1-log₂FC; 4.43-FC); b) myo-inositol transport (solute carrier family 5 member 3) (SLC5A3) (0.81-log₂FC; 1.76-FC); c) amino-acid transport (solute carrier family 38 member 2) (SLC38A2) (0.59-log₂FC; 1.50-FC); d) glycine transport (solute carrier family 6 member 9) (SLC6A9) (0.58-log₂FC; 1.50-FC); and e) acetyl-coenzyme A transporter (solute carrier family 33 member 1) (SLC33A1) (0.39-log₂FC; 1.32-FC). These proteins aid in ions, solutes, and metabolites across the biological membranes (Sahoo et al., 2014). Among the upregulated mRNA, three non-coding RNA—LOC284344 (1.32-log₂FC; 2.50-FC), LOC100379224 (0.56-log₂FC; 1.47-FC), and LOC646762 (0.47-log₂FC; 1.39-FC) —and two uncharacterized protein-coding RNA—LOC100130705 (1.23-log₂FC; 2.34-FC) (predicted protein, ATP6V1FNB, UniProtKB, A0A1B0GUX0) and C11orf94 (1.61-log₂FC; 3.06-FC) (UniProtKB - C9JXX5)—were found. The network map of all the upregulated genes in TM subjected to CMS is provided inSupplementary Figure S3.

Gene ontology enrichment analysis for upregulated genes in TM under CMS based on its molecular function and the biological process performed using ShinyGO v0.75 (Ge et al., 2020) are given in Table 2 and Supplementary Table S2, respectively. Based on molecular functions, most of the upregulated genes were associated with membrane transporter (eight genes), protein ubiquitination (nine genes), transcription repression (eight genes), stress-response (six genes), chaperone activity (five genes), and SMAD binding (four genes) (Table 2). Pathway enrichment analysis of genes based on their biological process (Supplementary Table S2) demonstrated genes that were upregulated and involved in transmembrane signaling, signal transduction, responders to stimuli such as stress, and the SMAD signaling pathway.

TABLE 2

TABLE 2. Pathway enrichment analysis for upregulated genes using ShinyGO based on molecular function.

On examining the 63 downregulated genes, we found that the coiled-coil domain containing 65 (CCDC65) gene (−1.629-log₂FC; −3.09-FC) was the most significantly downregulated. Table 3 shows the top 20 downregulated genes. We found that the only stress-response gene that was downregulated was DnaJ heat shock protein family (Hsp40) member B5 (DNAJB5) (−0.71-log₂FC; −1.64-FC). Other notable downregulated genes included the pro-apoptotic gene thioredoxin interacting protein (TXNIP) (−0.57-log₂FC; −1.49-FC), rho family interacting cell polarization regulator 2 (RIPOR2), and a negative regulator of small GTPase Rho (−1.23-log₂FC; −2.34-FC). Neural crest cells are progenitors of TM (Kaiser-Kupfer, 1989), and the expression of keratin 19 (KRT19), a neural crest marker (Lignell et al., 2017), was downregulated (−0.79-log₂FC; −1.72-FC) due to CMS. Hyaluronidase-1, an enzyme degrading hyaluronan in the ECM, was downregulated significantly (−0.85-log₂FC; −1.80-FC). Tight junction adhesion molecule claudin 11 (CLDN11) (−0.74-log₂FC; −1.67-FC) and cell adhesion molecule vascular cell adhesion molecule 1 (VCAM1) (−0.52-log₂FC; −1.44-FC) were decreased significantly indicating a potential loss of the cell–cell contact and communication in TM endothelial cells (Cerutti and Ridley, 2017).

TABLE 3

TABLE 3. Top 20 downregulated genes.

Interestingly, two lipid metabolism-related genes, including carnitine palmitoyltransferase 2 (CPT2) (−0.51-log₂FC; −1.42-FC) and short-chain specific acyl-CoA dehydrogenase (ACADS) (−0.62-log₂FC; −1.53-FC), which are involved in fatty acid beta-oxidation (Houten and Wanders, 2010), went down significantly. Pathway enrichment analysis for downregulated genes was performed using the ShinyGO v0.75 (Ge et al., 2020) tool based on their molecular function (Table 4). Based on the molecular functions, our results indicate that CMS on TM negatively regulated the cellular metabolic process at the transcription levels. A network map of downregulated genes in TM subjected to CMS is provided in Supplementary Figure S4. The original mRNA sequencing data files are available at GEO under the GEO accession: GSE195756.

TABLE 4

TABLE 4. Pathway enrichment analysis for downregulated genes using ShinyGO based on molecular function.

Differential expression of proteins in TM subjected to CMS

For proteomics analysis, TM cells were subjected to CMS for 24 h. The time provided for stretch and examination of protein changes was based on published data (Chen et al., 2016). The published study suggested a median half-life of a protein of around 9 h out of the 804 proteins traced for half-life. Nevertheless, roughly 272 proteins had a half-life ranging between 10 and 36 h (Chen et al., 2016). Therefore, we decided to examine protein changes after 24 h of CM using a TMT-based LC/MS-MS proteomics approach. In total, 4,758 proteins were detected. Further, based on FDR ≤5% and statistical significance (p ≤ 0.050), 147 proteins were significantly upregulated and 66 were significantly downregulated. Based on mean ± 2σ of log₂ ≥ 0.3 of confidence fold change limits (Supplementary Figure S2), 10 proteins were upregulated and 29 proteins were downregulated. Further, using the criteria mean ± 2σ of log₂ < 0.3–0.1, 82 proteins were upregulated and 17 proteins were downregulated in TM under CMS (Supplementary Table S3). The list of top 10 upregulated proteins is provided in Table 5.

TABLE 5

TABLE 5. Top 10 upregulated cellular proteins.

The most significantly upregulated protein was methyltransferase-like protein 13 (METTL13) (0.91-log₂FC; 1.879-FC), which catalyzes protein methylation of elongation factor 1-alpha (Jakobsson et al., 2018), indicating a potential increase in protein synthesis (Jakobsson et al., 2017; Jakobsson et al., 2018). Many proteins that bind to and/or are associated with actin cytoskeleton organization were upregulated. Among them, the significant ones were a) ataxin-3 (0.34-log₂FC; 1.26-FC), which is involved in cytoskeletal organization and cell survival (Rodrigues et al., 2010), and b) MKL/myocardin-like protein 2 (MKL2) which was upregulated (0.32-log₂FC; 1.25-FC), a coactivator that binds to transcriptional factor serum response factor (SRF), is regulated by activation of Rho GTPase, and induces genes involved in myogenic differentiation, actin cytoskeletal organization, focal adhesion assembly, and tissue fibrosis (Pipes et al., 2006; Pattabiraman and Rao, 2010; Pattabiraman et al., 2014; Xu et al., 2015). Other upregulated cytoskeleton-related and ECM proteins were cytoskeleton binding protein thymosin beta 4 (TMSB4X) (0.13-log₂FC; 1.10-FC) (Xue et al., 2014) and tropomodulin-3 (TMOD3) (0.14-log2FC; 1.11-FC) involved in F-actin stabilization (Parreno and Fowler, 2018), collagen 1A2 (COL1A2), which is a part of a large molecule type I collagen (0.14-log₂FC; 1.10-FC), pro-fibrotic mothers against decapentaplegic homolog 3 (SMAD3) (0.11-log₂FC; 1.08-FC) (Flanders, 2004), and laminin subunit gamma-1 (0.09-log₂FC; 0.09-FC), though it did not fall under our stringency criteria for consideration (log₂FC > 0.1). Mechanical stress on TM triggered upregulation of anti-apoptotic proteins such as DnaJ homolog subfamily A member 1 (DNAJA1) (0.46-log₂FC; 1.12-FC) significantly and the death domain-associated protein 6 (DAXX) (0.14-log₂FC; 1.10-FC) and tax1-binding protein 1 (TAX1BP1) (0.12-log₂FC; 1.08-FC). Additionally, nuclear receptor-binding factor 2 (NRBF2) (0.35-log₂FC; 1.27-FC) and maturation- gamma-aminobutyric acid receptor-associated protein-like 1 (GABARAPL2) (0.21-log₂FC; 1.16-FC)—proteins involved in autophagosome formation (Behrends et al., 2010; Birgisdottir et al., 2019)—were upregulated.

Interestingly, mechanical stretching of TM upregulated proteins was associated with cholesterol biosynthesis. Hydroxymethylglutaryl-CoA synthase (HMGCS1) (0.21-log₂FC; 1.16-FC), isopentenyl-diphosphate delta isomerase 1 (IDI1) (0.12-log₂FC; 1.08-FC), and squalene synthase (FDFT1) (0.33-log₂FC; 1.26-FC) were all upregulated in TM subjected to CMS, thus indicating the modulation of membrane rigidity, exocytosis, and mechanotransduction (Najafinobar et al., 2016; Lakk et al., 2021) under mechanical stress. Pathway enrichment of upregulated proteins by ShinyGO analysis (Ge et al., 2020) based on molecular function and biological process are given in Table 6 and Supplementary Table S4, respectively. The molecular functions of most of the upregulated proteins in TM subjected to CMS were strongly associated with RNA binding and/or involved in post-translational modifications. Network analysis by STRING shows a map of upregulated proteins whose biological functions are interlinked (Supplementary Figure S5).

TABLE 6

TABLE 6. Pathway enrichment analysis for upregulated cellular proteins using ShinyGO based on molecular function.

We further performed a confirmatory immunoblotting analysis of TM cells subjected to CMS for selected proteins upregulated in proteomics. Figure 2A shows significant increase in the expression of ECM proteins such as COL1A2 (n = 4, p = 0.0048) and total laminin (n = 4, p = 0.003) and MKL2 (n = 4, p = 0.049), a coactivator of SRF and a regulator of actin polymerization, in TM cells subjected to CMS (ST) compared to control (CTL). Though there was no significant change in plasminogen activator inhibitor-1 (PAI-1) expression in proteomic analysis, we found it to be significantly increasing (n = 4, p = 0.007) in immunoblotting (Figure 2A), suggesting an increase in fibrogenic response (Ghosh and Vaughan, 2012) in TM subjected to CMS. Additionally, investigating cholesterol biosynthesis, our data from immunoblotting support the proteomics data in demonstrating a significant increase in HMGCS1 (n = 4, p = 0.030) (Figure 2B). The production of HMGCS is controlled by the activation of a transcription factor sterol regulatory element-binding protein 2 (SREBP2), a master regulator of cholesterol biosynthesis (Madison, 2016). Upon activation, SREBP traffics from ER to Golgi and gets cleaved, and the cleaved N-terminal protein moves into the nucleus to initiate the transcription by binding to the sterol response element (SRE) (Madison, 2016). In TM subjected to CMS, on examining the active form of SREBP2 (nSREBP2) by immunoblotting, we identified a significant increase (n = 4, p = 0.048) (Figure 2B), indicating a direct relationship between mechanical stress on inducing cholesterol synthesis.

FIGURE 2

FIGURE 2. Immunoblot analysis for ECM and fibrosis-related proteins and proteins involved in cholesterol biosynthesis. (A) Confirmatory analysis for protein expression changes in HTM cells subjected to CMS for 24 h. COL1A2, total laminin, MKL2, and PAI-1 showed a significant upregulation at 24 h of CMS (ST) compared to unstretched control (CTL). (B) Confirmatory analysis for proteins involved in cholesterol biosynthesis in HTM cells subjected to CMS for 24 h. HMGCS1 and n-SREBP2 showed a significant upregulation at 24 h of CMS (ST) compared to unstretched control (CTL). The results were based on immunoblot analysis with subsequent densitometric analysis. GAPDH was immunoblotted as the loading control. Circles and boxes in the histograms represent sample numbers for CTL and ST, respectively. Values represent mean ± SEM, where n = 4 (biological replicates). *p ≤ 0.050 was considered significant.

Among the downregulated proteins, the coiled-coil domain-containing proteins, namely, CCDC113 (−0.95-log₂FC; 0.52-FC) and CCDC170 (−0.91-log₂FC; 0.53-FC), a component of centriole satellite and protein involved in microtubule organization and stabilization, respectively, were the most downregulated (Table 7). Intriguingly, CCDC113 was downregulated at mRNA and protein levels.

TABLE 7

TABLE 7. Top 20 downregulated cellular proteins.

Collagen alpha-2(IV) chain (COL4A2), one of the most abundant type IV collagens involved in the formation of the basement membrane, was downregulated (−0.15-log₂FC; 0.90-FC). Actin binding Nesprin-2 (−0.77-log₂FC; 0.59-FC) (Dawe et al., 2009) and the dedicator of cytokinesis protein 2 (DOCK2) (−0.76-log2FC; 0.59-FC), which is involved in activating Rho family of GTPases (Laurin and Côté, 2014), were downregulated significantly. Pro-apoptotic protein, CASP8-associated protein 2 (CASP8AP2), was downregulated (−0.77-log₂FC; 0.59-FC). Though proteins involved in autophagosome formation and maturation were upregulated, we found that autophagy-related protein 13 (ATG13), essential for autophagy induction and autophagosome formation (Behrends et al., 2010), was significantly downregulated (−0.39-log₂FC; 0.77-FC). Thus, indicating a potential activation and a balancing act on the autophagic mechanisms under mechanical stress in the TM cells. Pathway enrichment analysis performed using ShinyGO v0.75 (Ge et al., 2020) tool for downregulated proteins based on molecular function (Table 8) showed several proteins associated with biological quality, cellular localization, and stress response. Network analysis of downregulated proteins is shown in Supplementary Figure S6.

TABLE 8

TABLE 8. Pathway enrichment analysis for downregulated cellular proteins using ShinyGO based on molecular function.

Lipid changes in TM subjected to CMS

Lipidomics was performed in control and stretched HTM cells from cellular lipids and conditioned media.

The PCA for the following lipid classes yielded two distinct clusters between control and stretched HTM cells in total phospholipids (Supplementary Figure S7A), total ceramide (Cer) (Supplementary Figure S7B), total cholesteryl ester (Supplementary Figure S7C), and total diacylglycerol (DG) (Supplementary Figure S7D). The first and second principal components, PC1 and PC2, explained >50% of the variance, indicating that total phospholipids, Cer, cholesteryl esters, and DG were highly different between control and stretched HTM cells. The PCA for total TAG showed no distinct cluster between control and stretched HTM cells (Supplementary Figure S7E).

In the lipids from conditioned media, the PCA revealed two distinct clusters of total phospholipids in control and stretched HTM conditioned media (Supplementary Figure S7F). The first two principal components explained >50% of the variance, indicating that total phospholipids differed highly between control and stretched HTM conditioned media (Supplementary Figure S7F). However, in PCA for the lysophospholipids, there was no distinct grouping in control and stretched HTM conditioned media (Supplementary Figure S7G).

Further analysis of each class of lipids is shown as follows:

(1) Total 603 phospholipid classes, including phosphatidylcholine (PC), phosphatidylinositol (PI), phosphatidylethanolamine (PE), phosphatidylserine (PS), phosphatidylglycerol (PG), and sphingomyelin (SM), were identified. Cluster heatmap of the top 25 phospholipids showed that most were increased in stretched HTM cells compared to the control (Figure 3A). Paired t-test showed that there was a significant increase in total phospholipids (p = 0.0001) in stretched HTM cells compared to control (Figure 3F).

(2) Among all phospholipids, 101 PC classes were identified. PC is an abundant component in the cell membrane (Kanno et al., 2007). Therefore, we next measured the change of total PC level in phospholipids of control and stretched HTM cells. The paired t-test showed that total PC was significantly increased (p = 0.044) in stretched HTM cells compared to control (Figure 3G). In addition to PC, we also compared the changes of PI, PE, PG, and PS in control and stretched HTM cells. A total of 133 PI classes were identified. The paired t-test showed that total PI was significantly increased (p = 0.0001) in stretched HTM cells (Supplementary Figure S8A). A total of 141 PE classes were identified, and the paired t-test showed that total PE was significantly increased (p = 0.0002) in stretched HTM cells (Supplementary Figure S8B). A total of 134 PG classes were identified in lipidomics analysis, and the paired t-test showed that total PG was significantly increased (p = 0.007) in stretched HTM cells (Supplementary Figure S8C). A total of 131 PS classes were identified, and the paired t-test showed no difference in the total PS level between control and stretched HTM cells (Supplementary Figure S8D).

(3) Since sphingosine-1-phosphate (S1P) is known to regulate IOP (Mettu et al., 2004; Stamer et al., 2009; Sumida and Stamer, 2010, 2011; Honjo et al., 2018; Ho et al., 2020), we next checked the phospholipid sphingomyelin (SM) level, a source of S1P synthesis, in control and stretched HTM cells. Among all phospholipids, 27 SM classes were identified. We also found an increase in SM with no significance in stretched HTM cells compared to control (Figure 3H).

(4) In addition to SM, we also checked the levels of Cer, another lipid involved in S1P synthesis in the control and stretched HTM cells. A total of 80 Cer classes were identified. The cluster heatmap showed the top 25 upregulated Cer in stretched HTM cells compared to control (Figure 3B), and the paired t-test for the total Cer resulted in significantly increased Cer levels (p = 0.012) in stretched HTM cells (Figure 3I). Such increases in SM and Cer indicate a potential increase in the biogenesis of S1P in stretched HTM cells.

(5) A total of 25 cholesteryl ester classes were identified. The cluster heatmap showed that cholesteryl esters were decreased in two stretched HTM cell lines and increased in one line (Figure 3C). The paired t-test showed that the total cholesteryl ester significantly decreased (p = 0.02) in stretched HTM cells compared to the control (Figure 3J).

(6) A total of 24 diacylglycerol (DG) classes were identified. Heatmap from cluster analysis showed that total DG was increased in stretched HTM cell samples (Figure 3D), and in accordance with this observation, the paired t-test also showed that total DG was significantly increased (p = 0.025) in stretched HTM cells compared to the control (Figure 3K).

(7) A total of 328 triglyceride (TAG) classes were identified. The cluster heatmap of the top 25 TAGs showed that some were increased, and some were decreased in stretched HTM cells (Figure 3E). The paired t-test showed that total TAG levels were significantly decreased (p = 0.048) in stretched HTM cells compared to the controls (Figure 3L).

FIGURE 3

FIGURE 3. Cellular lipid analysis. (A) Heatmap of the total cellular phospholipids in control and stretched HTM cells. (B) Heatmap of the total cellular ceramide (Cer) in control and stretched HTM cells. (C) Heatmap of the total cellular cholesteryl ester in control and stretched HTM cells. (D) Heatmap of the total cellular diacylglycerol (DG) in control and stretched HTM cells. (E) Heatmap of the total cellular triglyceride (TAG) in control and stretched HTM cells. All heatmaps are presented by control (left 3 red columns) and stretched HTM cells (right green columns), total n = 3 pairs (biological replicates). Analysis parameters utilized the Ward clustering algorithm, Euclidean distance measure, and autoscaling based upon samples (performed using MetaboAnlayst 5.0). Comparison between control and stretched HTM cells. (F) Total cellular phospholipids were significantly increased in stretched HTM cells. (G) Total cellular phosphatidylcholine (PC) was significantly increased in stretched HTM cells. (H) Total cellular sphingomyelin (SM) was increased in stretched HTM cells but was nonsignificant. (I) Total cellular Cer was significantly increased in stretched HTM cells. (J) Total cellular cholesteryl ester was significantly decreased in stretched HTM cells. (K) Total cellular DG was significantly increased in stretched HTM cells. (L) Total cellular TAG was significantly decreased in stretched HTM cells. Values represent the mean ± SEM, where n = 3 (biological replicates). *p ≤ 0.050 was considered significant.

In addition to the analysis of cellular lipids, we also evaluated lipid changes in the conditioned media obtained from control and stretched HTM. Since lysophosphatidic acid (Liu et al., 2018) is an important extracellular lipid growth factor involved in the regulation of TM cell function and IOP (Mettu et al., 2004; Honjo et al., 2018; Ho et al., 2020), we analyzed the changes in total phospholipids (PC, PI, PE, PS, PG, and SM) and lysophospholipids (lysophosphatidylcholine, lysophosphatidylethanolamine, and lysophosphatidylglycerol) level in conditioned media, which are the resources of extracellular LPA. Lipidomics analysis identified 384 phospholipid and 43 lysophospholipid classes in control and stretched HTM conditioned media. Cluster heatmap of the top 25 phospholipids showed that most were increased in stretched HTM conditioned media compared to control (Figure 4A), and paired t-test showed a small increase in total phospholipids level in stretched HTM conditioned media (Figure 4C). Cluster heatmap of the top 25 lysophospholipids showed that they were decreased in stretched HTM conditioned media compared to the control (Figure 4B), and paired t-test also presented a significant decrease in total lysophospholipids (p = 0.029) in stretched HTM conditioned media (Figure 4D).

FIGURE 4

FIGURE 4. Secretory lipid analysis. (A) Heatmap of the total media phospholipids in control and stretched HTM cells. (B) Heatmap of the total media lysophospholipids in control and stretched HTM cells. Comparison between control and stretched HTM cell culture media. (C) Total media phospholipids were increased in stretched HTM cell culture media but were nonsignificant. (D) Total media lysophospholipids were significantly decreased in HTM cell culture media. Values represent the mean ± SEM, where n = 3 (biological replicates). *p ≤ 0.050 was considered significant.

The lipidomics analysis revealed interesting changes in lipid classes in TM experiencing CMS, indicating an active remodeling of lipid components in the TM. The intracellular and extracellular lipid changes, such as Cer, SM, and lysophospholipids, can promote S1P and LPA biosynthesis, further regulating HTM cell biomechanics and TM stiffness.

Discussion

To date, two studies have utilized the “omics” approach to identify changes in genes and pathways in TM subjected to CMS. Luna et al. (2009) used microarray analysis to find the differentially expressed genes, and a more recent study (Youngblood et al., 2020) utilized RNA sequencing to identify mRNAs, micro RNAs, and long non-coding RNAs. Interestingly, a comparative analysis of our study where we performed mRNA-seq with the study of Youngblood et al. (2020) that utilized total RNA-seq analysis identified a similar upregulation of genes such as MARCKS like 1 (MARCKSL1), leucine-rich adaptor protein 1 like (LURAPIL), SLC38A2, a regulator of chromosome condensation 1 (RCC1), SLC5A3, CEP295 N-terminal like (CEP295NL), and downregulation of thioredoxin interacting protein (TXNIP).

To better understand and integrate the stress-responsive mechanisms utilized by TM due to mechanical stress, in this exploratory study, we have performed an unbiased multiomics analysis to identify links between changes in mRNA, proteins, and lipids using RNA sequencing, proteomics, and lipidomics on TM cells subjected to CMS. Our study reveals that mechanical stress induced a potential protein quality control response to counter oxidative stress and induction of cellular contractility and cell adhesive interactions combined with pro-fibrogenic signaling via increased lipid signaling and metabolic response. In addition, we find that the cells turn on various homeostatic control signaling machinery to counteract such stress.

TM counters CMS-induced oxidative stress and apoptosis

Based on the RNA-seq analysis, our study found 211 differentially expressed genes. Many genes found to be differentially expressed in our study belonged to stress-response, protein ubiquitination, transcription repression, and chaperone activity. Understanding the role of oxidative stress and reactive oxygen species is of significant interest in the pathogenesis of ocular hypertension and glaucoma. Multiple lines of investigations pointed out the relationship between oxidative stress-related damage to TM and elevated IOP (Saccà and Izzotti, 2008; Saccà et al., 2016; Tanito et al., 2016; Wu et al., 2022). First, we propose that CMS induces oxidative stress on the TM and activates transcription and translation to negate oxidative stress. Our RNA-seq analysis revealed that downstream response genes under oxidative stress, including HSPA1B and HES1, were upregulated, whereas CCL2, TXNIP, and MEKK6/MAP3K6 were downregulated. An earlier study also indicated the induction of oxidative stress in TM upon stretching (Vittal et al., 2005). Our proteomics data support the potential negative regulation of oxidative stress by the downregulation of the Kelch-like ECH-associated protein 1 (KEAP1). Mechanistically, nuclear factor erythroid 2- (NFE2-) related factor 2 (NRF2) is a regulator of cellular resistance to oxidative stress, and under basal conditions, NRF2 is suppressed by binding to the repressor protein KEAP1. During oxidative stress, NRF2 becomes active by dissociating from KEAP1 (Deshmukh et al., 2017). In addition, we found upregulation of casein kinase 2α (CSNK2A2), which is required for the phosphorylation of NRF2, which aids in the translocation of NRF2 to the nucleus and activates several cytoprotective genes (Apopa et al., 2008). Thus, the negative regulation of KEAP1 in TM indicates a potential activation of NRF2 because of the likely oxidative stress damage of TM subjected to CMS. Our study also points out that the loss of counter-mechanisms to oxidative stress in TM can be an important step in ocular hypertension and glaucoma pathogenesis. Interestingly, TXNIP and MAP3K6 are pro-apoptotic genes, and the decrease in their transcription aid in protecting the cells from undergoing cell death (Takeda et al., 2007; Li J. et al., 2017). We observed that proteins involved in inhibition of apoptosis were increased in TM undergoing CMS. Upregulation of DnaJ homolog subfamily A member 1 (DNAJA1), a co-chaperone for HSPA8/Hsc70 prevents BAX translocation to mitochondria, can inhibit apoptosis (Gotoh et al., 2004); death domain-associated protein 6 (DAXX), a transcription repressor is a negative regulator of p53-mediated apoptosis (Zhao et al., 2004; McDonough et al., 2009); and tax1-binding protein 1 (TAX1BP1) is an inhibitor of TNF-induced apoptosis (Shembade et al., 2007). Incidentally, DNAJA1 was upregulated at transcript and protein levels, signifying its effects on regulating chaperone function and inhibiting apoptosis. In contrast, CASP8AP2, which is required for FAS-mediated apoptosis, was downregulated. Even in our mRNA sequencing data, we found that TXNIP, which is known to be involved in apoptosis (Li J. et al., 2017), was going down. We believe these apoptotic inhibitory aspects in TM cells could be a cell survival mechanism to overcome mechanical stress. The potential role of oxidative stress-induced DNA damage can upregulate the expression of DNA repair enzymes (Bridge et al., 2014). In our RNA-seq analyses, the top upregulated gene was MSH4, which is involved in meiotic recombination (Snowden et al., 2004) during DNA repair (Tock and Henderson, 2018), indicating potential protection against oxidative DNA damage. MSH4 also interacts with von Hippel-Lindau binding protein 1 (VBP1), which binds to VHL ubiquitin E3 ligase that regulates cellular response to oxidative stress, functioning in the protein degradation pathway (Xu and Her, 2013).

TM overcomes mechanical stress by upregulating transcriptional and translational control of autophagy and protein quality control

Autophagy is an adaptive response mechanism for the survival of cells experiencing stress. Autophagy induction in TM subjected to CMS was previously reported (Shim et al., 2020). In our study, we found the upregulation of HSPB8 mRNA, which promotes autophagosome formation (Li X. C. et al., 2017), and an increase in NRBF2 and GABARAPL2 proteins, which are essential for the formation and maturation of autophagosome, respectively (Behrends et al., 2010). Hence, the increase in proteins related to autophagy can potentially aid the TM cells to combat stress (Rao et al., 2021) and modulate the fibrogenic pathway in TM (Nettesheim et al., 2019), and an increase in autophagic process/response serves as an anti-apoptotic signal (Liu et al., 2015). Protein folding mechanisms can be compromised in cells under stress due to mechanical insults (Höhfeld et al., 2021). Such misfolded proteins are cleared by the process of ubiquitination. Our study showed the induction of various HSP genes and TNF receptor-associated factor 1(TRAF1) involved in the regulation of proteostasis. In addition, we found upregulation of TAB2 and downregulation of TNFAIP3-interacting protein 2 (TNIP2), which are involved in protein quality control mechanism via the NFκB pathway (Kanayama et al., 2004; Banks et al., 2016). The NFκB pathway plays an important role in protein clearance by regulating autophagy (Nivon et al., 2016) and is known to regulate TGF-β2-mediated ocular hypertension (Hernandez et al., 2020).

Lipids are involved in the autophagy pathway (Soto-Avellaneda and Morrison, 2020). Phosphoinositides, a class of phospholipids derived from PI, play a major role in the regulation of autophagy. They can control the pathway that directly activates or deactivates mTORC1, which controls the balance between growth and autophagy (Soto-Avellaneda and Morrison, 2020). Generally, mTORC1 is active in the presence of nutrients, driving the growth and anabolic process while suppressing autophagy. In the absence/depletion of nutrients, mTORC1 is inactivated, and autophagy is promoted. We speculate that CMS on TM promotes nutrient utilization to meet energy expenditure. Since the cells are stretched in the conditioned media for 24 h, nutrients are used up from the conditioned media and leading to the depletion of some nutrients. In addition to PI, in the absence/depletion of nutrients, PC and its product phosphatidic acid (PA), as well as DG, can promote autophagy by mediating signaling and membrane remodeling to support autophagy activation. The Cer and its product S1P can contribute to mTORC1 inhibition and promote autophagy (Dall'Armi et al., 2013). In our lipidomics analysis, we found significantly increased PI, PC, DG, and Cer in stretched HTM cells implying the balance between growth and autophagy under CMS. On the contrary, we observed increased SREBP2 activation and lipid biogenesis in stretched HTM cells, which can be downstream of mTORC1 signaling (Lewis et al., 2011; Peterson et al., 2011). All these lipid changes imply that the HTM cell is achieving a balance between growth and autophagy under mechanical stress.

Induction of actin-based contraction and ECM buildup by the interplay of the gene, protein, and lipid signaling in TM undergoing CMS

Cells respond to various signals, especially mechanical stimuli, through changes in genes and proteins associated with actin-cytoskeleton and ECM interactions (Martino et al., 2018). Not surprisingly, this study identified the upregulation of actin-cytoskeleton-associated genes, including SRPGAP1, SPTBN5, FMNL2, Ral A, and CHORDC1, and the downregulation of RIPOR2/RIPR2 under CMS. Interestingly, CHORDC1/CHP1, a negative regulator of Rho GTPase signaling by inhibiting Rho kinase (Ferretti et al., 2010), acts as a circuit breaker of Rho-Rho kinase signaling, indicating the potential homeostatic effect on Rho GTPase-mediated mechanotransduction. Looking at the regulation of protein, we found upregulation of cytoskeletal-associated proteins such as a) ataxin-3 involved in the maintenance of actin, b) MKL2, a transcriptional coregulator of SRF, c) thymosin beta 4, which is involved in sequestering of g-actin, and d) tropomodulin-3 that is involved in actin stabilization and downregulation of DOCK2, which can activate Rac GTPases (Sanui et al., 2003), and Nesprin-2, which modulates mechanotransduction between the cytoplasm and nuclear lamina (Stewart-Hutchinson et al., 2008). Of importance is the upregulation of MKL2, which is involved in the regulation of actin cytoskeletal organization and fibrosis in HTM cells through the transcriptional control of SRF-mediated actin polymerization (Pattabiraman and Rao, 2010; Pattabiraman et al., 2014). Thymosin beta 4 and tropomodulin-3 increase in TM under experimental IOP elevation (Gonzalez et al., 2000; Borrás, 2003; Vittitow and Borrás, 2004). It is interesting to note that the most upregulated protein is the methyltransferase-like protein 13, catalyzing the methylation of elongation factor 1-alpha (eEF1A) (Jakobsson et al., 2018). eEF1A has been reported to be colocalized with F-actin and associated with changes in the cytoskeletal organization (Liu et al., 1996). In gene and protein expression analysis, the most downregulated are coiled-coil domain-containing proteins, which are involved in microtubule binding, organization, and stabilization (Pal et al., 2010; Kumari and Panda, 2018). Microtubules are reported to self-repair in response to mechanical stress (Schaedel et al., 2015). It has been reported that microtubule disruption leads to a cellular contraction in human TM cells in response to stimuli (Gills et al., 1998). We speculate that the downregulation of proteins related to microtubule organization can potentially disrupt microtubule assembly or act as a response mechanism to mechanical stress on TM. The major lipid growth factors, including LPA and S1P, are known to regulate actin cytoskeleton and barrier function in endothelial cells and function in TM by binding their cognate receptors to regulate IOP (Mettu et al., 2004; Sumida and Stamer, 2011; Rao, 2014). The S1P is synthesized inside the cell, which is generated by phosphorylation of sphingosine through SphK (Fyrst and Saba, 2010). Sphingosine is derived by the hydrolysis of Cer during the sequential degradation of plasma membrane glycosphingolipids and SM (Proia and Hla, 2015). We found that SM and Cer increased in stretched HTM cells, and our mRNA sequencing analysis also found an increase in the expression of SphK in stretched HTM cells. These results indicate that, under CMS, HTM cells generate more substrates in the S1P biosynthesis pathway. Combined with increased SphK, more S1P will be produced inside the cells and exported into the extracellular environment. Outside the cell, S1P can bind to the S1P receptors (S1PR1-5) located on the cell membrane and cause a series of downstream effects (Proia and Hla, 2015). A speculative sequence of events in TM due to altered S1P and LPA identified in this study and previous studies are shown in Supplementary Figure S9. Additionally, in the extracellular space, phospholipids are converted to lysophospholipids by phospholipase A1/2, which is converted to LPA by the enzyme autotaxin (ATX) (Perrakis and Moolenaar, 2014). LPA can bind to their cognate LPA receptors, activate a series of downstream signaling cascades, and regulate HTM cell biomechanics (Supplementary Figure S9). Our lipidomics analysis demonstrated no significant difference in total phospholipid levels in conditioned media between control and stretched HTM. However, the total lysophospholipid levels were significantly decreased in stretched HTM conditioned media. We predict this is due to the dynamic conversion of lysophospholipid to LPA as previously reported that ATX and its activity are increased in stretched TM and the conditioned media (Iyer et al., 2012; Ho et al., 2018).

Mechanical stress positively regulates growth factor-induced fibrogenic activation in TM (Liton et al., 2005; Ashwinbalaji et al., 2020). Like earlier studies (Vittal et al., 2005; Luna et al., 2009), our analysis found regulation in the genes related to enzymes involved in ECM regulation—MMP16 (upregulated) and hyaluronidase-1(downregulated). The differences in the degree of changes in these studies could result from differences in experimental technique (microarray vs. mRNA-seq) and altered response by the donor tissue to prepare the HTM cells. We found a potential increase in fibrogenic activation with an increase in SMAD3 protein known to signal via the canonical TGF-β signaling pathway and PAI-1 levels in TM subjected to CMS (Flanders, 2004; Pattabiraman et al., 2014) and an increase in COL1A2 and total laminin proteins. Interestingly, the COL4A2 protein, one of the major structural components of the basement membrane, was downregulated, correlating with earlier microarray analysis (Luna et al., 2009).

The cellular phospholipids and cholesterol composition play an important role in the cell-ECM interactions (Márquez et al., 2008; Horn and Jaiswal, 2019), and their dysregulation can result in pathology as observed in lung fibrosis (Burgy et al., 2022) and renal pathology (Abrass, 2004). A previous study found that patients with glaucoma have higher cholesterol levels (Posch-Pertl et al., 2022), and patients with hyperlipidemia will have an increased risk of getting glaucoma (Wang and Bao, 2019). In addition, cholesterol-lowering statin, an inhibitor of HMG-CoA reductase, can decrease the risk of POAG and be an effective therapy for POAG (McGwin et al., 2004; Stein et al., 2012; Wu et al., 2020). Statins are proposed to decrease IOP by inhibiting the Rho and Rho-kinase activity in TM cells (Song et al., 2005; Von Zee et al., 2009; Cong et al., 2018). This inhibition can alter TM cell morphology and decrease actomyosin contractile activity and ECM assembly. The current understanding of how cholesterol changes in stretched HTM cells is still unclear. Due to poor ionization and the low amount of free cholesterol in the cell, detection of free cholesterol is difficult, and free cholesterol is usually detected by GC-MS (Edwards et al., 2011). We evaluated the cholesteryl ester level changes in stretched HTM cells. We found that cholesteryl ester was significantly decreased in stretched HTM cells. Interestingly, our proteomics and immunoblotting analysis found that HMGCS, a critical enzyme in the cholesterol biosynthesis pathway, was significantly increased in HTM cells. Moreover, immunoblotting of whole cell lysate post mechanical stretch, we found activation of SREBP2, a transcriptional factor promoting cholesterol biosynthesis (Madison, 2016). The decreased cholesteryl ester levels observed in stretched HTM cells are speculated to result from less stockpile of cholesterol inside the cell as it is used for rapid membrane changes and the release of free cholesterol from cholesteryl esters (van Meer, 2001; Walther and Farese, 2012). Low cholesterol levels act as a stimulus for SREBP activation (Madison, 2016), which can turn on more cholesterol synthesis. These results put together indicate an increased cholesterol biosynthesis in stretched HTM cells.

DG is an intracellular lipid second messenger involved in multiple processes and pathways. Cellular DG can be synthesized from TAG and phospholipid (Carrasco and Mérida, 2007; Fagone and Jackowski, 2009; Coleman and Mashek, 2011). Cellular DG can directly interact with membrane-bound PKC and lead to its activation. The role of PKC in TM has been studied before. The DG analog diC8 can stimulate membrane-bound PKC and further cause increased TM cell contraction (Thieme et al., 1999). PKC activators such as phorbol-12-myristate 13-acetate (PMA) and phorbol-12,13-dibutyrate (PDBu) increased the formation of actin stress fibers and focal adhesions and myosin light chain (MLC) phosphorylation in TM cells (Khurana et al., 2003). In contrast, PKC inhibitors such as H7 or chelerythrine can induce relaxation in TM (Wiederholt et al., 1998; Thieme et al., 1999). These inhibitors have been used successfully to lower IOP in the animal model (Tian et al., 1999). In our lipidomics analysis, we found that both DG and its substrate—phospholipid—are increased in stretched HTM cells, indicating that under CMS, more DG generated can result in the activation of PKC and further cause increased TM cell contraction.

Interplay between mRNA, protein, and lipid pathway changes in response to CMS

Though there is no direct method to derive the interplay between mRNA, protein, and lipid pathways due to CMS, to create a conceptual interactome, we utilized the STRING database and analyzed the interactions and interaction strength between molecules either upregulated or downregulated within pathways. Interestingly, we found strong and significant interactions in the genes and proteins upregulated in association with proteostasis (Figure 5), actin cytoskeleton and ECM-associated pathway (Figure 6), and apoptosis pathway (Supplementary Figure S10) and downregulated in association with lipid metabolism (Supplementary Figure S11). Other network analysis of the downregulated oxidative stress-related genes and proteins (Supplementary Figure S12), downregulated apoptosis-associated proteins (Supplementary Figure S13), upregulation and downregulation of autophagy-related proteins (Supplementary Figures S14, S15), and upregulated lipid-related genes and proteins (Supplementary Figure S16) showed minimal interactors.

FIGURE 5

FIGURE 5. STRING network analysis of upregulated proteostasis-associated genes and proteins. Colored lines between the nodes indicate their basis for interconnection. Table here gives a list of inclusive genes and proteins. Overall, the protein–protein interaction (PPI) enrichment p-value is shown for connected nodes and is considered significant if p < 0.05 with individual PPI interaction significance not shown.

FIGURE 6

FIGURE 6. STRING network analysis of actin cytoskeleton and ECM-associated genes and proteins. (A) Network analysis of all upregulated actin cytoskeleton and ECM-associated genes and proteins. (B) Network analysis of all downregulated actin cytoskeleton and ECM-associated genes and proteins. Table here gives a list of inclusive genes and proteins. Colored lines between the nodes indicate their basis for interconnection. Overall PPI enrichment p-value is shown for connected nodes and is considered significant if p < 0.05 with individual PPI interaction significance not shown.

One limitation to the study is the choice of a low cut-off for significance, which was not arbitrary. We believe that even a small fold change in proteins can yield functional significance. In signal transduction, stoichiometry is extremely important because enzymes, repressors, and activators, among others, directly induce signal transduction cascades (Lee and Yaffe, 2016). Also, because we have investigated the changes in RNA for 8 h and proteins and lipids for 24 h, we believe that such changes can also include the consequence of a homeostatic mechanism. Another limitation is only performing proteomics and lipidomics analysis on the cellular proteins and lipids, respectively, because the secreted biomolecules also play a significant role in the outside-in signaling. Further studies on the secreted proteins and the investigation into cell-based and animal model-based studies can reveal the significance of such interaction in TM physiology and their potential pathological relevance in ocular hypertension.

Disease relevance of RNA, proteins, and lipids modulated in TM under CMS

Some of the genes, proteins, and lipid pathways we report here are associated with systemic hypertension, elevated IOP, and POAG pathogenesis. Comparing the GWAS catalog (https://www.ebi.ac.uk/gwas/home) under open-angle glaucoma trait with the differentially regulated genes and proteins found in our analysis, we found nearly 15 genes and protein changes. The upregulated genes and proteins include FMNL2, DNAJB1, DNAJA1, DNAJB1, KPNA3, NAP1L4, DNAJA1, OXNAD1, NAA50, and GPALPP1. The downregulated genes and proteins include DNAJB5, NAT6, SNAPC2, RNASE4, and PMM2. Additionally, aldehyde dehydrogenase 3 family member A1 (ALDH3A1) mRNA downregulated in our study has been reported to decrease in TM from POAG patients (Liu et al., 2013). Aldehyde dehydrogenases play an important role in lipid peroxidation and in the maintenance of redox homeostasis (Choudhary et al., 2005). Interestingly, SLC38A2, a calcium transporter that shows a significant increase in both mRNA and protein and splicing factor, and the suppressor of white-apricot homolog (SFSWAP), which is increased in protein expression analysis, are implicated in systemic hypertension (Iwamoto et al., 2004; Zhang et al., 2010; de Las Fuentes et al., 2013). We found that SLC2A12 and VCAM1 genes decrease significantly under CMS, and, interestingly, SLC2A12 has been shown to influence the POAG risk (Gharahkhani et al., 2021). As SLC2A12 is known to code for a glucose transporter (Wood et al., 2003), further studies on the role of the SLC2A12 and the correlation with metabolic regulations in TM can provide deeper insights into the relationship between hyperglycemia and POAG (Zhou et al., 2014; Hymowitz et al., 2016). The cell adhesion protein VCAM1 has been predicted to be a marker for low flow regions in the TM (Staverosky et al., 2020), thus implicating their significant function in the regulation of AH humor outflow via the TM outflow pathway. Thus, further investigations on the role of VCAM1 will help better understand the role of cell adhesion proteins and the pathogenesis of ocular hypertension and POAG.

In general, fluctuations in arterial pressure induce changes in IOP (Bakke et al., 2009). Population-based studies have shown a positive association between blood pressure and IOP (McLeod et al., 1990; Chen and Lai, 2005; Wang et al., 2009). Significant downregulation of the phosphomannomutase 2 (PMM2) level was observed in our proteomic study. PMM2 has been shown by multiple GWAS studies to be a novel locus for POAG (Chen et al., 2014; Wiggs and Pasquale, 2017). PMM2 is known to be involved in glycosylation, enabling posttranslational modifications (addition of N-linked oligosaccharides) on proteins (López-Gálvez et al., 2020). Among the upregulated genes, we found induction of CAVIN4, which is involved in the formation of caveolae (Ogata et al., 2014). Interestingly, genetic polymorphism in caveolin-1 and -2 is associated with POAG (Thorleifsson et al., 2010), and caveolin knockout mice had elevated IOP (Elliott et al., 2016). Additionally, studies have shown oxidative DNA damage to TM cells in patients with POAG (Sacca et al., 2005), implicating oxidative stress as a cause of TM damage in glaucoma (Wu et al., 2022). Finally, the induction of actin cytoskeleton-based contractility and activation of pro-fibrotic pathways in TM have been implicated in glaucoma (Zhavoronkov et al., 2016), and our study identifies mechanical stretch-induced activation of the fibrogenic process, indicating CMS as an appropriate model to study pathological signaling.

Conclusion

Based on multiomics analysis of the different pathways activated and suppressed in TM due to mechanical stress (Figure 7), we propose that, in TM, CMS activated cholesterol biosynthesis and lipid messengers contribute to the maintenance of cell membrane fluidity and cell permeability as well as in the regulation of actin-based cell tension and ECM based cell stiffness. Additionally, TM cells turn on the protein quality control machinery, autophagy induction, oxidative damage response, and anti-apoptotic and cell survival signaling to maintain homeostasis against the mechanical insult.

FIGURE 7

FIGURE 7. Graphical abstract derived from the multiomics analysis representing the modulation in signaling mechanisms activated and suppressed in TM due to the mechanical (created using BioRender.com).

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

Conceptualization: PP. Methodology, cell cultures, immunoblotting: AS and TW. Formal analysis: AS, TW, and PP. Investigation: AS, TW, and PP. Proteome analysis sample preparation: AS and AW. Proteomics LC-MS/MS/data analysis/data evaluation/verification and methods: AW and AM. Lipidomics data collection and data analysis: FH, SB, and TW. RNA-seq data extraction and analysis integration, including volcano plot, heat maps, and PCA plots: RS. Write-up and data curation: AS, TW, RS, and PP. Writing—original draft preparation: AS, TW, and PP. Figure preparation: AS, TW, RS, and PP. Writing: review and editing, AS, TW, and PP. Visualization: AS, TW, and PP: Supervision: PP: Project administration: PP: Funding acquisition: PP. All authors have read and approved the final manuscript.

Funding

This project was supported by the National Institutes of Health/National Eye Institute (R01EY029320) and a grant from Research to Prevent Blindness to IU. The mass spectrometry work performed in this work was done by the Indiana University School of Medicine Center for Proteome Analysis. Acquisition of the IUSM Center for Proteome Analysis instrumentation used for this project was provided by the Indiana University Precision Health Initiative. The proteomics work was supported, in part, with support from the Indiana Clinical and Translational Sciences Institute, funded in part by Award Number UL1TR002529 from the National Institutes of Health, National Center for Advancing Translational Sciences, Clinical and Translational Sciences Award, and the Cancer Center Support Grant for the IU Simon Comprehensive Cancer Center (Award Number P30CA082709) from the National Cancer Institute.

Acknowledgments

We would like to acknowledge Dr. Emma Doud from the Center for Proteome Analysis, IUSM, for helping with LC-MS/MS, Data Analysis, Data evaluation, verification, and methods write-up. We thank Dr. Christina Ferreira from the Bindley Bioscience Center, Purdue University, for the discussions, planning, and help and training with the execution of the MRM-based Lipidomics analysis, data analysis, and manuscript writing. We thank Jennifer Arcuri from SB’s lab for the discussions on the methodology for lipid extraction. We would like to acknowledge the RNA sequencing facility at the Center for Medical Genomics, Indiana University School of Medicine, which is partially supported by the Indiana University Grand Challenges Precision Health Initiative, metabolite profiling facility at the Bindley Bioscience Center, Purdue University, for providing the support for lipid data acquisition and analysis, and Miami Integrative Metabolomics Research Center.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2022.874828/full#supplementary-material

References

Abrass, C. K. (2004). Cellular lipid metabolism and the role of lipids in progressive renal disease. Am. J. Nephrol. 24, 46–53. doi:10.1159/000075925

PubMed Abstract | CrossRef Full Text | Google Scholar

Acott, T. S., and Kelley, M. J. (2008). Extracellular matrix in the trabecular meshwork. Exp. Eye Res. 86, 543–561. doi:10.1016/j.exer.2008.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Acott, T. S., Vranka, J. A., Keller, K. E., Raghunathan, V., and Kelley, M. J. (2021). Normal and glaucomatous outflow regulation. Prog. Retin. Eye Res. 82, 100897. doi:10.1016/j.preteyeres.2020.100897

PubMed Abstract | CrossRef Full Text | Google Scholar

Allison, K., Patel, D., and Alabi, O. (2020). Epidemiology of glaucoma: The past, present, and predictions for the future. Cureus 12, e11686. doi:10.7759/cureus.11686

PubMed Abstract | CrossRef Full Text | Google Scholar

Apopa, P. L., He, X., and Ma, Q. (2008). Phosphorylation of Nrf2 in the transcription activation domain by casein kinase 2 (CK2) is critical for the nuclear translocation and transcription activation function of Nrf2 in IMR-32 neuroblastoma cells. J. Biochem. Mol. Toxicol. 22, 63–76. doi:10.1002/jbt.20212

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashwinbalaji, S., Haribalaganesh, R., Krishnadas, S., Muthukkaruppan, V., and Senthilkumari, S. (2020). SB772077B (SB77) alleviated the aqueous outflow resistance mediated by cyclic mechanical stress in perfused human cadaveric eyes. Sci. Rep. 10, 10202. doi:10.1038/s41598-020-67087-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bakke, E. F., Hisdal, J., and Semb, S. O. (2009). Intraocular pressure increases in parallel with systemic blood pressure during isometric exercise. Investig. Ophthalmol. Vis. Sci. 50, 760–764. doi:10.1167/iovs.08-2508

PubMed Abstract | CrossRef Full Text | Google Scholar

Banks, C. A., Boanca, G., Lee, Z. T., Eubanks, C. G., Hattem, G. L., Peak, A., et al. (2016). TNIP2 is a hub protein in the NF-κB network with both protein and RNA mediated interactions. Mol. Cell. Proteomics 15, 3435–3449. doi:10.1074/mcp.m116.060509

PubMed Abstract | CrossRef Full Text | Google Scholar

Behrends, C., Sowa, M. E., Gygi, S. P., and Harper, J. W. (2010). Network organization of the human autophagy system. Nature 466, 68–76. doi:10.1038/nature09204

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B Methodol. 57, 289–300. doi:10.1111/j.2517-6161.1995.tb02031.x

CrossRef Full Text | Google Scholar

Birgisdottir, Å. B., Mouilleron, S., Bhujabal, Z., Wirth, M., Sjøttem, E., Evjen, G., et al. (2019). Members of the autophagy class III phosphatidylinositol 3-kinase complex I interact with GABARAP and GABARAPL1 via LIR motifs. Autophagy 15, 1333–1355. doi:10.1080/15548627.2019.1581009

PubMed Abstract | CrossRef Full Text | Google Scholar

Bligh, E. G., and Dyer, W. J. (1959). A rapid method of total lipid extraction and purification. Can. J. Biochem. Physiol. 37, 911–917. doi:10.1139/o59-099

PubMed Abstract | CrossRef Full Text | Google Scholar

Borras, T. (2003). Gene expression in the trabecular meshwork and the influence of intraocular pressure. Prog. Retin. Eye Res. 22, 435–463. doi:10.1016/s1350-9462(03)00018-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Bradley, J. M., Kelley, M. J., Zhu, X., Anderssohn, A. M., Alexander, J. P., and Acott, T. S. (2001). Effects of mechanical stretching on trabecular matrix metalloproteinases. Investig. Ophthalmol. Vis. Sci. 42, 1505–1513.

PubMed Abstract | Google Scholar

Breese, M. R., and Liu, Y. (2013). NGSUtils: A software suite for analyzing and manipulating next-generation sequencing datasets. Bioinformatics 29, 494–496. doi:10.1093/bioinformatics/bts731

PubMed Abstract | CrossRef Full Text | Google Scholar

Bridge, G., Rashid, S., and Martin, S. A. (2014). DNA mismatch repair and oxidative DNA damage: Implications for cancer biology and treatment. Cancers 6, 1597–1614. doi:10.3390/cancers6031597

PubMed Abstract | CrossRef Full Text | Google Scholar

Burgy, O., Loriod, S., Beltramo, G., and Bonniaud, P. (2022). Extracellular lipids in the lung and their role in pulmonary fibrosis. Cells 11, 1209. doi:10.3390/cells11071209

PubMed Abstract | CrossRef Full Text | Google Scholar

Carrasco, S., and Mérida, I. (2007). Diacylglycerol, when simplicity becomes complex. Trends Biochem. Sci. 32, 27–36. doi:10.1016/j.tibs.2006.11.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cerutti, C., and Ridley, A. J. (2017). Endothelial cell-cell adhesion and signaling. Exp. Cell. Res. 358, 31–38. doi:10.1016/j.yexcr.2017.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H. Y., and Lai, S. W. (2005). Relation between intraocular pressure and systemic health parameters in Taiwan. South. Med. J. 98, 28–32. doi:10.1097/01.smj.0000145389.15201.7f

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Smeekens, J. M., and Wu, R. (2016). Systematic study of the dynamics and half-lives of newly synthesized proteins in human cells. Chem. Sci. 7, 1393–1400. doi:10.1039/c5sc03826j

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Lin, Y., Vithana, E. N., Jia, L., Zuo, X., Wong, T. Y., et al. (2014). Common variants near ABCA1 and in PMM2 are associated with primary open-angle glaucoma. Nat. Genet. 46, 1115–1119. doi:10.1038/ng.3078

PubMed Abstract | CrossRef Full Text | Google Scholar

Choudhary, S., Xiao, T., Vergara, L. A., Srivastava, S., Nees, D., Piatigorsky, J., et al. (2005). Role of aldehyde dehydrogenase isozymes in the defense of rat lens and human lens epithelial cells against oxidative stress. Investig. Ophthalmol. Vis. Sci. 46, 259–267. doi:10.1167/iovs.04-0120

PubMed Abstract | CrossRef Full Text | Google Scholar

Clyde-Brockway, C. E., Flaherty, E. A., and Paladino, F. (2021). Lipid profiles of Pacific green and hawksbill sea turtle plasma.

Google Scholar

Coleman, R. A., and Mashek, D. G. (2011). Mammalian triacylglycerol metabolism: Synthesis, lipolysis, and signaling. Chem. Rev. 111, 6359–6386. doi:10.1021/cr100404w

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, L., Fu, S., Zhang, J., Zhao, J., and Zhang, Y. (2018). Effects of atorvastatin on porcine aqueous humour outflow and trabecular meshwork cells. Exp. Ther. Med. 15, 210–216. doi:10.3892/etm.2017.5353

PubMed Abstract | CrossRef Full Text | Google Scholar

Dall'Armi, C., Devereaux, K. A., and Di Paolo, G. (2013). The role of lipids in the control of autophagy. Curr. Biol. 23, R33–R45. doi:10.1016/j.cub.2012.10.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Dawe, H. R., Adams, M., Wheway, G., Szymanska, K., Logan, C. V., Noegel, A. A., et al. (2009). Nesprin-2 interacts with meckelin and mediates ciliogenesis via remodelling of the actin cytoskeleton. J. Cell. Sci. 122, 2716–2726. doi:10.1242/jcs.043794

PubMed Abstract | CrossRef Full Text | Google Scholar

de Las Fuentes, L., Sung, Y. J., Schwander, K. L., Kalathiveetil, S., Hunt, S. C., Arnett, D. K., et al. (2013). The role of SNP-loop diuretic interactions in hypertension across ethnic groups in HyperGEN. Front. Genet. 4, 304. doi:10.3389/fgene.2013.00304

PubMed Abstract | CrossRef Full Text | Google Scholar

de Lima, C. B., Ferreira, C. R., Milazzotto, M. P., Sobreira, T. J. P., Vireque, A. A., and Cooks, R. G. (2018). Comprehensive lipid profiling of early stage oocytes and embryos by MRM profiling. J. Mass Spectrom. 53, 1247–1252. doi:10.1002/jms.4301

PubMed Abstract | CrossRef Full Text | Google Scholar

Deshmukh, P., Unni, S., Krishnappa, G., and Padmanabhan, B. (2017). The keap1-nrf2 pathway: Promising therapeutic target to counteract ROS-mediated damage in cancers and neurodegenerative diseases. Biophys. Rev. 9, 41–56. doi:10.1007/s12551-016-0244-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Dipali, S. S., Ferreira, C. R., Zhou, L. T., Pritchard, M. T., and Duncan, F. E. (2019). Histologic analysis and lipid profiling reveal reproductive age-associated changes in peri-ovarian adipose tissue. Reprod. Biol. Endocrinol. 17, 46. doi:10.1186/s12958-019-0487-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Dobin, A., Davis, C. A., Schlesinger, F., Drenkow, J., Zaleski, C., Jha, S., et al. (2013). Star: Ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21. doi:10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

Duffy, L., and O’Reilly, S. (2018). Functional implications of cross-linked actin networks in trabecular meshwork cells. Cell. Physiol. Biochem. 45, 783–794. doi:10.1159/000487170

PubMed Abstract | CrossRef Full Text | Google Scholar

Edwards, S. H., Kimberly, M. M., Pyatt, S. D., Stribling, S. L., Dobbin, K. D., and Myers, G. L. (2011). Proposed serum cholesterol reference measurement procedure by gas chromatography-isotope dilution mass spectrometry. Clin. Chem. 57, 614–622. doi:10.1373/clinchem.2010.158766

PubMed Abstract | CrossRef Full Text | Google Scholar

Elliott, M. H., Ashpole, N. E., Gu, X., Herrnberger, L., McClellan, M. E., Griffith, G. L., et al. (2016). Caveolin-1 modulates intraocular pressure: Implications for caveolae mechanoprotection in glaucoma. Sci. Rep. 6, 37127. doi:10.1038/srep37127

PubMed Abstract | CrossRef Full Text | Google Scholar

Fagone, P., and Jackowski, S. (2009). Membrane phospholipid synthesis and endoplasmic reticulum function. J. Lipid Res. 50 (1), S311–S316. doi:10.1194/jlr.r800049-jlr200

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferretti, R., Palumbo, V., Di Savino, A., Velasco, S., Sbroggiò, M., Sportoletti, P., et al. (2010). Morgana/chp-1, a ROCK inhibitor involved in centrosome duplication and tumorigenesis. Dev. Cell. 18, 486–495. doi:10.1016/j.devcel.2009.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Flanders, K. C. (2004). Smad3 as a mediator of the fibrotic response. Int. J. Exp. Pathol. 85, 47–64. doi:10.1111/j.0959-9673.2004.00377.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Fyrst, H., and Saba, J. D. (2010). An update on sphingosine-1-phosphate and other sphingolipid mediators. Nat. Chem. Biol. 6, 489–497. doi:10.1038/nchembio.392

PubMed Abstract | CrossRef Full Text | Google Scholar

Gao, K., Song, S., Johnstone, M. A., Zhang, Q., Xu, J., Zhang, X., et al. (2020). Reduced pulsatile trabecular meshwork motion in eyes with primary open angle glaucoma using phase-sensitive optical coherence tomography. Investig. Ophthalmol. Vis. Sci. 61, 21. doi:10.1167/iovs.61.14.21

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, S. X., Jung, D., and Yao, R. (2020). ShinyGO: A graphical gene-set enrichment tool for animals and plants. Bioinformatics 36, 2628–2629. doi:10.1093/bioinformatics/btz931

PubMed Abstract | CrossRef Full Text | Google Scholar

Gharahkhani, P., Jorgenson, E., Jorgenson, P., Hysi, A. P., Khawaja, S., Pendergrass, X., et al. (2021). Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries. Nat. Commun. 12, 1258. doi:10.1038/s41467-020-20851-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, A. K., and Vaughan, D. E. (2012). PAI-1 in tissue fibrosis. J. Cell. Physiol. 227, 493–507. doi:10.1002/jcp.22783

PubMed Abstract | CrossRef Full Text | Google Scholar

Gills, J. P., Roberts, B. C., and Epstein, D. L. (1998). Microtubule disruption leads to cellular contraction in human trabecular meshwork cells. Investig. Ophthalmol. Vis. Sci. 39, 653–658.

PubMed Abstract | Google Scholar

Gonzalez, P., Epstein, D. L., and Borrás, T. (2000). Genes upregulated in the human trabecular meshwork in response to elevated intraocular pressure. Investig. Ophthalmol. Vis. Sci. 41, 352–361.

PubMed Abstract | Google Scholar

Gotoh, T., Terada, K., Oyadomari, S., and Mori, M. (2004). hsp70-DnaJ chaperone pair prevents nitric oxide- and CHOP-induced apoptosis by inhibiting translocation of Bax to mitochondria. Cell. Death Differ. 11, 390–402. doi:10.1038/sj.cdd.4401369

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurney, A., Axelrod, F., Bond, C. J., Cain, J., Chartier, C., Donigan, L., et al. (2012). Wnt pathway inhibition via the targeting of Frizzled receptors results in decreased growth and tumorigenicity of human tumors. Proc. Natl. Acad. Sci. U.S.A. 109, 11717–11722. doi:10.1073/pnas.1120068109

PubMed Abstract | CrossRef Full Text | Google Scholar

Hann, C. R., and Fautsch, M. P. (2015). Recent developments in understanding the role of aqueous humor outflow in normal and primary open angle glaucoma. Curr. Ophthalmol. Rep. 3, 67–73. doi:10.1007/s40135-015-0072-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hernandez, H., Roberts, A. L., and McDowell, C. M. (2020). Nuclear factor-kappa beta signaling is required for transforming growth factor Beta-2 induced ocular hypertension. Exp. Eye Res. 191, 107920. doi:10.1016/j.exer.2020.107920

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, L. T. Y., Osterwald, A., Ruf, I., Hunziker, D., Mattei, P., Challa, P., et al. (2020). Role of the autotaxin-lysophosphatidic acid axis in glaucoma, aqueous humor drainage and fibrogenic activity. Biochimica Biophysica Acta (BBA) - Mol. Basis Dis. 1866, 165560. doi:10.1016/j.bbadis.2019.165560

PubMed Abstract | CrossRef Full Text | Google Scholar

Ho, L. T. Y., Skiba, N., Ullmer, C., and Rao, P. V. (2018). Lysophosphatidic acid induces ECM production via activation of the mechanosensitive YAP/TAZ transcriptional pathway in trabecular meshwork cells. Investig. Ophthalmol. Vis. Sci. 59, 1969–1984. doi:10.1167/iovs.17-23702

PubMed Abstract | CrossRef Full Text | Google Scholar

Höhfeld, J., Benzing, T., Bloch, W., Fürst, D. O., Gehlert, S., Hesse, M., et al. (2021). Maintaining proteostasis under mechanical stress. EMBO Rep. 22, e52507. doi:10.15252/embr.202152507

PubMed Abstract | CrossRef Full Text | Google Scholar

Honjo, M., Igarashi, N., Kurano, M., Yatomi, Y., Igarashi, K., Kano, K., et al. (2018). Autotaxin-lysophosphatidic acid pathway in intraocular pressure regulation and glaucoma subtypes. Investig. Ophthalmol. Vis. Sci. 59, 693–701. doi:10.1167/iovs.17-23218

PubMed Abstract | CrossRef Full Text | Google Scholar

Horn, A., and Jaiswal, J. K. (2019). Structural and signaling role of lipids in plasma membrane repair. Curr. Top. Membr. 84, 67–98. doi:10.1016/bs.ctm.2019.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Houten, S. M., and Wanders, R. J. A. (2010). A general introduction to the biochemistry of mitochondrial fatty acid β‐oxidation. J Inher Metab Disea 33, 469–477. doi:10.1007/s10545-010-9061-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hymowitz, M. B., Chang, D., Feinberg, E. B., and Roy, S. (2016). Increased intraocular pressure and hyperglycemic level in diabetic patients. PLOS ONE 11, e0151833. doi:10.1371/journal.pone.0151833

PubMed Abstract | CrossRef Full Text | Google Scholar

Iwamoto, T., Kita, S., Zhang, J., Blaustein, M. P., Arai, Y., Yoshida, S., et al. (2004). Salt-sensitive hypertension is triggered by Ca2+ entry via Na+/Ca2+ exchanger type-1 in vascular smooth muscle. Nat. Med. 10, 1193–1199. doi:10.1038/nm1118

PubMed Abstract | CrossRef Full Text | Google Scholar

Iyer, P., Lalane, R., Morris, C., Challa, P., Vann, R., and Rao, P. V. (2012). Autotaxin-lysophosphatidic acid axis is a novel molecular target for lowering intraocular pressure. PLoS One 7, e42627. doi:10.1371/journal.pone.0042627

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakobsson, M. E., Małecki, J. M., Halabelian, L., Nilges, B. S., Pinto, R., Kudithipudi, S., et al. (2018). The dual methyltransferase METTL13 targets N terminus and Lys55 of eEF1A and modulates codon-specific translation rates. Nat. Commun. 9, 3411. doi:10.1038/s41467-018-05646-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakobsson, M. E., Małecki, J., Nilges, B. S., Moen, A., Leidel, S. A., and Falnes, P. Ø. (2017). Methylation of human eukaryotic elongation factor alpha (eEF1A) by a member of a novel protein lysine methyltransferase family modulates mRNA translation. Nucleic Acids Res. 45, 8239–8254. doi:10.1093/nar/gkx432

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, M., and Erickson, K. (2000). Aqueous humor and the dynamics of its flow. Princ. Pract. Ophthalmol., 2577–2595.

Google Scholar

Johnstone, M. A. (2004a). The aqueous outflow system as a mechanical pump. J. glaucoma 13, 421–438. doi:10.1097/01.ijg.0000131757.63542.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnstone, M. A. (2004b). The aqueous outflow system as a mechanical pump. J. Glaucoma 13, 421–438. doi:10.1097/01.ijg.0000131757.63542.24

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnstone, M., Xin, C., Tan, J., Martin, E., Wen, J., and Wang, R. K. (2021). Aqueous outflow regulation - 21st century concepts. Prog. Retin. Eye Res. 83, 100917. doi:10.1016/j.preteyeres.2020.100917

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaiser-Kupfer, M. I. (1989). Neural crest origin of trabecular meshwork cells and other structures of the anterior chamber. Am. J. Ophthalmol. 107, 671–672. doi:10.1016/0002-9394(89)90266-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanayama, A., Seth, R. B., Sun, L., Ea, C. K., Hong, M., Shaito, A., et al. (2004). TAB2 and TAB3 activate the NF-κB pathway through binding to polyubiquitin chains. Mol. Cell. 15, 535–548. doi:10.1016/j.molcel.2004.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kanno, K., Wu, M. K., Scapa, E. F., Roderick, S. L., and Cohen, D. E. (2007). Structure and function of phosphatidylcholine transfer protein (PC-TP)/StarD2. Biochimica Biophysica Acta (BBA) - Mol. Cell. Biol. Lipids 1771, 654–662. doi:10.1016/j.bbalip.2007.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, K. E., Bhattacharya, S. K., Borrás, T., Brunner, T. M., Chansangpetch, S., Clark, A. F., et al. (2018). Consensus recommendations for trabecular meshwork cell isolation, characterization and culture. Exp. Eye Res. 171, 164–173. doi:10.1016/j.exer.2018.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, K. E., Kelley, M. J., and Acott, T. S. (2007). Extracellular matrix gene alternative splicing by trabecular meshwork cells in response to mechanical stretching. Investig. Ophthalmol. Vis. Sci. 48, 1164–1172. doi:10.1167/iovs.06-0875

PubMed Abstract | CrossRef Full Text | Google Scholar

Khurana, R. N., Deng, P. F., Epstein, D. L., and Vasantha Rao, P. (2003). The role of protein kinase C in modulation of aqueous humor outflow facility. Exp. Eye Res. 76, 39–47. doi:10.1016/s0014-4835(02)00255-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kneitz, B., Cohen, P. E., Avdievich, E., Zhu, L., Kane, M. F., Hou, H., et al. (2000). MutS homolog 4 localization to meiotic chromosomes is required for chromosome pairing during meiosis in male and female mice. Genes. Dev. 14, 1085–1097. doi:10.1101/gad.14.9.1085

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, A., and Panda, D. (2018). Regulation of microtubule stability by centrosomal proteins. IUBMB Life 70, 602–611. doi:10.1002/iub.1865

PubMed Abstract | CrossRef Full Text | Google Scholar

Lakk, M., Hoffmann, G. F., Gorusupudi, A., Enyong, E., Lin, A., Bernstein, P. S., et al. (2021). Membrane cholesterol regulates TRPV4 function, cytoskeletal expression, and the cellular response to tension. J. lipid Res. 62, 100145. doi:10.1016/j.jlr.2021.100145

PubMed Abstract | CrossRef Full Text | Google Scholar

Last, J. A., Pan, T., Ding, Y., Reilly, C. M., Keller, K., Acott, T. S., et al. (2011). Elastic modulus determination of normal and glaucomatous human trabecular meshwork. Investig. Ophthalmol. Vis. Sci. 52, 2147–2152. doi:10.1167/iovs.10-6342

PubMed Abstract | CrossRef Full Text | Google Scholar

Laurin, M., and Côté, J. F. (2014). Insights into the biological functions of Dock family guanine nucleotide exchange factors. Genes. Dev. 28, 533–547. doi:10.1101/gad.236349.113

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, M. J., and Yaffe, M. B. (2016). Protein regulation in signal transduction. Cold Spring Harb. Perspect. Biol. 8. doi:10.1101/cshperspect.a005918

PubMed Abstract | CrossRef Full Text | Google Scholar

LeMosy, E. K., Lightner, V. A., and Erickson, H. P. (1996). Structural analysis of a human glial variant laminin. Exp. Cell. Res. 227, 80–88. doi:10.1006/excr.1996.0252

PubMed Abstract | CrossRef Full Text | Google Scholar

Lewis, C. A., Griffiths, B., Santos, C. R., Pende, M., and Schulze, A. (2011). Regulation of the SREBP transcription factors by mTORC1. Biochem. Soc. Trans. 39, 495–499. doi:10.1042/bst0390495

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Yue, Z., Xiong, W., Sun, P., You, K., and Wang, J. (2017a). TXNIP overexpression suppresses proliferation and induces apoptosis in SMMC7221 cells through ROS generation and MAPK pathway activation. Oncol. Rep. 37, 3369–3376. doi:10.3892/or.2017.5577

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X. C., Hu, Q. K., Chen, L., Liu, S. Y., Su, S., Tao, H., et al. (2017b). HSPB8 promotes the fusion of autophagosome and lysosome during autophagy in diabetic neurons. Int. J. Med. Sci. 14, 1335–1341. doi:10.7150/ijms.20653

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, Y., Smyth, G. K., and Shi, W. (2014). featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930. doi:10.1093/bioinformatics/btt656

PubMed Abstract | CrossRef Full Text | Google Scholar

Lignell, A., Kerosuo, L., Streichan, S. J., Cai, L., and Bronner, M. E. (2017). Identification of a neural crest stem cell niche by Spatial Genomic Analysis. Nat. Commun. 8, 1830. doi:10.1038/s41467-017-01561-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Liton, P. B., Liu, X., Challa, P., Epstein, D. L., and Gonzalez, P. (2005). Induction of TGF-β1 in the trabecular meshwork under cyclic mechanical stress. J. Cell. Physiol. 205, 364–371. doi:10.1002/jcp.20404

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, B., McNally, S., Kilpatrick, J. I., Jarvis, S. P., and O'Brien, C. J. (2018). Aging and ocular tissue stiffness in glaucoma. Surv. Ophthalmol. 63, 56–74. doi:10.1016/j.survophthal.2017.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G., Edmonds, B. T., and Condeelis, J. (1996). pH, EF-1α and the cytoskeleton. Trends Cell. Biol. 6, 168–171. doi:10.1016/0962-8924(96)20013-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G. Y., Jiang, X. X., Zhu, X., He, W. Y., Kuang, Y. L., Ren, K., et al. (2015). ROS activates JNK-mediated autophagy to counteract apoptosis in mouse mesenchymal stem cells in vitro. Acta Pharmacol. Sin. 36, 1473–1479. doi:10.1038/aps.2015.101

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Allingham, R. R., Qin, X., Layfield, D., Dellinger, A. E., Gibson, J., et al. (2013). Gene expression profile in human trabecular meshwork from patients with primary open-angle glaucoma. Investig. Ophthalmol. Vis. Sci. 54, 6382–6389. doi:10.1167/iovs.13-12128

PubMed Abstract | CrossRef Full Text | Google Scholar

López-Gálvez, R., de la Morena-Barrio, M. E., López-Lera, A., Pathak, M., Miñano, A., Serrano, M., et al. (2020). Factor XII in PMM2-CDG patients: Role of N-glycosylation in the secretion and function of the first element of the contact pathway. Orphanet J. Rare Dis. 15, 280. doi:10.1186/s13023-020-01564-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Luna, C., Li, G., Liton, P. B., Epstein, D. L., and Gonzalez, P. (2009). Alterations in gene expression induced by cyclic mechanical stress in trabecular meshwork cells. Mol. Vis. 15, 534–544.

PubMed Abstract | Google Scholar

Maddala, R., Skiba, N. P., and Rao, P. V. (2017). Vertebrate lonesome kinase regulated extracellular matrix protein phosphorylation, cell shape, and adhesion in trabecular meshwork cells. J. Cell. Physiol. 232, 2447–2460. doi:10.1002/jcp.25582

PubMed Abstract | CrossRef Full Text | Google Scholar

Madison, B. B. (2016). Srebp2: A master regulator of sterol and fatty acid synthesis. J. Lipid Res. 57, 333–335. doi:10.1194/jlr.c066712

PubMed Abstract | CrossRef Full Text | Google Scholar

Márquez, M. G., Leocata Nieto, F. L., Fernández-Tome, M. C., Favale, N. O., and Sterin-Speziale, N. (2008). Membrane lipid composition plays a central role in the maintenance of epithelial cell adhesion to the extracellular matrix. Lipids 43, 343–352. doi:10.1007/s11745-008-3152-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Martino, F., Perestrelo, A. R., Vinarský, V., Pagliari, S., and Forte, G. (2018). Cellular mechanotransduction: From tension to function. Front. Physiol. 9, 824. doi:10.3389/fphys.2018.00824

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsuo, T., Uchida, H., and Matsuo, N. (1996). Bovine and porcine trabecular cells produce prostaglandin F2 alpha in response to cyclic mechanical stretching. Jpn. J. Ophthalmol. 40, 289–296.

PubMed Abstract | Google Scholar

McCarthy, D. J., Chen, Y., and Smyth, G. K. (2012). Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 40, 4288–4297. doi:10.1093/nar/gks042

PubMed Abstract | CrossRef Full Text | Google Scholar

McDonough, H., Charles, P. C., Hilliard, E. G., Qian, S. B., Min, J. N., Portbury, A., et al. (2009). Stress-dependent Daxx-CHIP interaction suppresses the p53 apoptotic program. J. Biol. Chem. 284, 20649–20659. doi:10.1074/jbc.m109.011767

PubMed Abstract | CrossRef Full Text | Google Scholar

McGwin, G., Mcneal, S., Owsley, C., Girkin, C., Epstein, D., and Lee, P. P. (2004). Statins and other cholesterol-lowering medications and the presenceof glaucoma. Arch. Ophthalmol. 122, 822–826. doi:10.1001/archopht.122.6.822

PubMed Abstract | CrossRef Full Text | Google Scholar

McLeod, S. D., West, S. K., Quigley, H. A., and Fozard, J. L. (1990). A longitudinal study of the relationship between intraocular and blood pressures. Investig. Ophthalmol. Vis. Sci. 31, 2361–2366.

PubMed Abstract | Google Scholar

Mettu, P. S., Deng, P. F., Misra, U. K., Gawdi, G., Epstein, D. L., and Rao, P. V. (2004). Role of lysophospholipid growth factors in the modulation of aqueous humor outflow facility. Investig. Ophthalmol. Vis. Sci. 45, 2263–2271. doi:10.1167/iovs.03-0960

PubMed Abstract | CrossRef Full Text | Google Scholar

Najafinobar, N., Mellander, L. J., Kurczy, M. E., Dunevall, J., Angerer, T. B., Fletcher, J. S., et al. (2016). Cholesterol alters the dynamics of release in protein independent cell models for exocytosis. Sci. Rep. 6, 33702. doi:10.1038/srep33702

PubMed Abstract | CrossRef Full Text | Google Scholar

Nettesheim, A., Shim, M. S., Hirt, J., and Liton, P. B. (2019). Transcriptome analysis reveals autophagy as regulator of TGFβ/Smad-induced fibrogenesis in trabecular meshwork cells. Sci. Rep. 9, 16092. doi:10.1038/s41598-019-52627-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Nivon, M., Fort, L., Muller, P., Richet, E., Simon, S., Guey, B., et al. (2016). NFκB is a central regulator of protein quality control in response to protein aggregation stresses via autophagy modulation. MBoC 27, 1712–1727. doi:10.1091/mbc.e15-12-0835

PubMed Abstract | CrossRef Full Text | Google Scholar

Ogata, T., Naito, D., Nakanishi, N., Hayashi, K., Taniguchi, T., Miyagawa, K., et al. (2014). MURC/Cavin-4 facilitates recruitment of ERK to caveolae and concentric cardiac hypertrophy induced by α1-adrenergic receptors. Proc. Natl. Acad. Sci. U.S.A. 111, 3811–3816. doi:10.1073/pnas.1315359111

PubMed Abstract | CrossRef Full Text | Google Scholar

Okada, Y., Matsuo, T., and Ohtsuki, H. (1998). Bovine trabecular cells produce TIMP-1 and MMP-2 in response to mechanical stretching. Jpn. J. Ophthalmol. 42, 90–94. doi:10.1016/s0021-5155(97)00129-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pal, D., Wu, D., Haruta, A., Matsumura, F., and Wei, Q. (2010). Role of a novel coiled-coil domain-containing protein CCDC69 in regulating central spindle assembly. Cell. Cycle 9, 4117–4129. doi:10.4161/cc.9.20.13387

PubMed Abstract | CrossRef Full Text | Google Scholar

Parreno, J., and Fowler, V. M. (2018). Multifunctional roles of tropomodulin-3 in regulating actin dynamics. Biophys. Rev. 10, 1605–1615. doi:10.1007/s12551-018-0481-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pattabiraman, P. P., Inoue, T., and Rao, P. V. (2015a). Elevated intraocular pressure induces Rho GTPase mediated contractile signaling in the trabecular meshwork. Exp. Eye Res. 136, 29–33. doi:10.1016/j.exer.2015.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Pattabiraman, P. P., Maddala, R., and Rao, P. V. (2014). Regulation of plasticity and fibrogenic activity of trabecular meshwork cells by Rho GTPase signaling. J. Cell. Physiol. 229, 927–942. doi:10.1002/jcp.24524

PubMed Abstract | CrossRef Full Text | Google Scholar

Pattabiraman, P. P., and Rao, P. V. (2010). Mechanistic basis of Rho GTPase-induced extracellular matrix synthesis in trabecular meshwork cells. Am. J. Physiology-Cell Physiology 298, C749–C763. doi:10.1152/ajpcell.00317.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Pattabiraman, P. P., Rinkoski, T., Poeschla, E., Proia, A., Challa, P., and Rao, P. V. (2015b). RhoA GTPase-induced ocular hypertension in a rodent model is associated with increased fibrogenic activity in the trabecular meshwork. Am. J. Pathology 185, 496–512. doi:10.1016/j.ajpath.2014.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Pattabiraman, P. P., and Toris, C. B. (2016). The exit strategy: Pharmacological modulation of extracellular matrix production and deposition for better aqueous humor drainage. Eur. J. Pharmacol. 787, 32–42. doi:10.1016/j.ejphar.2016.04.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Perrakis, A., and Moolenaar, W. H. (2014). Autotaxin: Structure-function and signaling. J. Lipid Res. 55, 1010–1018. doi:10.1194/jlr.r046391

PubMed Abstract | CrossRef Full Text | Google Scholar

Peterson, T. R., Sengupta, S. S., Harris, T. E., Carmack, A. E., Kang, S. A., Balderas, E., et al. (2011). mTOR complex 1 regulates lipin 1 localization to control the SREBP pathway. Cell. 146, 408–420. doi:10.1016/j.cell.2011.06.034

PubMed Abstract | CrossRef Full Text | Google Scholar

Pipes, G. C., Creemers, E. E., and Olson, E. N. (2006). The myocardin family of transcriptional coactivators: Versatile regulators of cell growth, migration, and myogenesis. Genes. Dev. 20, 1545–1556. doi:10.1101/gad.1428006

PubMed Abstract | CrossRef Full Text | Google Scholar

Posch-Pertl, L., Michelitsch, M., Wagner, G., Wildner, B., Silbernagel, G., Pregartner, G., et al. (2022). Cholesterol and glaucoma: A systematic review and meta-analysis. Acta Ophthalmol. 100, 148–158. doi:10.1111/aos.14769

PubMed Abstract | CrossRef Full Text | Google Scholar

Proia, R. L., and Hla, T. (2015). Emerging biology of sphingosine-1-phosphate: Its role in pathogenesis and therapy. J. Clin. Investig. 125, 1379–1387. doi:10.1172/jci76369

PubMed Abstract | CrossRef Full Text | Google Scholar

Radons, J. (2016). The human HSP70 family of chaperones: Where do we stand? Cell. Stress Chaperones 21, 379–404. doi:10.1007/s12192-016-0676-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramos, R. F., Sumida, G. M., and Stamer, W. D. (2009). Cyclic mechanical stress and trabecular meshwork cell contractility. Investig. Ophthalmol. Vis. Sci. 50, 3826–3832. doi:10.1167/iovs.08-2694

PubMed Abstract | CrossRef Full Text | Google Scholar

Rao, A., Sahay, P., Chakraborty, M., Prusty, B. K., Srinivasan, S., Jhingan, G. D., et al. (2021). Switch to autophagy the key mechanism for trabecular meshwork death in severe glaucoma. Opth Vol. 15, 3027–3039. doi:10.2147/opth.s292218

CrossRef Full Text | Google Scholar

Rao, P. V. (2014). Bioactive lysophospholipids: Role in regulation of aqueous humor outflow and intraocular pressure in the context of pathobiology and therapy of glaucoma. J. Ocular Pharmacol. Ther. 30, 181–190. doi:10.1089/jop.2013.0194

PubMed Abstract | CrossRef Full Text | Google Scholar

Raviola, G., and Raviola, E. (1981). Paracellular route of aqueous outflow in the trabecular meshwork and canal of Schlemm. A freeze-fracture study of the endothelial junctions in the sclerocorneal angel of the macaque monkey eye. Investigative Ophthalmol. Vis. Sci. 21, 52–72.

PubMed Abstract | Google Scholar

Reina-Torres, E., Boussommier-Calleja, A., Sherwood, J. M., and Overby, D. R. (2020). Aqueous humor outflow requires active cellular metabolism in mice. Investig. Ophthalmol. Vis. Sci. 61, 45. doi:10.1167/iovs.61.10.45

PubMed Abstract | CrossRef Full Text | Google Scholar

Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodrigues, A. J., do Carmo Costa, M., Silva, T. L., Ferreira, D., Bajanca, F., Logarinho, E., et al. (2010). Absence of ataxin-3 leads to cytoskeletal disorganization and increased cell death. Biochimica Biophysica Acta (BBA) - Mol. Cell. Res. 1803, 1154–1163. doi:10.1016/j.bbamcr.2010.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Saccà, S. C., Gandolfi, S., Bagnis, A., Manni, G., Damonte, G., Traverso, C. E., et al. (2016). From DNA damage to functional changes of the trabecular meshwork in aging and glaucoma. Ageing Res. Rev. 29, 26–41. doi:10.1016/j.arr.2016.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Saccà, S. C., and Izzotti, A. (2008). Oxidative stress and glaucoma: Injury in the anterior segment of the eye. Prog. Brain Res. 173, 385–407. doi:10.1016/s0079-6123(08)01127-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Saccà, S. C., Pascotto, A., Camicione, P., Capris, P., and Izzotti, A. (2005). Oxidative DNA damage in the human trabecular meshwork. Arch. Ophthalmol. 123, 458–463. doi:10.1001/archopht.123.4.458

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahoo, S., Aurich, M., Jonsson, J., and Thiele, I. (2014). Membrane transporters in a human genome-scale metabolic knowledgebase and their implications for disease. Front. Physiology 5. doi:10.3389/fphys.2014.00091

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanui, T., Inayoshi, A., Noda, M., Iwata, E., Stein, J. V., Sasazuki, T., et al. (2003). DOCK2 regulates Rac activation and cytoskeletal reorganization through interaction with ELMO1. Blood 102, 2948–2950. doi:10.1182/blood-2003-01-0173

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, Y., Matsuo, T., and Ohtsuki, H. (1999). A novel gene (oculomedin) induced by mechanical stretching in human trabecular cells of the eye. Biochem. Biophysical Res. Commun. 259, 349–351. doi:10.1006/bbrc.1999.0797

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaedel, L., John, K., Gaillard, J., Nachury, M. V., Blanchoin, L., and Théry, M. (2015). Microtubules self-repair in response to mechanical stress. Nat. Mater 14, 1156–1163. doi:10.1038/nmat4396

PubMed Abstract | CrossRef Full Text | Google Scholar

Semënov, M. V., Zhang, X., and He, X. (2008). DKK1 antagonizes Wnt signaling without promotion of LRP6 internalization and degradation. J. Biol. Chem. 283, 21427–21432. doi:10.1074/jbc.m800014200

PubMed Abstract | CrossRef Full Text | Google Scholar

Shembade, N., Harhaj, N. S., Liebl, D. J., and Harhaj, E. W. (2007). Essential role for TAX1BP1 in the termination of TNF-α-, IL-1- and LPS-mediated NF-κB and JNK signaling. Embo J. 26, 3910–3922. doi:10.1038/sj.emboj.7601823

PubMed Abstract | CrossRef Full Text | Google Scholar

Shim, M. S., Nettesheim, A., Dixon, A., and Liton, P. B. (2021). Primary cilia and the reciprocal activation of AKT and SMAD2/3 regulate stretch-induced autophagy in trabecular meshwork cells. Proc. Natl. Acad. Sci. U. S. A. 118. doi:10.1073/pnas.2021942118

CrossRef Full Text | Google Scholar

Shim, M. S., Nettesheim, A., Hirt, J., and Liton, P. B. (2020). The autophagic protein LC3 translocates to the nucleus and localizes in the nucleolus associated to NUFIP1 in response to cyclic mechanical stress. Autophagy 16, 1248–1261. doi:10.1080/15548627.2019.1662584

PubMed Abstract | CrossRef Full Text | Google Scholar

Snowden, T., Acharya, S., Butz, C., Berardini, M., and Fishel, R. (2004). hMSH4-hMSH5 recognizes Holliday Junctions and forms a meiosis-specific sliding clamp that embraces homologous chromosomes. Mol. Cell. 15, 437–451. doi:10.1016/j.molcel.2004.06.040

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, J., Deng, P. F., Stinnett, S. S., Epstein, D. L., and Rao, P. V. (2005). Effects of cholesterol-lowering statins on the aqueous humor outflow pathway. Investig. Ophthalmol. Vis. Sci. 46, 2424–2432. doi:10.1167/iovs.04-0776

PubMed Abstract | CrossRef Full Text | Google Scholar

Soto-Avellaneda, A., and Morrison, B. E. (2020). Signaling and other functions of lipids in autophagy: A review. Lipids Health Dis. 19, 214. doi:10.1186/s12944-020-01389-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Stamer, W. D., Read, A. T., Sumida, G. M., and Ethier, C. R. (2009). Sphingosine-1-phosphate effects on the inner wall of Schlemm's canal and outflow facility in perfused human eyes. Exp. Eye Res. 89, 980–988. doi:10.1016/j.exer.2009.08.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Staverosky, J., Dhamodaran, K., Acott, T., Raghunathan, V., and Vranka, J. (2020). Isolation and Characterization of primary human trabecular meshwork cells from segmental flow regions: New tools for understanding segmental flow. Exp. Eye Res. 197, 108046. doi:10.1016/j.exer.2020.108046

PubMed Abstract | CrossRef Full Text | Google Scholar

Stein, J. D., Newman-Casey, P. A., Talwar, N., Nan, B., Richards, J. E., and Musch, D. C. (2012). The relationship between statin use and open-angle glaucoma. Ophthalmology 119, 2074–2081. doi:10.1016/j.ophtha.2012.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Stewart-Hutchinson, P. J., Hale, C. M., Wirtz, D., and Hodzic, D. (2008). Structural requirements for the assembly of LINC complexes and their function in cellular mechanical stiffness. Exp. Cell. Res. 314, 1892–1905. doi:10.1016/j.yexcr.2008.02.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Stumpff, F., and Wiederholt, M. (2000). Regulation of trabecular meshwork contractility. Ophthalmologica 214, 33–53. doi:10.1159/000027471

PubMed Abstract | CrossRef Full Text | Google Scholar

Suárez-Trujillo, A., Senn, L. K., Teeple, K., Casey, T. M., and Stewart, K. R. (2020). A standardized model to study effects of varying 24-h colostrum dose on postnatal growth and development. Transl. Anim. Sci. 4, txaa212. doi:10.1093/tas/txaa212

PubMed Abstract | CrossRef Full Text | Google Scholar

Sumida, G. M., and Stamer, W. D. (2011). S1P2 receptor regulation of sphingosine-1-phosphate effects on conventional outflow physiology. Am. J. Physiology-Cell Physiology 300, C1164–C1171. doi:10.1152/ajpcell.00437.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Sumida, G. M., and Stamer, W. D. (2010). Sphingosine-1-phosphate enhancement of cortical actomyosin organization in cultured human Schlemm's canal endothelial cell monolayers. Investig. Ophthalmol. Vis. Sci. 51, 6633–6638. doi:10.1167/iovs.10-5391

PubMed Abstract | CrossRef Full Text | Google Scholar

Susanna, R., De Moraes, C. G., Cioffi, G. A., and Ritch, R. (2015). Why do people (still) go blind from glaucoma? Trans. Vis. Sci. Tech. 4, 1. doi:10.1167/tvst.4.2.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Takeda, K., Shimozono, R., Noguchi, T., Umeda, T., Morimoto, Y., Naguro, I., et al. (2007). Apoptosis signal-regulating kinase (ASK) 2 functions as a mitogen-activated protein kinase kinase kinase in a heteromeric complex with ASK1. J. Biol. Chem. 282, 7522–7531. doi:10.1074/jbc.m607177200

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanito, M., Kaidzu, S., Takai, Y., and Ohira, A. (2016). Association between systemic oxidative stress and visual field damage in open-angle glaucoma. Sci. Rep. 6, 25792. doi:10.1038/srep25792

PubMed Abstract | CrossRef Full Text | Google Scholar

Tham, Y. C., Li, X., Wong, T. Y., Quigley, H. A., Aung, T., and Cheng, C. Y. (2014). Global prevalence of glaucoma and projections of glaucoma burden through 2040. Ophthalmology 121, 2081–2090. doi:10.1016/j.ophtha.2014.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Thieme, H., Nass, J. U., Nuskovski, M., Bechrakis, N. E., Stumpff, F., Strauss, O., et al. (1999). The effects of protein kinase C on trabecular meshwork and ciliary muscle contractility. Investig. Ophthalmol. Vis. Sci. 40, 3254–3261.

PubMed Abstract | Google Scholar

Thorleifsson, G., Walters, G. B., Hewitt, A. W., Masson, G., Helgason, A., DeWan, A., et al. (2010). Common variants near CAV1 and CAV2 are associated with primary open-angle glaucoma. Nat. Genet. 42, 906–909. doi:10.1038/ng.661

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, B., Gabelt, B. T., Peterson, J. A., Kiland, J. A., and Kaufman, P. L. (1999). H-7 increases trabecular facility and facility after ciliary muscle disinsertion in monkeys. Investig. Ophthalmol. Vis. Sci. 40, 239–242.

PubMed Abstract | Google Scholar

Tock, A. J., and Henderson, I. R. (2018). Hotspots for initiation of meiotic recombination. Front. Genet. 9, 521. doi:10.3389/fgene.2018.00521

PubMed Abstract | CrossRef Full Text | Google Scholar

Vahabikashi, A., Gelman, A., Dong, B., Gong, L., Cha, E. D., Schimmel, M., et al. (2019). Increased stiffness and flow resistance of the inner wall of Schlemm’s canal in glaucomatous human eyes. Proc. Natl. Acad. Sci. 116, 26555–26563.doi:10.1073/pnas.1911837116

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Meer, G. (2001). Caveolin, cholesterol, and lipid droplets? J. Cell. Biol. 152, F29–F34. doi:10.1083/jcb.152.5.f29

PubMed Abstract | CrossRef Full Text | Google Scholar

Vittal, V., Rose, A., Gregory, K. E., Kelley, M. J., and Acott, T. S. (2005). Changes in gene expression by trabecular meshwork cells in response to mechanical stretching. Investig. Ophthalmol. Vis. Sci. 46, 2857–2868. doi:10.1167/iovs.05-0075

PubMed Abstract | CrossRef Full Text | Google Scholar

Vittitow, J., and Borrás, T. (2004). Genes expressed in the human trabecular meshwork during pressure-induced homeostatic response. J. Cell. Physiol. 201, 126–137. doi:10.1002/jcp.20030

PubMed Abstract | CrossRef Full Text | Google Scholar

Von Zee, C. L., Richards, M. P., Bu, P., Perlman, J. I., and Stubbs, E. B. (2009). Increased RhoA and RhoB protein accumulation in cultured human trabecular meshwork cells by lovastatin. Investig. Ophthalmol. Vis. Sci. 50, 2816–2823. doi:10.1167/iovs.08-2466

CrossRef Full Text | Google Scholar

Walther, T. C., and Farese, R. V. (2012). Lipid droplets and cellular lipid metabolism. Annu. Rev. Biochem. 81, 687–714. doi:10.1146/annurev-biochem-061009-102430

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Johnstone, M. A., Xin, C., Song, S., Padilla, S., Vranka, J. A., et al. (2017). Estimating human trabecular meshwork stiffness by numerical modeling and advanced OCT imaging. Investigative Ophthalmol. Vis. Sci. 58, 4809–4817. doi:10.1167/iovs.17-22175

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., and Bao, X. (2019). Hyperlipidemia, blood lipid level, and the risk of glaucoma: A meta-analysis. Investig. Ophthalmol. Vis. Sci. 60, 1028–1043. doi:10.1167/iovs.18-25845

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Xu, L., Jonas, J. B., Wong, T. Y., Cui, T., Li, Y., et al. (2009). Major eye diseases and risk factors associated with systemic hypertension in an adult Chinese population: The beijing eye study. Ophthalmology 116, 2373–2380. doi:10.1016/j.ophtha.2009.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Weinreb, R. N., Leung, C. K., Crowston, J. G., Medeiros, F. A., Friedman, D. S., Wiggs, J. L., et al. (2016). Primary open-angle glaucoma. Nat. Rev. Dis. Prim. 2, 16067. doi:10.1038/nrdp.2016.67

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiederholt, M., Groth, J., and Strauss, O. (1998). Role of protein tyrosine kinase on regulation of trabecular meshwork and ciliary muscle contractility. Investig. Ophthalmol. Vis. Sci. 39, 1012–1020.

PubMed Abstract | Google Scholar

Wiederholt, M., Thieme, H., and Stumpff, F. (2000). The regulation of trabecular meshwork and ciliary muscle contractility. Prog. Retin Eye Res. 19, 271–295. doi:10.1016/s1350-9462(99)00015-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Wiggs, J. L., and Pasquale, L. R. (2017). Genetics of glaucoma. Hum. Mol. Genet. 26, R21–R27. doi:10.1093/hmg/ddx184

PubMed Abstract | CrossRef Full Text | Google Scholar

Wood, I. S., Hunter, L., and Trayhurn, P. (2003). Expression of Class III facilitative glucose transporter genes (GLUT-10 and GLUT-12) in mouse and human adipose tissues. Biochem. Biophys. Res. Commun. 308, 43–49. doi:10.1016/s0006-291x(03)01322-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, A., Khawaja, A. P., Pasquale, L. R., and Stein, J. D. (2020). A review of systemic medications that may modulate the risk of glaucoma. Eye Lond. Engl. 34, 12–28. doi:10.1038/s41433-019-0603-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, H., Shui, Y.-B., Liu, Y., Liu, X., and Siegfried, C. J. (2022). Trabecular meshwork mitochondrial function and oxidative stress: Clues to racial disparities of glaucoma. Ophthalmol. Sci. 2, 100107. doi:10.1016/j.xops.2021.100107

CrossRef Full Text | Google Scholar

Xin, C., Johnstone, M., Wang, N., and Wang, R. K. (2016). OCT study of mechanical properties associated with trabecular meshwork and collector channel motion in human eyes. PLoS One 11, e0162048. doi:10.1371/journal.pone.0162048

PubMed Abstract | CrossRef Full Text | Google Scholar

Xin, C., Wang, R. K., Song, S., Shen, T., Wen, J., Martin, E., et al. (2017). Aqueous outflow regulation: Optical coherence tomography implicates pressure-dependent tissue motion. Exp. Eye Res. 158, 171–186. doi:10.1016/j.exer.2016.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, H., Wu, X., Qin, H., Tian, W., Chen, J., Sun, L., et al. (2015). Myocardin-related transcription factor A epigenetically regulates renal fibrosis in diabetic nephropathy. J. Am. Soc. Nephrol. 26, 1648–1660. doi:10.1681/asn.2014070678

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, Y., and Her, C. (2013). VBP1 facilitates proteasome and autophagy-mediated degradation of MutS homologue hMSH4. Faseb J. 27, 4799–4810. doi:10.1096/fj.13-235127

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, B., Leyrat, C., Grimes, J. M., and Robinson, R. C. (2014). Structural basis of thymosin-β4/profilin exchange leading to actin filament polymerization. Proc. Natl. Acad. Sci. U. S. A. 111, E4596–E4605. doi:10.1073/pnas.1412271111

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, E., Van Nimwegen, E., Zavolan, M., Rajewsky, N., Schroeder, M., Magnasco, M., et al. (2003). Decay rates of human mRNAs: Correlation with functional characteristics and sequence attributes. Genome Res. 13, 1863–1872. doi:10.1101/gr.1272403

PubMed Abstract | CrossRef Full Text | Google Scholar

Youngblood, H., Cai, J., Drewry, M. D., Helwa, I., Hu, E., Liu, S., et al. (2020). Expression of mRNAs, miRNAs, and lncRNAs in human trabecular meshwork cells upon mechanical stretch. Investig. Ophthalmol. Vis. Sci. 61, 2. doi:10.1167/iovs.61.5.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Ren, C., Chen, L., Navedo, M. F., Antos, L. K., Kinsey, S. P., et al. (2010). Knockout of Na+/Ca2+ exchanger in smooth muscle attenuates vasoconstriction and L-type Ca2+ channel current and lowers blood pressure. Am. J. Physiol. Heart Circ. Physiol. 298, H1472–H1483. doi:10.1152/ajpheart.00964.2009

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, L. Y., Liu, J., Sidhu, G. S., Niu, Y., Liu, Y., Wang, R., et al. (2004). Negative regulation of p53 functions by Daxx and the involvement of MDM2. J. Biol. Chem. 279, 50566–50579. doi:10.1074/jbc.m406743200

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhavoronkov, A., Izumchenko, E., Kanherkar, R. R., Teka, M., Cantor, C., Manaye, K., et al. (2016). Pro-fibrotic pathway activation in trabecular meshwork and lamina cribrosa is the main driving force of glaucoma. Cell. Cycle 15, 1643–1652. doi:10.1080/15384101.2016.1170261

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, M., Wang, W., Huang, W., and Zhang, X. (2014). Diabetes mellitus as a risk factor for open-angle glaucoma: A systematic review and meta-analysis. PloS one 9, e102972e102972. doi:10.1371/journal.pone.0102972

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, L., Chen, S., and Chen, Y. (2011). Unraveling the biological functions of Smad7 with mouse models. Cell. Biosci. 1, 44. doi:10.1186/2045-3701-1-44

PubMed Abstract | CrossRef Full Text | Google Scholar