|
|
||||||||
From the Department of Ophthalmology, Duke University Medical Center, Durham, North Carolina.
| Abstract |
|---|
|
|
|---|
METHODS. A polymerase chain reaction (PCR)amplified cDNA library was constructed using RNA from the TM of a 67-year-old normal, perfused human eye. A total of 1060 clones were randomly selected for sequencing of one end. These sequences were searched against nonredundant GenBank and dbEST databases for similarity comparison by using a FASTA file and the BLASTcl3 program. Relative expression patterns of those clones that matched other expressed sequence tags (ESTs) were determined using the National Center for Biotechnology Information (NCBI) Unique Human Gene Sequence Collection (UniGene) database.
RESULTS. Of the 1060 clones analyzed, 519 (48.9%) had sequences identical with known genes, 125 (11.8%) matched ESTs, and 189 (17.8%) did not match any database sequences. Of the remaining clones, 31 (3%) corresponded to mitochondrial transcripts and 196 (18.5%) to repetitive and noninformative sequences. It is notable that some of the genes highly represented in this library are not ubiquitously expressed in other tissues, which suggests a potentially important role in the HTM. As evidence for the presence of true novel genes in the library, one of the clones was fully sequenced. This clone comprised a complete open reading frame of 966 nucleotides, and its deduced amino acid sequence corresponded to a protein 33% similar to the MAS-related G-proteincoupled receptor.
CONCLUSIONS. The identification of the more highly expressed genes in HTM and the discovery of novel genes expressed in this tissue provides basic information for further research on the physiology of the TM and for the identification of glaucoma candidate genes.
| Introduction |
|---|
|
|
|---|
Despite its small size, the TM has a complex architecture, including several layers and cell types that are morphologically and functionally different. The more anterior portion, the nonfiltering meshwork, consists of three to five elongated beams covered by cells. The posterior portion, or filtering meshwork, consists of three layers: the uveal meshwork, formed by irregular strands; the corneoscleral meshwork, formed by 8 to 15 layers of interconnected flat beams; and the outermost layer, the juxtacanicular region, formed by two to five layers of loosely arranged cells just underneath SC.
Human TM (HTM) samples include three main cell types: (1) trabecular cells, which cover the beams and strands of the nonfiltering, uveal, and corneoscleral meshwork regions, originate from the neural crest, and may be predominately involved in phagocytosis and tissue repair; (2) juxtacanalicular cells, which are randomly distributed within the extracellular matrix of the juxtacanalicular meshwork, are not attached to basal membranes, and exhibit long, slender processes that can reach the SC cells; they are believed to originate from the perivascular cells of Rouget; (3) and SC inner wall endothelial cells, which constitute the last barrier for the aqueous humor before it reaches the lumen of SC. These cells appear to be of vascular origin and have been implicated in modulating the resistance to aqueous humor outflow.1 4 5 6 7 In addition to these three main cell types, electron microscopy of HTM sections often reveals the presence of other cells, such as macrophages.8 These cells are found within the intertrabecular spaces, are usually loaded with pigment granules, and probably leave the meshwork through paracellular routes of the inner wall.9
The TM is implicated in functions such as phagocytosis, extracellular matrix dynamics, secretion, cytoskeletal reorganization, and detoxification of aqueous humor.10 11 It is also believed to be involved in the maintenance of the immune privilege in the eye.12 13 14 However, most of the specific mechanisms by which the TM accomplishes these functions are not well understood.
The biologic features of TM, like those of any tissue, are determined largely at the level of gene expression. The identification of the genes expressed in TM constitutes a necessary step toward the understanding of its physiology. A powerful approach to the analysis of expression in a particular tissue or cell type involves single-pass partial sequencing of clones from a cDNA library together with the comparative analysis of the resultant expressed sequence tags (ESTs) with entries in public databases. A number of sequencing projects have provided a large collection of human ESTs from a variety of tissues.15 16 17 18 However, the gene expression profile in the HTM has not been studiedin part, because of difficulties in obtaining enough RNA from such a small tissue that in turn preclude the construction of cDNA libraries by conventional methods. Polymerase chain reaction (PCR)based methods for construction of cDNA libraries have been developed to circumvent these problems.19 20 21
However, libraries generated using standard PCR methods tend to include a number of artifacts because of the overcycling of the PCR reaction. In addition, these libraries result in an enrichment of small cDNA clones that do not possess the entire open reading frame (ORF) of the genes. Instead, a recently developed, exponential PCR-based method generates high-quality cDNA libraries, enriched for full-length clones. The new method has been used to make good cDNA libraries from mouse blastocyst22 and human CD34+ hematopoietic stemprogenitor cells.23 It also has been successfully used in our laboratory to generate HTM cDNA probes for the hybridization of gene arrays.24 25 Using this new approach, we constructed a cDNA library from the TM of a 67-year-old human donor with no history of glaucoma, and analyzed 1060 clones by single-pass sequencing.
For this first analysis of the HTM gene expression profile, we chose a TM sample that had been perfused under organ culture conditions. Perfusion of the anterior segment of the eye is one of the best and most commonly used models for studying outflow physiology.26 27 28 Although perfused HTMs may not exactly reproduce the pattern of gene expression of the in vivo tissue, they are one of the closest. The model helps recover TM cell activity and maintains the expression of genes responsive to flow and IOP. In comparison, gene expression of a nonperfused HTM sample from a postmortem donor diverges from the in vivo expression pattern because of ongoing cell death and absence of IOP and aqueous humor flow.
Our results provide the first profile of gene expression in a human TM. This information will undoubtedly help in understanding the normal physiology of the TM and will provide new candidate genes potentially involved in the pathophysiology and the susceptibility to glaucoma.
| Materials and Methods |
|---|
|
|
|---|
Light and Electron Microscopy
A second pair of eyes from a 32-year-old normal donor was perfused
in similar conditions and used for morphology analysis. After 48 hours
in culture, whole anterior segments were fixed by overnight perfusion
with 2% glutaraldehyde at a pressure of 15 mm Hg. The anterior
segments were then cut meridionally into quadrants and immersed in the
same fixative until processing. From each of two quadrants,
wedge-shaped specimens containing the angle region with the TM were cut
and embedded in glycol methacrylate using an embedding kit (JB-4;
Polysciences, Warrenton, PA). Blocks were counterstained with eosin and
hematoxylin and sectioned to 2 µm. At least two sections from each of
the two quadrants of the processed eyes were analyzed by light
microscopy to examine the morphology of the tissue.
For electron microscopy, wedges of meridional sections of the TM were postfixed in 1% osmium tetroxide, block stained in 2% uranyl acetate, dehydrated, and embedded (Spurr embedding medium; Polysciences). Ultrathin sections were stained with potassium manganate and lead and examined with a transmission electron microscope (1200Ex; JEOL, Tokyo, Japan).
Construction of an HTM cDNA Library
RNA extraction and exponential amplification of the reverse
transcribed cDNA were performed as previously
described.24
25
Briefly, the TM from a single eye was
dissected and its RNA extracted with an RNeasy kit (Qiagen, Chatsworth,
CA). One forth of the total RNA was lyophilized and used as template
for the generation of a reversed transcribed first cDNA using a DNA
synthesis kit (SMART-PCR; Clontech, Palo Alto, CA) and Superscript II,
RNase H- reverse transcriptase (Gibco,
Gaithersburg, MD) to a final volume of 50 µl. The number of cycles
needed for exponential phase amplification of this cDNA was determined
by running a series of 16 µl analytical PCR amplifications at 17, 20,
and 23 cycles using the same kit. To prevent overrepresentation of
clones with short inserts (which would ligate to the vector more
efficiently than larger ones), the resultant cDNAs were size
fractionated using agarose gel electrophoresis. For that, the
PCR-amplified sample was purified and concentrated in an elution column
(QIAquick; Qiagen) to 30 µl and loaded onto a single well of a 2%
low-melting agarose gel. Four separate slices corresponding
approximately to molecular weights of 0.6 to 1 kb (fraction I), 1 to 2
kb (fraction II), 2 to 3 kb (fraction III), and more than 3 kb
(fraction IV) were cut from the gel and melted at 65°C for 10
minutes. The cDNA from each slice was purified and concentrated to a
total of 30 µl with an elution column (QIAquick; Qiagen). One tenth
of each eluted cDNA was used for ligation into a cloning vector
(pCR2.1-TOPO; Invitrogen, Carlsbad, CA) followed by transformation of
Escherichia coli competent cells (TOP10F' One Shot;
Invitrogen) and plating onto ampicillin-agar plates with
5-bromo-4-chloro-3-indoyl-ß-galactosidase (X-gal). A similar number
of white cell colonies from each fraction were randomly selected for
DNA sequencing. Glycerol stocks from each of these clones have been
kept under identification numbers HTM1-0001 to HTM1-1060.
Template Preparation and Sequencing
Plasmid purification and automatic DNA sequencing were performed
by Midland Certified Reagent (Midland TX). Briefly, plasmid DNAs were
purified by the alkaline lysis method and 0.5 to 1 µg used as
template for cycle sequencing reactions (Thermosequenase kit; Amersham,
Arlington Heights, IL) with IRD800-labeled M13 as the forward primer.
Electrophoresis was then conducted (model LS4200; Li-cor, Lincoln, NE).
For clones in which the M13 primer produced a poor-quality sequence, a
second reaction was prepared using the M13 reverse primer. Because the
SMART PCR library cloned by the TA method is not directional, the ESTs
were generated from either the 5' or the 3' end of the cDNA clones. In
the case of full-length cDNA sequencing, custom-made oligonucleotides
were made inside the clone insert, and the full sequence was obtained
by overlapping readings of both strands of the cDNA.
Sequence Analysis
Sequences corresponding to the vector and to the PRC primer used
to construct the library were deleted manually from each clone. All
sequences were arranged in a unique FASTA file and searched against
databases using the network BLAST client program
(BLASTcl3)29
of the National Center for Biotechnology
Information (NCBI; URL: http://www.ncbi.nlm.nih.gov/blast/blast.cgi).
The searches were performed between September and November 1999.
The first sequence search was conducted using BlastN against all entries in the nonredundant GenBank database. Sequences with an E-value expected (expected number of matches to be found merely by chance)30 equaling zero were considered to identify known genes or to have partial similarity to known genes. Sequences with E-values close to zero, were manually analyzed for overall similarity or identical stretches that could indicate a different gene from the same family, small sequencing errors, alternative splicing, chimeric clones or presence of vector sequences. Sequences including fragments with similarity with repetitive sequences (e.g., ALU, LINE) were considered perfect matches only when BLAST searches performed after eliminating the nonrepetitive part of the sequence also provided an E-value equal to zero.
DNA sequences for which no match was found in the nonredundant GenBank database were subjected to a second search with the BlastN against the dbEST database. From those, DNA sequences that matched dbESTs were searched against the Unique Human Gene Sequence Collection (UniGene) database (URL: http://www.ncbi.nlm.nih.gov/UniGene/) to determine whether they were integrated in a gene cluster. Expression patterns and chromosomal localization of these DNA sequences found to belong to a cluster were thus obtained.
Last, DNA sequences found to contain no similarities to any database were translated in all reading frames and further searched against the nonredundant protein database using the BLASTX program. The same DNA sequences were also searched for the presence of open reading frames using the ORF Finder program (http://www.ncbi.nlm.nih.gov/gorf/gorf.html).
Hybrid Panel Analysis
To determine the chromosomal location of the gene represented by
HTM1-0025 cDNA, we used the hybrid panel G3 from Stanford University
(Stanford, CA). We designed primers for the specific PCR amplification
of a fragment of the gene. Positive amplification was seen only with
the human genomic DNA control but not with the hamster DNA. After
analyzing PCR reactions from the 83 hybrids in the G3 panel, we
submitted the results to the SHGC server
(rhserver@paxil.stanford.edu).
| Results |
|---|
|
|
|---|
|
|
|
Other potentially relevant genes were present in this
library but did not show in the 1060 clones analyzed by
single-pass sequencing. For instance, the presence of
TIGR/MYOC as well as that of other olfactomedins was
confirmed by PCR, but none of the 1060 clones corresponded to these
genes. Despite the importance of TIGR/MYOC in glaucoma,
this result is not surprising. In our laboratory, Northern
blot analysis consistently shows that the expression of
this gene in the HTM is lower than that of
B-crystallin (represented by only one clone in the library).
Gene Expression Profile in HTM
Group I: Clones Similar to Known Genes.
The 519 group I sequences matched to known genes with
E-values of zero represented 338 different genes. The
description and accession number of each of these genes, along with
their distribution into nine different functional categories are
summarized in Table 2
. The number of
clones and cDNA species from each category are represented in Figure 3
.
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The percentage of ESTs obtained in group I (known genes), group II (dbESTs), and group III (novel sequences) in the HTM library was not substantially different from those in other tissues.16 17 18 22 23 The classification of the known genes by function (Table 2 , Fig. 3 ) also revealed a general pattern similar to that previously described for an average of 37 human tissues.33 The functional group showing a relatively higher representation in the TM library was the one for cell signaling and cell communication. This finding may be an indication of the importance of sensing changes in IOP and sending signals to other cells to respond to these changes. The group with a higher level of gene redundancy was the one involved in extracellular matrix synthesis and degradation, which could be a reflection of the importance of extracellular matrix dynamics in the TM metabolism. Although some genes coding for proteins involved in vesicle dynamics were present in the library, genes involved directly in phagocytosis were not highly represented. This could be the result of the perfusion model used. Because the perfusion medium did not contain aqueous humor proteins, DNA, and/or debris, the mechanisms to eliminate them may have been underregulated.
In group I (known genes), some of the observed most abundant transcripts are usually not present at such high frequency in other tissues and may therefore provide useful clues about the biologic features of TM. The high expression of other genes, such as the glycolytic enzymes LDHA, GAPDH, and TPI is expected, given that TM obtains most of its energy though anaerobic metabolism.34 More surprising is the high abundance of transcripts coding for proteins such as matrix GLA (MGP), chitinase 3-like 1 (CHI3L1), apolipoprotein D, small inducible cytokine (SCYA20), regulator of G protein signal, stromelysin 1, and from the uncharacterized genes KIAA0258 and DKFZp586O0118.
Matrix GLA is the first protein clearly associated with protection against calcification of soft tissues.35 36 High levels of matrix GLA expression are observed in response calcification of human vascular cells in vitro37 and in atherosclerotic lesions.38 39 40 41 The presence of matrix GLA as one of the more abundant transcripts in our library may indicate that protection against calcification is important in TM. Calcification is associated with the age-related loss of elasticity and contractility in some tissues and could have a similar effect on TM, impairing the ability of TM to regulate outflow. Of note, dexamethasone treatment that can result in reversible glaucoma, is known to induce calcification in blood vessels.42
Chitinase 3-like or human cartilage protein gp-39 is another extracellular protein expressed in cartilage and strongly induced in macrophages from atherosclerotic lesions.43 This protein is also one of the antigens implicated in the autoimmune disease, rheumatoid arthritis.44 45 Although its function is not yet clear, its involvement in tissue remodeling and in the pathologic course of blood vessel disease makes this gene a good candidate to play a potential role in TM conductivity and the pathology of glaucoma.
It is more difficult to speculate about the possible role of such proteins as ApoD, SCYA20, and RGS5, because of our incomplete knowledge about their function in other tissues. However, data suggesting the involvement of ApoD in neuronal degeneration and regeneration,46 SCYA20 in inflammation and vascular permeability,47 48 49 and RGS5 in the modulation of rectifying potassium (GIRK) channels50 51 suggest possible links to TM physiology.
Stromelysin 1 already has been identified as an important protein in the regulation of outflow resistance.52 53 54 Stromelysin 1 can be induced by a number of factors.55 56 57 58 59 60 However, it is not usually found as such an abundant transcript in noninduced tissues, with only 21 ESTs reported in the UniGene cluster. Its relative abundance in the TM library (approximately 0.4%) is consistent with the concept that the TMs ECM is under constant remodeling.
The presence of four clones representing each of the uncharacterized genes KIAA0258 and cDNA DKFZp586O0118 is particularly interesting, because none of them appears to be highly expressed in other tissues, and therefore they may be important for some specific functions of TM. This finding further emphasizes the need for high-throughput sequencing of HTM cDNAs. Some of the key genes essential for the physiology of this tissue may not yet have been fully characterized.
The identification of ESTs from group II (previously reported ESTs from uncharacterized genes), and group III (novel ESTs not previously reported in other tissues) is particularly important for two reasons: First, given the incomplete penetrance, the influence of nongenetic factors, and the late onset of the disease, an efficient search for glaucoma susceptibility genes necessitates not only positional cloning but also a candidate gene approach. The identification of novel genes expressed in HTM together with their chromosome location may greatly contribute to the identification of glaucoma genes. Second, because glaucoma does not appear to be associated with problems in other tissues, proteins involved in this disease may have a TM-preferred tissue expression pattern.
ESTs from uncharacterized genes (group II) can be analyzed for chromosome location and pattern of tissue expression using the information available in databases such as UniGene (Table 4) . Some of the ESTs found in our TM library are located in chromosomes with regions linked to glaucoma.61 62 63 As new linkage analysis provides more precise information about all glaucoma loci and more extensive analysis of TM libraries provides more information about ESTs from this tissue, it will be possible to select candidate genes from this group of clones. Furthermore, other ESTs from this group belong to UniGene clusters that appear to have a pattern of expression restricted to very specific tissues. This is the case of clusters such as Hs.97431 which includes only four ESTs, all from a testis library; Hs.132591, which includes five ESTs, all from brain; and Hs.199612, which has only one EST from kidney. This observation suggests that these genes may have some specific functions associated with the particular physiology of these tissues.
The identification of novel ESTs (group III) is particularly important in highly specialized and poorly studied tissues such as TM. However, analysis of novel, previously unreported ESTs has to be undertaken with some caution. cDNA libraries are never totally free of artifacts, which can include chimeric clones and/or fragments of genomic DNA. The presence of four chimeric clones in our library indicates that some of the novel clones could also be artifacts. However, four chimeras in 519 clones in the known gene pool represents a small percentage, probably not higher than that normally found in conventional non-PCRamplified libraries. More important, the characterization of clone HTM1-0025 shows that true novel genes can be discovered within the subset of novel ESTs. The features of this clone (complete ORF, similarity with a known family, and presence of the polyadenylation signal) confirm its identity as a novel gene and not as an artifact.
Finally, when evaluating single-pass sequencing analysis of cDNA libraries such as the one presented in this study, it is important to remember that any single library has some bias associated with the specific source of RNA. For instance, HTM samples acquired from postmortem donors may not accurately reflect the in vivo expression profile because of selective cell death (uveal layer), absence of IOP, and no aqueous flow. Using perfused TM solves some of these problems by recovering cell activity and maintaining IOP and aqueous flow. However, a limitation of this system is that the TM cells are exposed to tissue culture medium, not to physiological aqueous humor. They therefore do not have the exposure to some components, such as growth factors and reactive oxygen species, that is likely to affect their gene expression. In addition, it is important to consider the high degree of individual variability associated with the analysis of human samples. The particular history of the donor, medication, age, sex, and race also affect gene expression. Further analysis of more individuals is needed for a definitive characterization of the HTM gene expression profile.
In summary, our study describes the first expression profile of known and novel genes in a human TM. We identified 338 known genes, 120 ESTs previously reported in other tissues, and 189 new ESTs. The relative higher abundance of some of the gene species suggests new potential functions of the TM tissue, such as the prevention of cell calcification. The identification of novel ESTs should provide new tools for understanding TM physiology and for identifying genes involved in the genetic susceptibility to glaucoma.
| Acknowledgements |
|---|
| Footnotes |
|---|
Supported by National Eye Institute Grants EY11906 and EY01894, and Research to Prevent Blindness. TB is a Jules and Doris Stein Research to Prevent Blindness Professor Awardee.
Submitted for publication March 22, 2000; revised June 23, 2000; accepted July 5, 2000.
Commercial relationships policy: N.
Corresponding author: Teresa Borrás, Duke University Medical Center, Wadsworth Building, Erwin Road, Box 3802, Durham, NC 27710. borra001{at}mc.duke.edu
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
W. Xue, N. Comes, and T. Borras Presence of an Established Calcification Marker in Trabecular Meshwork Tissue of Glaucoma Donors Invest. Ophthalmol. Vis. Sci., July 1, 2007; 48(7): 3184 - 3194. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. W. Choy, C. C. Wang, A. Ogura, T. K. Lau, M. S. Rogers, K. Ikeo, T. Gojobori, D. S. C. Lam, and C. P. Pang Genomic annotation of 15,809 ESTs identified from pooled early gestation human eyes Physiol Genomics, March 13, 2006; 25(1): 9 - 15. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Xue, R. Wallin, E. A. Olmsted-Davis, and T. Borras Matrix GLA Protein Function in Human Trabecular Meshwork Cells: Inhibition of BMP2-Induced Calcification Process. Invest. Ophthalmol. Vis. Sci., March 1, 2006; 47(3): 997 - 1007. [Abstract] [Full Text] [PDF] |
||||
![]() |
J C H Tan, F B Kalapesi, and M T Coroneo Mechanosensitivity and the eye: cells coping with the pressure. Br. J. Ophthalmol., March 1, 2006; 90(3): 383 - 388. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. B. Liton, X. Liu, W. D. Stamer, P. Challa, D. L. Epstein, and P. Gonzalez Specific Targeting of Gene Expression to a Subset of Human Trabecular Meshwork Cells Using the Chitinase 3-Like 1 Promoter Invest. Ophthalmol. Vis. Sci., January 1, 2005; 46(1): 183 - 190. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Zhao, K. E. Ramsey, D. A. Stephan, and P. Russell Gene and Protein Expression Changes in Human Trabecular Meshwork Cells Treated with Transforming Growth Factor-{beta} Invest. Ophthalmol. Vis. Sci., November 1, 2004; 45(11): 4023 - 4034. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Wajih, T. Borras, W. Xue, S. M. Hutson, and R. Wallin Processing and Transport of Matrix {gamma}-Carboxyglutamic Acid Protein and Bone Morphogenetic Protein-2 in Cultured Human Vascular Smooth Muscle Cells: EVIDENCE FOR AN UPTAKE MECHANISM FOR SERUM FETUIN J. Biol. Chem., October 8, 2004; 279(41): 43052 - 43060. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Ahmed, M. Torrado, R. D. Zinovieva, V. V. Senatorov, G. Wistow, and S. I. Tomarev Gene Expression Profile of the Rat Eye Iridocorneal Angle: NEIBank Expressed Sequence Tag Analysis Invest. Ophthalmol. Vis. Sci., September 1, 2004; 45(9): 3081 - 3090. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Gonzalez, M. Caballero, P. B. Liton, W. D. Stamer, and D. L. Epstein Expression Analysis of the Matrix GLA Protein and VE-Cadherin Gene Promoters in the Outflow Pathway Invest. Ophthalmol. Vis. Sci., May 1, 2004; 45(5): 1389 - 1395. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. I. Tomarev, G. Wistow, V. Raymond, S. Dubois, and I. Malyukova Gene Expression Profile of the Human Trabecular Meshwork: NEIBank Sequence Tag Analysis Invest. Ophthalmol. Vis. Sci., June 1, 2003; 44(6): 2588 - 2596. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. R. Lo, L. L. Rowlette, M. Caballero, P. Yang, M. R. Hernandez, and T. Borras Tissue Differential Microarray Analysis of Dexamethasone Induction Reveals Potential Mechanisms of Steroid Glaucoma Invest. Ophthalmol. Vis. Sci., February 1, 2003; 44(2): 473 - 485. [Abstract] [Full Text] [PDF] |
||||