PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of iovsIOVSARVO
 
Invest Ophthalmol Vis Sci. Apr 2011; 52(5): 2727–2740.
Published online Apr 25, 2011. doi:  10.1167/iovs.10-6482
PMCID: PMC3088560
Changes in Gene Expression in Human Meibomian Gland Dysfunction
Shaohui Liu,1,2 Stephen M. Richards,1,2 Kristine Lo,2,3 Mark Hatton,1,2,4 Aaron Fay,2,3 and David A. Sullivancorresponding author1,2
From the 1Schepens Eye Research Institute,
the 2Department of Ophthalmology, and
the 3Massachusetts Eye and Ear Infirmary, Harvard Medical School, Boston, Massachusetts; and
4Ophthalmic Consultants of Boston, Boston, Massachusetts.
corresponding authorCorresponding author.
Corresponding author; David A. Sullivan, Schepens Eye Research Institute, 20 Staniford Street, Boston, MA, 02114; david.sullivan/at/schepens.harvard.edu.
Received August 27, 2010; Revised December 15, 2010; Revised February 13, 2011; Accepted March 1, 2011.
Purpose.
Meibomian gland dysfunction (MGD) may be the leading cause of dry eye syndrome throughout the world. However, the precise mechanism(s) underlying the pathogenesis of this disease is unclear. This study was conducted to identify meibomian gland genes that may promote the development and/or progression of human MGD.
Methods.
Lid tissues were obtained from male and female MGD patients and age-matched controls after eyelid surgeries (e.g., to correct entropion or ectropion). Meibomian glands were isolated and processed for RNA extraction and the analysis of gene expression.
Results.
The results show that MGD is associated with significant alterations in the expression of almost 400 genes in the human meibomian gland. The levels of 197 transcripts, including those encoding various small proline-rich proteins and S100 calcium-binding proteins, are significantly increased, whereas the expression of 194 genes, such as claudin 3 and cell adhesion molecule 1, is significantly decreased. These changes, which cannot be accounted for by sex differences, are accompanied by alterations in many gene ontologies (e.g., keratinization, cell cycle, and DNA repair). The findings also show that the human meibomian gland contains several highly expressed genes that are distinct from those in an adjacent tissue (i.e., conjunctival epithelium).
Conclusions.
The results demonstrate that MGD is accompanied by multiple changes in gene expression in the meibomian gland. The nature of these alterations, including the upregulation of genes encoding small proline-rich proteins and S100 calcium-binding proteins, suggest that keratinization plays an important role in the pathogenesis of MGD.
Meibomian glands play a critical role in the health and well-being of the ocular surface.16 These glands secrete a lipid and protein mixture that provides a clear optical surface for the cornea, interferes with bacterial colonization, and retards tear overflow.18 The glandular secretions also enhance the stability and reduce the evaporation of the tear film.17 Conversely, meibomian gland dysfunction (MGD) destabilizes the tear film, increases its evaporation and osmolarity and is believed to be the key trigger for the induction of evaporative dry eye syndrome.16,814
The major cause of MGD appears to be excretory duct obstruction, due to hyperkeratinization of the ductal epithelium and an increased viscosity of meibum.1522 This obstruction may lead to cystic dilatation of glandular ducts, acinar cell atrophy, and a loss of secretory meibocytes.15,23 The MGD may also facilitate bacterial growth on the lid margin2426 and promote inflammation in the adjacent conjunctiva.27,28 The development of MGD has been linked to several risk factors, including aging,2934 androgen deficiency,3542 isotretinoin treatment,4345 and possibly postmenopausal estrogen therapy.4649 However, the molecular mechanisms that underlie the pathogenesis of MGD are unknown. This lack of information, in turn, has hampered the generation of safe and effective therapies for the treatment of MGD.
We hypothesize that alterations in gene expression promote the development and/or progression of human MGD. We also hypothesize that identification of such genetic changes will provide unique insight into the pathogenesis of MGD and will reveal novel glandular targets for possible therapeutic intervention. The purpose of this investigation was to begin to test these hypotheses by taking advantage of new advances in cDNA microarray technology, computational biology, and bioinformatics.
Human Subjects, Clinical Evaluation, and Tissue Collection
Human eyelid tissues were obtained from adult individuals with MGD (mean age, 73 ± 4 years) and from age-matched controls (mean age, 65 ± 9 years), after lid resection surgery at the Massachusetts Eye and Ear Infirmary (Table 1). Before surgery, normal or dysfunctional meibomian glands were diagnosed in the participants by two ophthalmologists (i.e., one evaluated two controls and two patients; the other four controls and four patients), according to a published classification system.44 In brief, the grading scheme was 0 for clear excreta with or without small particles, 1 for opaque cloudy excreta with normal viscosity, 2 for opaque excreta with increased viscosity, and 3 for secreta that retained shape after digital expression. Subjects were excluded from the study if they had active infection, used topical antiglaucoma or anti-inflammatory medications, or wore contact lenses. Lid segments were placed in RNA stabilizer (RNAlater; Ambion, Austin, TX) and stored at −80°C. Meibomian glands (2–5 glands/tarsal plate) were then isolated under a dissecting microscope (Bausch & Lomb, Rochester, NY), by the removal of skin, subcutaneous tissue, muscle, and palpebral conjunctiva, and were processed in entirety for molecular biological procedures. The use of human eyelid samples that would otherwise have been discarded was approved by the Institutional Review Boards of Massachusetts Eye and Ear Infirmary and Schepens Eye Research Institute and adhered to the tenets of the Declaration of Helsinki.
Table 1.
Table 1.
Age, Sex, Meibomian Gland Secretion Quality, and Medication History of Subjects
Molecular Biological Procedures
Total RNA was extracted (RNAqueous Kit; Ambion) and analyzed (RNA Nano 6000 Series II Chip with a 2100 Bioanalyzer; Agilent Technologies, Palo Alto, CA) to verify RNA integrity. The RNA concentrations and 260/280-nm ratios were determined with spectrophotometer (a NanoDrop 1000; Thermo Scientific, Waltham MA).
The RNA (100 ng) samples were processed by Asuragen (Austin, TX) for the quantitation of mRNA levels by using microarray expression analysis (HumanHT-12 ver. 3 Expression BeadChips; Illumina, San Diego, CA). The chips target greater than 25,000 annotated genes with more than 48,000 probes derived from NCBI (National Center for Biotechnology Information, Bethesda, MD) reference sequences and the UniGene databases (http://www.ncbi.nlm.nih.gov/UniGene; NCBI). Briefly, biotin-labeled cRNA samples were prepared (MessageAmp II; Ambion, Inc.), quantitated by UV spectrophotometry and examined with a capillary electrophoresis system (2100 Bioanalyzer; Agilent). The labeled cRNAs were used to probe the chips. Hybridization, washing, and scanning of the arrays were performed according to the manufacturer's instructions. Data were processed with appropriate software (BeadStudio software ver. 3; Illumina), by using background subtraction and cubic spline normalization. Standardized hybridization intensity values were adjusted by adding a constant, such that the lowest intensity value for any sample was 16.51 Hybridization intensity measurements are generally relied on in microarray experiments to infer transcript abundance levels for specific genes.51
Normalized data were evaluated with commercial software (GeneSifter.Net; Geospiza, Seattle, WA). This comprehensive program also produced gene ontology and z-score reports. Ontologies encompassed biological processes, molecular functions, and cellular components and were organized according to the guidelines of the Gene Ontology Consortium (http://www.geneontology.org/).52 Gene expression data were analyzed, with and without log transformation, and statistical evaluation of these data was performed with Student's t-test (two-tailed, unpaired). Genes that were expressed in similar directions in different experiments were identified by using the intersector program (GeneSifter.Net; Geospiza; www.public.genesifter.net). All data from the microarray chips (BeadChips; Illumina) are accessible for download through the National Center for Biotechnology Information's Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo) via series accession number (GSE17822).
Real-Time PCR Procedures
The differential expression of selected genes was confirmed by using quantitative real-time PCR (qPCR) procedures. The cDNAs were transcribed by using a transcriptase (SuperScript III Reverse Transcriptase; Invitrogen, Grand Island, NY) and random hexamer primers (Invitrogen). Real-time PCR reactions (n = 3 control and 5 MGD separate samples) were performed in triplicate (TaqMan Gene Assays; Applied Biosystems, Inc. [ABI], Foster City, CA) and specific primers and probes for small proline-rich protein 2A (Hs03046643_s1), keratin 10 (Hs00166289_m1), S100 calcium-binding protein A9 (S100A9; Hs00610058_m1), S100 calcium-binding protein A8 (Hs00374264_g1), and β-actin endogenous control (4326315E). Differential gene expression was calculated according to the comparative CT method (outlined in ABI User Bulletin 2; updated 2001).
Gene Expression
Analysis of gene expression in human meibomian glands (n = 12 samples) led to the identification of 11,173 genes that exceeded the manufacturer's quality value of 0.95 (Illumina). Of these genes, 423 were associated with lipid metabolic processes, 185 with lipid biosynthesis, 116 with fatty acid metabolic processes, 89 with phospholipid metabolic processes, and 66 genes with lipid transport.
The intensity range for gene expression spanned from 29 to 29,219 units, with the mean values of 7 genes averaging above 20,000 units, 47 genes over 10,000 units, 204 genes above 3,000 units, and 689 genes between 1,000 to 10,000 units. Examples of high-intensity genes are shown in Table 2. These include transcripts such as stearoyl-CoA desaturase, mucin-like 1, keratinocyte differentiation-associated protein, fatty acid synthase, and S100 calcium-binding proteins A8 and A9 (S100A8, S100A9). Several of these glandular transcripts have recently been discovered as proteins in human meibomian gland secretions,8 including adenylate cyclase 7; dystrophin, Ellis-van Creveld syndrome protein; keratin-1, -5, -6, -7, -9, -10, -13, -14, and -16; the lacrimal proline-rich protein proline-rich 4 (similar to 1); lipocalin 1; seven transmembrane helix receptor mucin-like 1 (also termed small breast epithelial mucin); synaptonemal complex protein 3; transcription elongation factor A (SII)-like 1; UDP-N-acetylglucosamine pyrophosphorylase 1; and zinc metalloproteinase.
Table 2.
Table 2.
Gene Expression in Human Meibomian Glands
Our data also demonstrate that human meibomian gland gene expression is not the same as that of another lid tissue, human conjunctival epithelium. We compared 300 sequences with the highest average intensity in our samples' raw data to those from human conjunctival epithelia (Gene Expression Omnibus Series accession number GSE5543).53 By this process we found that 42.3% of the sequences matched and most were for ribosomal genes. If ribosomal genes were excluded, then only 24.1% of the remaining sequences were similar between the meibomian gland and conjunctiva. In effect, the most highly expressed genes in the human meibomian gland are unlike those in human conjunctival epithelia.
Effect of MGD on Gene Expression
Evaluation of non- and log-transformed data from controls subjects (n = 6) and MGD patients (n = 6) showed that significant differences exist in meibomian gland gene expression between these groups. As demonstrated in Table 3, MGD was associated with a significant increase in the expression of 197 genes and a significant decrease in the activity of 194 genes, compared with that of the control population. Particularly notable was the MGD-related upregulation of genes for small proline-rich (SPRR) proteins 1A, 2A, 2E, 2F, and 3; cornulin; cystatin A; S100A8; and S100A9 (Table 4).
Table 3.
Table 3.
Influence of MGD on Gene Expression
Table 4.
Table 4.
Effect of MGD on Gene Expression Ratios within the Human Meibomian Gland
Many of the genes upregulated in MGD, as shown in Table 4, are located on chromosome 1. Significantly increased gene expression in MGD was also found on chromosomes 12 and 14, whereas higher gene activity was observed on chromosomes 13 and 17 in control tissues (Table 5).
Table 5.
Table 5.
Impact of MGD on Chromosomal Gene Expression in Human Meibomian Glands
To confirm in part the microarray results, we analyzed selected genes by qPCR. This experimental approach verified the significant (P < 0.05) upregulation of MGD on all tested genes, including small proline-rich protein 2A (2.6-fold ↑), S100A8 (12.8-fold ↑), S100A9 (5.9-fold ↑), and keratin 10 (4.6-fold ↑).
Influence of MGD on the Expression of Gene Ontologies
The influence of MGD was quite striking in its stimulation of various gene ontologies (Table 6), such as keratinization (Table 7); cell cycle (Table 8); DNA repair (Table 9); and energy generation, translation, structural molecule activity, cytoplasmic vesicles, and cornified envelope components (Table 10). In contrast, meibomian glands of control subjects presented with a significant increase in ontologies associated with cellular nitrogen compound metabolic processes, transport, transferase activity, iron and manganese binding, and intracellular membrane-bound organelles (Table 6).
Table 6.
Table 6.
Influence of MGD on the Expression of Gene Ontologies in the Human Meibomian Gland
Table 7.
Table 7.
Effect of MGD on Genes Related to Keratinization
Table 8.
Table 8.
Influence of MGD on Genes Associated with the Cell Cycle
Table 9.
Table 9.
Impact of MGD on Genes Associated with DNA Repair
Table 10.
Table 10.
Effect of MGD on Genes Associated With Energy Generation, Translation, Structural Molecule Activity, Cytoplasmic Vesicles and Cornified Envelope Components in the Human Meibomian Gland
Examination of the data also showed that certain gene ontologies were overrepresented in both control and dysfunctional meibomian glands, including nucleotide metabolic processes, proteasomal ubiquitin-dependent protein catabolic processes, DNase activity and cytoplasm (Table 6). Other gene ontologies were underrepresented in both groups, such as cell communication and signal transducer activity (Table 6).
Impact of Sex on Gene Expression
Given that the control and patient groups were not sex matched (Table 1), it is possible that sex differences in meibomian gland gene expression may have contributed to the MGD-related results. More specifically, because more females were present in the MGD group, it may be that our data reflect female- rather than disease-associated genes. To explore this possibility, we compared meibomian gland gene expression between men (n = 5) and women (n = 7) and then compared these data to those between controls and patients.
Our analyses of non– and log-transformed data demonstrate that the sex of the individual has a significant impact on gene expression in human meibomian glands. We identified 257 genes that show sex-related differences, with 137 genes more highly expressed in men and 120 in women. Examples of these genes, which include prolactin-induced protein (male > female) and S100 calcium-binding protein A7 (female > male), are depicted in the Table 11. Sex also has a significant effect on the expression of gene ontologies in the human meibomian gland, such as processes involving Golgi vesicle transport, translation, amino acid ligase activity, protein binding, endoplasmic reticulum, and ribosomes (Table 12).
Table 11.
Table 11.
Influence of Sex on Gene Expression Ratios within the Human Meibomian Gland
Table 12.
Table 12.
Effect of Sex on the Expression of Gene Ontologies in the Human Meibomian Gland
However, these sex-related differences do not account for the MGD-associated changes in gene expression. Comparison of these data sets showed that only four genes (e.g., zinc finger protein 658 and S100 calcium-binding protein A7) were upregulated, and only four genes were downregulated (e.g., der1-like domain family, member 2, and esterase D/formylglutathione hydrolase), in meibomian glands of both women (versus men) and patients (versus controls). Moreover, comparison of the ontologies listed in Tables 6 and and1212 demonstrate that, of the 40 biological process (36 different, 4 the same), 19 molecular function (16 different, 3 the same), and 31 cellular component (25 different, 6 the same) ontologies, the majority (n = 77) were not the same (n = 13) in both women and patients.
Our study shows that MGD exerts MGD exerted a significant effect on gene expression in the human meibomian gland. Thus, MGD was found to be associated with significant alterations in the expression of almost 400 genes, including those encoding SPRR, S100A8,and S100A9. These changes, which cannot be accounted for by sex-associated differences between control and patient populations, are accompanied by the up- and downregulation of many ontologies (e.g., keratinization, cell cycle, DNA repair, translation, cytoplasmic vesicles, and cornified envelope components). Our study also demonstrates that the human meibomian gland contains several highly expressed genes that are distinct from those in an adjacent tissue (i.e., conjunctival epithelium). Overall, the findings support our hypothesis that alterations in gene expression promote the development and/or progression of human MGD. Our data also reveal novel glandular targets for possible therapeutic intervention in MGD.
We identified several SPRR, S100A8, and S100A9 genes that are significantly upregulated in human MGD. The SPRRs are precursor proteins of the cornified cell envelope of terminally differentiating stratified squamous epithelia and influence the flexibility, rigidity, or toughness of this envelope.5559 The increased expression of these SPRR genes reflect the stimulation of meibomian gland ontologies associated with keratinization; keratinocyte and epidermal cell differentiation; ectoderm, epidermis, and tissue development; and cornified envelope. Of particular interest, the SPRRs, whether alone55 or in combination with S100A8 and -A9,60 are known to promote keratinization, such as in the conjunctiva61,62 and elsewhere in the body.63,64 Further, increased S100A8 and S100A9 expression induces proliferation of human keratinocytes.65 If translated, these SPRR, S100A8, and S100A9 proteins may be responsible, at least in part, for the hyperkeratinization of meibomian gland ductal epithelium. This process is believed to be a primary cause of MGD.16,814
It is possible that SPRRs and S100A8 and S100A9 represent novel and unique targets for the potential topical treatment of MGD. First, downregulation of SPRRs may represent a mechanism to inhibit the keratinization process. Second, attenuation of S100A8 and S100A9 may permit suppression of both the keratinization and the possible conjunctival inflammation associated with MGD. The proteins S100A8 and S100A9 form a heterodimeric protein complex called calprotectin, which is a Toll-like receptor 4 agonist and a damage associated molecular pattern molecule (DAMP).66 A DAMP can initiate and perpetuate an immune response in noninfectious inflammation.66 Typically, calprotectin may play an important role in regulating cellular homeostasis (e.g., cytoskeleton), but after secretion may become a danger signal during inflammatory processes.66 The secreted S100A8 and -A9 amplifies inflammatory responses66,67 and could induce extravasation of leukocytes into the adjacent conjunctiva. Such conjunctival inflammation may occur during MGD.27,28 For comparison, the S100A8 and -A9 proteins may be involved in the pathophysiology of pterygium,68,69 as well as in Sjögren's syndrome, rheumatoid arthritis, inflammatory bowel disease, vasculitis, lichen planus, and lupus erythematosus.7073
However, although calprotectin is commonly identified as a promoter and/or marker of inflammation,67,7476 it may also act, depending on its concentration (e.g., high), in an anti-inflammatory manner.70,77 In fact, calprotectin may exert a variety of protective antioxidant, anti-invasive, and antimicrobial functions.70,77,78 An anti-inflammatory action of calprotectin would be especially intriguing, given that the most highly expressed gene in the human meibomian gland is leukocyte-associated immunoglobulin-like receptor (LAIR)-1. The LAIR-1 protein is an inhibitory receptor that attenuates immune cell activation and suppresses proinflammatory cytokine production.79,80 Working in tandem, calprotectin and LAIR-1 could theoretically be responsible for the typical absence of inflammation20,23,81 within the meibomian gland in MGD. In addition, the ability of calprotectin to make epithelial cells more resistant to bacterial invasion82 could be one reason why such bacterial infection is minimal or absent in obstructive MGD.20 Overall, calprotectin could possibly protect the meibomian gland against inflammation and infection in MGD. Whether it exerts either or both of these activities or is instead proinflammatory remains to be determined.
Additional meibomian gland genes that are modulated in MGD may play a role in these keratinization, antimicrobial, and structural processes. For example, the upregulated genes cornulin and cystatin A may promote epidermal differentiation,83,84 and the increased peptidase inhibitor 3, skin-derived, transcripts, if translated, may inhibit bacterial infection.84 In addition, the genes encoding claudin 3, cell adhesion molecule 1, and CD151 molecule are downregulated in MGD, a response that may disrupt cell adhesion.84 Our data also show that many of the most highly upregulated genes in MGD are from chromosome 1. In fact, the SPRR, S100A8, and S100A9 genes colocalize on chromosome 1 at q21, the region called the epidermal differentiation complex.60 There are more than 850 diseases related to this chromosome, including deafness, Alzheimer's, glaucoma, and breast cancer.85
The influence of MGD is striking in its stimulation of various gene ontologies, such as cell cycle, DNA repair, generation of precursor metabolites, and translation. These responses may be analogous to those in other sebaceous glands, wherein stem cell niches are mobilized to repair injured or diseased tissue.86,87 Of interest, a deficiency in the capacity to repair DNA may lead to neoplasia.88,89
A notable finding was our discovery that significant sex-related differences exist in gene expression in the human meibomian gland. Genes encoding such proteins as lysozyme and neurocalcin δ were upregulated in males (M), whereas leptin gene expression was increased in females (F). Analogous sex-associated differences have been observed for lysozyme in stimulated submandibular/sublingual saliva (M > F),90 neurocalcin δ in the olfactory bulb (M > F),91 and leptin in the plasma (F > M).92 More importantly, sex-associated variations are also known to occur in the morphologic appearance, neutral and polar lipid profiles, and secretory output of the meibomian gland, as well as in the gene expression of murine meibomian tissues.31,3337,39,40,42,9398 Many of these differences appear to be due to the effects of sex steroid hormones and sex-specific aromatase activities (Liu S, et al. IOVS 2007;43:ARVO E-Abstract 5657).49,98 It is possible that these sex-related differences in the meibomian gland contribute to the increased prevalence of dry eye disease in females.11
Our data demonstrate that the human meibomian gland contains several highly expressed genes that are associated with lipid dynamics, glandular structure, immunologic factors, and hormone action. For instance, these genes encode proteins involved in cholesterol transport (Niemann-Pick disease, type C2); cholesterol (farnesyl-diphosphate farnesyltransferase 1; sterol-C4-methyl oxidase-like; sterol-C5-desaturase-like), fatty acid (stearoyl-CoA desaturase; thyroid hormone responsive SPOT14; fatty acid synthase), and lipid (acyl-CoA synthetase short-chain family member 2; transaldolase 1) biosynthesis; fatty acyl-CoA reduction (fatty acyl CoA reductase 2) and elongation (hydroxysteroid (17β) dehydrogenase 12); and glycosphingolipid catabolism (prosaposin).84,99 The primary product of stearoyl-CoA desaturase is oleic acid, which has been identified in meibomian gland secretions100,101 and shown to be decreased in meibomianitis.102 Of interest, androgens upregulate stearoyl-CoA desaturase mRNA levels in mouse male and female meibomian glands,49,103 and the targeted disruption or this rate-limiting enzyme causes meibomian gland atrophy.104 Additional highly expressed genes are those related to extracellular matrix stabilization (inter-α (globulin) inhibitor H5), cell proliferation, migration and/or differentiation (clusterin associated protein 1; thymosin β4, X-linked; keratinocyte differentiation-associated protein; annexin A2), and proteasomal degradation (ubiquitin B and C).84,105,106 One of these structure-related genes encodes biglycan, which, if overexpressed, leads to failure of meibomian gland formation in mice.107 Immune-associated genes that present with high activities in the human meibomian gland are associated with antigen presentation (β2-microglobulin), dendritic and T-cell migration (chemokine (C-C motif) receptor 6), bacteriolytic action (lysozyme), and the inflammatory response (CCAAT/enhancer binding protein β, S100A8, S100A9, and LAIR-1).84 Highly expressed genes related to the endocrine system include those that may bind androgens (secretoglobin, family 2A, member 1) and catalyze the transformation of estrone into estradiol (hydroxysteroid [17β] dehydrogenase 12).84,108,109 Androgens are known to upregulate the prostatic levels of secretoglobin 2A1, a lipophilin that is also translated and strongly expressed in human meibomian gland epithelial cells.110
Several proteins encoded by meibomian gland genes have been identified in meibomian gland secretions.8 In this regard, it is noteworthy that S100A8, S100A9, prolactin-inducible protein, lipocalin-1, lysozyme, β2-microglobulin, and the lacrimal proline-rich protein proline-rich 4 have all been termed tear film biomarkers in dry eye,111,112 and some of these are proposed to be indicative of lacrimal gland aqueous deficiency.112 However, given that the meibomian gland contains transcripts for all these proteins, it may be that alterations in their translation and secretion influence the tear film levels.
In this study, a number of the subjects were taking medications, including lisinopril (antihypertensive), atorvastatin (lowers blood cholesterol levels), atenolol (β1 receptor selective antagonist for cardiovascular disease), levothyroxine (hypothyroid treatment), alendronate (osteoporosis therapy), escitalopram (antidepressant), hydrochlorothiazide (diuretic), and clonazepam (antiepileptic). It is possible that some of these pharmaceutical agents influenced the secretion quality and/or gene expression in the meibomian gland. For example, atorvastatin inhibits a rate-limiting enzyme involved in the cholesterol synthetic pathway113 and could theoretically alter the melting point and viscosity of meibomian gland secretions.5 Similarly, thyroxine could possibly influence meibomian gland activity, given that this hormone has been reported to regulate the function of other sebaceous glands,114 and antidepressant use is a risk factor for the development of dry eye syndrome.42 However, none of these pharmaceutical compounds has been shown to have any effect on the meibomian gland, and neither statins nor antihypertensive medications promote dry eye.42 Consequently, it remains to be determined whether any of these medications represented a confounding factor in this investigation.
Last, some limitations in this study should be considered. First, it is possible that very small pieces of adjacent lid tissue adhered to the meibomian glands during isolation under the dissecting microscope. However, such possible contamination would be similar for all tissues, given that the same individual removed glands from every lid segment in the same manner. Second, the number of tissue samples (i.e., five male and seven female) is relatively low for the analysis of sex-related differences in gene expression. Increasing the sample number is difficult, however, unless tissues are available from multiple lid surgeries. The tissue collection process for this study took more than 2 years to complete. Third, our data reflect mRNA levels. We have yet to determine whether these transcripts are translated into functional proteins.
Overall, our findings demonstrate that MGD is accompanied by multiple changes in gene expression in the meibomian gland. The nature of these alterations, including the upregulation of genes encoding SPRR and S100 calcium-binding proteins, suggest that keratinization plays an important role in the pathogenesis and/or progression of MGD.
Footnotes
Supported by Grant EY05612 from the National Institutes of Health and by Alcon Research, Ltd.
Disclosure: S. Liu, P; S.M. Richards, None; K. Lo, None; M. Hatton, None; A.M. Fay, None; D.A. Sullivan, Alcon (F), P
1. McCulley JP, Shine WE. Meibomian gland function and the tear lipid layer. Ocul Surf. 2003;1:97–106. [PubMed]
2. Foulks GN, Bron AJ. Meibomian gland dysfunction: a clinical scheme for description, diagnosis, classification, and grading. Ocul Surf. 2003;1:107–126. [PubMed]
3. Bron AJ, Tiffany JM, Gouveia SM, Yokoi N, Voon LW. Functional aspects of the tear film lipid layer. Exp Eye Res. 2004;78:347–360. [PubMed]
4. Bron AJ, Sci FM, Tiffany JM. The contribution of meibomian disease to dry eye. Ocul Surf. 2004;2:149–165. [PubMed]
5. Driver PJ, Lemp MA. Meibomian gland dysfunction. Surv Ophthalmol. 1996;40:343–367. [PubMed]
6. Tiffany J. Physiological functions of the meibomian glands. Prog Retin Eye Res. 1995;14:47–74.
7. Knop E, Knop N, Schirra F. Meibomian glands; Part II, physiology, characteristics, distribution and function of meibomian oil (in German). Ophthalmologe. 2009;106:884–892. [PubMed]
8. Tsai P, Evans JE, Green KM, et al. Proteomic analysis of human meibomian gland secretions. Br J Ophthalmol. 2006;90:372–377. [PMC free article] [PubMed]
9. Shimazaki J, Sakata M, Tsubota K. Ocular surface changes and discomfort in patients with meibomian gland dysfunction. Arch Ophthalmol. 1995;113:1266–1270. [PubMed]
10. Lemp MA. Report of the National Eye Institute/Industry Workshop on clinical trials in dry eyes. CLAO J. 1995;21:221–232. [PubMed]
11. 2007 Report of the International Dry Eye WorkShop (DEWS). Ocul Surf. 2007;5:65–204.
12. Lemp MA, Nichols KK. Blepharitis in the United States 2009: a survey-based perspective on prevalence and treatment. Ocul Surf. 2009;7:S1–S14. [PubMed]
13. Shine WE, McCulley JP. Keratoconjunctivitis sicca associated with meibomian secretion polar lipid abnormality. Arch Ophthalmol. 1998;116:849–852. [PubMed]
14. Mathers W. Evaporation from the ocular surface. Exp Eye Res. 2004;78:389–394. [PubMed]
15. Knop E, Knop N. Meibomian glands: part IV. Functional interactions in the pathogenesis of meibomian gland dysfunction (MGD) (in German). Ophthalmologe. 2009;106:980–987. [PubMed]
16. Knop E, Knop N, Brewitt H, et al. Meibomian glands, Part III: Dysfunction—argument for a discrete disease entity and as an important cause of dry eye (in German). Ophthalmologe. 2009;106:966–979. [PubMed]
17. Jester JV, Rajagopalan S, Rodrigues M. Meibomian gland changes in the rhino (hrrhhrrh) mouse. Invest Ophthalmol Vis Sci. 1988;29:1190–1194. [PubMed]
18. Jester JV, Nicolaides N, Smith RE. Meibomian gland studies: histologic and ultrastructural investigations. Invest Ophthalmol Vis Sci. 1981;20:537–547. [PubMed]
19. Jester JV, Nicolaides N, Kiss-Palvolgyi I, Smith RE. Meibomian gland dysfunction, II: the role of keratinization in a rabbit model of MGD. Invest Ophthalmol Vis Sci. 1989;30:936–945. [PubMed]
20. Gutgesell VJ, Stern GA, Hood CI. Histopathology of meibomian gland dysfunction. Am J Ophthalmol. 1982;94:383–387. [PubMed]
21. Jester JV, Nicolaides N, Smith RE. Meibomian gland dysfunction. I. Keratin protein expression in normal human and rabbit meibomian glands. Invest Ophthalmol Vis Sci. 1989;30:927–935. [PubMed]
22. Obata H. Anatomy and histopathology of human meibomian gland. Cornea. 2002;21:S70–S74. [PubMed]
23. Obata H, Horiuchi H, Miyata K, Tsuru T, Machinami R. Histopathological study of the meibomian glands in 72 autopsy cases (in Japanese). Nippon Ganka Gakkai Zasshi. 1994;98:765–771. [PubMed]
24. Dougherty JM, McCulley JP. Comparative bacteriology of chronic blepharitis. Br J Ophthalmol. 1984;68:524–528. [PMC free article] [PubMed]
25. Dougherty JM, McCulley JP. Bacterial lipases and chronic blepharitis. Invest Ophthalmol Vis Sci. 1986;27:486–491. [PubMed]
26. McCulley JP, Dougherty JM. Bacterial aspects of chronic blepharitis. Trans Ophthalmol Soc UK. 1986;105:314–318. [PubMed]
27. Messmer EM, Torres SE, Mackert MI, Zapp DM, Kampik A. In vivo confocal microscopy in blepharitis (in German). Klin Monatsbl Augenheilkd. 2005;222:894–900. [PubMed]
28. Matsumoto Y, Shigeno Y, Sato E, et al. The evaluation of the treatment response in obstructive meibomian gland disease by in vivo laser confocal microscopy. Graefe's Arch Clin Exp Ophthalmol. 2009;247:821–829. [PubMed]
29. Hom MM, Martinson JR, Knapp LL, Paugh JR. Prevalence of meibomian gland dysfunction. Optom Vis Sci. 1990;67:710–712. [PubMed]
30. Hykin PG, Bron AJ. Age-related morphological changes in lid margin and meibomian gland anatomy. Cornea. 1992;11:334–342. [PubMed]
31. Lin PY, Tsai SY, Cheng CY, Liu JH, Chou P, Hsu WM. Prevalence of dry eye among an elderly Chinese population in Taiwan: the Shihpai Eye Study. Ophthalmology. 2003;110:1096–1101. [PubMed]
32. Lin PY, Cheng CY, Hsu WM, et al. Association between symptoms and signs of dry eye among an elderly Chinese population in Taiwan: the Shihpai Eye Study. Invest Ophthalmol Vis Sci. 2005;46:1593–1598. [PubMed]
33. Sullivan BD, Evans JE, Dana MR, Sullivan DA. Influence of aging on the polar and neutral lipid profiles in human meibomian gland secretions. Arch Ophthalmol. 2006;124:1286–1292. [PubMed]
34. Den S, Shimizu K, Ikeda T, Tsubota K, Shimmura S, Shimazaki J. Association between meibomian gland changes and aging, sex, or tear function. Cornea. 2006;25:651–655. [PubMed]
35. Sullivan DA, Sullivan BD, Ullman MD, et al. Androgen influence on the meibomian gland. Invest Ophthalmol Vis Sci. 2000;41:3732–3742. [PubMed]
36. Krenzer KL, Dana MR, Ullman MD, et al. Effect of androgen deficiency on the human meibomian gland and ocular surface. J Clin Endocr Metab. 2000;85:4874–4882. [PubMed]
37. Sullivan BD, Evans JE, Krenzer KL, Dana MR, Sullivan DA. Impact of antiandrogen treatment on the fatty acid profile of neutral lipids in human meibomian gland secretions. J Clin Endocr Metab. 2000;85:4866–4873. [PubMed]
38. Cermak JM, Krenzer KL, Sullivan RM, Dana MR, Sullivan DA. Is complete androgen insensitivity syndrome associated with alterations in the meibomian gland and ocular surface? Cornea. 2003;22:516–521. [PubMed]
39. Sullivan BD, Evans JE, Cermak JM, Krenzer KL, Dana MR, Sullivan DA. Complete androgen insensitivity syndrome: effect on human meibomian gland secretions. Arch Ophthalmol. 2002;120:1689–1699. [PubMed]
40. Sullivan DA, Sullivan BD, Evans JE, et al. Androgen deficiency, meibomian gland dysfunction, and evaporative dry eye. Ann N Y Acad Sci. 2002;966:211–222. [PubMed]
41. Tamer C, Oksuz H, Sogut S. Androgen status of the nonautoimmune dry eye subtypes. Ophthalmic Res. 2006;38:280–286. [PubMed]
42. Schaumberg DA, Dana R, Buring JE, Sullivan DA. Prevalence and risk factors for dry eye disease among US men. Arch Ophthalmol. 2009;127:763–768. [PMC free article] [PubMed]
43. Caffery BE, Josephson JE. Ocular side effects of isotretinoin therapy. J Am Optom Assoc. 1988;59:221–224. [PubMed]
44. Mathers WD, Shields WJ, Sachdev MS, Petroll WM, Jester JV. Meibomian gland morphology and tear osmolarity: changes with Accutane therapy. Cornea. 1991;10:286–290. [PubMed]
45. Denis P, Nordmann JP, Saiag P, Liotet S, Laroche L, Saraux H. Chronic blepharoconjunctivitis during a treatment with acitretin (Soriatane) (in French). J Fr Ophtalmol. 1993;16:191–194. [PubMed]
46. Sullivan DA, Yamagami H, Liu M, et al. Sex steroids, the meibomian gland and evaporative dry eye. Adv Exp Med Biol. 2002;506:389–400. [PubMed]
47. Schaumberg DA, Buring JE, Sullivan DA, Dana MR. Hormone replacement therapy and dry eye syndrome. JAMA. 2001;286:2114–2119. [PubMed]
48. Suzuki T, Schirra F, Richards SM, Jensen RV, Sullivan DA. Estrogen and progesterone control of gene expression in the mouse meibomian gland. Invest Ophthalmol Vis Sci. 2008;49:1797–1808. [PubMed]
49. Sullivan DA, Jensen RV, Suzuki T, Richards SM. Do sex steroids exert sex-specific and/or opposite effects on gene expression in lacrimal and meibomian glands? Mol Vis. 2009;10:15:1553–1572. [PMC free article] [PubMed]
51. Shi L, Reid LH, Jones WD, et al. The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotech. 2006;24:1151–1161. [PMC free article] [PubMed]
52. Ashburner M, Ball CA, Blake JA, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–29. [PMC free article] [PubMed]
53. Turner HC, Budak MT, Akinci MA, Wolosin JM. Comparative analysis of human conjunctival and corneal epithelial gene expression with oligonucleotide microarrays. Invest Ophthalmol Vis Sci. 2007;48:2050–2061. [PMC free article] [PubMed]
54. Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, Conklin BR. MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol. 2003;4:R7. [PMC free article] [PubMed]
55. Hohl D, de Viragh PA, Amiguet-Barras F, Gibbs S, Backendorf C, Huber M. The small proline-rich proteins constitute a multigene family of differentially regulated cornified cell envelope precursor proteins. J Invest Dermatol. 1995;104:902–909. [PubMed]
56. Tesfaigzi J, Carlson DM. Expression, regulation, and function of the SPR family of proteins: a review. Cell Biochem Biophys. 1999;30:243–265. [PubMed]
57. Fischer DF, Sark MW, Lehtola MM, Gibbs S, van de Putte P, Backendorf C. Structure and evolution of the human SPRR3 gene: implications for function and regulation. Genomics. 1999;55:88–99. [PubMed]
58. Song HJ, Poy G, Darwiche N, et al. Mouse Sprr2 genes: a clustered family of genes showing differential expression in epithelial tissues. Genomics. 1999;55:28–42. [PubMed]
59. Candi E, Schmidt R, Melino G. The cornified envelope: a model of cell death in the skin. Nat Rev. 2005;6:328–340. [PubMed]
60. Mischke D, Korge BP, Marenholz I, Volz A, Ziegler A. Genes encoding structural proteins of epidermal cornification and S100 calcium-binding proteins form a gene complex (“epidermal differentiation complex”) on human chromosome 1q21. J Invest Dermatol. 1996;106:989–992. [PubMed]
61. Kawasaki S, Kawamoto S, Yokoi N, et al. Up-regulated gene expression in the conjunctival epithelium of patients with Sjögren's syndrome. Exp Eye Res. 2003;77:17–26. [PubMed]
62. Li S, Gallup M, Chen YT, McNamara NA. Molecular mechanism of proinflammatory cytokine-mediated squamous metaplasia in human corneal epithelial cells. Invest Ophthalmol Vis Sci. 2010;51(5):2466–2475. [PMC free article] [PubMed]
63. Iizuka H, Takahashi H, Honma M, Ishida-Yamamoto A. Unique keratinization process in psoriasis: late differentiation markers are abolished because of the premature cell death. J Dermatol. 2004;31(4):271–276. [PubMed]
64. Ishida-Yamamoto A, Iizuka H, Manabe M, et al. Altered distribution of keratinization markers in epidermolytic hyperkeratosis. Arch Dermatol Res. 1995;287(8):705–711. [PubMed]
65. Nukui T, Ehama R, Sakaguchi M, et al. S100A8/A9, a key mediator for positive feedback growth stimulation of normal human keratinocytes. J Cell Biochem. 2008;104:453–464. [PubMed]
66. Ehrchen JM, Sunderkötter C, Foell D, Vogl T, Roth J. The endogenous Toll-like receptor 4 agonist S100A8/S100A9 (calprotectin) as innate amplifier of infection, autoimmunity, and cancer. J Leukoc Biol. 2009;86:557–566. [PubMed]
67. Foell D, Wittkowski H, Vogl T, Roth J. S100 proteins expressed in phagocytes: a novel group of damage-associated molecular pattern molecules. J Leukocyte Biol. 2007;81:28–37. [PubMed]
68. Zhou L, Beuerman RW, Ang LP, et al. Elevation of human alpha-defensins and S100 calcium-binding proteins A8 and A9 in tear fluid of patients with pterygium. Invest Ophthalmol Vis Sci. 2009;50(5):2077–2086. [PubMed]
69. Riau AK, Wong TT, Beuerman RW, Tong L. Calcium-binding S100 protein expression in pterygium. Mol Vis. 2009;15:335–342. [PMC free article] [PubMed]
70. Perera C, McNeil HP, Geczy CL. S100 calgranulins in inflammatory arthritis. Immunol Cell Biol. 2010;88(1):41–49. [PubMed]
71. Nacken W, Roth J, Sorg C, Kerkhoff C. S100A9/S100A8: Myeloid representatives of the S100 protein family as prominent players in innate immunity. Microsc Res Tech. 2003;60:569–580. [PubMed]
72. Kunz M, Roth J, Sorg C, Kolde G. Epidermal expression of the calcium binding surface antigen 27E10 in inflammatory skin diseases. Arch Dermatol Res. 1992;284:386–390. [PubMed]
73. Gabrielsen TO, Dale I, Brandtzaeg P, et al. Epidermal and dermal distribution of a myelomonocytic antigen (L1) shared by epithelial cells in various inflammatory skin diseases. J Am Acad Dermatol. 1986;15:173–179. [PubMed]
74. van Rheenen PF, Van de Vijver E, Fidler V. Faecal calprotectin for screening of patients with suspected inflammatory bowel disease: diagnostic meta-analysis. BMJ. 2010;15;341:c3369. [PMC free article] [PubMed]
75. Thuijls G, Derikx JP, van Wijck K, et al. Non-invasive markers for early diagnosis and determination of the severity of necrotizing enterocolitis. Ann Surg. 2010;251(6):1174–1180. [PubMed]
76. Kostakis ID, Cholidou KG, Kallianidis K, Perrea D, Antsaklis A. The role of calprotectin in obstetrics and gynecology. Eur J Obstet Gynecol Reprod Biol. 2010;151(1):3–9. [PubMed]
77. Lim SY, Raftery M, Cai H, et al. S-nitrosylated S100A8: novel anti-inflammatory properties. J Immunol. 2008;181(8):5627–5636. [PubMed]
78. Hsu K, Champaiboon C, Guenther BD, Sorenson BS, et al. Anti-infective protective properties of s100 calgranulins. Antiinflamm Antiallergy Agents Med Chem. 2009;8(4):290–305. [PMC free article] [PubMed]
79. Meyaard L. The inhibitory collagen receptor LAIR-1 (CD305). J Leukoc Biol. 2008;83(4):799–803. [PubMed]
80. Omiya R, Tsushima F, Narazaki H, et al. Leucocyte-associated immunoglobulin-like receptor-1 is an inhibitory regulator of contact hypersensitivity. Immunology. 2009;128(4):543–555. [PubMed]
81. Straatsma BR. Cystic degeneration of the meibomian glands. Arch Ophthalmol. 1959;61:918–927. [PubMed]
82. Champaiboon C, Sappington KJ, Guenther BD, Ross KF, Herzberg MC. Calprotectin S100A9 calcium-binding loops I and II are essential for keratinocyte resistance to bacterial invasion. J Biol Chem. 2009;284(11):7078–7090. [PubMed]
83. Contzler R, Favre B, Huber M, Hohl D. Cornulin, a new member of the “fused gene” family, is expressed during epidermal differentiation. J Invest Dermatol. 2005;124(5):990–997. [PubMed]
84. [Accessed August 11, 2010]. http://www.genecards.org/
86. Fuchs E, Horsley V. More than one way to skin. Genes Dev. 2008;22(8):976–985. [PubMed]
87. Knaggs HE, Holland DB, Morris C, Wood EJ, Cunliffe WJ. Quantification of cellular proliferation in acne using the monoclonal antibodyKi-67. J Invest Dermatol. 1994;102(1):89–92. [PubMed]
88. Singh RS, Grayson W, Redston M, et al. Site and tumor type predicts DNA mismatch repair status in cutaneous sebaceous neoplasia. Am J Surg Pathol. 2008;32(6):936–942. [PubMed]
89. Popnikolov NK, Gatalica Z, Colome-Grimmer MI, Sánchez RL. Loss of mismatch repair proteins in sebaceous gland tumors. J Cutan Pathol. 2003;30(3):178–184. [PubMed]
90. Yeh CK, Dodds MW, Zuo P, Johnson DA. A population-based study of salivary lysozyme concentrations and candidal counts. Arch Oral Biol. 1997;42(1):25–31. [PubMed]
91. Murias AR, Weruaga E, Recio JS, Airado C, Diaz D, Alonso JR. Distribution of neurocalcin-containing neurons reveals sexual dimorphism in the mouse olfactory bulb. Chem Senses. 2007;32(7):673–680. [PubMed]
92. Zhang S, Liu X, Brickman WJ, et al. Association of plasma leptin concentrations with adiposity measurements in rural Chinese adolescents. J Clin Endocrinol Metab. 2009;94(9):3497–3504. [PubMed]
93. Suzuki T, Sullivan BD, Liu M, et al. Estrogen and progesterone effects on the morphology of the mouse meibomian gland. Adv Exp Med Biol. 2002;506:483–488. [PubMed]
94. Chew CKS, Hykin PG, Janswijer C, Dikstein S, Tiffany JM, Bron AJ. The casual level of meibomian lipids in humans. Curr Eye Res. 1993;12:255–259. [PubMed]
95. Sullivan BD, Cermak JM, Sullivan RM, et al. Correlations between nutrient intake and the polar lipid profiles of meibomian gland secretions in women with Sjögren's syndrome. Adv Exp Med Biol. 2002;506:441–448. [PubMed]
96. Schaumberg DA, Sullivan DA, Buring JE, Dana MR. Prevalence of dry eye syndrome among US women. Am J Ophthalmol. 2003;136:318–326. [PubMed]
97. Chia EM, Mitchell P, Rochtchina E, Lee AJ, Maroun R, Wang JJ. Prevalence and associations of dry eye syndrome in an older population: the Blue Mountains Eye Study. Clin Experiment Ophthalmol. 2003;31:229–232. [PubMed]
98. Richards SM, Yamagami H, Schirra F, Suzuki T, Jensen RV, Sullivan DA. Sex-related effect on gene expression in the mouse meibomian gland. Curr Eye Res. 2006;31:119–128. [PubMed]
99. Cunningham BA, Moncur JT, Huntington JT, Kinlaw WB. “Spot 14” protein: a metabolic integrator in normal and neoplastic cells. Thyroid. 1998;8(9):815–825. [PubMed]
100. Nichols KK, Ham BM, Nichols JJ, Ziegler C, Green-Church KB. Identification of fatty acids and fatty acid amides in human meibomian gland secretions. Invest Ophthalmol Vis Sci. 2007;48(1):34–39. [PubMed]
101. Butovich IA, Uchiyama E, McCulley JP. Lipids of human meibum: mass-spectrometric analysis and structural elucidation. J Lipid Res. 2007;48(10):2220–2235. [PubMed]
102. Shine WE, McCulley JP. Association of meibum oleic acid with meibomian seborrhea. Cornea. 2000;19(1):72–74. [PubMed]
103. Schirra F, Suzuki T, Richards SM, et al. Androgen control of gene expression in the mouse meibomian gland. Invest Ophthalmol Vis Sci. 2005;46:3666–3675. [PubMed]
104. Miyazaki M, Man WC, Ntambi JM. Targeted disruption of stearoyl-CoA desaturase1 gene in mice causes atrophy of sebaceous and meibomian glands and depletion of wax esters in the eyelid. J Nutr. 2001;131(9):2260–2268. [PubMed]
105. Tsuchida S, Bonkobara M, McMillan JR, et al. Ariizumi K. Characterization of Kdap, a protein secreted by keratinocytes. J Invest Dermatol. 2004;122(5):1225–1234. [PubMed]
106. Takahashi M, Lin YM, Nakamura Y, Furukawa Y. Isolation and characterization of a novel gene CLUAP1 whose expression is frequently upregulated in colon cancer. Oncogene. 2004;23(57):9289–9294. [PubMed]
107. Hayashi Y, Liu CY, Jester JJ, et al. Excess biglycan causes eyelid malformation by perturbing muscle development and TGF-alpha signaling. Dev Biol. 2005;277(1):222–234. [PMC free article] [PubMed]
108. Remington SG, Nelson JD. Secretoglobins: sexually dimorphic expression of androgen-binding protein mRNA in mouse lacrimal glands. Invest Ophthalmol Vis Sci. 2005;46(1):31–38. [PubMed]
109. Luu-The V, Labrie F. The intracrine sex steroid biosynthesis pathway. Prog Brain Res. 2010;181:177–192. [PubMed]
110. Stoeckelhuber M, Messmer EM, Schmidt C, Xiao F, Schubert C, Klug J. Immunohistochemical analysis of secretoglobin SCGB 2A1 expression in human ocular glands and tissues. Histochem Cell Biol. 2006;126(1):103–109. [PubMed]
111. Zhou L, Beuerman RW, Chan CM, et al. Identification of tear fluid biomarkers in dry eye syndrome using iTRAQ quantitative proteomics. J Proteome Res. 2009;8(11):4889–4905. [PubMed]
112. Nichols JJ, Green-Church KB. Mass spectrometry-based proteomic analyses in contact lens-related dry eye. Cornea. 2009;28(10):1109–1117. [PubMed]
113. Naoumova RP, Dunn S, Rallidis L, et al. Prolonged inhibition of cholesterol synthesis explains the efficacy of atorvastatin. J Lipid Res. 1997;38:1496–1500. [PubMed]
114. Thody AJ, Shuster S. Control and function of sebaceous glands. Physiol Rev. 1989;69:383–416. [PubMed]
Articles from Investigative Ophthalmology & Visual Science are provided here courtesy of
Association for Research in Vision and Ophthalmology