Next Article in Journal
Multimodality Imaging in Tumor Angiogenesis: Present Status and Perspectives
Next Article in Special Issue
Regenerative Therapies in Dry Eye Disease: From Growth Factors to Cell Therapy
Previous Article in Journal
REP1 Modulates Autophagy and Macropinocytosis to Enhance Cancer Cell Survival
Previous Article in Special Issue
Analysis of the Pathogenic Factors and Management of Dry Eye in Ocular Surface Disorders
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 18
Issue 9
10.3390/ijms18091862
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Human Meibum Age, Lipid–Lipid Interactions and Lipid Saturation in Meibum from Infants

by
Samiyyah Sledge
1,
Collin Henry
1,
Douglas Borchman
1,*,
Marta C. Yappert
2,
Rahul Bhola
1,3,
Aparna Ramasubramanian
1,
Ryan Blackburn
1,
Jonathan Austin
1,
Kayla Massey
1,
Shanzeh Sayied
1,
Aliza Williams
1,
Georgi Georgiev
4 and
Kenneth N. Schikler
5
1
Department of Ophthalmology and Visual Sciences, University of Louisville, Louisville, KY 40202, USA
2
Department of Chemistry, University of Louisville, Louisville, KY 40292, USA
3
Division of Ophthalmology, Children’s Hospital of Orange County, Orange, CA 92868, USA
4
Model Membranes Lab, Department of Biochemistry, Faculty of Biology, St. Kliment Ohridski University of Sofia, Sofia 1164, Bulgaria
5
Department of Pediatrics, University of Louisville, Louisville, KY 40202, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(9), 1862; https://doi.org/10.3390/ijms18091862
Submission received: 28 April 2017 / Revised: 25 July 2017 / Accepted: 17 August 2017 / Published: 28 August 2017
(This article belongs to the Special Issue Dry Eye and Ocular Surface Disorders)
Browse Figures
Versions Notes

Abstract

:
Tear stability decreases with increasing age and the same signs of instability are exacerbated with dry eye. Meibum lipid compositional changes with age provide insights into the biomolecules responsible for tear film instability. Meibum was collected from 69 normal donors ranging in age from 0.6 to 68 years of age. Infrared spectroscopy was used to measure meibum lipid phase transition parameters. Nuclear magnetic resonance spectroscopy was used to measure lipid saturation. Increasing human meibum lipid hydrocarbon chain unsaturation with age was related to a decrease in hydrocarbon chain order, cooperativity, and in the phase transition temperature. The change in these parameters was most dramatic between 1 and 20 years of age. Meibum was catalytically saturated to determine the effect of saturation on meibum lipid phase transition parameters. Hydrocarbon chain saturation was directly related to lipid order, phase transition temperature, cooperativity, changes in enthalpy and entropy, and could account for the changes in the lipid phase transition parameters observed with age. Unsaturation could contribute to decreased tear film stability with age.

1. Introduction

Tear lipids, mostly from the Meibomian gland and a minor amount from sebaceous glands [1], may be important for tear stability [2,3,4]. Changes in tear film lipid composition with age could give us insights into lipid compositional-functional relationships with dry eye. For instance, the signs of dry eye such as decreased breakup time and increased blink rate, are exacerbations of the same signs observed with aging [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15]. The spontaneous blink rate of adults is as much as 20 times per minute, much higher than that of infants which blink less than one time a minute [15]. The spontaneous blink rate is related to the tear break-up time. Tear break-up time is as high as 35 s in infants and decreases to 8–16 s in adults. Tear break-up time is even lower (5 s) in adults with Meibomian gland dysfunction [16,17,18,19,20,21,22].
This project is an extension of previous nuclear magnetic resonance (NMR) [23,24] and Fourier transform infrared (FTIR) [1,24,25,26,27,28] spectral studies relating age with meibum composition, structure, and function. Using deuterated chloroform as a solvent rather than deuterated cyclohexane used previously [23,24], in the current study, the double bond resonance assigned to cholesterol was resolved and quantified from the double bond resonance associated with hydrocarbons using a 700 MHz NMR spectrometer. The 700 MHz NMR spectrometer is more powerful than the 500 MHz NMR spectrometer used previously [23,24]. Furthermore, catalytic hydrogenation was used to examine the relationships between hydrocarbon chain order and the level of saturation. This was an improvement over the previous study where only native meibum was compared with meibum that was 100% saturated [25], a level that is not physiological. This study provides insights into how the increase in meibum lipid saturation can be related to the observed decrease in tear film stability with age.

2. Results

2.1. Donor Demographics

Meibum was obtained from donors that did not have observable signs of Meibomian gland blockage or any general complaint of dry eye. Meibum for the catalytic saturation study was pooled from six Caucasian females 23, 24, 28, 32, 66, and 68 years old; a 29-year-old Caucasian male; and a 32-year-old African American female. Meibum for the saturation-age NMR study were obtained from the following donors: CM0.6, CM01, CM01, BM04, CM01, BM01, CM01, CF1.4, BM1.6, CM1.8, HF2, CM02a, CF02, CM02, CM02, AM03, CM03, CM03, CM04, CM04, CM04, CF04, CM04, CF05, CM05, CM05, CF05, XM08, XM08, CF8.5, BM10, CF11, CF11, CF13, CM14, CM17, BF19, CF20, CM22, CF23, CF24, CF24, CM25, CF26, CF27, CF28, HF29, CM29, CM31, CM31b, B32F, CF36, HF37, BF38, BM39, CF43, CM55, CM59, CM62, CF66, and CF68. The sample key is: C = Caucasian; H = Hispanic; A = Asian; B = African American; X = unknown race; F = female; M = male; and numbers are age. Meibum was collected only once from each donor.

2.2. H-NMR Spectroscopy

Average proton NMR (1H-NMR) spectra of human meibum were typical of meibum collected from a 700 MHz spectrometer (Figure 1). Band assignments were made based on previous 1H and carbon 13 NMR studies [23,29]. The largest resonance in this region was observed at 5.32 ppm with a shoulder at 5.35 ppm assigned to protons of the cis =CH moieties from hydrocarbon chains and to the proton attached to carbon #6 of cholesterol esters, respectively. The resonance at 4.6 ppm is from cholesteryl esters and the resonance at 4.0 ppm is from wax esters (Figure 1). The resonances near 5.1 ppm have been assigned to squalene [29].
The total level of double bonds from cholesteryl and wax =CH resonances (5.32 and 5.35 ppm) relative to the sum of wax (4.0 ppm) and cholesteryl ester (4.6 ppm) resonances increased significantly (p = 0.03) from 1.0 ± 0.1 for infants to 1.4 ± 0.1 in children. Unsaturation of human meibum increased with age and the relative level of hydrocarbon cis =CH unsaturation of infants was significantly lower (p < 0.0001) compared with that of adults (Figure 2a). Samples above 20 years of age were grouped together based on the developmental Tanner stage and changes in blink rate, free fatty acids and meibum lipid phase transition parameters (see Discussion).

2.3. Infrared Spectroscopy

Infrared spectroscopy was used to study lipid–lipid interactions and composition. The CH2 stretching and bending bands are predominant in the infrared spectra of lipids due to the large number of CH2 groups in their hydrocarbon chains. The CH stretching region of meibum is composed of five major bands (Figure 3) [26]. Note the catalytically hydrogenated sample has no =CH stretching band (Figure 3b). In this study, we used the frequency of the symmetric CH2 stretching band near 2850 cm−1 (ṽsym) to estimate the trans to gauche rotamer content of the hydrocarbon chains. The ṽsym increased with an increase in temperature and the number of gauche rotamers, concurrent with a decrease in intensity (Figure 4) [26,30]. The peak height of the CH2 symmetric stretching band at 9 °C was approximately 0.23 absorbance units. The absolute intensity of the CH stretching region decreased by about 20% with an increase in temperature from 9 to 65 °C which was attributed partially to a 50% decrease in the CH2 symmetric stretching band [30]. A sigmoidal equation was used to fit and quantify the lipid phase transitions [26]. Two of the fitted parameters, the minimum and maximum ṽsym, correspond to the most ordered and disordered states of hydrocarbon chains, respectively. Another fitted parameter was the phase transition temperature, which is the temperature at which half of the lipid molecules undergo a phase change. The fourth fitted parameter was the relative cooperativity of the phase transition that describes how the order of a lipid influences that of neighboring lipids. Broad phase transitions have a relatively smaller absolute value of the cooperativity. Lipid phase transition parameters for a pool of human meibum used in the saturation study are listed in Table 1. Lipid order was measured close to the surface temperature of the human eye, 33.4 °C, by extrapolating the ṽsym at 33.4 °C from the fit of the phase transition and then converting ṽsym to the percentage of trans rotamers [26]. The lipid order measured in this study, (31 ± 2) % trans rotamers), reinforced the correlation between a decrease in lipid order with increasing age (Figure 2b, r = 0.963, p < 0.01).
Meibum lipid was catalytically saturated and the lipid phase transition parameters were measured. Lipid order at 33.4 °C increased significantly (p < 0.0001) from 39 ± 3% to 82 ± 1% between 0% and 25% saturation (Figure 5a). Above 25% saturation, lipid order reached a maximum. The lipid phase transition temperature for meibum lipids increased significantly (p < 0.01, r = 0.963) with saturation, as expected, from about 30 to 51 °C (Figure 5b). The change in enthalpy (∆H) (Figure 5c), and entropy (∆S) (Figure 5d), and cooperativity (Figure 5e) of the lipid phase transition increased substantially with saturation, p < 0.01, 0.05, and 0.1; and r = 0.626, 0.544, and 0.950, respectively.
Arrhenius plots used to calculate the ∆H and ∆S values from the lipid phase transitions were linear, with correlation coefficients greater than 0.998 (Figure 6). For comparison of phase transition parameters of catalytically saturated meibum with age related changes, we refitted the phase transition curves from previous studies and recalculated the percent trans rotamers because in previous publications [1,25,26,27], the Equation used to curve fit the phase transitions was a general equation for sigmoidal curves. Equation (1) used in the current study is more physiologically relevant as it is related to the “Hill” Equation used to measure enzyme kinetics. Another reason to recalculate the previously measured phase transitions is that the minimum and maximum ṽsym used in the older studies were less accurate. In studies before our 2007 study [26], the maximum ṽsym of 2854.5 cm−1 was estimated from phosphatidylcholine in CHCl3. In this study, we used a maximum ṽsym of 2855.36 cm−1 calculated from an isomeric distribution of hexanes [26]. In addition, in previous studies, the minimum ṽsym of 2849 cm−1 was estimated from dipalmitoylphosphatidylcholine at −20 °C. In this study, we used a minimum ṽsym of 2848.00 cm−1 calculated from distearoylphosphatidylcholine at −50 °C [26]. Data using the parameters in citation 26 are plotted in Figure 2a and Figure 7c,d,f.

3. Discussion

A major finding of this study is that human meibum lipid hydrocarbon chain unsaturation increases with age in agreement with previous FTIR [27], and NMR [23] spectroscopic studies. However, a greater number of samples, 69, were measured in the current study and, in previous studies, a 500 MHz NMR spectrometer was used and the resonance from the =CH of cholesterol was not resolved from the =CH due to the hydrocarbon cis =CH resonance. The contribution of the =CH resonance of cholesterol was significant, 20% of the total intensity of =CH resonances. In the current study, we used a 700 MHz NMR spectrometer that allowed the resolution of the two resonances thus circumventing this shortfall. In the current study, unsaturation was related to the amount of wax and cholesteryl esters. This is more meaningful and accurate than relating saturation to the intensity of all the resonances as in the previous study. Furthermore, in the previous study [24], the resonance at 1.39 ppm was the major resonance in the NMR spectra and was from protonated h-cyclohexane, a contaminant of the d-cyclohexane that was incorrectly assigned to the meibum lipid CH2 resonance. In this study, we used CDCl3 to circumvent this issue.
Increasing human meibum hydrocarbon chain unsaturation with age (Figure 2a) was related to hydrocarbon chain order (fluidity, Figure 7c), and a significant decrease in cooperativity (Figure 7f, r = 0.940, p < 0.01), and the phase transition temperature (Figure 7d, r = 0.982, p < 0.01). The change in these parameters was most dramatic between 1 and 20 years of age. The significant decrease in the phase transition parameters between 1 and 20 years of age can be explained by the observation that the phase transition temperature is linearly related to meibum lipid order [54]. The largest decline in the meibum phase transition temperature and hence the largest decline in lipid order occurred between 1 and 20 years of age (Figure 2b and Figure 7d). The change in blink rate with age (Figure 7a) was closely related with the increase in hydrocarbon chain disorder (Figure 7c), decrease in the plasma levels of free fatty acids (Figure 7b), phase transition temperature (Figure 7d) and cooperativity (Figure 7f). Correlation does not necessitate cause, but it is interesting that the breaks in the curves in Figure 7 occur around 20 years of age, at a Tanner stage V and adult level of development [55]. It is reasonable to speculate that endocrine changes with adolescence could be responsible for the observed break in the curves since the metabolism of lipids is under hormonal control.
Our catalytic saturation study showed that meibum hydrocarbon chain saturation was directly related to lipid order, phase transition temperature, cooperativity, ∆H and ∆S. Saturated hydrocarbon chains contain more trans rotamers and pack much more tightly than unsaturated hydrocarbon chains due to bends introduced into the hydrocarbon chains from the cis C=C bond. As the saturated hydrocarbon chains containing more trans rotamers pack more tightly together compared with hydrocarbon chains containing cis C=C bonds, it takes more enthalpy to break the van der Waal’s interactions between the saturated hydrocarbon chains, thus the ∆H of the lipid phase transition is greater for saturated hydrocarbon chains compared with unsaturated chains. As expected, a 40% increase in saturation from adult meibum to meibum from infants (Figure 2a) would be expected to increase the phase transition of adult meibum from 28 to 40 °C, similar to the observed increase from 28 to 36 °C, (Figure 7d). Our catalytic hydrogenation study also showed that the saturation driven increase in the phase transition temperature (Figure 5b) could account for an increase in lipid hydrocarbon chain order (Figure 5a) from about 30% in adults to about 80% for infants, a little more than the 60% order observed for infants (Figure 2b). Other factors such as hydrocarbon chain branching and hydroxyl groups could contribute to disordering meibum [34], whereas protein (Figure 7e) could contribute to the ordering of meibum [32,56]. Saturation correlated with the phase transition temperature of pure and native membranes and may contribute to lipid order more than phospholipid, wax, cholesteryl ester content or hydrocarbon chain length or branching, [54].
The lipid phase transition temperature and cooperativity measured by FTIR in this work were reasonably close to those measured in our previous FTIR study and those of others using different techniques (Table 1), especially considering that the age, race and gender of the samples were not exact. The value we obtained for the ∆H of the meibum lipid phase transition is much larger than that reported using differential scanning calorimetry (DSC, Table 1). The reason for this difference may be due to technical differences or because DSC measures the total ∆H of the phase transition which includes the ∆H of hydrocarbon and interface interactions whereas the ∆H reported in the current study measures the ∆H for the transition of a mole of trans rotamers to a mole of gauche rotamers. There may be about eight trans rotamers per hydrocarbon chain. From the maximum and minimum infrared ṽsym of the phase transition, we calculate that in the ordered “gel phase” at low temperature, 72% of the rotamers are trans and 18% of the rotamers are trans in the disordered “liquid crystal phase” at higher temperature. Therefore, we estimate that DSC measures the ∆H for only 53% of the total isomers, the ones that undergo a trans to gauche change. Because the hydrocarbon chains are not completely ordered (solid) below the phase transition temperature and not completely disordered (liquid) above the phase transition temperature, the transition is called a gel to liquid crystalline phase transition and not a melting. Meibum compositional differences in hydrocarbon chain saturation can account for meibum structural differences with age as suggested in the current study. Lipid saturation [56], order and phase transition temperature [54] are higher in donors with Meibomian gland dysfunction compared with adults. Intuitively, meibum should be ordered enough to flow out of the Meibomian glands and fluid enough to spread on the surface of the tears. The relationships between meibum lipid structure and tear film function are less clear with dry eye that they are with age. The hydrocarbon chain order and phase transition temperature of meibum from donors with dry eye and unstable tears is 49% trans and 28 °C, respectively, comparable to that of donors younger than 10 years old with extremely stable tears, 50% to 60% trans and 35 °C, respectively. Therefore, other factors in addition to meibum lipid structure such as elevated levels of protein [32,56], cooperative unit size [33], loss of squalene [29], inflammation [57], sebum [1], interactions between meibum and moieties in tears [58,59,60], differences between the lipid composition of tears and meibum [1,26,58,60], and aqueous deficiency could all contribute to functional derangements with dry eye. Future studies focused on the role of meibum structure and tear film function are needed.
The infrared spectroscopic parameters discussed above are relevant to bulk meibum in the Meibomian gland and on the surface of the eye lid. The change in structural order of meibum with age could also be related to the structural order of meibum on the surface of tears since most (94%) of the lipid on the tear film surface is not in contact with the aqueous interface. However, we used Langmuir trough technology to measure how saturation influenced the surface properties of meibum [25] and compared native meibum with meibum that was 100% saturated, a level that was not physiological. We have completed a study comparing the rheology of meibum at physiological saturation levels for comparison with the composition, structure and functional data from the current study. We may speculate that the functional consequence of a more ordered, more elastic saturated meibum as observed for infants and the higher maximum surface pressure observed in pressure area curves of saturated meibum compared to native meibum suggests that more saturated meibum films could be more stable, especially under the high shear stress of a blink [25].

4. Experimental Section

4.1. Materials

Silver chloride windows for infrared spectroscopy were obtained from Crystran Limited, Poole, UK. Platinum (IV) oxide was obtained from the Sigma Chemical Company (St. Louis, MO, USA).

4.2. Diagnosis of Normal Status

Normal status was assigned when the patient’s Meibomian gland orifices showed no evidence of keratinization or plugging with turbid or thickened secretions and no dilated blood vessels were observed on the eyelid margin. Normal donors did not recall having dry eye symptoms. Written informed consent was obtained from all donors and protocols and procedures were approved by the University of Louisville Institutional Review Board # 11.0319, August 2016. All procedures were in accord with the Declaration of Helsinki.

4.3. Collection and Extraction of Lipid from Meibum

Meibum lipid was expressed from the eye lids [61] and was collected with a platinum spatula with attention to avoiding scraping of the eyelid margin. Donors had no signs or symptoms of dry eye. Expressed meibum was dissolved in 1.5 mL CDCl3. The samples were pooled for catalytic hydrogenation.

4.4. Catalytic Hydrogenation

Half the pooled meibum was decanted to be catalytically hydrogenated. Saturated meibum was prepared as we did for sphingomyelin [26,62,63]. Platinum (IV) oxide (7.4 mg) was used as a catalyst to reduce the samples with hydrogen at room temperature and atmospheric pressure for approximately 4 h with stirring. Centrifugation was used to separate the catalyst from the solution. Catalytically saturated samples were quantitatively mixed with sample that was not catalytically saturated to provide mixtures containing 1%, 2%, 3%, 4%, 5%, 10%, 25%, 50%, and 67% of catalytically saturated meibum.

4.5. Saturation Analysis Using H-NMR

On the day of NMR measurement, the sample was sonicated under an atmosphere of argon gas in an ultrasonic bath (Branson 1510, Branson Ultrasonics, Danbury, CT, USA) for 10 min and placed into a NMR tube for spectral measurement. Meibum-CDCl3 samples were transferred from the microvial to a NMR tube using a glass pipet. Spectral data were acquired using a Varian VNMRS 700 MHz NMR spectrometer (Varian, Lexington, MA, USA) equipped with a 5-mm 1H{13C/15N}13C enhanced pulse-field gradient cold probe (Palo Alto, CA, USA). Spectra were acquired with a minimum of 250 scans, 45° pulse width, and a relaxation delay of 1.000 s. All spectra were obtained at 25 °C. Spectra were processed and integration of spectral bands was performed with GRAMS/386 software (Galactic Industries, Salem, NH, USA).
To quantify the relative level of cis hydrocarbon =CH (5.32 ppm) bonds, the intensity of the =CH resonance from cholesteryl esters (5.35 ppm) was subtracted from the total area of the 5.32 and 5.35 ppm resonances then divided by the sum of the resonances from cholesteryl and wax esters at 4.6 and 4.1 ppm, respectively.

4.6. Measurement of Lipid Phase Transitions Using FTIR Spectroscopy

Lipid phase transitions were measured as described previously [25]. About 500 µL of sample in CDCl3 was applied to a AgCl infrared window. The solvent was evaporated under a stream of Argon gas and the window was placed in a lyophilizer for 4 h to remove all traces of solvent. A Fourier transform infrared spectrometer (Nicolet 5000 Magna Series; Thermo Fisher Scientific, Inc., Waltham MA, USA) was used to measure the infrared spectra of the lipid on a AgCl window. The window was placed in a temperature-controlled infrared cell. The sample temperature was adjusted by an insulated water coil connected to a circulating water bath (model R-134A; Neslab Instruments, Newton, MA, USA) surrounding the cell. A thermistor touching the sample cell window was used to measure the sample temperature. The sample was cooled or heated at a rate of 1 °C/15 min. Temperatures were maintained within ±0.01 °C. Exactly 100 interferograms were recorded and averaged. Spectral resolution was set to 1.0 cm−1. Infrared data analysis was then performed (GRAMS/386 software; Galactic Industries, Salem, NH, USA). ṽsym was used to estimate the content of trans and gauche rotamers in the hydrocarbon chains. The OH–CH2 stretching region of the spectra was baselined between 3500 and 2700 cm−1. ṽsym was calculated from the center of mass of the CH symmetric stretching band by integrating the top 10% of the intensity of the band. The baseline for integrating the top 10% of the intensity of the band was parallel to the OH–CH region baseline. The change in ṽsym versus temperature was used to characterize lipid phase transitions as described previously [25]. Since rotamers are in either trans or gauche conformations, phase transitions were fit to a two-state sigmoidal equation using Sigma plot 10 software (Systat Software, Inc. Chicago, IL, USA):
ṽsym = (ṽsym)minimum + ((ṽsym)maximum − (ṽsym)minimum)/(1 + (temperature/Tc)hillslope)
ṽsym is the frequency of the symmetric CH2 stretching band near 2850 cm−1. Tc is the phase transition temperature.
ṽsym at 33.4 °C was extrapolated from the fit of the phase transition and then converting to lipid order which is the percentage of trans rotamers [25]. ∆H and ∆S were calculated from the slopes of Arrhenius plots [25].

4.7. Statistics

Curves were fit using Sigma plot 10 software (Systat Software, Inc., Chicago, IL, USA) and the confidence levels, p, were obtained from a critical value table of the Pearson product–moment correlation coefficient. A value of p < 0.05 was considered statistically significant. Error bars are the standard error of the mean.

5. Conclusions

Hydrocarbon chain saturation was directly related to lipid order, phase transition temperature, cooperativity, changes in enthalpy (∆H) and entropy (∆S) and could account for the changes in the lipid phase transition parameters observed with age. Unsaturation could contribute to decreased tear film stability with age.

Acknowledgments

Major support was obtained from the National Institute of Health EYO RO126180 (DB) and an unrestricted grant from Research to Prevent Blindness, Inc., New York, NY, USA, Aliza Williams received a fellowship from an Institutional Development Award from the National Institute of General Medical Sciences of the National Institutes of Health under grant number P20GM103436. Kayla Massay received a fellowship from the Physiology and Cardiovascular Minority Program and Collin Henry received a Medical School Student fellowship from the “Summer Research Scholar Program”, the University of Louisville, Louisville, KY, USA.

Author Contributions

Study Design/Grant: Douglas Borchman, Marta C. Yappert, Rahul Bhola and Kenneth N. Schikler; Contributed Samples: Rahul Bhola, Aparna Ramasubramanian, Ryan Blackburn, Samiyyah Sledge, Jonathan Austin and Collin Henry; Collection and Analysis of Data: Samiyyah Sledge, Collin Henry, Douglas Borchman, Kayla Massey, Shanzeh Sayied, and Aliza Williams; all authors have contributed to Edit the Manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mudgil, P.; Borchman, D.; Gerlach, D.; Yappert, M.C. Sebum/Meibum surface film interactions and phase transitional differences. Investig. Ophthalmol. Vis. Sci. 2016, 57, 2401–2411. [Google Scholar] [CrossRef]
  2. Green-Church, K.B.; Butovich, I.; Willcox, M.; Borchman, D.; Paulsen, F.; Barabino, S.; Glasgow, B.J. The international workshop on meibomian gland dysfunction: Report of the subcommittee on tear film lipids and lipid-protein interactions in health and disease. Investig. Ophthalmol. Vis. Sci. 2011, 52, 1979–1993. [Google Scholar] [CrossRef] [PubMed]
  3. Foulks, G.N. The correlation between the tear film lipid layer and dry eye disease. Surv. Ophthalmol. 2007, 52, 369–374. [Google Scholar] [CrossRef] [PubMed]
  4. Bron, A.J.; Tiffany, J.M.; Gouveia, S.M.; Yokoi, N.; Voon, L.W. Functional aspects of the tear film lipid Layer. Exp. Eye Res. 2004, 78, 347–360. [Google Scholar] [CrossRef] [PubMed]
  5. Mantelli, F.; Tiberi, E.; Micera, A.; Lambiase, A.; Visintini, F.; Bonini, S. MUC5AC over expression in tear film of neonates. Graefes Arch. Clin. Exp. Ophthalmol. 2007, 245, 1377–1381. [Google Scholar] [CrossRef] [PubMed]
  6. Isenberg, S.J.; del Signore, M.; Chen, A.; Wei, J.; Guillon, J. The lipid layer and stability of the preocular tear film in newborns and infants. Ophthalmology 2003, 110, 1408–1411. [Google Scholar] [CrossRef]
  7. Bacher, L.F. Factors regulating eye blink rate in young infants. Optom. Vis. Sci. 2010, 87, 337–343. [Google Scholar] [CrossRef] [PubMed]
  8. Lawrenson, J.G.; Birhah, R.; Murphy, P.J. Tear-film lipid layer morphology and corneal sensation in the development of blinking in neonates and infants. J. Anat. 2005, 206, 265–270. [Google Scholar] [CrossRef] [PubMed]
  9. Sforza, C.; Rango, M.; Galante, D.; Bresolin, N.; Ferrario, V.F. Spontaneous blinking in healthy persons: An optoelectronic study of eyelid motion. Ophthalmic Physiol. Opt. 2008, 28, 345–353. [Google Scholar] [CrossRef]
  10. Cho, P.; Yap, M. Age, gender, and tear break-up time. Optom. Vis. Sci. 1993, 70, 828–831. [Google Scholar] [CrossRef] [PubMed]
  11. Mohidin, N.; Bay, T.C.; Yapm, M. Non-invasive tear break-up time in normal Malays. Clin. Exp. Optom. 2002, 85, 37–41. [Google Scholar] [CrossRef] [PubMed]
  12. Ozdemir, M.; Temizdemir, H. Age- and sex- related tear function changes in normal population. Eye 2010, 24, 79–83. [Google Scholar] [CrossRef]
  13. Maïssa, C.; Guillon, M. Tear film dynamics and lipid layer characteristics-effect of age and sex. Contact Lens Anterior Eye 2010, 33, 176–182. [Google Scholar] [CrossRef] [PubMed]
  14. Sun, W.S.; Baker, R.S.; Chuke, J.C.; Rouholiman, B.R.; Hasan, S.A.; Gaza, W.; Stava, M.W.; Porter, J.D. Age-related changes in human blinks. Passive and active changes in eyelid kinematics. Investig. Ophthalmol. Vis. Sci. 1997, 38, 92–99. [Google Scholar]
  15. Lavezzo, M.M.; Schellini, S.A.; Padovani, C.R. Eye blink in newborn and preschool-age children. Acta Ophthalmol. 2008, 86, 275–278. [Google Scholar] [CrossRef] [PubMed]
  16. Craig, J.; Tomlinson, A. Age and sex effects on the normal tear film. Adv. Exp. Med. Biol. 1998, 438, 411–415. [Google Scholar] [PubMed]
  17. Henderson, J.W.; Prough, W.A. Influence of age, sex on flow of tears. Arch. Ophthalmol. 1950, 43, 224–231. [Google Scholar] [CrossRef]
  18. Guillon, M.; Maïssa, C. Tear film evaporation: Effect of age and sex. Contact Lens Anterior Eye 2010, 33, 171–175. [Google Scholar] [CrossRef] [PubMed]
  19. Cruz, A.A.; Garcia, D.M.; Pinto, C.T.; Cehetti, S.P. Spontaneous eyeblink activity. Ocul. Surf. 2011, 9, 29–41. [Google Scholar] [CrossRef]
  20. Tsubota, K.; Hata, S.; Okusawa, Y.; Egami, F.; Ohtsuki, T.; Nakamori, K. Quantitative videographic analysis of blinking in normal subjects and patients with dry eye. Arch. Ophthalmol. 1996, 114, 715–720. [Google Scholar] [CrossRef] [PubMed]
  21. Tsubota, K. Tear dynamics and dry eye. Prog. Retin. Eye Res. 1998, 17, 565–596. [Google Scholar] [CrossRef]
  22. Tomlinson, A.; Khanal, S. Assessment of tear film dynamics: Quantification approach. Ocul. Surf. 2005, 3, 81–95. [Google Scholar] [CrossRef]
  23. Borchman, D.; Yappert, M.C.; Milliner, S.; Duran, D.; Cox, G.W.; Smith, R.J.; Bhola, R. 13C and 1H-NMR ester region resonance assignments and the composition of human infant and child meibum. Exp. Eye Res. 2013, 112, 151–159. [Google Scholar] [CrossRef] [PubMed]
  24. Borchman, D.; Foulks, G.N.; Yappert, M.C.; Milliner, S.E. Changes in human meibum lipid composition with age using NMR spectroscopy. Investig. Ophthalmol. Vis. Sci. 2012, 53, 475–482. [Google Scholar] [CrossRef] [PubMed]
  25. Mudgil, P.; Borchman, D.; Yappert, M.C.; Duran, D.; Cox, G.W.; Smith, R.J.; Bhola, R.; Dennis, G.R.; Whitehall, J.S. Human meibum saturation and lipid order. Exp. Eye Res. 2013, 116, 79–85. [Google Scholar] [CrossRef] [PubMed]
  26. Borchman, D.; Foulks, G.N.; Yappert, M.C.; Ho, D.V. Temperature-induced conformational changes in human tear lipids hydrocarbon chains. Biopolymers 2007, 87, 124–133. [Google Scholar] [CrossRef] [PubMed]
  27. Borchman, D.; Foulks, G.N.; Yappert, M.C.; Kakar, S.; Podoll, N.; Rychwalski, P.; Schwietz, E. Physical changes in human meibum with age as measured by infrared spectroscopy. Ophthalmic Res. 2010, 44, 34–42. [Google Scholar] [CrossRef] [PubMed]
  28. Borchman, D.; Foulks, G.N.; Yappert, M.C. Confirmation of changes in human meibum lipid infrared spectra with age using principal component analysis. Curr. Eye Res. 2010, 35, 778–786. [Google Scholar] [CrossRef] [PubMed]
  29. Borchman, D.; Yappert, M.C.; Milliner, S.; Bhola, R. Confirmation of squalene in human eye lid lipid by heteronuclear single quantum correlation spectroscopy. Lipids 2013, 48, 1269–1277. [Google Scholar] [CrossRef] [PubMed]
  30. Kóta, Z.; Debreczeny, M.; Szalontai, B. Separable contributions of ordered and disordered lipid fatty acyl chain segments to ṽCH2 bands in model and biological membranes: A fourier transform infrared spectroscopic study. Biospectroscopy 1999, 5, 169–178. [Google Scholar] [CrossRef]
  31. Butovich, I.A.; Arciniega, J.C.; Wojtowicz, J.C. Meibomian lipid films and the impact of temperature. Investig. Ophthalmol. Vis. Sci. 2010, 51, 5508–5518. [Google Scholar] [CrossRef]
  32. Butovich, I.A.; Lu, H.; McMahon, A.; Ketelson, H.; Senchyna, M.; Meadows, D.; Campbell, E.; Molai, M.; Linsenbardt, E. Biophysical and morphological evaluation of human normal and dry eye meibum using hot stage polarized light microscopy. Investig. Ophthalmol. Vis. Sci. 2014, 55, 87–101. [Google Scholar] [CrossRef] [PubMed]
  33. Lu, H.; Wojtowicz, J.C.; Butovich, I.A. Differential scanning calorimetric evaluation of human meibomian gland secretions and model lipid mixtures: Transition temperatures and cooperativity of melting. Chem. Phys. Lipids 2013, 170–171, 55–64. [Google Scholar] [CrossRef]
  34. Butovich, I.A. On the lipid composition of human meibum and tears: Comparative analysis of nonpolar lipids. Investig. Ophthalmol. Vis. Sci. 2008, 49, 3779–3789. [Google Scholar] [CrossRef] [PubMed]
  35. Petrikovsky, B.M.; Kaplan, G.; Holsten, N. Eyelid movements in normal human fetuses. J. Clin. Ultrasound 2003, 31, 299–301. [Google Scholar] [CrossRef] [PubMed]
  36. Bacher, L.F.; Allen, K.J. Sensitivity of the rate of spontaneous eye blinking to type of stimuli in young infants. Dev. Psychobiol. 2009, 51, 186–197. [Google Scholar] [CrossRef] [PubMed]
  37. Bacher, L.F.; Smotherman, W.P. Systematic temporal variation in the rate of spontaneous eye blinking in human infants. Dev. Psychobiol. 2004, 44, 140–145. [Google Scholar] [CrossRef] [PubMed]
  38. Zametkin, A.J.; Stevens, J.R.; Pittman, R. Ontogeny of spontaneous blinking and of habituation of the blink reflex. Ann. Neurol. 1979, 5, 453–457. [Google Scholar] [CrossRef] [PubMed]
  39. Bentivoglio, A.R.; Bressman, S.B.; Cassetta, E.; Carretta, D.; Tonali, P.; Albanese, A. Analysis of blink rate patterns in normal subjects. Mov. Disord. 1997, 12, 1028–1034. [Google Scholar] [CrossRef] [PubMed]
  40. Sahlin, S.; Chen, E. Gravity, blink rate, and lacrimal drainage capacity. Am. J. Ophthalmol. 1997, 124, 758–764. [Google Scholar] [CrossRef]
  41. Zaman, M.L.; Doughty, M.J. Some methodological issues in the assessment of the spontaneous eyeblink frequency in man. Ophthalmic Physiol. Opt. 1997, 17, 421–432. [Google Scholar] [CrossRef]
  42. Nakamori, K.; Odawara, M.; Nakajima, T.; Mizutani, T.; Tsubota, K. Blinking is controlled primarily by ocular surface conditions. Am. J. Ophthalmol. 1997, 124, 24–30. [Google Scholar] [CrossRef]
  43. Biousse, V.; Skibell, B.C.; Watts, R.L.; Loupe, D.N.; Drews-Botsch, C.; Newman, N.J. Ophthalmologic features of Parkinson’s disease. Neurology 2004, 62, 177–180. [Google Scholar] [CrossRef] [PubMed]
  44. Doughty, M.J. Effects of background lighting and retinal illuminance on spontaneous eyeblink activity of human subjects in primary eye gaze. Eye Contact Lens Sci. Clin. Pract. 2013, 39, 138–146. [Google Scholar] [CrossRef]
  45. Begley, C.G.; Port, N.; Bradley, A.; Braun, R. The effects of increasing ocular surface stimulation on blinking and tear secretion. Investig. Ophthalmol. Vis. Sci. 2015, 56, 4211–4220. [Google Scholar]
  46. Cardona, G.; Perez-Cabre, E.; Rodriguez, M. Blink rate and incomplete blinks in six different controlled hard-copy and electronic reading conditions. Investig. Ophthalmol. Vis. Sci. 2015, 56, 6679–6685. [Google Scholar]
  47. Michael, J.; Doughty, M.J. Assessment of short-term variability in human spontaneous blink rate during video observation with or without head/chin support. Clin. Exp. Optom. 2016, 99, 135–141. [Google Scholar]
  48. Kruis, A.; Slagter, H.A.; Bachhuber, D.R.W.; Davidson, R.J.; Lutz, A. Effects of meditation practice on spontaneous eyeblink rate. Psychophysiology 2016, 53, 749–758. [Google Scholar] [CrossRef] [PubMed]
  49. Nosch, D.S.; Pult, H.; Albon, J.; Purslow, C.; Murphy, P.J. Relationship between corneal sensation, blinking, and tear film quality. Optom. Vis. Sci. 2016, 93, 471–481. [Google Scholar] [CrossRef]
  50. Lane, K.J.; Ousler, G.W., III; Angjeli, E.; Smith, L.M.; Rodriguez, J.D.; Bateman, K.M.; Abelson, M.B. Diurnal tracking of blink and relationship to signs and symptoms of dry eye. Cornea 2016, 35, 1104–1111. [Google Scholar]
  51. Nakano, T. The right angular gyrus controls spontaneous eyeblink rate: A combined structural MRI and TMS study. Cortex 2017, 88, 186–191. [Google Scholar] [CrossRef] [PubMed]
  52. Frohnert, B.; Jacobs, D.R., Jr.; Steinberger, J.; Moran, A.; Steffen, L.M.; Sinaiko, A.R. Relation between serum free fatty acids and adiposity, insulin resistance, and cardiovascular risk factors from adolescence to adulthood. Diabetes 2013, 62, 3163–3169. [Google Scholar] [CrossRef] [PubMed]
  53. Felix, P.; Heald, G.A.; Seakold, W.; Morrison, D. Plasma levels of free fatty acids in adolescents. Am. J. Clin. Nutr. 1967, 20, 1010–1014. [Google Scholar]
  54. Borchman, D.; Foulks, G.N.; Yappert, M.C.; Bell, J.; Wells, E.; Neravetla, S.; Greenstone, V. Human meibum lipid conformation and thermodynamic changes with meibomian-gland dysfunction. Investig. Ophthalmol. Vis. Sci. 2011, 52, 3805–3817. [Google Scholar] [CrossRef] [PubMed]
  55. Marshall, W.A.; Tanner, J.M. Variations in the pattern of pubertal changes in boys. Arch. Dis. Child. 1970, 45, 3–23. [Google Scholar] [CrossRef]
  56. Borchman, D.; Foulks, G.N.; Yappert, M.C. Changes in human meibum lipid with meibomian gland dysfunction using principal component analysis. Exp. Eye Res. 2010, 91, 246–256. [Google Scholar] [CrossRef] [PubMed]
  57. Lim, A.; Wenk, M.R.; Tong, L. Lipid-based therapy for ocular surface inflammation and disease. Trends Mol. Med. 2015, 21, 736–748. [Google Scholar] [CrossRef] [PubMed]
  58. Borchman, D.; Foulks, G.N.; Yappert, M.C.; Tang, D.; Ho, D.V. Spectroscopic evaluation of human tear lipids. Chem. Phys. Lipids 2007, 147, 87–102. [Google Scholar] [CrossRef] [PubMed]
  59. Svitova, T.F.; Lin, M.C. Dynamic interfacial properties of human tear-lipid films and their interactions with model-tear proteins in vitro. Adv. Colloid Interface Sci. 2016, 233, 4–24. [Google Scholar] [CrossRef]
  60. Pucker, A.D.; Haworth, K.M. The presence and significance of polar meibum and tear lipids. Ocul. Surf. 2015, 13, 26–42. [Google Scholar] [CrossRef] [PubMed]
  61. Kilp, H.; Schmid, E.; Kirchner, L.; Pohl, A. Tear Film Observation by Reflecting Microscopy and Differential Interference Contrast Microscopy. In The Preocular Tear Film in Health, Disease and Contact Lens Wear; Holly, F.J., Ed.; The Dry Eye Institute: Lubbock, TX, USA, 1986; pp. 564–569. [Google Scholar]
  62. Ferguson, S.R.; Borchman, D.; Yappert, M.C. Confirmation of the identity of the major phospholipid in human lens membranes. Investig. Ophthalmol. Vis. Sci. 1996, 37, 1703–1706. [Google Scholar]
  63. Ferguson-Yankey, S.; Borchman, D.; Taylor, K.G.; DuPre, D.B.; Yappert, M.C. Conformational studies of sphingolipids by NMR spectroscopy. I. dihydrosphingomyelin. Biochim. Biophys. Acta 2000, 1467, 307–325. [Google Scholar] [CrossRef]
Ijms 18 01862 g001 550
Figure 1. Typical proton nuclear magnetic resonance (1H-NMR) spectrum of meibum from a five-year-old Caucasian male.
Figure 1. Typical proton nuclear magnetic resonance (1H-NMR) spectrum of meibum from a five-year-old Caucasian male.
Ijms 18 01862 g001
Ijms 18 01862 g002 550
Figure 2. (a) Average hydrocarbon chain double bond content of meibum calculated from 1H-NMR spectra. Numbers in parenthesis are the number of samples; and (b) Lipid order (gel) to disorder (liquid crystal) phase transition parameters of meibum. Recalculated from (see text): (●) Reference [27]; (■) Reference [1]; (▲) this study; (▼) and Reference [26]. (▬) Curve fit to data using the parameter, hyperbola, hyperbolic decay equation: f = y0 + (a × b)/(b + x). All donors were normal and did not have signs or symptoms of dry eye. Data are average ± the standard error of the mean.
Figure 2. (a) Average hydrocarbon chain double bond content of meibum calculated from 1H-NMR spectra. Numbers in parenthesis are the number of samples; and (b) Lipid order (gel) to disorder (liquid crystal) phase transition parameters of meibum. Recalculated from (see text): (●) Reference [27]; (■) Reference [1]; (▲) this study; (▼) and Reference [26]. (▬) Curve fit to data using the parameter, hyperbola, hyperbolic decay equation: f = y0 + (a × b)/(b + x). All donors were normal and did not have signs or symptoms of dry eye. Data are average ± the standard error of the mean.
Ijms 18 01862 g002
Ijms 18 01862 g003 550
Figure 3. Infrared CH stretching region for: (a) a pool of meibum lipids from four normal adult donors; (b) meibum lipids from the same pool that was catalytically hydrogenated and (c) AgCl window.
Figure 3. Infrared CH stretching region for: (a) a pool of meibum lipids from four normal adult donors; (b) meibum lipids from the same pool that was catalytically hydrogenated and (c) AgCl window.
Ijms 18 01862 g003
Ijms 18 01862 g004 550
Figure 4. Lipid phase transition of meibum pooled from six Caucasian females 23, 24, 28, 32, 66, and 68 years old; a 29-year-old Caucasian male; and an 32-year-old African American female. The CH2 symmetric stretching frequency is related to lipid structural order. The higher the value for the frequency, the more disordered the lipid. Symbols are different trials. Symbols are for different trials.
Figure 4. Lipid phase transition of meibum pooled from six Caucasian females 23, 24, 28, 32, 66, and 68 years old; a 29-year-old Caucasian male; and an 32-year-old African American female. The CH2 symmetric stretching frequency is related to lipid structural order. The higher the value for the frequency, the more disordered the lipid. Symbols are different trials. Symbols are for different trials.
Ijms 18 01862 g004
Ijms 18 01862 g005 550
Figure 5. Lipid phase transition parameters for human meibum from Figure 4 that was mixed with catalytically hydrogenated meibum from the same pool. The percent saturation is the percentage of catalytically saturated meibum mixed with normal meibum. Data are average ± the standard error of the mean and reflect the experimental error for one sample. (▬) Best fit, linear regression (a,c–e) with a line order of 1 and (b) linear regression with a line order of 3.
Figure 5. Lipid phase transition parameters for human meibum from Figure 4 that was mixed with catalytically hydrogenated meibum from the same pool. The percent saturation is the percentage of catalytically saturated meibum mixed with normal meibum. Data are average ± the standard error of the mean and reflect the experimental error for one sample. (▬) Best fit, linear regression (a,c–e) with a line order of 1 and (b) linear regression with a line order of 3.
Ijms 18 01862 g005
Ijms 18 01862 g006 550
Figure 6. A typical Arrhenius plot for a pool of human meibum used for the saturation study. The slope of the line is used to calculate the changes in enthalpy (∆H) and entropy (∆S) of the lipid phase transition.
Figure 6. A typical Arrhenius plot for a pool of human meibum used for the saturation study. The slope of the line is used to calculate the changes in enthalpy (∆H) and entropy (∆S) of the lipid phase transition.
Ijms 18 01862 g006
Ijms 18 01862 g007 550
Figure 7. (a) From References [5,8,14,15,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51]; (b) (▲) from Reference [52] and (●) Reference [53]; (c,d,f) hydrocarbon chain order at 33.4 °C and phase transition parameters for the lipid order (gel) to disorder (liquid crystal) phase transition of human meibum. (●) Reference [26]; (■) Reference [1]; (▲) this study and (▼) Reference [25]. (▬) Curve fit to data using the parameter, hyperbola, hyperbolic decay equation: f = y0 + (a × b)/(b + x); (e) Data from Reference [54]. All donors were normal and did not have signs or symptoms of dry eye. Data are average ± the standard error of the mean.
Figure 7. (a) From References [5,8,14,15,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51]; (b) (▲) from Reference [52] and (●) Reference [53]; (c,d,f) hydrocarbon chain order at 33.4 °C and phase transition parameters for the lipid order (gel) to disorder (liquid crystal) phase transition of human meibum. (●) Reference [26]; (■) Reference [1]; (▲) this study and (▼) Reference [25]. (▬) Curve fit to data using the parameter, hyperbola, hyperbolic decay equation: f = y0 + (a × b)/(b + x); (e) Data from Reference [54]. All donors were normal and did not have signs or symptoms of dry eye. Data are average ± the standard error of the mean.
Ijms 18 01862 g007
Table
Table 1. Meibum Lipid Phase Transition Parameters.
Table 1. Meibum Lipid Phase Transition Parameters.
ParameterThis Study FTIRFTIR [1]Microscopy [31]Birefringence [32]DSC [33]Light Reflectance [34]
Hydrocarbon Chain Order (% trans)31 ± 236 ± 2
Phase Transition Temperature (°C)29.6 ± 0.428.9 ± 0.632.1 ± 0.132 ± 230.2 ± 0.132 ± 1
Cooperativity10 ± 18 ± 119.3 ± 0.17 ± 1, 34 ± 1016, 40, 100
Change in enthalpy (Kcal/mol)205 ± 5196 ± 2 6 ± 1
Change in entropy (Kcal/mol/deg)0.68 ± 0.020.650 ± 0.008
Age (years)39 ± 731 ± 629, 5053, 38 ± 450 ± 247
Sample demographics1 pool of 81 pool of 71 pool of 21 donor and 1 pool of 41 pool of 21 pool of 5

Share and Cite

MDPI and ACS Style

Sledge, S.; Henry, C.; Borchman, D.; Yappert, M.C.; Bhola, R.; Ramasubramanian, A.; Blackburn, R.; Austin, J.; Massey, K.; Sayied, S.; et al. Human Meibum Age, Lipid–Lipid Interactions and Lipid Saturation in Meibum from Infants. Int. J. Mol. Sci. 2017, 18, 1862. https://doi.org/10.3390/ijms18091862

AMA Style

Sledge S, Henry C, Borchman D, Yappert MC, Bhola R, Ramasubramanian A, Blackburn R, Austin J, Massey K, Sayied S, et al. Human Meibum Age, Lipid–Lipid Interactions and Lipid Saturation in Meibum from Infants. International Journal of Molecular Sciences. 2017; 18(9):1862. https://doi.org/10.3390/ijms18091862

Chicago/Turabian Style

Sledge, Samiyyah, Collin Henry, Douglas Borchman, Marta C. Yappert, Rahul Bhola, Aparna Ramasubramanian, Ryan Blackburn, Jonathan Austin, Kayla Massey, Shanzeh Sayied, and et al. 2017. "Human Meibum Age, Lipid–Lipid Interactions and Lipid Saturation in Meibum from Infants" International Journal of Molecular Sciences 18, no. 9: 1862. https://doi.org/10.3390/ijms18091862

APA Style

Sledge, S., Henry, C., Borchman, D., Yappert, M. C., Bhola, R., Ramasubramanian, A., Blackburn, R., Austin, J., Massey, K., Sayied, S., Williams, A., Georgiev, G., & Schikler, K. N. (2017). Human Meibum Age, Lipid–Lipid Interactions and Lipid Saturation in Meibum from Infants. International Journal of Molecular Sciences, 18(9), 1862. https://doi.org/10.3390/ijms18091862

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop