near - infrared ( nir ) radiation can penetrate the skin and the sclera of the eye .
the high permeability of nir radiation also allows it to affect tissues deeper within the eye , such as muscles , the lens , and retina .
nir radiation can induce various biological effects,19 and intensive or long - term exposure to nir radiation is a factor in premature aging . despite the wide prevalence of a variety of ultraviolet ( uv ) blocking materials , such as sunblock
, sunglasses , glasses , films , and umbrellas , that are useful in protecting our tissue against uv exposure , nir can not be blocked sufficiently.8 consequently , in the absence of suitable protection , nir radiation can induce various kinds of tissue damage and diseases , such as cataracts and photoaging.7,8 the human cornea plays a critical role in refracting light onto the retina and also protects the eye against external agents .
since the epithelial layer of the cornea provides the first line of defense against environmental insults , the structural integrity of this layer is a key component of corneal function.10 although uv - induced corneal damage has been described in many previous studies,1115 the effects induced by nir radiation on the cornea have not been thoroughly investigated .
we hypothesized that nir irradiation simulating solar nir radiation that reaches human tissues can induce changes in gene expression . to test this hypothesis , a three - dimensional reconstructed human corneal epithelial model with multilayered ,
corneal epithelium - like structure was used to simulate the human eye , and we evaluated dna microarray and real - time polymerase chain reaction ( pcr ) analysis results from normal human corneal epithelial cells exposed to water - filtered broad - spectrum nir irradiation to simulate solar nir radiation that reaches the eye .
nir irradiation was performed with a broadband nir source ( titan ; cutera , brisbane , ca , usa ) .
the nir device emits an nir spectrum between 1,100 nm and 1,800 nm , with water filtering to remove wavelengths between 1,400 nm and 1,500 nm , and simulates solar nir radiation that reaches the skin of humans on the earth s surface . to avoid thermal effects ,
the sapphire contact cooling tip was set to a fixed temperature of 20c . in our previous in vitro study , ten rounds at 10 j / cm using continuous energy single irradiation pulses of 4.3 seconds achieved drastic reduction in cell count .
the three - dimensional reconstructed human corneal epithelial model ( labcyte cornea - model ) prepared from enzymatically digested normal human corneal epithelial tissues was purchased from japan tissue engineering corporation , aichi , japan as an in vitro model of corneal tissue.16 cells were cultured in media ( assay medium ; japan tissue engineering corporation ) , which was changed every 2 days until the cultures reached subconfluence.16 the subconfluent corneal cells were then subcultured with trypsin and seeded on a cell culture insert containing a microporous membrane with a 0.4 m pore size .
corneal cells were cultivated between the air and liquid interface to form a multilayered , corneal epithelium - like structure .
total rna from the corneal epithelial cells was extracted in qiazol reagent ( qiagen nv , venlo , the netherlands ) and spin - column purified using a rneasy mini spin column ( qiagen nv ) .
the rna was extracted from control and irradiated groups , which consisted of six whole individual three - dimensional reconstructed human corneal epithelial models .
the quantity and quality of the rna samples were then measured using uv absorbance ( nanodrop technologies , wilmington , de , usa ) , and the quality was also assessed using an agilent 2100 bioanalyzer series ii ( agilent technologies , santa clara , ca , usa ) .
cy3-labeled crna samples were synthesized using the low input quick amp labeling kit ( agilent technologies ) according to the manufacturer s instructions . for each time point , 50 ng of total rna was used to generate first - strand cdna . after the denaturation ( 10 minutes at 65c ) and crna synthesis ( 2 hours at 40c ) steps , the reactions were incubated at 70c for 15 minutes to inactivate the affinity - script enzyme ( agilent technologies ) . for labeling reactions ,
crna samples were each mixed with 6 l of transcription master mix cocktail containing cy3-ctp and incubated at 40c for 2 hours .
labeled crna was quantified on a nanodrop nd-1000 ( nanodrop technologies ) , and quality was also checked by agilent 2100 bioanalyzer series ii .
slides were scanned using the microarray scanner ( agilent technologies ) , and images were processed by the feature extraction software ( agilent technologies ) with background correction .
data were imported into genespring , and normalization was performed with a 75 percentile shift .
after normalization and filtering of the raw data output files , log - fold changes in gene expression in the nir - irradiated group compared with the control group were analyzed using the genespring gx 11.0 software ( agilent technologies ) to determine the upregulated and downregulated genes , using a twofold change as the cutoff .
upregulated and downregulated genes were further clustered into functional gene groups using the ontology analysis function with the genespring gx 11.0 software .
microarray data were deposited into the minimum information about a microarray experiment ( miame)-compliant gene expression omnibus ( geo ) database ( accession number : gse76500 , http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=gse76500 ) .
quantitative real - time pcr analysis was achieved using the dna engine opticon 2 system ( bio - rad laboratories inc .
all six whole individual three - dimensional reconstructed human corneal epithelial tissue culture models were used for quantitative real - time pcr analysis .
the primer sequences of epidermal growth factor receptor ( egfr ) , tubulin tyrosine ligase - like family member 5 ( ttll5 ) , and pq loop repeat containing 2 ( pqlc2 ) are given in table 1 .
since pqlc was ubiqitously and highly expressed in all samples , pqlc2 was used for normalization .
the variant - specific primer sequences for the four egfr variants are given in table 2 .
cdna templates were prepared for each rna sample by reverse transcription reaction using superscript iii ( thermo fisher scientific , waltham , ma , usa ) .
each quantitative pcr ( 25 l ) contained 1 l of cdna template , 1.5 l of forward and reverse primer ( 10 m each ) , 7.5 l of 2 master mix , and 3.5 l of nuclease - free water .
the thermal cycling conditions included one cycle at 95c for 5 minutes and 40 cycles at 95c for 10 seconds , 60c for 30 seconds , and 68c for 30 seconds .
all statistical data are presented as the mean  sem , and a p - value of < 0.05 was considered to be statistically significant .
nir irradiation was performed with a broadband nir source ( titan ; cutera , brisbane , ca , usa ) .
the nir device emits an nir spectrum between 1,100 nm and 1,800 nm , with water filtering to remove wavelengths between 1,400 nm and 1,500 nm , and simulates solar nir radiation that reaches the skin of humans on the earth s surface . to avoid thermal effects ,
the sapphire contact cooling tip was set to a fixed temperature of 20c . in our previous in vitro study , ten rounds at 10 j / cm using continuous energy single irradiation pulses of 4.3 seconds achieved drastic reduction in cell count .
the three - dimensional reconstructed human corneal epithelial model ( labcyte cornea - model ) prepared from enzymatically digested normal human corneal epithelial tissues was purchased from japan tissue engineering corporation , aichi , japan as an in vitro model of corneal tissue.16 cells were cultured in media ( assay medium ; japan tissue engineering corporation ) , which was changed every 2 days until the cultures reached subconfluence.16 the subconfluent corneal cells were then subcultured with trypsin and seeded on a cell culture insert containing a microporous membrane with a 0.4 m pore size .
corneal cells were cultivated between the air and liquid interface to form a multilayered , corneal epithelium - like structure .
total rna from the corneal epithelial cells was extracted in qiazol reagent ( qiagen nv , venlo , the netherlands ) and spin - column purified using a rneasy mini spin column ( qiagen nv ) . the rna was extracted from control and irradiated groups , which consisted of six whole individual three - dimensional reconstructed human corneal epithelial models .
the quantity and quality of the rna samples were then measured using uv absorbance ( nanodrop technologies , wilmington , de , usa ) , and the quality was also assessed using an agilent 2100 bioanalyzer series ii ( agilent technologies , santa clara , ca , usa ) .
cy3-labeled crna samples were synthesized using the low input quick amp labeling kit ( agilent technologies ) according to the manufacturer s instructions . for each time point , 50 ng of total rna was used to generate first - strand cdna . after the denaturation ( 10 minutes at 65c ) and crna synthesis ( 2 hours at 40c ) steps , the reactions were incubated at 70c for 15 minutes to inactivate the affinity - script enzyme ( agilent technologies ) . for labeling reactions ,
crna samples were each mixed with 6 l of transcription master mix cocktail containing cy3-ctp and incubated at 40c for 2 hours .
labeled crna was quantified on a nanodrop nd-1000 ( nanodrop technologies ) , and quality was also checked by agilent 2100 bioanalyzer series ii .
slides were scanned using the microarray scanner ( agilent technologies ) , and images were processed by the feature extraction software ( agilent technologies ) with background correction .
data were imported into genespring , and normalization was performed with a 75 percentile shift .
after normalization and filtering of the raw data output files , log - fold changes in gene expression in the nir - irradiated group compared with the control group were analyzed using the genespring gx 11.0 software ( agilent technologies ) to determine the upregulated and downregulated genes , using a twofold change as the cutoff .
upregulated and downregulated genes were further clustered into functional gene groups using the ontology analysis function with the genespring gx 11.0 software .
microarray data were deposited into the minimum information about a microarray experiment ( miame)-compliant gene expression omnibus ( geo ) database ( accession number : gse76500 , http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=gse76500 ) .
quantitative real - time pcr analysis was achieved using the dna engine opticon 2 system ( bio - rad laboratories inc .
all six whole individual three - dimensional reconstructed human corneal epithelial tissue culture models were used for quantitative real - time pcr analysis .
the primer sequences of epidermal growth factor receptor ( egfr ) , tubulin tyrosine ligase - like family member 5 ( ttll5 ) , and pq loop repeat containing 2 ( pqlc2 ) are given in table 1 .
since pqlc was ubiqitously and highly expressed in all samples , pqlc2 was used for normalization .
the variant - specific primer sequences for the four egfr variants are given in table 2 .
cdna templates were prepared for each rna sample by reverse transcription reaction using superscript iii ( thermo fisher scientific , waltham , ma , usa ) .
each quantitative pcr ( 25 l ) contained 1 l of cdna template , 1.5 l of forward and reverse primer ( 10 m each ) , 7.5 l of 2 master mix , and 3.5 l of nuclease - free water .
the thermal cycling conditions included one cycle at 95c for 5 minutes and 40 cycles at 95c for 10 seconds , 60c for 30 seconds , and 68c for 30 seconds .
all statistical data are presented as the mean  sem , and a p - value of < 0.05 was considered to be statistically significant .
to elucidate changes in genetic expression induced by exposure to nir irradiation that simulates solar nir radiation that reaches human eye tissues , dna microarray analysis was carried out using rna isolated from control and nir - irradiated corneal epithelial cells . using ~6.310 probes , we identified 52 upregulated genes and 98 downregulated genes in the nir - irradiated group compared with the control group .
gene ontology analysis of these 150 genes identified 41 upregulated genes that encode proteins with receptor - binding properties ( table 3 ) .
quantitative real - time pcr analysis was then performed for the two genes showing the highest degree of upregulation in the dna microarray analysis : egfr and ttll5 .
egfr showed a high degree of upregulation ( 19.4-fold ) relative to control cells , although statistically significant upregulation of egfr and ttll5 was not observed ( p=0.5127 and p=0.8273 , respectively ) ( figure 1 ) using the primer sequences listed in table 1 .
meanwhile , when the expression of egfr variants was assessed in a dna microarray analysis , expression increased by 1.0-fold and 1.3-fold for egfr variants 1 and 3 , respectively . upon inclusion of variant - specific primer sequences for egfr ( table 2 ) , quantitative real - time pcr analysis revealed that statistically significant upregulation of variants 1 and 3 was observed ( p<0.05 ) , whereas there was no statistically significant upregulation of variants 2 and 4 ( p=0.1266 and p=0.5127 , respectively ) ( figure 2 ) .
to the best of our knowledge , this is the first study that examines changes in genetic expression in three - dimensional reconstructed corneal epithelial tissue after nir irradiation to simulate solar nir radiation that reaches human eye tissues .
the biological effects of sun and uv exposure have been extensively investigated . however , over half of the solar energy is in the nir range , and most of the uv - blocking materials can not block nir radiation .
moreover , humans are continuously exposed to artificial nir radiation from sources , such as electrical appliances.17,18 despite the prevalence of nir radiation , the need for protective measures against this type of radiation has received less attention than that of uv . because solar nir radiation is filtered by atmospheric water,19,20 water - filtering in experimental settings is indispensable to provide an accurate simulation of solar nir radiation that reaches human tissues.5 cooling is also recommended when investigating the effects of nir radiation , as nir radiation can increase surface temperatures and induce thermal effects.5 indeed , without water filter or contact cooling , nir irradiation immediately increases the temperature of the superficial layer of culture fluid in a laboratory dish or on skin , as nir radiation is predominantly absorbed by hydrogen bond - containing molecules , such as water and hemoglobin . meanwhile , as the nir radiation penetrates deeper into the target , its energy can dissipate to the point where insufficient amounts of energy reach target cells in the base of culture dishes or deeper tissues in skin samples .
thus , to accurately investigate the biological effects of solar nir radiation that reaches human tissues , a water - filter that excludes wavelengths between 1,400 nm and 1,500 nm and a cooling system were used in this study.5 a three - dimensional reconstructed human corneal epithelial model was used as an in vitro model of corneal tissue .
the corneal epithelial model was proven to have a fully differentiated corneal epidermal tissue including all major corneal epithelial cell layers , such as the basal cell layer , the wing layer , and the superficial cell layer .
cytokeratin-3 is a specific marker of corneal epithelium , and it was expressed in all layers of the corneal epithelial model.16 muc-1 and muc-16 , which compose the transmembrane glycoproteins in the surface of the corneal epithelium , were well expressed in the superficial layer of the corneal epithelial model.16 cells in the corneal epithelium are connected by desmosomes , tight junctions , and adherence junctions.16 claudin-1 , which is a marker of a tight junction , desmograin-3 , which is a marker of desmosome , and e - cadherin , which is a marker of an adherence junction , are localized in the interface between cells at all cell layers , including the superficial layers.16 lamin is an important constituent in the basement membrane at the basal corneal epithelium junction , which is expressed continuously in basal cells of the corneal epithelial model basal layer.16 the basement membrane was smooth like that of the human eye.16 therefore , the corneal epithelial model reproduced many of the characteristics of the native human corneal tissue , and it provides a morphologically relevant means to assess nir effects as an alternative to animal testing .
the corneal epithelium , like other epithelial barriers in the human body , is continuously subjected to physical , chemical , and biological insults , which can produce a wound and/or loss of barrier functions .
proper healing of corneal wounds is vital for maintaining a clear , healthy cornea and preserving vision .
corneal epithelium responds rapidly to injury , wherein wound healing occurs by cells migrating as a sheet to cover the defect and reestablish barrier functions.21 in wounded cornea , the epithelium plays a central role , not only as a key cell type in corneal repair , but also as a growth factor source.22 as in other tissues , several growth factors are suggested to play a role in regulating corneal epithelial function and wound healing.21,23 following nir irradiation that simulates solar radiation , microarray analysis of corneal epithelial cell dna in this study showed that egfr was highly expressed , and statistically significant upregulation of egfr variants 1 and 3 was observed ( p<0.05 ) by quantitative real - time pcr .
egfr is a prototypic tyrosine kinase receptor that is part of a larger family of erbb receptors .
erbb2 , erbb3 , and erbb4 share with egfr characteristics that include extracellular ligand - binding sites , intracellular kinase domains , and tyrosine - rich regions .
egfr activated in response to injury serves as a powerful mediator of corneal epithelial wound healing10 by coordinating multiple extracellular signals generated in response to cell injury.2427 in turn , the levels of egfr and tyrosine kinase activity are major determinant factors for epithelial cell function in tissues and organs.22 echoing this concept is the finding that in some cancer patients , treatment with egfr - specific monoclonal antibodies cetuximab and gefitinib ( an egfr kinase inhibitor ) resulted in ocular abnormalities , including diffuse punctate keratitis and corneal erosion.2830 thus , maintaining a proper level of egfr signaling is critical for corneal homeostasis.22 ttll5 gene has 32 exons with high expression in heart and skeletal muscle and lower expression in many other tissues , including the eye and brain.31,32 it encodes a 1,281-amino acid protein that is localized to the cytoplasm and nucleus.31 ttll5 is thought to be a key initiator of polyglutamylation in -tubulin33 and has been previously reported to be essential for the correct function of sperm flagella.34 ttll5 plays a role in polyglutamylation of primary cilia in vitro.35 notably , genes involved in the polyglutamylation and deglutamylation of -tubulin have been associated with photoreceptor degeneration in mice.35 in this study , ttll5 showed the high degree of upregulation in the dna microarray analysis , whereas statistically significant upregulation of ttll5 was not observed in quantitative real - time pcr analysis .
further studies are needed to determine if a higher power output or increased treatment frequency may promote significant upregulation of ttll5 .
the findings here that egfr expression is significantly upregulated following exposure of corneal epithelial cells to nir irradiation support that solar nir radiation that reaches human tissues can induce corneal cell injury .
although significant egfr expression upregulation occurred after just one treatment at the tested power output , further studies are needed to determine if a higher power output or increased treatment frequency may promote even larger changes in expression .
furthermore , this was a preliminary study based on a fairly small number of samples , so a large - scale study will be needed .
nir irradiation simulating solar nir radiation that reaches human tissues induced upregulated expression of egfr in human corneal cells .
nir radiation represents more than half of solar energy and can not be sufficiently blocked by eyewear . moreover ,
thus , the results of this study indicate that enhanced protection from and avoidance of both uv and nir radiation should be considered to prevent eye damage .
 background and objectivehumans are increasingly exposed to near - infrared ( nir ) radiation from both natural ( eg , solar ) and artificial ( eg , electrical appliances ) sources . 
 although the biological effects of sun and ultraviolet ( uv ) exposure have been extensively investigated , the biological effect of nir radiation is still unclear . 
 we previously reported that nir as well as uv induces photoaging and standard uv - blocking materials , such as sunglasses , do not sufficiently block nir . 
 the objective of this study was to investigate changes in gene expression in three - dimensional reconstructed corneal epithelial tissue culture exposed to broad - spectrum nir irradiation to simulate solar nir radiation that reaches human tissues.materials and methodsdna microarray and quantitative real - time polymerase chain reaction analysis were used to assess gene expression levels in a three - dimensional reconstructed corneal epithelial model composed of normal human corneal epithelial cells exposed to water - filtered broad - spectrum nir irradiation with a contact cooling ( 20c ) . 
 the water - filter allowed 1,0001,800 nm wavelengths and excluded 1,4001,500 nm wavelengths.resultsa dna microarray with > 62,000 different probes showed 25 and 150 genes that were up- or downregulated by at least fourfold and twofold , respectively , after nir irradiation . in particular , epidermal growth factor receptor ( egfr ) was upregulated by 19.4-fold relative to control cells . 
 quantitative real - time polymerase chain reaction analysis revealed that two variants of egfr in human corneal epithelial tissue were also significantly upregulated after five rounds of 10 j / cm2 irradiation ( p<0.05).conclusionwe found that nir irradiation induced the upregulated expression of egfr in human corneal cells . since over half of the solar energy reaching the earth is in the nir region , which can not be adequately blocked by eyewear and thus can induce eye damage with intensive or long - term exposure , protection from both uv and nir radiation may prevent changes in gene expression and in turn eye damage . 
Introduction
Materials and methods
NIR irradiation
Corneal epithelial model
RNA extraction
cRNA synthesis, labeling, and purification
Microarray analysis
Quantitative real-time PCR analysis
Statistical analysis
Results
Discussion
Conclusion