Influence of dexamethasone-induced matrices on the TM transcriptome

Pathologic bidirectional interactions between the extracellular matrix (ECM) and cells within the human trabecular meshwork (hTM) contribute to ocular hypertension. An in vitro model is needed to study these cell-matrix interactions and their effect on outflow homeostasis. This study aimed to determine whether pathogenic ECM derived from dexamethasone (DEX)-treated hTM cultures induces clinically relevant glaucoma-like changes in healthy hTM cells at the transcriptional level. Corneoscleral rims from non-glaucoma donors were used to isolate primary hTM cells after validation according to the consensus recommendations for TM culture. Normal hTM cells (n = 5) were plated on a coverslip and treated with 100 nM DEX or ethanol for four weeks. These cultures were then decellularized, plated with primary hTM cells, and allowed to grow for another 72 h. RNA was extracted from these hTM cells for stranded total RNA-Seq. Sequencing libraries prepared using the Zymo-Seq RiboFree Total RNA library kit were pooled and sequenced using Illumina NovaSeq 6000. After quality control, sequence reads were aligned to the human genome build hg19. Differential expression (DE) analyses were performed using paired multi-factorial ANOVA. The expression of several DE genes associated with glaucoma (ANGPTL2, PDE7B, C22orf23, COL4A1, ADAM12, IFT122, SEMA6C) was validated using EvaGreen-based Droplet Digital PCR (ddPCR) assays. Gene ontology analyses of the DE genes were performed using the PANTHER and NDEx IQA databases, and functional analyses were performed with the DAVID Bioinformatics software. Using a cutoff of p-value <0.05 and fold change ≥2.0, our differential analysis identified 267 up- and 135 down-regulated genes in DEX-induced ECM-treated cells compared to the control. These differentially expressed genes were found to play a significant role in pathways such as cytokine and oxidative stress-induced inflammation, integrin signaling, matrix remodeling, and angiogenesis. These findings were further supported by previously performed proteomics studies using the same model. Using ddPCR, we validated the expression of seven genes associated with the risk of primary open-angle glaucoma. These results not only provide support for the pathogenic ECM model of steroid-induced glaucoma, but also demonstrate that the pathologic changes induced by this model are indeed found at the transcriptional level. These findings further demonstrate that matrix changes significantly influence cell expression profiles, which enable further understanding of the molecular mechanisms underlying glaucomatous changes in the TM. However, future studies with a larger and more diverse set of samples and longer time points are needed to confirm the utility of this model for mechanistic studies.