Skip to main content
Download PDF

In silicoanalysis of the molecular machinery underlying aqueous humor production: potential implications for glaucoma

Journal of Clinical Bioinformatics volume 3, Article number: 21 (2013) Cite this article

Abstract

Background

The ciliary body epithelia (CBE) of the eye produce the aqueous humor (AH). The equilibrium between the AH production by the CBE and the outflow through the trabecular meshwork ultimately determines the intraocular pressure (IOP). An increased IOP is a major risk factor for primary open angle glaucoma (POAG). This study aims to elucidate the molecular machinery of the most important function of the CBE: the AH production and composition, and aims to find possible new molecular clues for POAG and AH production-lowering drugs.

Methods

We performed a gene expression analysis of the non-pigmented (NPE) and pigmented epithelia (PE) of the human CBE of post mortem eyes. We used 44 k Agilent microarrays against a common reference design. Functional annotations were performed with the Ingenuity knowledge database.

Results

We built a molecular model of AH production by combining previously published physiological data with our current genomic expression data. Next, we investigated molecular CBE transport features which might influence AH composition. These features included caveolin- and clathrin vesicle-mediated transport of large biomolecules, as well as a range of substrate specific transporters. The presence of these transporters implies that, for example, immunoglobins, thyroid hormone, prostaglandins, cholesterol and vitamins can be secreted by the CBE along with the AH. In silico, we predicted some of the molecular apical interactions between the NPE and PE, the side where the two folded epithelia face each other. Finally, we found high expression of seven POAG disease genes in the plasma membrane of extracellular space of the CBE, namely APOE, CAV1, COL8A2, EDNRA, FBN1, RFTN1 and TLR4 and we found possible new targets for AH lowering drugs in the AH.

Conclusions

The CBE expresses many transporters, which account for AH production and/or composition. Some of these entries have also been associated with POAG. We hypothesize that the CBE may play a more prominent role than currently thought in the pathogenesis of POAG, for example by changing the composition of AH.

Background

The ciliary body epithelia (CBE) in the eye, consisting of the non-pigmented (NPE) and pigmented epithelia (PE), are responsible for the production of the aqueous humor (AH). This production includes passive diffusion, active transport and production of molecules from the CBE into the AH. Via the anterior chamber, the AH exits the eye trough the trabecular meshwork (TM) and the canal of Schlemm into the venous blood system. The balance between the production and outflow of AH ultimately determines the intraocular pressure (IOP).

In many physiological in vitro and in vivo animal studies, the major ion-channel and transporter proteins involved in AH production have already been identified. In short, the chain of events leading to AH production starts with the PE cells that actively take up ions at their stromal (basolateral) surface (major channels Na+-K+-2Cl- symporters, Cl--HCO3- and Na+-H+ antiporters, K+ channels). These ions flow passively through gap junctions between the PE and NPE cell. The NPE cells actively secrete ions into the AH (major channels Na+-K+-activated ATPase, H+-ATPase, Na+-K+-2Cl- symport, K+ and Cl- channels). Water follows by passive diffusion through water channels (aquaporins). Together, the studies on ion-channels led to a physiological model of AH production in the CBE (reviewed in [1, 2]). However, the specific genes that code for these ion channels are still unknown for many of these channels. The first aim of our current study was to investigate the human molecular background of the physiological AH production model.

Besides the ion-channels directly involved in AH production, the CBE expresses also other transporters [3, 4]. These transporters may play a role in the composition of the AH, which determines a range of other functionalities previously assigned to the CBE. These are, for example, nourishment of avascular tissues in the eye (lens, cornea), neuro-endocrine signaling to the TM and other tissues and maintenance of the immune privilege of the eye. Our second aim of this study was to investigate the genetic identity of these transporters in the CBE.

The NPE and PE are folded, polarized neuroepithelial cell layers that face each other at their apical sides. At this side of the NPE and PE, secreted and plasma membrane biomolecules of both epithelia can interact. These biomolecules include signaling molecules, growth factors and extracellular matrix molecules. We are interested in the molecular characteristics of this apical interaction between NPE and PE.

In disease state, like primary open angle glaucoma (POAG), AH composition may be disturbed [5–7]. POAG is a neurodegenerative disease characterized by abnormal cupping of the optic nerve head and corresponding visual field defects that cannot be explained by other diseases, and with a normal open anterior chamber of the eye [8]. A major risk factor for POAG is an increased intraocular pressure (IOP) [9, 10]. Proteomic analysis of POAG AH showed an increased level of total protein content and also the presence of larger and potentially harmful molecules compared to control AH [5–7]. This changed AH composition in POAG may ultimately affect the outflow facility of the TM via intrinsic molecular and structural changes, which will result in a decreased outflow and an increased IOP. The fourth aim of our study was to look for expression of POAG disease genes in the CBE, and we hypothesize about their potential role in the CBE in POAG when these genes are mutated.

The only proven therapy of POAG is lowering of the IOP. Drugs that lower the IOP target the CBE to decrease AH production or target the TM to increase AH outflow. With our newly formed hypothetical molecular model of AH production in the CBE (first aim of this study), we finally looked for potentially new targets to change the AH production with already existing drugs.

Methods

Ethics statement

This study was performed in agreement with the declaration of Helsinki on the use of human material for research. The human post-mortem donor eyes were obtained from the Corneabank, Beverwijk, the Netherlands. The Corneabank obtained permission (informed consent) from the donors for enucleation of the eyes and to use the eye for scientific purposes after removal of the cornea. All data were analyzed anonymously.

Tissue sampling, RNA processing and microarray

We performed gene expression analysis on separately laser dissected NPE (n = 7) and PE (n = 7) cells (only pars plicata of the CBE) with microarray studies. For detailed description of donor eye selection, tissue sampling, RNA processing and microarray performance, we refer to our previously published paper [11]. In short, we selected seven healthy human donor eyes that were snap frozen. From 20 μm cryosections, the NPE and PE were separately cut out with laser dissection microscopy (PALM Carl Zeiss, Microlaser Technologies AG, Germany). Cresyl Violet staining was used to distinguish NPE (Merck, Frankfurt, Germany, art.5235). After selection of NPE and PE cells, RNA was isolated, amplified and labeled. Microarrays were performed against common reference sample (RPE/choroid RNA), in order to compare NPE and PE and normalize the data. The total gene expression datasets can be found in Gene Expression Omnibus database (GSE37957).

Data analysis I: construction of sub-datasets

From the gene expression studies on the (N)PE, we ranked the genes by expression level and assigned percentile ranks (P) [12]. Next, we formed four sub-datasets: high expression (expression >90th P), moderate expression (50-90th P), low expression (10-50th P) and very low expression (<10th P). This means, for example, that a gene in the high expression sub-dataset has an expression intensity that falls into the highest 10% intensity values of the microarray. We also performed a statistical comparison between NPE and PE and stated gene expression statistically significant different between NPE and PE when p-value < 0.01 (correction for multiple testing; Significant Different sub-dataset). Finally, we determined the signature genes of NPE and PE, based on the selection criteria of fold-change >2.5 and p-value < 0.01 between NPE and PE.

From these different gene expression analyses, we studied the molecular background of AH production. In previous functional and physiological studies of in vitro and in vivo animal models of the CBE, ion channels that are involved in the AH production were widely studied. In a review of Civan and coworkers [2], these findings are put together in a model of AH production. We now looked for the gene expression level (high, moderate, low and very low) of the genes that code for these ion-channels. We also checked whether genes were a signature gene for one epithelial layer. We listed all the gene expression information of these ion-channel coding genes in Additional file 1: Table S1 and in Table 1 we presented genes with moderate or high expression level in at least one epithelial layer.

Table 1 Human CBE gene expression profiles for ion channels implicated in AH production
Full size table

Data analysis II: molecular machinery of AH production, CBE transport and (N)PE interaction

We further studied the molecular machinery of AH production, CBE transport and interaction of NPE and PE. We therefore created new sub-datasets. First, we combined the genes from the Highly Expressed NPE and PE sub-datasets (>90th P) and the genes from the Significantly Different sub-dataset (with mean gene expression value higher than 50th percentile), since we assumed that at least these genes are of biological importance for the (N)PE [12, 13]. We took all Significantly Different expressed genes together, so without subdivision between NPE and PE, since it is difficult to attribute specific biological meaning in a set of up- or down regulated genes, without a base-line, to one of two specific (N)PE layers that are both healthy and from the same person (p-value<0.01 and fold-change >2.5). Still, we believe that these differences, without attributing them to a specific epithelial layer, are of biological meaning. We therefore choose to include the highly and significant different expressed genes together for future analysis. The cell specific genes were determined with much more strict criteria and stated as signature genes. The signature genes were than specifically attributed to PE or NPE. This approach was also previously used [11]. From this pool we selected those genes that coded, according to the Ingenuity knowledge database, for either plasma membrane or extracellular secreted proteins. Obviously, the entries from this selection are most likely to be involved in transmembrane transport mechanisms and interactions between two cell types compared with genes localized in the cytoplasm or nucleus. We named this new dataset, which contained highly expressed genes in the plasma membrane or extracellular space of the CBE, the Transport sub-dataset (flow diagram in Figure 1) (total Transport sub-dataset in Additional file 2: Table S2).

Figure 1
figure 1

Flow diagram. Schematic overview of the subdivisions made in the gene expression data of the NPE and PE into different sub-datasets.

Full size image

Data analysis III: detailed analysis of the transport sub-dataset

We analyzed the Transport sub-dataset from five different points of view:

  1. 1)

    In silico hypothetical molecular model of the AH production: We translated the model of AH production from Civan and coworkers [2] into a hypothetical in silico molecular model of AH production based on our Transport sub-dataset. We specifically indicated the signature genes in the figure.

  2. 2)

    Transport mechanisms in the neuro-epithelia of the CBE: We studied all transporters and ion channels present in the Transport sub-dataset to infer possible specialized transport functions of the CBE. The signature genes of NPE or PE were specifically indicated. For conformation, we performed immunohistochemistry of two proteins of interest, caveolin-1 (CAV1) and clathrin heavy chain (CLTR). The methods of immunohistochemistry were described in detail elsewhere [11]. We used primary antibodies of Abcam (Cambridge, UK), for CAV1 a rabbit polyclonal (ab2910) and for CLTR also a rabbit polyclonal (ab21679).

  3. 3)

    Apical PE-NPE interaction and transport mechanisms: In order to study the possible molecular interaction PE and NPE, we selected from the Transport sub-dataset those genes which encoded proteins that are apically present in the plasma membrane or apically secreted into the extracellular space. We derived the protein localizations from the Ingenuity knowledge database and from Pubmed. We used the search term 'gene name AND (apical OR basolateral OR subcellular location OR immuno location)’. We formed two new sub-datasets, named the Apical Plasma Membrane sub-dataset and Apical Extracellular secreted sub-dataset (flow-diagram in Figure 1). With these two newly formed sub-datasets, we analyzed functional interaction by exploring possible direct connections between the entries using Ingenuity’s tool 'Connect’, resulting in two functional molecular networks. The first one indicating the functional interaction between apical secreted molecules in the extracellular space between NPE and PE and the second one indicating the transport molecules in the apical plasma membrane that are connected to the apical secreted molecules from the former network. Entries connected to themselves have a circle on top; entries that did not have a connection with another gene were removed (detailed description of Ingenuity see http://www.ingenuity.com and [11]). Yellow symbols indicate PE signature genes.

  4. 4)

    CBE transport and possible implications for POAG: Based on the list of (candidate) POAG disease genes that we recently selected (reviewed by [14]), we looked in our Transport sub-dataset for presence of these (candidate) POAG disease genes.

  5. 5)

    Possible pharmacological targets for AH production altering drugs within the CBE. Within the Transport sub-dataset we looked for genes that code for targets of known AH production lowering drugs. We also looked in the hypothetical molecular model for AH production that we built (first aim) for drugs that target on proteins in this model. These drugs might be of interest for future POAG medication strategies to alter/decreased AH production.

Results

1. In silicohypothetical molecular model of the AH production

Twenty four highly expressed genes in our Transport sub-dataset coded for ion-channels involved in AH production. Furthermore, 32 additional ion-channel coding genes were moderately expressed in the NPE and/or PE (Table 1). Of all these genes, four were signature genes of the PE compared the NPE, namely ATP1B1, GJA4, KCNB2, and KCND2. We translated the model of AH production of Civan and coworkers [2] into a hypothetical molecular model of AH production in which we included the genes of the Transport sub-dataset (Figure 2).

Figure 2
figure 2

Hypothetical molecular model of aqueous humor production. On the left, there is the physiological model of aqueous humor (AH) production of Civan and coworkers [2]. This model includes the different ion-channels involved in the AH production. On the right, there is hypothetical in silico molecular model of AH production. This model is based on our gene expression data of the ciliary body epithelia (CBE) and our strict selection criteria of the 10% highest and significant different expressed genes of the non-pigmented (NPE) and pigmented epithelia (PE) (all together in the Transport sub-dataset; see Methods). For the Na+/K+/2Cl- symporter (coded by SLC12A1 and SLC12A2) we did not found a coding gene in the Transport sub-dataset (indicated with * in the model). Genes coding for ion-channels with moderate expression levels (for example SLC12A2) can be found in Table 1. The underlined genes (ATP1B1, GJA4, KCNB2, KCND2) are signature genes of the PE compared to the NPE (fold change >2.5 and p-value < 0.01; see Methods section).

Full size image

2. Transport mechanisms in the neuro-epithelia of the CBE

To determine which CBE highly expressed genes could ultimately be involved in the composition of the AH, we assessed the expression of all transporters and (cat-)ion channels and their specific substrate in the Transport sub-dataset. Total list is outlined in Table 2. Different groups of transporters can be distinguished.

Table 2 All transporters in the Transport sub-dataset of the ciliary body epithelia
Full size table

Vesicle mediated transporter gene expression in CBE: caveolin and clathrin

We identified many different entries in our Transport sub-dataset involved in endo- and exocytosis as well as vesicle mediated transport (ANXA2, -4, -5, -7, CAV1, CLTA, CLTC, EXOC7, SORL1, STXBP3, SV2B, SYPL1, SYT13, VAT1, and VTI1B).

Using immunohistochemistry, we confirmed the presence of caveolin-1 (CAV1) and clathrin heavy chain proteins (CLTR) in the CBE (Figure 3). More specifically, we found strong staining of CAV1 in the NPE and possible light staining in the PE. Furthermore, we observed strong staining of CLTR in both the NPE and PE of the CBE, especially in the apical plasma membrane.

Figure 3
figure 3

Immunohistochemistry of caveolin 1 and clathrin in the ciliary body. Immunofluorescence staining for caveolin 1 (CAV1) and clathrin (CLTR) in human ciliary body sections. Our gene expression data show that CAV1 is significantly different expressed between NPE and PE (p-value < 0.01), with higher expression in NPE and all CLTR coding genes (CLTA, CLTB and CLTC) are highly (>90th P) expressed in both NPE and PE. In this immunofluorescence picture, CAV1 is clearly present in the non-pigmented epithelium (NPE) and possibly also in the pigmented epithelium (PE). We found CLTR in both NPE and PE, especially in the apical plasma membrane. Both CAV1 and CLTR were found in the vascular endothelium of the stroma. Negative controls were for both proteins negative (not shown). Legend: blue = dapi = cell nucleus; red = cy3 = protein of interest.

Full size image

Specific transporter activity in CBE: hormones, vitamins, metabolism, metal ions

We found a number of genes coding for neuro-endocrine transporters, like glutamate transporters (SLC1A3, SLC1A4, SLC1A7), glutamate-gated ion channels (GRIN2C, GRIN3A), GABA transporter (SLC6A13) and transporters of thyroid hormone (SLC16A2, TTR), thiamine (SLC19A2), vitamin C (SLC23A2; signature gene of PE), lipids (SLC27A1, APOE, LDLR, SORL1, and SCARB1) and glucose (SLC2A1, SLC2A4). We also identified genes coding for transporters of zinc (SLC39A14, SLC39A8), copper (SLC31A2) and iron (SLC40A1 (signature gene of PE), TF, HEPH) in the Transport sub-dataset.

Transporter activity in CBE: organic anions and ions (other than in the hypothetical molecular model of AH production)

We identified high expression of several genes coding for ABC transporters (ABCA1, ABCA4, ABCA5, ABCA7, ABCB1, ABCC1, ABCC5, ABCC9 and ABCG2), responsible for organic anion transport. ABCA1 and ABCC1 are signature genes of the PE when compared to the NPE. Moreover, we found genes coding for phosphate transporters (SLC20A1, SLC20A2), carnitine transporters (SLC22A4, SLC22A5), and amino acid transporters (SLC38A1, SLC6A15 (signature gene of PE), SCL7A5, SLC7A7, SLC7A8). Next, we found genes coding for other ATP-transporters than that were already attributed for the AH production, namely ATP2B3, ATP2B4, ATP13A1, -3, and ATP9A, and genes coding for K+/Cl- transporters (SLC12A4, SLC12A7) and the Na+/K+/Ca2+ exchanger (SLC26A7, signature gene of PE).

3. Apical PE-NPE interaction and transport mechanisms

Molecular interactions in the extracellular space between the (N)PE apical sides

We explored the potential functional relationships between the predicted secreted entries of the Apical Extracellular Secreted sub-dataset (Figure 4). This network represents potential molecular interactions taking place in the extracellular space between the NPE and PE. The top functions assigned to this network by Ingenuity were “Neurological disease, psychological disorders and cellular movement”. The network contained 34 genes of a total of 67 (51%) secreted genes in the (N)PE apical extracellular space. Interestingly, this network contained many entries involved in matrix assembly and cytoskeleton organization (COL18A1, CTGF, DCN, FBN1, FBLN1, LTBP1, SPARC, and SPP1). Indeed, these proteins may partly fill the extracellular space between the apical membranes of the NPE and PE, as was previously shown for FBN1 by Gabriel and coworkers [15]. We also found genes coding for components of the classical complement cascade (C1R, C1S, C3, and CFB), the coagulation pathway (F10, PLAT and PROS1) and immune and/or inflammatory responses (CCL2 and IL6) and we identified many genes coding for a range of growth factors (CTGF, EGF, IGF1, IGF2, IGFBP2-6, LTBP1, PTN, TGFB2 and VEGFA). Finally, we identified several entries apparently involved in amyloid beta metabolism and plaque formation (A2M, APOE, CLU, SERPINA3, and TTR).

Figure 4
figure 4

Molecular interactions in the apical extracellular space of the pigmented (PE) and non-pigmented epithelium (NPE). Molecular network build in Ingenuity (http://www.Ingenuity.com), representing entries that are highly expressed and secreted in the apical extracellular space between NPE and PE. The top functions assigned to this network by Ingenuity were “Neurological disease, psychological disorders and cellular movement”. This network is contains many entries involved in matrix assembly and cytoskeleton organization (COL18A1, CTGF, DCN, FBN1, FBLN1, LTBP1, SPARC, and SPP1) and molecules involved in growth and proliferation (CTGF, EGF, IGF1, IGF2, IGFBP2-6, LTBP1, PTN, TGFB2 and VEGFA). IGF1 is a signature gene of the PE (yellow symbol). Together, these molecules are likely to be involved in the local turnover of the ECM between the PE-NPE layers of the CBE. We also identified several entries within the PE-NPE apical extracellular space of the CBE that are involved in amyloid-beta metabolism and plaque formation (A2M, APOE, CLU, SERPINA3, and TTR). Legend of different lines: Circles on top indicate direct relationships of entries with itself. Solid lines indicate that binding of the two proteins occur, arrows indicate that the first protein interfere with the expression or activity of the second protein.

Full size image

Molecular interactions in apical transport of the CB

We were interested in the transport mechanisms of the entries in the apical extracellular space of the NPE and PE. Therefore, we explored in silico potential functional relationships between the genes coding for transporters in the Apical Plasma Membrane sub-dataset and the network described above. The resulting network is displayed in Figure 5. The top functions assigned by Ingenuity to this network were “Neurological disease, psychological disorders and cell-to-cell signaling and interaction”. This network showed direct interactions of 15 transporter genes (from a total of 50 genes; 30%) in the apical plasma membrane with the apical extracellular secreted molecules. Clearly, several entries involved in vesicle-mediated transport were connected to the extracellular secreted entries, namely ANXA7, CAV1, CLTA, CLTC, and SORL1. Furthermore, we identified several genes in the apical plasma membrane that codes for proteins involved in lipid transport (ABCA1 (signature gene of PE), LDLR and SCARB1), iron transport (SLC40A1; signature gene of PE) and xenobiotic/drug transport (ABCC5 and ABCG2) that are connected to molecules in the apical extracellular space of the CBE.

Figure 5
figure 5

Molecular interactions in apical transport of the ciliary body epithelia. Molecular network build in Ingenuity (http://www.Ingenuity.com), representing possible transporters in the apical plasma membranes of the pigmented and non-pigmented epithelia that transport the apical extracellular excreted entries outlined in Figure 4. The top functions assigned by Ingenuity to this network were “Neurological disease, psychological disorders and cell-to-cell signaling and interaction”. The major transport mechanisms to which these genes were connected were vesicle-mediated transporters (ANXA7, CAV1, CLTA, CLTC, and SORL1) and lipid (ABCA1, LDLR and SCARB1) and xenobiotic (ABCC5 and ABCG2) transporters. Two PE signature genes, ABCA1 and SLC40A1 (yellow symbols) are involved in this network. Legends of different lines: see Figure 4.

Full size image

4. CBE transport and possible implications for POAG

In our Transport sub-dataset we found seven genes that were previously associated with POAG, namely APOE, CAV1, COL8A2, EDNRA, FBN1, RFTN1 and TLR4. None of these genes were signature genes of the NPE or PE.

5. Possible pharmacological targets for AH production altering drugs within the CBE

The three classes of drugs that can lower AH production are adrenergic agonists, carbonic anhydrase (CA) inhibitors and beta blockers. Within our Transport sub-dataset we found high expression of CA4 and CA14 in the plasma membrane of the CBE.

We also looked closer to the ion-channels in the hypothetical model of AH production (Figure 2) for potentially new drugs that might target these channels. For the Na+/K + ATPase (coded by ATP1A1 and ATP1A2) there are several inhibitors, namely digoxin, omeprazole, ethacrynic acid and perphenazine. Two sodium channels coded by SCN4B and SCN8A genes are inhibited by riluzole, a drug used in ALS. Finally, the potassium channels coded by KCNA5 and KCNQ1 can be blocked by quinidine, amiodarone, tedisamil and flecainide (antianginal and antiarrhythmic drugs) or opened by nicorandil (frequently used against hypertension).

Discussion

In order to shed light on the molecular mechanisms in the CBE related to AH production and composition, we analyzed our expression data of the NPE and PE from five different angles. First, we studied the expression profiles of the genes coding for all ion-channels involved in AH production and we built in silico a hypothetical molecular model for AH production. Next, we assessed putative transport mechanisms in the NPE and PE, including vesicle mediated and specialized transport. Third, we analyzed, in silico, several apical molecular networks present between the NPE and PE, amongst others with respect to transport and extracellular matrix molecules. Fourth, we identified a number of clues which suggest that the CBE plays an important (transport) role in POAG. And finally, we looked for existing and possible new pharmacological targets in the CBE which could potentially modify the AH production.

Our microarray approach and corresponding functional annotations have several limitations, which are extensively discussed elsewhere [12, 13]. The gene expression datasets used in the current study were first roughly analyzed in a previous paper, which describes the general properties of the CBE [11]. The value of the current study is that we now analyzed these datasets in detail, primarily focused on the molecular machinery underlying transport mechanisms. Fortunately, the physiology of the transport mechanisms in the CBE has previously already been studied extensively in in vitro and in vivo models, and was leading for our molecular analysis and interpretation.

Obviously, the amount of (combination of old and) new data presented in this manuscript prohibits verification in vitro and in vivo models in one single study, let alone in this manuscript.

Finally, it must be taken into account that we used very strict selection criteria for our Transport sub-dataset. Therefore, the genes that we selected for further future analysis in this manuscript may not cover all genes that are involved in the transport and AH production functions of the CBE. It is not unlikely that there is functional redundancy, and that several similar genes may code for more or less effective additional transport proteins, or that more – yet to be identified – proteins may be involved in AH production. On the other hand, because of our strict selection criteria, the genes that are present in this model, are likely to be truly involved in the molecular machinery of AH production and CBE transport mechanisms.

1. In silicohypothetical molecular model of the AH production

We build a hypothetical molecular model of AH production based on the gene expression data of (the ion-channels in) the CBE (Figure 2). We selected genes under strict selection criteria in the Transport sub-dataset.

Due to these strict selection criteria, a number of additional genes or proteins, previously implicated in AH production, did not pop up in our present hypothetical molecular model (Figure 2). For example, for the Na+/K+/2Cl- symporter (coded by SLC12A1 or SLC12A2) we did not found a coding gene in our Transport sub-dataset, but these genes were moderately (SLC12A2) and lowly (SLC12A1) expressed in the (N)PE (Table 1).

2. Transport mechanisms in the neuro-epithelia of the CBE

The CBE expresses many genes coding for ion-channels and transporters at high levels. Due to their known substrate specificity of their corresponding proteins, we can also predict the biomolecules transported by these CBE entries. Consequently, our inventory, in combination with data from the literature, showed that the CBE most likely transports a wide variety of molecules, including cations, organic cations and anions, amino acids, glutamate, GABA, iron, copper, zinc, thyroid hormone, vitamin A, B and C, glucose, and fatty acids. Indeed, the majority of these ions and biomolecules were also found in previous proteomic studies of the AH, which suggests that the CBE is an important ocular entry point of these molecules [16–20].

Interestingly, we also observed high expression of many entries involved in vesicle mediated trans- or endocytosis. For example, we found high expression of CAV1 and CLTC in the CBE and we confirmed their presence by immunohistochemistry (Figure 3). Caveolin and clathrin are involved in vesicle mediated endo-, exo- and transcytosis of large molecules through the cell [21–25]. Caveolin-mediated vesicles transport for example iron, transferrin, insulin, lipids, albumin, chemokines and pathogens [21]. Clathrin-mediated vesicles can transport nutrients, viruses, toxins, plasma membrane proteins and signaling receptors, including ß2-adrenergic receptors, CD4, insulin receptors, T-cell and B-cell receptors [24, 25]. While caveolin transport is more or less a-specific, clathrin-mediated vesicle transport is only possible with adaptor proteins (receptors). We also found expression of several genes coding for other proteins involved in less specific vesicle transport, like ANXA2, EXOC7, SORL1, STXBP3, SYT13, VAT1, and VTI1B

In summary, our data, in combination with the literature, suggest that both general and specific vesicle mediated transport play a role in the transport of a range of biomolecules over the CBE neuroepithelia. Indeed, this was already suggested by a histochemical electron microscope study [26]. We hypothesize that these transport mechanisms may (also) be involved in AH composition, in particular with regard to the presence of large molecules.

3. Apical PE-NPE interaction and transport mechanisms

In silico, we assigned several interacting proteins to the extracellular space between the PE and NPE, where transported molecules from the blood to the AH are likely to pass. We found several proteins involved matrix assembly and cytoskeleton organization, a network of protease (inhibitors) and a group of proteins involved in cell growth and proliferation pathways. These molecules are likely to be involved in the local turnover of the ECM between the PE-NPE layers of the CBE. We also identified several entries within the PE-NPE apical extracellular space of the CB that are involved in amyloid-beta metabolism and plaque formation. The major transport mechanisms to which these entries were connected were vesicle-mediated transporters and lipid and xenobiotic transporters. The molecules within the PE-NPE extracellular space and their transporters may influence intrinsic filtration or secretion capabilities of the CBE. If disrupted by a genetic mutation, or altered by age-related oxidative stress, these molecules change and may affect the AH secretion and composition.

4. CBE transport and possible implications for POAG

We identified seven POAG candidate disease genes in our Transport sub-dataset. These genes were APOE, CAV1, COL8A2, EDNRA, FBN1, RFTN1 and TLR4. We hypothesized that mutations in these genes might modify the AH production and composition via the CBE. Although the production rate of AH is hardly changed in POAG [27], the composition of AH does, with overall an increase of (large and harmful) proteins in POAG AH compared to controls [5–7]. The POAG disease genes that we found in our Transport sub-dataset code for structural components (FBN1, COL8A2), transporters (CAV1, APOE, RFTN1) and signaling molecules (CAV1, EDNRA), all potentially involved in AH production and composition. Further research is warranted to prove or reject their involvement in disturbed AH dynamics of the CBE during POAG. It would be of interest to study mouse models with mutations in these genes, measure IOP, optic nerve head characteristics and AH composition. It would also be of interest to study the effect of IOP-lowering drugs in animals with mutations in these target genes in comparison with controls.

5. Possible pharmacological targets for AH production altering drugs within the CBE

Last but not least, we annotated our molecular model (Figure 2, discussed above) with possible new therapeutic targets that might alter AH production. First, we found several drugs that block Na+/K + ATPases: digoxin, ethacrynic acid, omeprazol, and perphenazine. Digoxin has antiarrhythmic effects via myocytes and some older literature sources describe the lowering effect of digoxin on the IOP in human [28–30] and cats [31], probably due to blocking of Na+/K+-ATPase in the CBE. Ethacrynic is a diuretic drug and multiple animal studies already showed the IOP lowering effect of intracameral or topical administered ethacrynic acid in the eye, probably through increased AH outflow, but possibly also by decreased AH production [32–39]. The major disadvantage of athacrynic acid is the risk for (severe) corneal edema. Omeprazole is a proton pump inhibitor that inhibits gastric acid secretion in the stomach. In the literature, we found one clinical study in which treatment with omeprazole indeed lowered the IOP [40]. The drug perphenazine is an antagonist of dopamine receptors and has antiemetic and antipsychotic effects. We found no studies of perhanzine in relation to IOP altering effects.

Secondly, we found one drug that block sodium channels, namely riluzole. This drug is used in ALS. In the literature, we couldn’t find studies describing an effect on IOP with riluzole use.

Finally, we found several drugs that influence the activity of potassium channels. Quinidine, amiodarone, tedisamil and flecainide block potassium channels, while nicorandil opens them. Only for the latter drug, nicorandil, two studies exist on the IOP effects. Chiang and coworkers [41] found increased IOPs after nicorandil use in rabbits, whereas Chowdhury and coworkers [42] found an IOP lowering effect of nicorandil in rats.

In conclusion, we found possible new therapeutic targets for AH production modification in the CBE. These targets are of interest, since IOP lowering is currently the only proven therapy for POAG and current medical treatment modalities often fail to reduce the IOP sufficiently. More research into these new targets might result in new treatment options for POAG, which may make surgery, with its side effects and risks, less often needed. Extensive studies in animal models are warranted to study potential beneficial effects of these drugs on the IOP and their potential harmful side effects.

Conclusions

The transport mechanisms in the CBE are numerous, ranging from extensively studied specific ion channels and pumps that form the basis of AH production, to more generalized transport functions like vesicle mediated transport. The CBE transports (neuro-) endocrine and metabolic particles and is involved in specific apical transport of ECM proteins and growth factors. Taken together, these transport mechanisms are important in determining the AH production rate and composition. In POAG, the composition of the AH is changed and several POAG disease gene are highly expressed, and their corresponding proteins localized in the CBE plasma membrane or extracellular space. We propose that mutations in these POAG genes cause disturbances in the AH composition by the CBE, resulting in a changed AH content and disturbed outflow facility by the TM. Finally, we found possible new therapeutic targets for AH production modification in the CBE.

Abbreviations

AH:

Aqueous humor

CBE:

Ciliary body epithelia

IOP:

Intraocular pressure

NPE:

Non-pigmented epithelium

PE:

Pigmented epithelium

POAG:

Primary open angle glaucoma

TM:

Trabecular meshwork.

References

  1. Civan MM: Transport components and net secretion of the aqueous humor and their integrated regulation. Current Topics in Membranes: The eye's aqueous humor. From secretion to glaucoma, Volume 45. Edited by: Civan MM. 1998, San Diego, London, Boston, New York, Sydney, Tokyo, Toronto: Academic Press, 1-24.

    Google Scholar 

  2. Civan MM, Macknight AD: The ins and outs of aqueous humour secretion. Exp Eye Res. 2004, 78: 625-631. 10.1016/j.exer.2003.09.021.

    Article  CAS  PubMed  Google Scholar 

  3. Kraft ME, Glaeser H, Mandery K, Konig J, Auge D, Fromm MF, et al: The prostaglandin transporter OATP2A1 is expressed in human ocular tissues and transports the antiglaucoma prostanoid latanoprost. Invest Ophthalmol Vis Sci. 2010, 51: 2504-2511. 10.1167/iovs.09-4290.

    Article  PubMed  Google Scholar 

  4. Shin BC, Suzuki T, Tanaka S, Kuraoka A, Shibata Y, Takata K: Connexin 43 and the glucose transporter, GLUT1, in the ciliary body of the rat. Histochem Cell Biol. 1996, 106: 209-214. 10.1007/BF02484402.

    Article  CAS  PubMed  Google Scholar 

  5. Duan XM, Xue P, Wang NL, Dong Z, Lu QJ, Yang FQ: Proteomic analysis of aqueous humor from patients with primary open angle glaucoma. Mol Vis. 2010, 16: 2839-2846.

    PubMed Central  CAS  PubMed  Google Scholar 

  6. Grus FH, Joachim SC, Sandmann S, Thiel U, Bruns K, Lackner KJ, et al: Transthyretin and complex protein pattern in aqueous humor of patients with primary open-angle glaucoma. Mol Vis. 2008, 14: 1437-1445.

    PubMed Central  CAS  PubMed  Google Scholar 

  7. Izzotti A, Longobardi M, Cartiglia C, Sacca SC: Proteome alterations in primary open angle glaucoma aqueous humor. J Proteome Res. 2010, 9: 4831-4838. 10.1021/pr1005372.

    Article  CAS  PubMed  Google Scholar 

  8. Casson RJ, Chidlow G, Wood JP, Crowston JG, Goldberg I: Definition of glaucoma: clinical and experimental concepts. Clin Experiment Ophthalmol. 2012, 40: 341-349. 10.1111/j.1442-9071.2012.02773.x.

    Article  PubMed  Google Scholar 

  9. Leske MC, Connell AM, Wu SY, Hyman LG, Schachat AP: Risk factors for open-angle glaucoma. The Barbados Eye Study. Arch Ophthalmol. 1995, 113: 918-924. 10.1001/archopht.1995.01100070092031.

    Article  CAS  PubMed  Google Scholar 

  10. Weinreb RN, Khaw PT: Primary open-angle glaucoma. Lancet. 2004, 363: 1711-1720. 10.1016/S0140-6736(04)16257-0.

    Article  PubMed  Google Scholar 

  11. Janssen SF, Gorgels TG, Bossers K, ten Brink JB, Essing AH, Nagtegaal M, et al: Gene expression and functional annotation of the human ciliary body epithelia. Plos One. 2012, 7: e44973-10.1371/journal.pone.0044973.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  12. Booij JC, van Soest S, Swagemakers SMA, Essing AHW, Verkerk AJMH, van der Spek PJ, et al: Functional annotation of the human retinal pigment epithelium transcriptome. Bmc Genomics. 2009, 10: 164-10.1186/1471-2164-10-164.

    Article  PubMed Central  PubMed  Google Scholar 

  13. Booij JC, ten Brink JB, Swagemakers SM, Verkerk AJ, Essing AH, van der Spek PJ, et al: A new strategy to identify and annotate human RPE-specific gene expression. Plos One. 2010, 5: e9341-10.1371/journal.pone.0009341.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  14. Janssen SF, Gorgels TGNF, Ramdas WD, Klaver CCW, Van Duijn CM, Jansonius NM, et al: The Vast Complexity of Primary Open Angle Glaucoma: Disease Genes, Risks, Molecular Mechanisms and Pathobiology. Prog Retin Eye Res. 2013, in press

    Google Scholar 

  15. Gabriel LA, Wang LW, Bader H, Ho JC, Majors AK, Hollyfield JG, et al: ADAMTSL4, a secreted glycoprotein widely distributed in the eye, binds fibrillin-1 microfibrils and accelerates microfibril biogenesis. Invest Ophthalmol Vis Sci. 2011, 1: 461-469.

    Google Scholar 

  16. Chowdhury UR, Madden BJ, Charlesworth MC, Fautsch MP: Proteome Analysis of Human Aqueous Humor. Invest Ophthalmol Vis Sci. 2010, 51: 4921-4931. 10.1167/iovs.10-5531.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Cumurcu T, Mendil D, Etikan I: Levels of zinc, iron, and copper in patients with pseudoexfoliative cataract. Eur J Ophthalmol. 2006, 16: 548-553.

    CAS  PubMed  Google Scholar 

  18. Kallberg ME, Brooks DE, Gelatt KN, Garcia-Sanchez GA, Szabo NJ, Lambrou GN: Endothelin-1, nitric oxide, and glutamate in the normal and glaucomatous dog eye. Vet Ophthalmol. 2007, 10: 46-52. 10.1111/j.1463-5224.2007.00529.x.

    Article  CAS  PubMed  Google Scholar 

  19. Leite MT, Prata TS, Kera CZ, Miranda DV, Barros SBD, Melo LAS: Ascorbic acid concentration is reduced in the secondary aqueous humour of glaucomatous patients. Clin Exp Ophthalmol. 2009, 37: 402-406. 10.1111/j.1442-9071.2009.02046.x.

    Article  Google Scholar 

  20. Pelletier CC, Lambert JL, Borchert M: Determination of glucose in human aqueous humor using Raman spectroscopy and designed-solution calibration. Appl Spectrosc. 2005, 59: 1024-1031. 10.1366/0003702054615133.

    Article  CAS  PubMed  Google Scholar 

  21. Frank PG, Pavlides S, Lisanti MP: Caveolae and transcytosis in endothelial cells: role in atherosclerosis. Cell Tissue Res. 2009, 335: 41-47. 10.1007/s00441-008-0659-8.

    Article  CAS  PubMed  Google Scholar 

  22. Jung N, Haucke V: Clathrin-mediated endocytosis at synapses. Traffic. 2007, 8: 1129-1136. 10.1111/j.1600-0854.2007.00595.x.

    Article  CAS  PubMed  Google Scholar 

  23. Parkar NS, Akpa BS, Nitsche LC, Wedgewood LE, Place AT, Sverdlov MS, et al: Vesicle Formation and Endocytosis: Function, Machinery, Mechanisms, and Modeling. Antioxid Redox Signal. 2009, 11: 1301-1312. 10.1089/ars.2008.2397.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Puertollano R: Clathrin-mediated transport: assembly required - Workshop on Molecular Mechanisms of Vesicle Selectivity. Embo Reports. 2004, 5: 942-946. 10.1038/sj.embor.7400249.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Tuma PL, Hubbard AL: Transcytosis: Crossing cellular barriers. Physiol Rev. 2003, 83: 871-932.

    Article  CAS  PubMed  Google Scholar 

  26. Eggli PS, Graber W, Vanderzypen E: Formation and Release of Vesicles from the Basal Surfaces of Rat Eye Nonpigmented Ciliary Epithelial-Cells - A Novel Secretory Mechanism. Anat Rec. 1991, 231: 156-166. 10.1002/ar.1092310203.

    Article  CAS  PubMed  Google Scholar 

  27. Larsson LI, Rettig ES, Brubaker RF: Aqueous flow in open-angle glaucoma. Arch Ophthalmol. 1995, 113: 283-286. 10.1001/archopht.1995.01100030037018.

    Article  CAS  PubMed  Google Scholar 

  28. Hardt BW, Johnen R, Fahle M: The influence of systemic digitalis application on intraocular pressure. Graefes Arch Clin Exp Ophthalmol. 1982, 219: 76-79. 10.1007/BF02173445.

    Article  CAS  PubMed  Google Scholar 

  29. Hegazy MA, Labib MA, Ghaleb H, Barradah A: Effect of digoxin on intraocular pressure. Bull Ophthalmol Soc Egypt. 1967, 60: 83-91.

    CAS  PubMed  Google Scholar 

  30. Peczon JD: Clinical evaluation of digitalization in glaucoma. Arch Ophthalmol. 1964, 71: 500-504. 10.1001/archopht.1964.00970010516012.

    Article  CAS  PubMed  Google Scholar 

  31. Ferraiolo BL, Pace DG: Digoxin-induced decrease in intraocular pressure in the cat. Eur J Pharmacol. 1979, 55: 19-22. 10.1016/0014-2999(79)90143-2.

    Article  CAS  PubMed  Google Scholar 

  32. Arnold JJ, Choksi Y, Chen X, Shimazaki A, Hatten J, Toone EJ, et al: Eyedrops containing SA9000 prodrugs result in sustained reductions in intraocular pressure in rabbits. J Ocul Pharmacol Ther. 2009, 25: 179-186. 10.1089/jop.2008.0107.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Croft MA, Hubbard WC, Kaufman PL: Effect of ethacrynic acid on aqueous outflow dynamics in monkeys. Invest Ophthalmol Vis Sci. 1994, 35: 1167-1175.

    CAS  PubMed  Google Scholar 

  34. Green K, Mayberry L: Drug effects on the hydraulic conductivity of the isolated rabbit ciliary epithelium. Q J Exp Physiol. 1985, 70: 271-281.

    Article  CAS  PubMed  Google Scholar 

  35. Searles RV, Johnson M, Shikher V, Balaban CD, Severs WB: Effects of ethacrynic acid on intraocular pressure of anesthetized rats. Proc Soc Exp Biol Med. 1999, 220: 184-188. 10.3181/00379727-220-44362.

    Article  CAS  PubMed  Google Scholar 

  36. Shimazaki A, Ichikawa M, Rao PV, Kirihara T, Konomi K, Epstein DL, et al: Effects of the new ethacrynic acid derivative SA9000 on intraocular pressure in cats and monkeys. Biol Pharm Bull. 2004, 27: 1019-1024. 10.1248/bpb.27.1019.

    Article  CAS  PubMed  Google Scholar 

  37. Tingey DP, Schroeder A, Epstein MP, Epstein DL: Effects of topical ethacrynic acid adducts on intraocular pressure in rabbits and monkeys. Arch Ophthalmol. 1992, 110: 699-702. 10.1001/archopht.1992.01080170121038.

    Article  CAS  PubMed  Google Scholar 

  38. Tingey DP, Ozment RR, Schroeder A, Epstein DL: The effect of intracameral ethacrynic acid on the intraocular pressure of living monkeys. Am J Ophthalmol. 1992, 113: 706-711.

    Article  CAS  PubMed  Google Scholar 

  39. Wang RF, Podos SM, Serle JB, Lee PY, Neufeld AH, Deschenes R: Effects of topical ethacrynic acid ointment vs timolol on intraocular pressure in glaucomatous monkey eyes. Arch Ophthalmol. 1994, 112: 390-394. 10.1001/archopht.1994.01090150120033.

    Article  CAS  PubMed  Google Scholar 

  40. Kountouras J, Mylopoulos N, Chatzopoulos D, Zavos C, Boura P, Konstas AG, et al: Eradication of Helicobacter pylori may be beneficial in the management of chronic open-angle glaucoma. Arch Intern Med. 2002, 162: 1237-1244. 10.1001/archinte.162.11.1237.

    Article  PubMed  Google Scholar 

  41. Chiang CH, Lin CH: Effects of cromakalim and nicorandil on intraocular pressure after topical administration in rabbit eyes. J Ocul Pharmacol Ther. 1995, 11: 195-201. 10.1089/jop.1995.11.195.

    Article  CAS  PubMed  Google Scholar 

  42. Chowdhury UR, Bahler CK, Hann CR, Chang M, Resch ZT, Romero MF, et al: ATP-sensitive potassium (KATP) channel activation decreases intraocular pressure in the anterior chamber of the eye. Invest Ophthalmol Vis Sci. 2011, 52: 6435-6442. 10.1167/iovs.11-7523.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgement

We thank the Corneabank, Beverwijk, the Netherlands, for providing donor eyes.

This study was supported by grants from the General Dutch Foundation Preventing Blindness, the National Foundation for Blinds and Low Vision, the Foundation Blinden-Penning, the Foundation Glaucoomfonds, the Rotterdam Foundation of Blinds and the Professor Mulder Foundation (all coordinated by UitZicht, project # UitZicht2008-7). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

  1. Department of Clinical and Molecular Ophthalmogenetics the Netherlands Institute for Neuroscience (NIN), Royal Netherlands Academy of Arts and Sciences (KNAW), Amsterdam, the Netherlands

    Sarah F Janssen, Theo GMF Gorgels & Arthur AB Bergen

  2. Department of Bioinformatics, Erasmus University Medical Center, Rotterdam, the Netherlands

    Peter J van der Spek

  3. Department of Ophthalmology, University Medical Center Groningen, University of Groningen, Groningen, the Netherlands

    Nomdo M Jansonius

  4. Department of Ophthalmology Academic Medical Center, Amsterdam, the Netherlands

    Arthur AB Bergen

  5. Department of Clinical Genetics Academic Medical Center, Amsterdam, the Netherlands

    Arthur AB Bergen

Authors
  1. Sarah F Janssen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Theo GMF Gorgels
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Peter J van der Spek
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Nomdo M Jansonius
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Arthur AB Bergen
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Arthur AB Bergen.

Additional information

Competing interests

The authors declare that they have no competing interests

Authors' contributions

SFJ carried out the laboratory work, performed the bioinformatic analysis and drafted the manuscript. TGG participated in the design of the study and the manuscript. PJS provided bioinformatics software and background knowledge. NMJ participated in the design of the study and the manuscript. AAB performed bioinformatic analysis and drafted the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Janssen, S.F., Gorgels, T.G., van der Spek, P.J. et al. In silicoanalysis of the molecular machinery underlying aqueous humor production: potential implications for glaucoma. J Clin Bioinform 3, 21 (2013). https://doi.org/10.1186/2043-9113-3-21

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2043-9113-3-21

Keywords

Journal of Clinical Bioinformatics

ISSN: 2043-9113