The TGFβ Induced MicroRNAome of the Trabecular Meshwork

Primary open-angle glaucoma (POAG) is a progressive optic neuropathy with a complex, multifactorial aetiology. Raised intraocular pressure (IOP) is the most important clinically modifiable risk factor for POAG. All current pharmacological agents target aqueous humour dynamics to lower IOP. Newer therapeutic agents are required as some patients with POAG show a limited therapeutic response or develop ocular and systemic side effects to topical medication. Elevated IOP in POAG results from cellular and molecular changes in the trabecular meshwork driven by increased levels of transforming growth factor β (TGFβ) in the anterior segment of the eye. Understanding how TGFβ affects both the structural and functional changes in the outflow pathway and IOP is required to develop new glaucoma therapies that target the molecular pathology in the trabecular meshwork. In this study, we evaluated the effects of TGF-β1 and -β2 treatment on miRNA expression in cultured human primary trabecular meshwork cells. Our findings are presented in terms of specific miRNAs (miRNA-centric), but given miRNAs work in networks to control cellular pathways and processes, a pathway-centric view of miRNA action is also reported. Evaluating TGFβ-responsive miRNA expression in trabecular meshwork cells will further our understanding of the important pathways and changes involved in the pathogenesis of glaucoma and could lead to the development of miRNAs as new therapeutic modalities in glaucoma.


Introduction
Glaucoma is a progressive optic neuropathy with loss of retinal ganglion cells and characteristic visual field and optic nerve changes [1].Glaucoma affects over 60 million people worldwide and the prevalence is projected to increase to over 118 million people by 2040 [2].Primary open-angle glaucoma (POAG) forms the most prevalent subtype of glaucoma and has a complex, multifactorial aetiology.Raised intraocular pressure (IOP) is the most important clinically modifiable risk factor for POAG.IOP is generated by the conventional aqueous humour outflow pathway via the trabecular meshwork and all current pharmacological agents target aqueous humour dynamics to lower IOP [3,4].
However, some patients with POAG show a limited therapeutic response or are refractory to these pharmacological agents [5][6][7], while others have clinically significant adverse physiological impacts, such as bradycardia or bronchospasm [8].Until the development of Rho kinase (ROCK) inhibitors, none of the topical drug treatments for POAG targeted the underlying cellular pathophysiology in the trabecular meshwork [9,10].ROCK inhibitors pharmacologically manipulate the cytoskeleton of trabecular meshwork and Schlemm's canal cells reducing outflow resistance and lowering IOP [9,10].Relatively common ocular side effects like conjunctival hyperaemia have limited tolerability of ROCK inhibitors in the pharmacological management of POAG [3,4,10].Therefore, there is a clinical need to develop novel therapies for POAG specifically targeting the molecular pathology in the trabecular meshwork to lower IOP [3,4,11].
Elevated IOP in POAG results from cellular and molecular changes in the trabecular meshwork (TM) driven by increased levels of transforming growth factor β (TGFβ) in the anterior segment of the eye [12].TGF-β1 and -β2 result in pathogenic changes in the human trabecular meshwork (TM) cell population and phenotype, which contribute to increased IOP [13].Various studies have identified elevated TGFβ2 levels in the aqueous humour of POAG patients [14][15][16][17][18], but the cause of these elevated levels is unclear.Elevated TGFβ2 levels have also been reported in the glaucomatous TM, indicating that increased levels of total and mature TGFβ may play an important role in the pathogenesis of POAG [19].TGFβ2 perfusion in an anterior eye segment organ culture model resulted in elevated intra-ocular pressure (IOP) and fibrillary material accumulation in the trabecular meshwork [15,20].In pseudoexfoliation glaucoma (XFG), an aggressive form of secondary open-angle glaucoma, both the latent and active forms of TGFβ1 are increased in the aqueous humour and extracellular matrix (ECM) [21][22][23].
The TGFβ receptor complex transmits signals via both canonical and non-canonical pathways [24][25][26].Canonical TGFβ signalling results in Smad2/3-Smad4 complexes translocating to the nucleus where they function as transcription factors that can regulate gene expression including microRNAs (miRNAs) [27][28][29][30].MiRNAs are small, singlestranded, noncoding RNAs which are important regulators of eukaryotic gene expression in health and disease [31].Normally, miRNAs bind to the 3 ′ untranslated region (UTR) of their target mRNAs resulting in mRNA degradation or translation repression [32].Most TGFβ signalling pathway components are known to be targeted by one or more miRNAs, and miRNA regulation of TGFβ signalling molecules influences the pathogenesis of fibrotic diseases [33].TGFβ signalling can also regulate miRNA expression at transcriptional and post-transcriptional levels [25].TGFβ is both cell-and context-specific in terms of its functionality.Therefore, research must focus on building our understanding of how TGFβ affects both the structural and functional changes in the outflow pathway and therefore IOP to develop new glaucoma therapies that target the molecular pathology in the TM.
Identifying the expression profiles of miRNAs in POAG and XFG could help us to better understand the changes in gene expression in the TM as well as gain insight into potential pathways that may be involved in the pathogenesis of the disease.The ability of miRNA manipulations to alter gene expression has raised the possibility of miRNAbased therapeutics [34][35][36][37].Using miRNA-Seq the miRNA expression profile in the normal human trabecular meshwork has been established [38].However, no previous studies have looked at miRNA expression changes in the TM following TGFβ treatment.Therefore, our study is the first miRNA-Seq study to evaluate the effects of TGF-β1 and -β2 treatment on miRNA expression in normal human primary TM cells.Evaluating TGFβ-responsive miRNA expression in the TM will further our understanding of the important pathways and changes involved in the pathogenesis of POAG and XFG and could lead to the development of miRNAs as new therapeutic modalities in glaucoma.

Human Trabecular Meshwork Culture and Characterisation
Cadaveric eyes (n = 4) were provided by the Liverpool Research Eye Bank, approved by the local ethics review board (RETH000833), and handled following the tenets of the Declaration of Helsinki.Donor eyes were obtained from the Royal Liverpool University Hospital mortuary.Medical history for the donor eyes was unknown, however, no donors had previous ocular surgery or a known glaucoma diagnosis.Donor eyes were excluded if the maximum post-mortem time exceeded 48 h.Primary normal human trabecular meshwork (TM) cells were isolated using the blunt dissection method as reported previously [39].Cells were maintained in Dulbecco's Modified Eagle Media (DMEM)-low glucose (Sigma, Gillingham, UK) supplemented with 10% foetal calf serum (Biosera, Heathfield, UK), 2 mM L-glutamine (Sigma, Gillingham, UK), Pen/Step (Sigma, Gillingham, UK), and 2.5 µg/mL Fungizone (amphotericin B, Sigma, Gillingham, UK).Samples were incubated at 37 • C (5% CO 2 and 95% humidity).Human TM characterisation was performed as previously described [39] including dexamethasone upregulation of myocilin protein expression (polyclonal rabbit anti-myocilin primary antibody was a kind gift from Dr. W. Daniel Stamer) as previously described by our group [40].

TGFβ Treatment and RNA Extraction
Human TM cells between passages 4 and 6 were grown to 70-80% confluency and growth arrested overnight using serum free medium before TGF-β1 or -β2 stimulation.with either recombinant human TGFβ1 (240-B-010, R&D Systems, Abingdon, UK) or TGFβ2 (302-B2-010, R&D Systems, UK) at a concentration of 5 ng/mL for 24 h.Total RNA was extracted from cells using the miRNeasy Mini Kit (Qiagen, Manchester, UK) following the manufacturer's instructions.The RNA concentration was measured using the NanoDrop 2000 (Thermofisher Scientific, Horsham, UK) and RNA quality was determined by the Bioanalyser 2100 (Agilent Technologies, Stockport, UK) using an RNA 6000 Nano Kit (Agilent, Santa Clara, CA, USA).

Small RNA Sequencing and Data Analysis
Small RNA sequencing was performed at the Genomics Core Technology Unit (Queen's University, Belfast, UK) using the Qiagen QIAseq miRNA Library Kit (Qiagen, Manchester, UK) to construct small RNA sequencing libraries.Sequencing was performed using a high-output sequencing kit (76 cycles) on a NextSeq500 (Illumina Inc., San Diego, CA, USA) according to the manufacturer's instructions.The raw sequence data reads were evaluated using FastQC v0.11.9 [41] to check the quality of the reads.Cutadapt 3.0 [42] was applied to remove low-quality reads (reads with a Q score/Phred score of less than 30; 99% accuracy) and adapters to ensure only the highest-quality sequences were included for further analysis.Following quality control, the Subread and RSubread packages [43] were used to align the remaining sequences to the human genome reference (GRCh38) from the Human GENCODE Gene Set [44].The reference database miRbase (v22.1) was used for miRNA alignment [45].A count matrix of mapped reads for each miRNA was generated using featureCount [46].These raw counts were used for downstream differential expression analysis.Counts data were generated for each paired treatment group and the significance of miRNA expression was compared using the generalised linear model (GLM) approach implemented by Bioconductor edgeR (v3.18), following the Trimmed Mean of M (TMM) method of normalisation [47].miRNAs were considered to be differentially expressed if the p-value was less than 0.05.Volcano plots of differentially expressed miRNAs (DEmiRs) were generated using the Enhanced Volcano package in R Version 1.10.0[41].

Functional and Pathway Enrichment Analysis
Following miRNA-Seq analysis, functional enrichment analysis was performed through the Database for Annotation, Visualization and Integrated Discovery (DAVID) bioinformatic package (https://david.ncifcrf.gov/;accessed on 1 June 2023) [48,49].Firstly, to gather a list of target genes for significant miRNAs, TargetScan (https://www.targetscan.org/vert_71/; accessed on 1 June 2023) [50] and miRTarbase (https://dianalab.e-ce.uth.gr/tarbasev9; accessed on 1 June 2023) [51] were used to identify strong experimentally validated genes for significant miRNAs with a threshold of FDR < 0.05 to ensure a high confidence.The list of gene targets was added to the DAVID software (v2023q4) to assess significantly enriched pathways.Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was carried out to investigate relationships between the significantly expressed genes and their related pathways.Significance was calculated in DAVID using Fisher's exact test to obtain p-values.Only pathways with a p-value < 0.05 were considered significant.

miRNA RT-qPCR Validation of miRNA-Seq Data
For miRNA analyses, 5 ng/µL of total RNA was reverse transcribed into cDNA using the miRCURY LNA RT Kit (Qiagen, UK) following manufacturer specifications.RT-qPCR was performed using miRCURY LNA SYBR Green PCR Kit (Qiagen, UK) on a LightCycler ® 480 real-time PCR system (Roche Diagnostics, Lucerne, Switzerland), following manufacturer instructions.miRNA primers were obtained from Qiagen (miRCURY LNA primer assays (Qiagen, UK) (Table S1).All miRNA expression was measured in duplicate at CT threshold levels and normalised with the average CT values of a housekeeping control; U6.Values were expressed as fold increase over the corresponding values for control by the 2-∆∆CT method.Two independent experiments were performed, and the average (±SD) results were calculated using GraphPad Prism software (Version Prism 10.1.1)(GraphPad Software, San Diego, CA, USA).Data were expressed as the mean values ± SD and graphed using log scale.Statistical significance was analysed using a student t-test.Differences in the mean were considered statistically significant if p < 0.05.

Descriptive Features of Small RNA-Seq Data
Three small RNA libraries (donor control, donor TGFβ1-treated and donor TGFβ2treated) were sequenced.In total, 10+ million 75 bp reads were obtained with an average of 13.2 million reads per sample ranging from 10 to 15 million reads per sample.Only data with a Q score greater than 30 (>99.9% correct) were utilised in the mapping step of the analysis pipeline.The mapping for each sample was on average 86% and the uniformity of the mapping statistics suggests that the samples are comparable.

Differential Expression of TGFβ1-Responsive miRNAs in TM Cells
To determine changes in miRNA expression in the TM in response to TGFβ1, normal primary human TM cells were treated with 5 ng/mL TGFβ1 and analysed using small RNA-Seq.In total, 107 miRNAs were statistically significantly altered (p < 0.05), with 58 significantly up-regulated and 49 significantly down-regulated.All differentially expressed TGFβ1-induced miRNAs were graphed on a volcano plot (Figure 1) which allows for easy visual identification of miRNAs with large and statistically significant fold changes.The top 30 up-regulated and 30 down-regulated DEmiRs were ranked by fold change and are shown in Tables 1 and 2, respectively.

Functional Enrichment Analysis of the TGFβ1 Differentially Expressed Genes
To determine the function of the TGFβ1-induced DEmiRs and the potential role they play within the TM, KEGG enrichment analysis was performed using DAVID software (https://david.ncifcrf.gov;v2023q4) on the target genes of the significant DEmiRs (p < 0.05).To gather a list of target genes for the significant miRNAs, TargetScan and miRTarbase were used to identify strong experimentally validated genes for significant miRNAs with a threshold of FDR < 0.05 to ensure high confidence.The list of gene targets was added to the DAVID software to assess significantly enriched pathways.In total, 135 unique pathways were found to be associated with TGFβ1-regulated miRNAs.The top 30 KEGG pathways for gene targets of the up-and down-regulated TGFβ1 responsive miRNAs are shown in Figure 2 and Figure 3, respectively.Pathways relating to TGFβ signalling, MAPK signalling, PI3K-Akt signalling, Hippo signalling, Wnt signalling, focal adhesion, and regulation of the actin cytoskeleton are relevant to glaucoma/PXFG pathogenesis and are listed in Table 3 (significantly up-regulated miRNAs) and Table 4 (significantly down-regulated miRNAs) along with the associated miRNAs and gene targets.Pathways relating to TGFβ signalling, MAPK signalling, PI3K-Akt signalling, Hippo signalling, Wnt signalling, focal adhesion, and regulation of the actin cytoskeleton are relevant to glaucoma/PXFG pathogenesis and are listed in Table S2 (significantly up-regulated miRNAs) and Table S3 (significantly down-regulated miRNAs) along with the associated miRNAs and gene targets.
to TGFβ signalling, MAPK signalling, PI3K-Akt signalling, Hippo signalling, Wnt si ling, focal adhesion, and regulation of the actin cytoskeleton are relevant to coma/PXFG pathogenesis and are listed in Table S2 (significantly up-regulated miRN and Table S3 (significantly down-regulated miRNAs) along with the associated miR and gene targets.

Validation of TGFβ1 Responsive DEmiRs by RT-qPCR
To validate the TGFβ1 miRNA-Seq results, the expression of a selected panel of significant TGFβ1 responsive miRNAs was analysed by RT-qPCR (Figure 4).

Validation of TGFβ1 Responsive DEmiRs by RT-qPCR
To validate the TGFβ1 miRNA-Seq results, the expression of a selected panel of significant TGFβ1 responsive miRNAs was analysed by RT-qPCR (Figure 4).

Differential Expression of TGFβ2-Responsive miRNAs in TM Cells
To determine changes in miRNA expression in the TM in response to TGFβ2, normal primary human TM cells were treated with 5 ng/mL TGFβ2 and analysed using small RNA-Seq.In total, 67 miRNAs were statistically significantly altered (p < 0.05), with 38 significantly up-regulated and 28 significantly down-regulated.All differentially expressed TGFβ2-induced miRNAs were graphed on a volcano plot (Figure 5).The upregulated and down-regulated DEmiRs were ranked by fold change and are shown in Tables 3 and 4, respectively.
Cells 2024, 13, x FOR PEER REVIEW 11 of

Differential Expression of TGFβ2-Responsive miRNAs in TM Cells
To determine changes in miRNA expression in the TM in response to TGFβ2, norm primary human TM cells were treated with 5 ng/mL TGFβ2 and analysed using sma RNA-Seq.In total, 67 miRNAs were statistically significantly altered (p < 0.05), with 3 significantly up-regulated and 28 significantly down-regulated.All differentially ex pressed TGFβ2-induced miRNAs were graphed on a volcano plot (Figure 5).The up-reg ulated and down-regulated DEmiRs were ranked by fold change and are shown in Table 3 and 4

Functional Enrichment Analysis of the TGFβ2 Responsive Differentially Expressed Genes
To determine the function of the TGFβ2-induced DEmiRs and the potential role the play within the TM, KEGG enrichment analysis was performed using DAVID software o the target genes of the significant DEmiRs adjusting for multiple corrections (FDR < 0.05 To gather a list of target genes for the significant miRNAs, TargetScan and miRTarbas were used to identify experimentally validated genes for significant miRNAs with threshold of FDR < 0.05 to ensure a high confidence.The genes were added to the DAVI software to assess significant pathways.In total, 58 pathways were found to be associate with TGFβ2-regulated miRNAs.The top 30 KEGG pathways for gene targets of the up and down-regulated TGFβ2 responsive miRNAs are shown in Figure 6 and Figure 7, r spectively.Pathways relating to MAPK signalling, PI3K-Akt signalling, Hippo signallin Wnt signalling, focal adhesion, apoptosis, cellular senescence, and regulation of the acti cytoskeleton are relevant to POAG pathogenesis and are listed in Table S4 (significantl up-regulated miRNAs) and Table S5 (significantly down-regulated miRNAs) along wit the associated miRNAs and gene targets.

Functional Enrichment Analysis of the TGFβ2 Responsive Differentially Expressed Genes
To determine the function of the TGFβ2-induced DEmiRs and the potential role they play within the TM, KEGG enrichment analysis was performed using DAVID software on the target genes of the significant DEmiRs adjusting for multiple corrections (FDR < 0.05).To gather a list of target genes for the significant miRNAs, TargetScan and miRTarbase were used to identify experimentally validated genes for significant miRNAs with a threshold of FDR < 0.05 to ensure a high confidence.The genes were added to the DAVID software to assess significant pathways.In total, 58 pathways were found to be associated with TGFβ2-regulated miRNAs.The top 30 KEGG pathways for gene targets of the upand down-regulated TGFβ2 responsive miRNAs are shown in Figure 6 and Figure 7, respectively.Pathways relating to MAPK signalling, PI3K-Akt signalling, Hippo signalling, Wnt signalling, focal adhesion, apoptosis, cellular senescence, and regulation of the actin cytoskeleton are relevant to POAG pathogenesis and are listed in Table S4 (significantly up-regulated miRNAs) and Table S5 (significantly down-regulated miRNAs) along with the associated miRNAs and gene targets.

Validation of TGFβ2 Responsive DEmiRs by RT-qPCR
To validate the TGFβ2 miRNA-Seq results, the expression of a selected panel of significant TGF2-responsive miRNAs was analysed by RT-qPCR (Figure 8).

Validation of TGFβ2 Responsive DEmiRs by RT-qPCR
To validate the TGFβ2 miRNA-Seq results, the expression of a selected panel of significant TGF2-responsive miRNAs was analysed by RT-qPCR (Figure 8).

Altered miRNA Expression Grouped by miRNA Families
Several miRNAs form a miRNA family which is defined as two or more miRNAs with high sequence similarity and are derived from molecular ancestors during evolution [52].MiRNA family members can target inter-related genes and pathways [52] and two well-studied examples in cancer and fibrosis are the miR-17-92 [53,54] and miR-29 families [55].

miR-17-92 Family
The miR-17-92 family is encoded by the MIR17HG gene, which is transcribed as a polycistronic transcript that produces six mature miRNAs: miR-17, miR-18, miR-19a, miR-19b, miR-20, and miR-92 (Figure 9) [56].Interestingly, we identified the miR-17-92 family as differentially expressed in TM cells following both TGFβ1 and -β2 stimulation.One striking difference we found between the TGFβ1 and TGFβ2 datasets was that TGFβ1 treatment appeared to have a larger effect on the family members.To confirm the differential expression of this family in TM cells, we performed RT-qPCR, treating the cells with either TGFβ1 or TGFβ2.In TGFβ1-treated TM cells, the RT-qPCR results followed the same trend as the miRNA-Seq data, apart from miR-20a-5p which showed donor variability.The RT-qPCR results in TGFβ2-treated TM cells were less significant, which followed the pattern found in the miRNA-Seq data.miR-17-5p and miR-18a-5p were both significantly up-regulated, like the miRNA-Seq data, while the rest of the family were not significantly altered (Table 5 and Figures 9-11).

Altered miRNA Expression Grouped by miRNA Families
Several miRNAs form a miRNA family which is defined as two or more miRNAs with high sequence similarity and are derived from molecular ancestors during evolution [52].MiRNA family members can target inter-related genes and pathways [52] and two wellstudied examples in cancer and fibrosis are the miR-17-92 [53,54] and miR-29 families [55].

miRNA Strands
During miRNA biogenesis, one strand of the miRNA duplex is selectively loaded into the AGO protein where the specificity of the miRISC is determined based on the complementarity between the miRNA and the 3′UTR of the target mRNA [60].The miRNA can originate from the 5′ side of the pre-miRNA to which it is referred to as the "-5p" strand or the 3′ side of the pre-miRNA, to which it is referred to as the "-3p" strand.Typically, one strand is loaded into the AGO protein becoming the functional miRNA or guide strand while the other strand is discarded.This strand is known as the miRNA* or passenger strand [61,62,63].The fact that both strands may be loaded into the AGO protein is becoming more widely studied and either strand may form fully functioning mature miRNAs [64,65].Therefore, it is important to know whether one or both of -5p/-3p strands

miRNA Strands
During miRNA biogenesis, one strand of the miRNA duplex is selectively load into the AGO protein where the specificity of the miRISC is determined based on the co plementarity between the miRNA and the 3′UTR of the target mRNA [60].The miRN can originate from the 5′ side of the pre-miRNA to which it is referred to as the "-5 strand or the 3′ side of the pre-miRNA, to which it is referred to as the "-3p" strand.T ically, one strand is loaded into the AGO protein becoming the functional miRNA or gu strand while the other strand is discarded.This strand is known as the miRNA* or p senger strand [61][62][63].The fact that both strands may be loaded into the AGO protein becoming more widely studied and either strand may form fully functioning mature m NAs [64,65].Therefore, it is important to know whether one or both of -5p/-3p strands a miRNA are playing a functional role in glaucoma pathology.We were interested in derstanding the effects of TGFβ on the expression of the paired -5p and -3p strands of

miRNA Strands
During miRNA biogenesis, one strand of the miRNA duplex is selectively loaded into the AGO protein where the specificity of the miRISC is determined based on the complementarity between the miRNA and the 3 ′ UTR of the target mRNA [60].The miRNA can originate from the 5 ′ side of the pre-miRNA to which it is referred to as the "-5p" strand or the 3 ′ side of the pre-miRNA, to which it is referred to as the "-3p" strand.Typically, one strand is loaded into the AGO protein becoming the functional miRNA or guide strand while the other strand is discarded.This strand is known as the miRNA* or passenger strand [61][62][63].The fact that both strands may be loaded into the AGO protein is becoming more widely studied and either strand may form fully functioning mature miRNAs [64,65].Therefore, it is important to know whether one or both of -5p/-3p strands of a miRNA are playing a functional role in glaucoma pathology.We were interested in understanding the effects of TGFβ on the expression of the paired -5p and -3p strands of the same miRNA in the TM cells.
In our datasets, we discovered both the guide and passenger strands of several miR-NAs were altered by TGFβ treatment.miR-708-3p and -5p are significantly up-regulated in TGFβ1 and TGFβ2 stimulated TM cells (Figure 14).miR-21-3p, also referred to as miR-21*, is significantly up-regulated in both TGFβ1 and TGFβ2 treated cells (Figure 15), while miR-21-5p, referred to as the miR-21 guide, was significantly down-regulated in TGFβ1 treated cells and slightly up-regulated, although not significant, in TGFβ2-treated TM cells (Figure 15).Cells 2024, 13, x FOR PEER REVIEW 18 of 33 in TGFβ1 and TGFβ2 stimulated TM cells (Figure 14).miR-21-3p, also referred to as miR-21*, is significantly up-regulated in both TGFβ1 and TGFβ2 treated cells (Figure 15), while miR-21-5p, referred to as the miR-21 guide, was significantly down-regulated in TGFβ1 treated cells and slightly up-regulated, although not significant, in TGFβ2-treated TM cells (Figure 15).

Discussion
Understanding the role of the TGFβ-induced microRNAome in the trabecular meshwork (TM) provides insight into the molecular pathology of primary open-angle glaucoma (POAG) and pseudo-exfoliation glaucoma (XFG).Determining the specific genes and pathways impacted by dysregulated miRNA expression will support the development of miRNA-based therapeutics for glaucoma [66][67][68].In this study, we have treated human primary TM cells with either TGFβ1 or TGFβ2, as cellular models of POAG (TGFβ2) and XFG (TGFβ1), and using small RNA sequencing, we identified 186 TGFβ1regulated miRNAs in TM cells and 72 TGFβ2-regulated miRNAs suggesting TGFβ1 stimulation had a stronger effect on the miRNA expression profile in TM cells.Increased levels of TGFβ in the anterior segment of the eye induce fibrotic changes in the TM in glaucoma (POAG and XFG) including altered turnover of extracellular matrix (ECM) components, formation of cross-linked actin networks (CLANS), upregulation of alpha-smooth muscle actin (αSMA), aberrant formation of actin stress fibres and epithelial to mesenchymal transition (EMT) [12,17,18,69,70].miR-122-5p was one of the highest upregulated miRNAs in response to TGFβ1 in the TM in our study.This response was also previously reported in the TM, and miR-122 was associated with the regulation of the TGFβ/Smad pathway [71].Within the TGFβ signalling pathway, miR-122-5p has predicted target interactions with TGFβR1, TGFβR2, LTBP1, and SMAD2.The levels of miR-122-5p are significantly elevated in the aqueous humour in XFG [72] in addition to TGFβ1 [21][22][23].miR-122 is one of the most abundant

Discussion
Understanding the role of the TGFβ-induced microRNAome in the trabecular meshwork (TM) provides insight into the molecular pathology of primary open-angle glaucoma (POAG) and pseudo-exfoliation glaucoma (XFG).Determining the specific genes and pathways impacted by dysregulated miRNA expression will support the development of miRNA-based therapeutics for glaucoma [66][67][68].In this study, we have treated human primary TM cells with either TGFβ1 or TGFβ2, as cellular models of POAG (TGFβ2) and XFG (TGFβ1), and using small RNA sequencing, we identified 186 TGFβ1-regulated miR-NAs in TM cells and 72 TGFβ2-regulated miRNAs suggesting TGFβ1 stimulation had a stronger effect on the miRNA expression profile in TM cells.Increased levels of TGFβ in the anterior segment of the eye induce fibrotic changes in the TM in glaucoma (POAG and XFG) including altered turnover of extracellular matrix (ECM) components, formation of cross-linked actin networks (CLANS), upregulation of alpha-smooth muscle actin (αSMA), aberrant formation of actin stress fibres and epithelial to mesenchymal transition (EMT) [12,17,18,69,70].miR-122-5p was one of the highest upregulated miRNAs in response to TGFβ1 in the TM in our study.This response was also previously reported in the TM, and miR-122 was associated with the regulation of the TGFβ/Smad pathway [71].Within the TGFβ signalling pathway, miR-122-5p has predicted target interactions with TGFβR1, TGFβR2, LTBP1, and SMAD2.The levels of miR-122-5p are significantly elevated in the aqueous humour in XFG [72] in addition to TGFβ1 [21][22][23].miR-122 is one of the most abundant miRNAs in the liver, playing an important role in liver fibrosis [73] and other organ fibrosis targeting the TGFβ signalling pathway [74].
There was significant upregulation of miR-182-5p in TGFβ1-treated primary human TM cells.Elevated expression of miR-182-5p was reported during stress-induced premature senescence in cultured TM cells [75].The absolute expression of miR-182-5p in the aqueous humour samples from glaucoma patients was elevated 2-fold [76].An SNP (rs76481776) in the MIR182 gene was associated with POAG in the NEIGHBORHOOD GWAS dataset although the mechanism linking this SNP with POAG has not been elucidated [76].Originally described as a sensory organ-specific miRNA, miR-182 is also involved in immunity, cancer and regulation of TGFβ signalling [77][78][79].By targeting SMAD7, a negative regulator of the TGFβ signalling pathway, miR-182-5p amplifies TGFβ induced epithelial to mesenchymal transition (EMT) and metastasis of cancer cells while inhibition of miR-182-5p reduces pulmonary fibrosis [80,81].
miR-145-5p is abundantly expressed in the TM and smooth muscle in the eye and is highly expressed in the aqueous humour in POAG patients [67,82,83].It is a member of the miR-143/145 cluster which plays a role in the regulation of IOP [67].miR-143/145 double knockout mice resulted in a 19% decrease in IOP [67].miR-143/145 increases IOP by modulating actin dynamics enhancing the contractility of TM cells [67].Manipulation of miR-143/145 levels in the TM may offer therapeutic potential in glaucoma [67].
The expression of miR-146b-5p was downregulated with TGFβ treatment.miR-146b-5p is a member of the miR-146 family of miRNAs, consisting of miR-146a-5p and miR-146b-5p.These two miRNAs only differ by two nucleotides on the 3 ′ end of their mature strand, sharing the same seed region [84].During replicative senescence in human TM cells, miR-146a upregulation limited inflammatory responses [85].miR-146a regulates the pro-inflammatory NF-κB signalling pathway by inhibiting interleukin-1 receptor-associated kinase 1 (IRAK1) to inhibit inflammation [86].Following lentiviral intracameral delivery of miR-146a in rats, there was a sustained reduction of IOP of 4.4 ± 2.9 mmHg over 8 months [66].The mechanism of IOP lowering was postulated as likely complex and potentially involving alterations in TGFβ signalling, ROCK inhibition, and/or NF-κB signalling [66].While miR-146b-5p has not been studied in the TM, it also inhibited NF-κB-induced interleukin 6 (IL-6) expression in breast cancer cells [87].This suggests that miR-146b-5p may play a similar role to miR-146a-5p in the TM.

miRNA Clusters and Families
The action of miRNAs can be synergistic, and this is exemplified by miRNA clusters which consist of multiple miRNAs with a common promoter resulting in co-expression and coordinated action [62,88,89].The miR-17-92 cluster was enriched in both the TGFβ1 and TGFβ2 datasets.This cluster consists of six mature miRNAs: miR-17, miR-18a, miR-19a, miR-19b, miR-20a, and miR-92a [53].Through gene duplication, this cluster has evolved to form two paralogs: the miR-106a-363 cluster and the miR-106b-25 cluster, shown in Figure 9 [54].As some of the miRNAs share a seed sequence, they have been divided into four main miRNA families: the miR-17, miR-18, miR-19, and miR-92 families [54].An important pathway targeted by members of the miR-17-92 family is the TGFβ signalling pathway [90,91].miR-17-5p is down-regulated in response to TGFβ1 treatment in TM cells.TM cells under oxidative stress down-regulate miR-17-5p, which may regulate the proliferation and apoptosis of TM cells through its direct targeting of tumour suppressor PTEN [92], which is up-regulated in TM cells following TGFβ treatment [93].Previous research from our group has shown miR-18a-5p expression increased in TM cells following TGFβ2, consistent with our miRNA-Seq results [94].miR-18a-5p targets connective tissue growth factor (CTGF), which is a fibrotic gene elevated in the TM of glaucoma patients [95].CTGF induces actin stress fibres and increases TM cell contractility by activating RhoA [96].Lentiviral-mediated overexpression of miR-18a reduced TGFβ2-induced CTGF expression in TM cells and showed a reduction in TGFβ2-induced contraction of collagen gels [94].miR-18a-5p is a potential miRNA therapeutic in glaucoma because of its ability to inhibit CTGF-associated increased TM cell contractility [94].
The synergistic action of miRNAs is also supported by miRNA families [89,97].A miRNA family consists of two or more miRNAs with high sequence similarity and can be located in one or more distinct clusters [89].The miR-29 family regulates a plethora of fibrosis-associated genes in various cell types including lungs, liver, heart, eye, and other organs [55,98].The family is a known downstream target in the TGFβ/Smad pathway, and the phosphorylation of Smad3 by TGFβ causes miR-29 to be downregulated [55].Our study detected down-regulation of miR-29b-3p in both TGFβ1 and TGFβ2 treated TM cells, as seen in previous reports [99,100].Transfection of human TM cells with miR-29b-3p mimic down-regulated ECM proteins including collagens, laminin subunit gamma 1 (LAMC1), and fibrillin 1 (FBN1), and secreted protein acidic rich in cysteine (SPARC), a gene involved in ECM remodelling [101].Alterations in SPARC, and collagens I and IV, cause significant changes in IOP in transgenic mice [102].Our KEGG analysis identified the PI3K-Akt signalling pathway to be over-represented with miR-29b-3p expression.Overexpression of miR-29b-3p represses the PI3K-Akt pathway reducing collagen I expression in human Tenon's ocular fibroblasts [103].miR-29b-3p down-regulation in the TM may contribute to increased TGFβ-induced ECM components.miR-29 plays a critical role in regulating ECM production and is an anti-fibrotic miRNA [104] and miRNA-29b mimics attenuate pulmonary fibrosis in vivo [105].

miRNA Strands
The miRNA biogenesis pathway involves the sequential processing of the pri-miRNA into pre-miRNA and finally into a mature miRNA [106].Following Drosha processing of pri-miRNAs in the nucleus, pre-miRNAs are exported into the cytoplasm and cleaved by Dicer, resulting in the miRNA duplex [63].Transcription produces equal amounts of both strands of miRNA duplexes; however, their accumulation is mostly asymmetric at steady state [107].As proposed in the oncology field [108], the results from our study highlight that both 3p-and 5p-arms from a miRNA warrant independent study.miR-21-5p is known to be one of the most overexpressed miRNAs in response to tissue injury and to play an important role in fibrosis [109,110].During miRNA biogenesis, pre-miR-21 is exported by Exportin 5 and processed by Dicer to release mature hsa-miR-21 (also known as hsa-miR-21-5p, the biologically dominant arm) and hsa-miR-21-3p (formerly named hsa-miR-21*), previously considered the less abundant or active miRNA strand [111][112][113][114]. Pro-fibrotic miR-21-5p binds to Smad7, an inhibitory Smad, and thus amplifies the TGFβ signalling pathway, causing fibrotic responses [115,116].There is crosstalk between miR-21-5p and a variety of signalling pathways: TGFβ/SMAD, PI3K/AKT and ERK/MAPK signalling pathways, in the regulation of fibrotic processes [117].A role of miR-21-5p in regulating IOP and outflow facility has been reported [118].Topical administration of a synthetic miR-21-5p mimic increased miR-21-5p expression in the TM while reducing IOP by 17% [118].Using RNA-sequencing and pathway analysis, with the predicted downstream target genes of miR-21-5p identified, they found a pathway involving FGF18, SMAD7, and MMP9 based on protein-protein interaction networks [118].RT-qPCR confirmed the downregulation of SMAD7 and FGF18 by a miR-21-5p mimic suggesting that miR-21-5p targets SMAD7 and FGF18 to encourage ECM degradation by MMP9 in the TM [118].In our data miR-21-5p expression was unaltered by TGFβ2 treatment in the TM but significantly downregulated with TGFβ1 treatment.Unexpectedly, our data show an upregulation of miR-21-3p with both TGFβ1 and -β2 treatment.There is emerging evidence for a biological role for miR-21-3p in malignancy [108,119], vascular biology [120,121] and in the regulation of TGFβ signalling [111].In hepatocellular carcinoma, miR-21-3p regulates both TGFβ and Hippo signalling via SMAD7 and YAP1 [111].Overexpression of miR-21-3p directly silences SMAD7 expression and reduces the stability of the SMAD7/YAP1 complex allowing YAP1 translocation to the nucleus and resultant profibrotic gene expression [111].Further work is required to understand the role of the -5p and -3p miR-21 strands in TM pathophysiology.miR-708 is not a widely studied miRNA in ocular tissues; however, its expression has been reported in retinal ganglion cells [122].miR-708-5p is the more abundant strand and is involved in oncogenesis [123,124].miR-708-5p was upregulated in the TM in response to TGFβ2 but the passenger strand (miR-708-3p) was significantly upregulated in response to TGFβ1 and -β2.A disintegrin and metalloproteinase 17 (ADAM17), which is overexpressed in fibrotic disorders [125][126][127] and is expressed in TM cells [128,129], is a direct target of miR-708-3p.By targeting ADAM17, miR-708-3p represses the GATA/STAT3 signalling pathway in idiopathic pulmonary fibrosis (IPF) reducing fibrosis [130].GATA6 promotes fibroblast differentiation into myofibroblasts in IPF by mediating the α-SMA-inducing signal of TGFβ1 [131,132], and STAT3, is abundantly expressed in multiple fibrotic disorders [133].In breast cancer cells, miR-708-3p inhibits EMT by directly targeting ZEB1, cadherin 2, and vimentin [134].Therefore, miR-708-3p could present a new therapeutic target for TM fibrosis by targeting the ADAM17-GATA/STAT pathway and EMT.

miRNA Regulation of Signalling Pathways
A cell can simultaneously express multiple miRNAs to regulate gene expression in a holistic, intricate network with a single miRNA targeting multiple mRNAs and a single mRNA targeted by multiple miRNAs [31,62,88].Our data can be considered in terms of specific miRNAs (miRNA-centric) but miRNAs work in networks to control cellular pathways and processes and a pathway-centric view of miRNA action is also required [135].Several of the DEmiRs altered by TGFβ in the TM have been implicated in the regulation of the TGFβ signalling pathway [90,91,136].Interestingly, the TGFβ signalling pathway was only significantly enriched in TGFβ1 down-regulated miRNAs in the TM.The transforming growth factor beta receptor 2 (TGFβR2) is targeted by miR-18a-5p [137], miR-20a-5p, miR-29b-3p, and miR-204-5p which were all significantly downregulated in the TM in response to TGFβ1.Smad2/3, the regulatory Smads in the canonical TGFβ signalling pathway, are also targeted by miR-18a-5p [138][139][140][141]. Lentiviral expression of miR-18a-5p in bleomycin mice presented lowered levels of phosphorylated Smad2/3 and a reduction in pulmonary fibrosis [137].There are several overlapping miRNAs involved in the regulation of both the TGFβ signalling pathway and the Hippo signalling pathway.
The Hippo signalling pathway is considered a tumour suppressor pathway, playing important roles in cell differentiation and cell proliferation [142,143].Our pathway enrichment analysis identified the Hippo signalling pathway as significantly enriched in both TGFβ1 and -β2 miRNA-Seq datasets in the TM, associated with both up-and downregulated miRNAs.TGFβ2 upregulated the expression of miR-181c-5p in TM cells and miR-181c-5p inhibits Hippo signalling through its target large tumour suppressor 2 (LATS2) and Salvador Family WW Domain Containing Protein 1 (SAV1) [144][145][146].SAV1 bind to the Macrophage-Stimulating 1/2 (MST1/2) kinases forming an enzyme complex which phosphorylates LATS1/2 [147,148].The LATS1/2 kinases phosphorylate Yes-Associated Protein 1 (YAP), which in turn prevents nuclear translocation and signals its proteasomal degradation [149].Therefore, miR-181c-5p disrupts the negative regulation of Hippo signalling and promotes the activation of YAP by targeting LATS2 and SAV1 [144,145].YAP/TAZ also stimulates the nuclear accumulation of SMAD complexes to increase their transcriptional activity [150][151][152].YAP/TAZ nuclear levels are elevated by TGFβ2 in both normal and glaucomatous human TM cells while inhibition of YAP/TAZ resulted in reduced focal adhesions, ECM remodelling and cell contractility [153].When Hippo signalling is inhibited the nuclear translocation of YAP/TAZ induces the transcriptional activity of TEA domain (TEAD) family members, which increases the expression of CTGF which is associated with glaucoma and TM fibrosis [94][95][96].
Inhibition of Wnt signalling has previously been shown in glaucomatous TM cells and is associated with TM cell stiffening [191,192].Secreted frizzled-related protein-1 (SFRP1), targeted by miR-582-3p [193] which was down-regulated in TGFβ2-treated TM cells, inhibits the Wnt signalling pathway and is associated with increased IOP [194].Activating Wnt signalling through miRNA targets could restore the normal phenotype caused by Wnt inactivation, through repression of ECM genes (SPARC and CTGF), crosslinking genes (LOX), and inhibitors of MMPs (TIMP1 and PAI-1) [195,196].

miRNA-Based Therapeutics
miRNA-based therapeutics offer translational benefits compared to other treatment modalities [197][198][199][200]. Mature miRNAs are highly conserved across vertebrate species, facilitating the evaluation of identical miRNA-based therapeutics in preclinical efficacy, safety, and pharmacodynamics, and subsequent human clinical trials [198].Intracameral delivery of miR-200c and miR-146a and topical application of miR-21-5p in rodent eyes lowers IOP and raises the possibility for the clinical translation of miRNA-based therapeutics for glaucoma [66,68,118].Mature miRNA sequences are naturally occurring, endogenous short molecules (approximately 19-25 nucleotides) allowing several approaches to modulate miRNA levels as a therapeutic intervention [197,198].MicroRNA mimics can restore miRNA levels and activity while mature miRNAs can be inhibited using miRNA sponges or antisense oligonucleotides (antimiRs) [197,198].Target site masking antisense oligonucleotides can mask a specific miRNA binding site in a target mRNA/gene [197].Chemical modification of miRNA-based therapeutics improves their pharmacokinetics and pharmacodynamics [197].
TGFβ signalling can be targeted using small molecule inhibitors of the TGF-β receptor kinases, neutralising antibodies to disrupt ligand-receptor interactions and antisense oligonucleotides (e.g., siRNA) [201,202].TGF-β is involved in many physiological processes including cell survival, metabolism, growth, proliferation, differentiation, adhesion, migration, and death and as such is central to tissue and immune homeostasis [203].Therefore, the critical function of TGFβ in maintaining tissue homeostasis makes targeting TGFβ a challenge and systemic inhibition of TGF-β signalling might evoke serious on-target side effects [201,[203][204][205]. Gene targeting by miRNAs involves a seed regionmediated miRNA::mRNA interaction [88].This is an important attribute that can be exploited in miRNA-based therapeutics as a single miRNA can target multiple genes with relatively weak gene suppression of 30-60% [206][207][208].This aggregate effect on multiple genes can influence gene networks and pathways [200] making miRNA-based therapeutics superior to other approaches targeting single genes or proteins (e.g., siRNAs or small molecule inhibitors) as miRNAs can regulate complex biological processes and networks [199].Rather than targeting a single gene or inhibiting a pathway like TGFβ with deleterious effects, miRNA-based therapeutics can fine-tune signalling pathways dampening responses [198][199][200].In our data, we have shown the enriched pathways and the key miRNAs and their mRNA targets in TM cells stimulated with either TGF-β1 or -β2.This is important for determining which gene or pathway to target and with which miRNA or miRNAs because miRNA action is cell and context-specific.These data are also important for determining on-target specificity and undesired off-target and on-target effects [197,200].Future studies assessing the miRNA::mRNA interactome using in silico and in vitro approaches can help characterise the "miRNA targetome" to support the development of miRNA-based therapeutics [200].
The delivery of miRNA therapeutics to the TM also poses challenges including penetration of the corneal barrier, prevention of premature degradation, and specific targeting of TM cells.A solution to these issues already exists in the use of nanoparticle delivery systems [209,210].Nanoparticles and nanocarriers can increase the stability of the miRNA molecule before being absorbed by the cell, essentially acting like a vesicle.Furthermore, nanoparticles can be coated to allow more specific targeting of a particular tissue [209].A TM-targeting coating for nanoparticles has already been developed using hyaluronan, which has been used to deliver a CTGF-targeting siRNA to the TM of ex vivo mouse, porcine, and human eyes [209].This specific targeting of the TM would also help prevent off-target effects in the eye by using cell-specific miRNA modulation [197].

Conclusions
To our knowledge, this is the first study reporting a genome-wide miRNA expression in primary normal human TM cell samples following either TGFβ1 or -β2 stimulation using miRNA-Seq.We have identified differentially expressed miRNAs that could target genes associated with the pathogenesis of POAG and XFG, and the identified enriched pathways associated with these miRNAs.Access to diseased TM disease from glaucoma patients is limited but would add to our understanding of miRNA dysregulation in the TM and disease pathogenesis.Further work to identify the miRNA::mRNA interactome is required to fully understand the TGFβ induced microRNAome of the TM.This knowledge is also required to develop miRNA-based therapeutics for POAG and XFG which will require consideration of miRNA-centric and pathway-centric effects.There are several challenges to overcome to develop miRNA therapeutics to reduce IOP in glaucoma [41], but early pre-clinical studies are encouraging [66,68,118,210].

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/cells13121060/s1,Table S1.RT-qPCR miRCURY LNA miRNA primer assays sequences.Table S2.Significantly up-regulated TGFβ1-responsive miRNAs associated with glaucoma-related signalling pathways and their gene targets.Table S3.Significantly downregulated TGFβ1-responsive miRNAs associated with glaucoma-related signalling pathways and their gene targets.Table S4.Significantly up-regulated TGFβ2-reponsive miRNAs associated with glaucoma-related signalling pathways and their gene targets.Table S5.Significantly down-regulated TGFβ2-responsive miRNAs associated with glaucoma-related signalling pathways and their gene targets.Table S6.Significantly down-regulated TGFβ2-responsive miRNAs associated with glaucomarelated signalling pathways and their gene targets.

Figure 1 .
Figure 1.Significant differential expression of miRNAs in response to TGFβ1 treatment (5 ng/mL 24 h).Volcano plot identifying differentially expressed miRNAs in response to TGFβ1 treatment in TM cells.A threshold was applied using a statistical significance of p < 0.05 and log2FC > 1. Grey and green dots represent miRNAs failing to meet the significance cut-off.Blue dots represent miRNAs that have met the p-value but are below the log2FC criteria.Red dots represent those miRNAs that have met the threshold with a p < 0.05 and log2FC > 1.Those on the left corner of the plot are downregulated whilst those on the right area of the plot are the up-regulated miRNAs.Grey dots = Not Significant (NS), Red dots = p-value < 0.05 & Log2FC ≥ 1, Blue dots = p-value < 0.05.p-value is represented by the dashed horizontal line and Log2FC is represented by the vertical dashed lines.

Figure 1 .
Figure 1.Significant differential expression of miRNAs in response to TGFβ1 treatment (5 ng/mL 24 h).Volcano plot identifying differentially expressed miRNAs in response to TGFβ1 treatment in TM cells.A threshold was applied using a statistical significance of p < 0.05 and log 2 FC > 1. Grey and green dots represent miRNAs failing to meet the significance cut-off.Blue dots represent miRNAs that have met the p-value but are below the log 2 FC criteria.Red dots represent those miRNAs that have met the threshold with a p < 0.05 and log 2 FC > 1.Those on the left corner of the plot are downregulated whilst those on the right area of the plot are the up-regulated miRNAs.Grey dots = Not Significant (NS), Red dots = p-value < 0.05 & Log 2 FC ≥ 1, Blue dots = p-value < 0.05.p-value is represented by the dashed horizontal line and Log 2 FC is represented by the vertical dashed lines.

Figure 2 .
Figure 2. KEGG pathway analysis of significantly up-regulated TGFβ1-responsive miRNAs i cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) wa formed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differen expressed miRNAs up-regulated in response to TGFβ1 in TM cells.The X-axis label represen richment Score, the amount of differentially expressed genes enriched in the pathway, and Y label represents the enriched pathway.The size and colour of the bubble represent the amou differentially expressed genes enriched in the pathway and the enrichment significance, re tively.

Figure 2 .
Figure 2. KEGG pathway analysis of significantly up-regulated TGFβ1-responsive miRNAs in TM cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) was performed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differentially expressed miRNAs up-regulated in response to TGFβ1 in TM cells.The X-axis label represents Enrichment Score, the amount of differentially expressed genes enriched in the pathway, and Y-axis label represents the enriched pathway.The size and colour of the bubble represent the amount of differentially expressed genes enriched in the pathway and the enrichment significance, respectively.

Figure 3 .
Figure 3. KEGG pathway analysis of significantly down-regulated TGFβ1-responsive miRNAs in TM cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) was performed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differentially expressed miRNAs down-regulated in response to TGFβ1 in TM cells.The X-axis label represents Enrichment Score, the amount of differentially expressed genes enriched in the pathway, and Y-axis label represents the enriched pathway.The size and colour of the bubble represent the amount of differentially expressed genes enriched in the pathway and the enrichment significance, respectively.

Figure 3 .
Figure 3. KEGG pathway analysis of significantly down-regulated TGFβ1-responsive miRNAs in TM cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) was performed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differentially expressed miRNAs down-regulated in response to TGFβ1 in TM cells.The X-axis label represents Enrichment Score, the amount of differentially expressed genes enriched in the pathway, and Y-axis label represents the enriched pathway.The size and colour of the bubble represent the amount of differentially expressed genes enriched in the pathway and the enrichment significance, respectively.

Figure 4 .
Figure 4. Expression levels of four candidate miRNAs in TGFβ1-stimulated TM cells.Candidate miRNAs were identified as significantly differentially expressed through miRNA-Seq.(A) hsa-miR-122-5p (B) hsa-miR-146b-5p (C) hsa-miR-182-5p (D) hsa-miR-204-5p.Vehicle controls are denoted as "Control" on graphs and TGFβ1 treated cells as "TGFβ1".Individual values for donors 1, 2, and 3 are shown.Data were normalised to U6 control and analysed using the ΔΔCt method.Circles represent individual donor gene without TGFβ1 treatment and squares represent donors with TGFβ1 treatment.The colour of circle or square represents the individual donor human primary TM cells.An asterisk denotes significant differential gene expression after treatment (* = p < 0.05, ** = p < 0.005).

Figure 4 .
Figure 4. Expression levels of four candidate miRNAs in TGFβ1-stimulated TM cells.Candidate miRNAs were identified as significantly differentially expressed through miRNA-Seq.(A) hsa-miR-122-5p (B) hsa-miR-146b-5p (C) hsa-miR-182-5p (D) hsa-miR-204-5p.Vehicle controls are denoted as "Control" on graphs and TGFβ1 treated cells as "TGFβ1".Individual values for donors 1, 2, and 3 are shown.Data were normalised to U6 control and analysed using the ∆∆Ct method.Circles represent individual donor gene without TGFβ1 treatment and squares represent donors with TGFβ1 treatment.The colour of circle or square represents the individual donor human primary TM cells.An asterisk denotes significant differential gene expression after treatment (* = p < 0.05, ** = p < 0.005).

Figure 5 .
Figure 5. Significant differential expression of miRNAs in response to TGFβ2 treatment (5 ng/mL 2 h).Volcano Plot identifying differentially expressed miRNAs in response to TGFβ2 treatment in TM cells.A threshold was applied using a statistical significance of p < 0.05 and log2FC > 1. Grey an green dots represent miRNAs failing to meet the significance cut-off.Blue dots represent miRNA that have met the p-value but are below the log2FC criteria.Red dots represent those miRNAs th have met the threshold with a p < 0.05 and log2FC > 1.Those on the left corner of the plot are down regulated whilst those on the right area of the plot are the up-regulated miRNAs.Grey dots-N Significant (NS), Red dots = p-value < 0.05 & Log2FC ≥ 1, Blue dots = p-value < 0.05.p-value is repr sented by the dashed horizontal line and Log2FC is represented by the vertical dashed lines.

Figure 5 .
Figure 5. Significant differential expression of miRNAs in response to TGFβ2 treatment (5 ng/mL 24 h).Volcano Plot identifying differentially expressed miRNAs in response to TGFβ2 treatment in TM cells.A threshold was applied using a statistical significance of p < 0.05 and log 2 FC > 1. Grey and green dots represent miRNAs failing to meet the significance cut-off.Blue dots represent miRNAs that have met the p-value but are below the log 2 FC criteria.Red dots represent those miRNAs that have met the threshold with a p < 0.05 and log 2 FC > 1.Those on the left corner of the plot are downregulated whilst those on the right area of the plot are the up-regulated miRNAs.Grey dots-Not Significant (NS), Red dots = p-value < 0.05 & Log 2 FC ≥ 1, Blue dots = p-value < 0.05.p-value is represented by the dashed horizontal line and Log 2 FC is represented by the vertical dashed lines.

24 ,Figure 6 .
Figure 6.KEGG pathway analysis of significantly up-regulated TGFβ2-responsive miRNAs i cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) wa formed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differen expressed miRNAs up-regulated in response to TGFβ2 in TM cells.The X-axis label represen richment Score, the amount of differentially expressed genes enriched in the pathway, and Y label represents the enriched pathway.The size and colour of the bubble represent the amou differentially expressed genes enriched in the pathway and the enrichment significance.

Figure 6 .
Figure 6.KEGG pathway analysis of significantly up-regulated TGFβ2-responsive miRNAs in TM cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) was performed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differentially expressed miRNAs up-regulated in response to TGFβ2 in TM cells.The X-axis label represents Enrichment Score, the amount of differentially expressed genes enriched in the pathway, and Y-axis label represents the enriched pathway.The size and colour of the bubble represent the amount of differentially expressed genes enriched in the pathway and the enrichment significance.

Figure 7 .
Figure 7. KEGG pathway analysis of significantly down-regulated TGFβ2-responsive miRNAs in TM cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) was performed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differentially expressed miRNAs down-regulated in response to TGFβ2 in TM cells.The X-axis label represents Enrichment Score, the amount of differentially expressed genes enriched in the pathway, and Y-axis label represents the enriched pathway.The size and colour of the bubble represent the amount of differentially expressed genes enriched in the pathway and the enrichment significance.

Figure 7 .
Figure 7. KEGG pathway analysis of significantly down-regulated TGFβ2-responsive miRNAs in TM cells.KEGG pathway analysis of the significant differentially expressed miRNAs (p < 0.05) was performed.The enrichment plot presents the top 30 KEGG pathway enrichment terms for differentially expressed miRNAs down-regulated in response to TGFβ2 in TM cells.The X-axis label represents Enrichment Score, the amount of differentially expressed genes enriched in the pathway, and Y-axis label represents the enriched pathway.The size and colour of the bubble represent the amount of differentially expressed genes enriched in the pathway and the enrichment significance.

Figure 8 .
Figure 8. Expression levels of five candidate miRNAs in TGFβ2-stimulated TM cells.Candidate miRNAs were identified as significantly differentially expressed through miRNA-Seq.(A) hsa-miR-21-3p (B) hsa-miR-29b-3p (C) hsa-miR-145-5p (D) hsa-miR-204-5p.Vehicle controls are denoted as "Control" on graphs and TGFβ2 treated cells as "TGFβ2".Individual values for donors 1, 2, and 3 are shown.Data were normalised to U6 control and analysed using the ΔΔCt method.Circles represent individual donor gene without TGFβ2 treatment and squares represent donors with TGFβ2 treatment.The colour of circle or square represents the individual donor human primary TM cells.An asterisk denotes significant differential gene expression after treatment (* = p < 0.05, ** = p < 0.005).

Figure 8 .
Figure 8. Expression levels of five candidate miRNAs in TGFβ2-stimulated TM cells.Candidate miRNAs were identified as significantly differentially expressed through miRNA-Seq.(A) hsa-miR-21-3p (B) hsa-miR-29b-3p (C) hsa-miR-145-5p (D) hsa-miR-204-5p.Vehicle controls are denoted as "Control" on graphs and TGFβ2 treated cells as "TGFβ2".Individual values for donors 1, 2, and 3 are shown.Data were normalised to U6 control and analysed using the ∆∆Ct method.Circles represent individual donor gene without TGFβ2 treatment and squares represent donors with TGFβ2 treatment.The colour of circle or square represents the individual donor human primary TM cells.An asterisk denotes significant differential gene expression after treatment (* = p < 0.05, ** = p < 0.005).

Figure 13 .
Figure 13.Expression levels of miR-29 family in TGFβ2-stimulated TM cells.(A) hsa-miR-29a-3p (B) hsa-miR-29b-3p (C) hsa-miR-29c-3p.Vehicle controls are denoted as "Control" on graphs and TGFβ2 treated cells as "TGFβ2".Individual values for donors 1, 2, and 3 are shown.Data were normalised to U6 control and analysed using the ΔΔCt method.Circles represent individual donor gene without TGFβ2 treatment and squares represent donors with TGFβ2 treatment.The colour of circle or square represents the individual donor human primary TM cells.An asterisk denotes significant differential gene expression after treatment (** = p < 0.005).

Figure 13 .
Figure 13.Expression levels of miR-29 family in TGFβ2-stimulated TM cells.(A) hsa-miR-29a-3p hsa-miR-29b-3p (C) hsa-miR-29c-3p.Vehicle controls are denoted as "Control" on graphs TGFβ2 treated cells as "TGFβ2".Individual values for donors 1, 2, and 3 are shown.Data were n malised to U6 control and analysed using the ΔΔCt method.Circles represent individual donor g without TGFβ2 treatment and squares represent donors with TGFβ2 treatment.The colour of ci or square represents the individual donor human primary TM cells.An asterisk denotes signific differential gene expression after treatment (** = p < 0.005).

Figure 13 .
Figure 13.Expression levels of miR-29 family in TGFβ2-stimulated TM cells.(A) hsa-miR-29a-3p (B) hsa-miR-29b-3p (C) hsa-miR-29c-3p.Vehicle controls are denoted as "Control" on graphs and TGFβ2 treated cells as "TGFβ2".Individual values for donors 1, 2, and 3 are shown.Data were normalised to U6 control and analysed using the ∆∆Ct method.Circles represent individual donor gene without TGFβ2 treatment and squares represent donors with TGFβ2 treatment.The colour of circle or square represents the individual donor human primary TM cells.An asterisk denotes significant differential gene expression after treatment (** = p < 0.005).

Table 5 .
Fold Change and p-value of the miR-17-92 family members in the miRNA-Seq datasets.

Table 5 .
Fold Change and p-value of the miR-17-92 family members in the miRNA-Seq datasets.

Table 5 .
Fold Change and p-value of the miR-17-92 family members in the miRNA-Seq datasets.