|
|
||||||||
1From the Schepens Eye Research Institute, and the 2Department of Ophthalmology, Harvard Medical School, Boston, Massachusetts; the 3Department of Physics, University of Massachusetts, Boston, Massachusetts; and the 4Department of Oral Medicine, Infection, and Immunity, Harvard School of Dental Medicine, Boston, Massachusetts.
| Abstract |
|---|
|
|
|---|
METHODS. Meibomian glands were obtained from orchiectomized mice that were treated with placebo or testosterone for 14 days. Tissues were processed for the analysis of differentially expressed mRNAs by using gene bioarrays, gene chips, and real-time PCR procedures. Bioarray data were analyzed with GeneSifter software (VizX Labs LLC, Seattle, WA).
RESULTS. The results show that testosterone influenced the expression of more than 1590 genes in the mouse meibomian gland. This hormone action involved a significant upregulation of 1080 genes (e.g., neuromedin B), and a significant downregulation of 518 genes (e.g., small proline-rich protein 2A). Some of the most significant androgen effects were directed toward stimulation of genes associated with lipid metabolism, sterol biosynthesis, fatty acid metabolism, protein transport, oxidoreductase activity, and peroxisomes.
CONCLUSIONS. These findings demonstrate that testosterone regulates the expression of numerous genes in the mouse meibomian gland and that many of these genes are involved in lipid metabolic pathways.
However, the mechanism(s) underlying this androgen influence on meibomian gland lipogenesis and function is unknown. It is possible that the hormonal regulation of this tissue is analogous to that of sebaceous glands, given that the meibomian gland is a large sebaceous gland and that androgens control the development, differentiation, and lipid production of these glands throughout the body.16 17 Androgen effects on sebaceous glands are mediated primarily through hormone binding to androgen receptors within acinar cell nuclei.17 18 19 This receptor interaction leads to increased gene transcription and the elaboration of proteins that stimulate the synthesis and secretion of lipids.17 18 19 20 In many sebaceous glands, androgen activity is also enhanced by, or dependent on, the presence of 5
-reductase, an enzyme that converts testosterone into the potent androgen, 5
-dihydrotestosterone.17
Consistent with this possibility are the findings that meibomian glands of males and females contain androgen receptor mRNA, androgen receptor protein within acinar epithelial cell nuclei, and the mRNAs for both types 1 and 2 5
-reductase.21 22 Given these observations, we hypothesized that the androgen control of the meibomian gland, as with other sebaceous glands, involves the regulation of gene expression. The purpose of the present study was to test this hypothesis.
| Materials and Methods |
|---|
|
|
|---|
Molecular Biological Procedures
To determine the effect of testosterone on meibomian gland gene expression, total RNA was isolated from tissues by using TRIzol reagent (Invitrogen Corp., Carlsbad, CA). When indicated, samples were also exposed to RNase-free DNase (Invitrogen), analyzed spectrophotometrically at 260 nm to determine concentration and examined on 6.7% formaldehyde/1.3% agarose (Gibco/BRL, Grand Island, NY) gels to verify RNA integrity. The RNA samples were then processed by utilizing several different technical approaches.
The principle method to evaluate differential gene expression involved the use of CodeLink Uniset Mouse I Bioarrays (
10,000 genes; Amersham, Piscataway, NJ). Before array studies, the integrity of glandular RNA preparations was further assessed with a RNA 6000 Nano LabChip with an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). The RNA samples were then processed for CodeLink Bioarray hybridization, as previously described.23 Briefly, cDNA was synthesized from RNA (2 µg) with a CodeLink Expression Assay Reagent Kit (Amersham) and isolated with a QIAquick purification kit (Qiagen, Valencia, CA). After sample drying, cRNA was generated with a CodeLink Expression Assay Reagent Kit (Amersham), recovered with an RNeasy kit (Qiagen) and quantified with an UV spectrophotometer. Fragmented, biotin-labeled cRNA was incubated and agitated (300 rpm shaker) on a CodeLink Bioarray at 37°C for 18 hours. The Bioarray was then washed and exposed to streptavidin-Alexa 647. Bioarrays were scanned by using ScanArray Express software and a ScanArray Express HT scanner (Packard BioScience, Meriden, CT) with the laser set at 635 nm, laser power at 100%, and photomultiplier tube voltage at 60%. Scanned image files were examined by using CodeLink image and data analysis software (Amersham), which generated both raw and normalized hybridization signal intensities for each array spot. The
10,000 spot intensities on the microarray image were standardized to a median of 1. Normalized data, with signal intensities exceeding 0.50, were analyzed with GeneSifter software (VizX Labs LLC, Seattle, WA; vizxlabs.com). This program also produced gene ontology and z-score reports. These ontologies, which were organized according to the guidelines of the Gene Ontology Consortium (http://www.geneontology.org/GO.doc.html),24 included biological processes, molecular functions and cellular components. Statistical analysis of individual gene expression data was conducted with Students t-test (two-tailed, unpaired). Data were evaluated with and without log transformation.
The data from the individual Bioarrays (n = 6) are accessible for downloading through the National Center for Biotechnology Informations Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) via series accession number GSE1582. The data will also be available for evaluation through GeneSifter (http://genesifter.net/datacenter/).
Differentially expressed mRNAs were also analyzed by using GEM 1 (>8,000 genes) and GEM 2 (>9,500 genes) gene chips (Incyte Genomics, Inc., St. Louis, MO). Poly(A) mRNA was isolated from meibomian gland RNA samples by using the MicroPoly(A) Pure mRNA Isolation Kit (Ambion, Inc., Austin, TX). The mRNA concentration was determined with a RiboGreen RNA Quantitation Kit (Molecular Probes, Eugene, OR), according to Incytes protocol. After designating mRNA samples (800 ng) for use with either cy3 or cy5 probes, preparations were suspended in TE buffer, placed in siliconized RNase-Free Microfuge Tubes (Ambion) and shipped on dry ice to Incyte for hybridization. Microarray data were sent electronically to the Harvard Center for Genomic Research (Cambridge, MA) and results were downloaded into the Resolver Gene Expression Data Analysis System, version 3.1 (Rosetta Inpharmatics, Kirkland, WA). This system displayed the sequence identification and description of all chip nucleotides, the signal strength of the treatment (i.e., testosterone) and control (i.e., placebo) channels, the relationship between the two channels in terms of ratio and fold change, the comparative P-value and information concerning various quality control fields. In addition, this system determined the error-weighted average ratios for each chip, and normalized data across chips, thereby permitting the combination of GEM 1 and 2 microarray results to achieve a stronger analysis of gene expression. The error model applied by Rosetta Resolver on Incytes microarrays has been described in the addenda of recent literature reports.25 26
To verify the differential expression of selected mRNAs, quantitative real-time PCR (qPCR) was utilized. cDNAs were transcribed from mRNA samples by employing SuperScript II Reverse Transcriptase (Invitrogen) and oligo dT priming (Promega, Madison, WI). Primers were designed by using Primer Express Software, version 1.5 (Applied Biosystems, Inc., Foster City, CA). Specificity of the primers was verified by performing BLASTn searches on all relevant NCBI nucleotide databases. Particular focus was placed on identifying primers with a 16- to 40-bp length, 20% to 80% GC content, and a melting temperature between 58°C and 60°C, that would generate amplicons between 140 bp and 160 bp. The qPCR was performed by utilizing the specific primers at optimal concentrations (Table 1) and Applied Biosystems SYBR Green PCR Master Mix, MicroAmp Optical 96-Well Reaction Plate, ABI PRISM Optical Adhesive Covers and GeneAmp 5700 Sequence Detection System, according to the manufacturers protocol. The instruments dissociation protocol did not show any secondary PCR products in any of the amplifications. Gene expression was determined by using either the Relative Standard Curve Method or the Comparative CT Method27 and standardizing levels to that of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) or tubulin,
1 mRNA.
|
| Results |
|---|
|
|
|---|
Evaluation of non- and log-transformed data from three separate experiments demonstrated that testosterone influenced the expression of more than 1590 genes in the mouse meibomian gland. This hormonal action involved a significant upregulation of 1080 genes (e.g., neuromedin B) and a significant downregulation of 518 genes (e.g., small proline-rich protein 2A; Table 2 ). Of particular interest was the finding that androgen treatment increased the activity of genes encoding various steroid receptors (e.g., types of estrogen, progesterone, and retinoic acid-binding sites), steroidogenic enzymes (e.g., 17ß-hydroxysteroid dehydrogenase 7), and endocrine factors (e.g., insulin-like growth factor 1; Table 3 ). Moreover, testosterone altered the expression of several immune-associated genes (e.g., caspase 7; Table 3 ).
|
|
)- and down (
) regulation of numerous genes associated with lipid metabolism (46
, 14
), lipid transport (8
, 1
), sterol biosynthesis (5
, 1
), fatty acid metabolism (16
, 4
), intracellular protein transport (31
, 10
), oxidoreductase activity (69
, 16
), peroxisomes (15
, 0
), mitochondria (76
, 13
), and early endosomes (5
, 0
; Tables 4 5 ). For comparison, testosterones effects were least directed toward processes such as cell adhesion, actin binding, and cytoskeleton (Table 5) .
|
|
|
To verify in part the bioarray (CodeLink; GE Healthcare) results, additional meibomian gland mRNA samples (n = 22 mice/group/experiment) were processed for gene chip (GEM 1 and 2; Incyte Genomics, Inc.) analyses. The gene chips and bioarrays have 4717 sequences in common. This approach showed that 474 genes were up (n = 319)- or downregulated (n = 155) by testosterone on both the bioarrays (P < 0.05) and gene chips (androgen/placebo ratio = > or <0; data not shown). If comparisons were restricted to those gene-chip genes that had expression ratios >1.5 (
or
), then 83 genes were identified as being similarly influenced by androgen on both platforms (e.g., Table 6 ).
|
|
| Discussion |
|---|
|
|
|---|
The mechanism by which testosterone influences meibomian gland gene expression undoubtedly involves an association with saturable, high-affinity, and steroid-specific receptors in acinar epithelial cell nuclei. Androgen receptors are members of the steroid/thyroid hormone/retinoic acid family of ligand-activated transcription factors and appear to mediate the classic actions of androgens throughout the body.28 29 After androgen binds to the receptor, the monomeric, activated hormonereceptor complex invariably associates with an androgen response element in the regulatory region of specific target genes; typically dimerizes with another androgen-bound complex; and, in combination with appropriate coactivators and promoter elements, controls gene transcription.28 29 In support of this hypothesis, it has been shown that androgen receptors exist in sebaceous gland epithelial cells17 18 19 and androgen activity in these cells may be compromised by androgen receptor defects or antagonists.30 31 32 Similarly, androgen receptors exist in meibomian gland epithelial cells,21 22 and androgen receptor disruption or the use of antiandrogen medications is associated with significant meibomian gland dysfunction and striking alterations in the neutral and polar lipid profiles of meibomian gland secretions.2 3 12 13 14 In addition, we have recently found that many androgen-regulated genes in the meibomian gland appear to depend on the presence of functional androgen receptors.33
It is important to note, though, that other processes may also be involved in, or mediate, androgen influence on meibomian gland gene expression. For example, the apparent androgen control of transcriptional activity may actually reflect hormone-induced alterations in mRNA stability,34 a possibility that remains to be explored. Another possibility is that testosterones impact on the meibomian gland is not direct, but rather is mediated through estrogen activity. The meibomian gland contains the mRNA for aromatase cytochrome P-450 (Schirra F, Suzuki T, Dickinson DP, Townsend DJ, Gipson IK, Sullivan DA, manuscript submitted) an enzyme that transforms testosterone into 17ß-estradiol.35 Moreover, the meibomian gland harbors estrogen receptor mRNA and protein.21 36 However, an estrogen mediation of androgen effects in the meibomian gland is highly unlikely. Recent research has shown that 17ß-estradiol treatment of ovariectomized mice elicits a pattern of gene expression in meibomian tissue that is dissimilar from that induced by testosterone (Suzuki T, Schirra F, Jensen RV, Richards SM, Sullivan DA, manuscript submitted).37 Furthermore, unlike androgens, estrogens appear to decrease sebaceous gland function,32 38 and this effect has been proposed to be due to an antagonism of androgen action.39 40
Our finding that androgens modulate gene expression in the mouse meibomian gland is consistent with our earlier preliminary observations in rabbits.41 Thus, by using RNA arbitrarily primed polymerase chain reactions, sequencing gels, and autoradiography, we were able to identify 58 differentially expressed mRNAs in the meibomian glands of orchiectomized rabbits treated topically with testosterone- or vehicle. However, analysis of 22 of the corresponding cDNA bands demonstrated that the majority had no significant homology to sequences in the GenBank database (http://www.ncbi.nlm.nih.gov/Genbank; provided in the public domain by the National Center for Biotechnology Information, Bethesda, MD), presumably due to limited data for rabbit sequences. Consequently, to identify genes regulated by androgens in the meibomian gland, we selected mice as an experimental model because of the extensive genetic information available for this species.
Considering that androgens modulate meibomian gland lipogenesis, the genes upregulated in this tissue by testosterone and the proteins they encode are particularly intriguing. Fatty acid synthase is a critical lipogenic enzyme that is known to be regulated by androgens in other tissues42 43 and is expressed in meibomian gland epithelial cells (Richards SM, et al. IOVS 2002;43:ARVO E-Abstract 3150). Fatty acid transport protein 4 facilitates the cellular uptake and metabolism of long- and very-long-chain fatty acids,44 whereas elongation of very-long-chain fatty acids-like 1 and 3 promote the tissue-specific synthesis of very-long-chain fatty acids and sphingolipids.45 46 These proteins could be involved in the androgen-induced increase of long-chain fatty acids in the total lipid fraction of rabbit meibomian glands.1 Monoglyceride lipase hydrolyzes tri- and monoglycerides to fatty acids and glycerol.47 Abca1 and Abcd3, which are members of the adenosine triphosphate (ATP)-binding cassette family, transport various molecules across extra- and intracellular membranes. Abca1 functions as a cholesterol efflux pump in the lipid-removal pathway48 and thereby serves as a key regulator of cholesterol distribution.49 50 Abcd3 modulates the importation of fatty acids and/or fatty acyl-CoAs into peroxisomes.51 3-Hydroxy-3-methylglutaryl-coenzyme A reductase is the rate-limiting enzyme of sterol biosynthesis.46 Oxysterol binding protein-like 1A, sterol carrier protein 2, liver, lipocalin 3, and phosphatidylcholine transfer protein are involved in the binding and/or transfer of phospholipids.46 However, whether these proteins play a definitive role in androgen-meibomian gland interactions has yet to be determined.
Androgens also were shown to control a series of genes that may be very important in the endocrine regulation of the meibomian gland. Thus, testosterone increased the mRNA levels of 17ß-hydroxysteroid dehydrogenase 7, a member of the enzyme family that regulates the interconversion of 17-ketosteroids with their corresponding 17ß-hydroxysteroids.52 This enzymatic activity is essential for the metabolism of all active androgens and estrogens in peripheral sites52 and may mediate the local, intracrine synthesis of androgens from adrenal precursors in the meibomian gland. Testosterone also enhanced the mRNA content of insulin-like growth factor 1, a pleiotropic protein that stimulates DNA synthesis and differentiation in sebaceous cells.53 Insulin-like growth factor 1 may also promote steroidogenesis54 and has been shown to be regulated by androgens in other tissues.55 Moreover, androgen treatment increased the expression of the gene encoding estrogen receptor 2 (ß). This receptor, which is upregulated by androgen in the prostate,56 may inhibit the activity of estrogen receptor 1 (
).57 Testosterone also elevated the mRNA levels of 11ß-hydroxysteroid dehydrogenase 1, an enzyme that catalyzes the conversion of cortisol to the inactive metabolite cortisone.46 Of interest, testosterone downregulated the gene expression of aldehyde dehydrogenase family 1, subfamily A3 (also called retinaldehyde dehydrogenase 3), an enzyme that stimulates retinoic acid biosynthesis.46
In addition to these actions, androgens promoted the expression of genes involved in the sorting (e.g., adaptor protein complexes, RAB9), trafficking (e.g., ADP-ribosylation factor 5, sorting nexin 2), and hydrolysis (e.g., cathepsin B) of proteins in various cellular locations, including the endosome, Golgi apparatus, endoplasmic reticulum, lysosome, proteasome, nucleus, and mitochondrion.46 Androgens also increased the mRNA levels of epimorphin (an extracellular protein that directs epithelial cell morphogenesis), neuromedin B (a bombesin-like peptide that stimulates epithelial cell proliferation),58 and phospholipases C-ß3 and -ß4 (mediators of the production of the second-messenger molecules diacylglycerol and inositol 1,4,5-trisphosphate).46 Testosterone decreased the mRNA amounts of Ia-associated invariant chain, which plays an essential role in major histocompatibility (MHC) class II antigen processing.46
In summary, our results show that testosterone regulates the expression of a number of genes in the mouse meibomian gland and that many of these genes are involved in the production, metabolism, transport, and release of lipids, as well as in steroidogenic pathways. We are currently attempting to determine the meibomian gland distribution of these genes (e.g., by in situ hybridization), in order to identify the cellular targets for androgen action. This procedure will also verify mRNA location within the gland, compared with the conjunctiva, given that very small parts of this latter tissue adhered to the meibomian gland during isolation. In concert with these studies, we are endeavoring to determine whether a variety of these hormone-regulated genes are translated (e.g., by immunohistochemistry and Western blot analysis). This combined information may help to explain, at least in part, the mechanism by which topical androgens reportedly stimulate the synthesis and secretion of meibomian gland lipids, prolong the tear film breakup time, and alleviate dry eye.7 8
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication April 4, 2005; revised May 25, 2005; accepted August 15, 2005.
Disclosure: F. Schirra, None; T. Suzuki, Allergan, Inc. (F); S.M. Richards, None; R.V. Jensen, None; M. Liu, Allergan, Inc. (F); M.J. Lombardi, None; P. Rowley, None; N.S. Treister, None; D.A. Sullivan, Allergan, Inc. (C, F)
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: David A. Sullivan, Schepens Eye Research Institute, 20 Staniford Street, Boston, MA 02114; sullivan{at}vision.eri.harvard.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. D. Sullivan, J. E. Evans, M. R. Dana, and D. A. Sullivan Influence of aging on the polar and neutral lipid profiles in human meibomian gland secretions. Arch Ophthalmol, September 1, 2006; 124(9): 1286 - 1292. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |