Abstract
Conjunctivochalasis is a common cause of tear dysfunction due to the conjunctiva becoming loose and wrinkly with age. The current solutions to this disease include either surgical excision in the operating room, or thermoreduction of the loose tissue with hot wire in the clinic. We developed a near-infrared laser thermal conjunctivoplasty system. The system utilizes a 1460-nm programmable laser diode system as the light source. At this wavelength, a water absorption peak exists and the blood absorption is minimal, so the heating of redundant conjunctiva is even and there is no bleeding. A miniaturized handheld probe delivers the laser light and reshapes the laser into a 10 × 1 mm2 line on the working plane. A foot pedal is used to deliver a preset number of calibrated laser pulses. A fold of loose conjunctiva is grasped by a pair of forceps. The NIR laser light is delivered through an optical fiber and a laser line is aimed exactly on the conjunctival fold by a cylindrical lens. Ex vivo experiments using porcine eye was performed to investigate the induced shrinkage of conjunctiva and decide the optimal laser parameters. It was found that up to 45% of conjunctiva shrinkage could be achieved.
Similar content being viewed by others
Introduction
The sclera, white part of the eye, is covered by a thin, clear membrane called the conjunctiva. Like skin, the conjunctiva becomes loose and wrinkly with age; this degenerative condition is called conjunctivochalasis1. The loose folds of conjunctiva often disrupt the uniform distribution of tears and can cause constant eye irritation and blurred vision2,3. In severe cases, the conjunctival folds protrude onto the inferior eye lid margin and are traumatized by the lid during blinking. Furthermore, the lid’s skin is also irritated and altered by the displaced tear.
Conjunctivochalasis is a common cause of tear dysfunction (“dry eyes”); however, it does not respond to the usual dry eye treatments such as artificial tears, punctal plugs and anti-inflammatory drops. Effective treatment requires surgical reduction or excision of the redundant conjunctiva to reestablish the inferior tear meniscus and normal tear dynamics. Conjunctivochalasis is typically diagnosed by evaluating the conjunctiva for redundant folds that prolapse onto the lower eyelid and obliterate the tear meniscus in that region. Cross-sectional optical coherence tomography (OCT) allows us for evaluating the severity of the condition and the effectiveness of surgery to remove the redundant tissue4,5,6.
Surgical means used to remove redundant conjunctival tissue is an effective way to treat conjunctivochalasis. Several techniques currently are employed in the surgery, such as thermocautery, electrocautery, and paste-pinch-cut7,8,9,10. Thermocautery and electrocautery reach temperatures exceeding the boiling point of water and can burn the conjunctival epithelium, underlying stroma, and surrounding tissue. This creates a full thickness burn wound that is generally painful, takes up to one month to fully heal, and occasionally induces excessive inflammation with scarring. Furthermore, a chronic inflammatory conjunctival mass called pyogenic granuloma could result, which would necessitate long-term anti-inflammatory eye drops and possibly further surgery. Poor cosmetic appearance (i.e. red blots due to bleeding in surface tissue) or scarring during the long healing period also deters patients from these techniques. Some surgeons have reservations about using thermo- or electrocautery because of the risk of causing excessive damage to the conjunctiva or cornea if the patient’s eye moves during the procedure. Surgical conjunctival excision with the addition of an amniotic membrane transplant (attached by fibrin glue or suture) can improve the healing course, but must be performed in the operating room, which markedly increases cost. The paste-pinch-cut technique uses forceps to hold the conjunctival and cuts the redundant part directly. The wound will be polymerized with pre-pasted fibrin sealant. Compared with thermocautery and electrocautery, this approach could bring less pain to the patient but must be performed in an operating room and still needs several weeks to resolve the inflamed tissus surface11. Laser conjunctivoplasty, on the other hand, can well constrain the thermal energy in the target area and depth and precisely control the amount of energy thus avoiding burning the conjunctival epithelium and underlying stroma.
Argon lasers have been used for laser conjunctivoplasty12,13. However, they work at a wavelength of 532 nm and the laser energy is primarily absorbed by blood and can lead to vessel rupture and hemorrhage. Because the water content of conjunctiva is around 70%14, heating the conjunctiva using water absorption is a more effective and safe way to achieve shrinkage. In this paper, we proposed to use the near-infrared (NIR) light for LTC, because the blood absorption at this waveband is much weaker than that at visible region15,16 while two strong water absorption peaks exist at 1460 nm and 1940 nm17. Compared with the expensive and inefficient lasers at 1940 nm, the diode laser at 1460 nm has an electro-optical efficiency of >40% and the price of the commercial products is reasonable. Thus we developed a miniature-probed-based LTC system using a laser diode system at this wavelength. Ex vivo experiments on porcine eye were conducted to demonstrate the effectiveness of the NIR LTC.
Results
Temperature changes in time domain
During the laser heating procedure, the temperature continues to increase until the laser is turned off. However, to optimize the shrinking, the temperature should be kept in a range to avoid burning and boiling. The collagen shrinkage temperate is around 62 °C18 and to avoid boiling, the temperature should be kept below 100 °C. So idea temperature is between 62 °C and 100 °C. In order to achieve this, we use pulse laser for the heating and a pulse repetition rate of 1 Hz and a working duration of 4 seconds (a total of 4 pulses) were selected.
Figure 1 shows the temperature changes during the LTC. The highest temperatures in each of recorded thermal camera image sequences were obtained. It can be seen that the peak temperature increases very fast during the laser pulse illumination period and decreases relatively slowly during the laser-off period. The effects of peak laser power and pulse duty cycle were investigated. With the laser peak power of 1.51 W and duty cycle of 10% (100 ms pulse duration), the temperature increase is very small and the temperature is always below 40 °C [Fig. 1(a)]. When the duty cycle was increased to 20% (200 ms pulse duration) and 30% (300 ms pulse duration), the cumulative effect of multiple pulses can be clearly seen. The temperature will increase during the laser-illumination period and the temperature will drop during the laser-off period. This temperature increase-decrease process will continue with the increase of the pulse numbers. For a peak power of 3 W, the maximum temperature for the 30% duty cycle exceeds 100 °C while the maximum temperature for the 10% duty cycle is around 65 °C [Fig. 1(b)]. When the peak power was further increased to 5.77 W, the maximum temperature exceeds 100 °C for all three different duty cycles [Fig. 1(c)].
OCT imaging
Real-time OCT imaging was used to monitor the LTC-induced shrinkage process. This dynamic process can be clearly visualized in a video of the time-sequence OCT B-scans (Media 1), but is difficult to see from the OCT images before [Fig. 2 (a)] and after [Fig. 2(b)] the LTC. To further visualize and quantify the shrinkage from these images, particle image velocimetry (PIV) was employed19. The green boxes in Fig. 2(a) and (b) were employed to calculate the PIV using ImageJ20. The results are shown in Fig. 2(c). The different colors correspond to different speeds and the arrows indicate the moving direction. The red arrow groups at both sides point to the central region, where was heated by the laser. This is a clear indication of the shrinkage process.
LTC induced shrinkage and optimized laser parameters
We investigate the influence of different laser parameters on the shrinkage. Figure 3 shows the microscope images of the porcine eyes before and after LTC with the laser powers of 2.33 W, 3.1 W, and 3.66 W. The pulse duty cycle was set to 20% here. As demonstrated, the shrinkage increases as the laser peak power increases. For the peak power of 2.33 W, the shrinkage is 21%. As the peak power increases to 3 W and 3.66 W, the shrinkage increased to 36% and 45%, respectively.
From the experiment above, 3 W peak power and 20% duty cycle is optimum for both the shrinkage and the temperature (the measured temperature is 88 °C). So we used the 3-W peak power to investigate the influence of the pulse duty cycle on the shrinkage. Figure 4 shows the microscope images of the shrinkage at the pulse duty cycles of 10%, 20% and 30%. It can be seen from the images, at 10% of the duty cycle, the shrinkage is minimal (∼8%). When the duty cycle increases to 30%, the shrinkage is 45%. However, the tissue surface was damaged due to a high temperature of more than 100 °C [Fig. 1(b)].
To determine the optimal laser parameters for the LTC, it is necessary to analyze the relationship between the shrinkage and the tissue temperature at different laser parameters. For a specific laser parameter combination, the shrinkage and temperature was measured 6–8 times on different samples and the averaged values were obtained. Figure 5(a) shows the shrinkage and temperature as a function of laser peak power. The default duty cycle of 20% was used. It can be seen the temperature keeps increasing as the peak power increases while the shrinkage starts to saturate when the peak power is higher than 3 W. A peak power of 3 W is ideal for the LTC application. The shrinkage and temperature as a function of pulse duty cycle is shown in Fig. 5(b). The peak power was set at 3 W. When the duty cycle is 30%, the temperature is above 110 °C while the shrinkage is comparable with that of the 20% duty cycle.
Shrinkage using angled forceps
The probe is designed to work with the angled forceps so that the working distance is guaranteed, the shining angle is normal and only the tissue fold held by the forceps is heated. Due to the force from the forceps, the tissue may not be able to recover to its original shape and this may introduce a false shrinkage. To eliminate this, we test the distance changes between the parallel maker lines due to only the force from the angled forceps. As shown in the left panels of Fig. 6(b), a slight change (10%) between the marker lines was found. After the LTC, the shrinkage was measured to be 45% [right panels of Fig. 6(b)].
Discussion and Conclusion
The ex vivo porcine eye experiments of the NIR LTC is a proof-of-concept of this new technique. The texture and contents of the conjunctiva of the porcine eye is very close to those of the human eye, so we can expect to achieve similar shrinkage from conjunctivoplasty surgery of human subjects by using this NIR LTC system. However, before testing it on human subjects, further in vivo animal experiments are still needed based on the following considerations. (1) The initial temperature of the porcine eye was about 20 °C, which is lower than that of the human patients. Therefore, the optimal laser parameters decided here may not be appropriate for the actual treatment of redundant conjunctiva. The laser parameters need to be calibrated by using in vivo animals with a similar body temperature as human being. (2) The ex vivo eye does not have blood circulations, therefore a living eye is needed to assesses conjunctival hyperemia (a measure of inflammation), hemorrhage (bleeding inside tissue), and blanching (a measure of ischemia). (3) The ex vivo eye does not have a healing response, therefore in vivo animal experiment are needed to assess epithelial healing, fibroblast repopulation of the stroma, any scarring formation, and whether the shrinkage is maintained in the long term.
In summary, we have developed a NIR LTC using a 1460-nm wavelength laser, which can be efficiently absorbed by water. We developed a compact handheld probe with optimized features for LTC, including a line-shape focus spot and angled forceps. The laser parameters were optimized and ex vivo experiments of on porcine eye was performed. Up to 45% of conjunctiva shrinkage was observed.
Materials and Methods
1460-nm programmable laser diode system
Figure 7(a) shows the schematic of the programmable pulsed laser diode system for the LTC. A fiber-coupled high-power laser diode module with a maximum continuous-wave (CW) power of 12 W (M1F2S22-1470.10-12C-SS5.x, DILAS, Tucson, AZ, USA) was employed as the light source. The light was output through a multimode fiber with a core diameter of 200 μm and a numerical aperture (NA) of 0.22. The other end of the fiber was connected to a handheld probe. A 650-nm laser diode was integrated into the source laser for aiming purpose. The power for the aiming light was 200 μW.
Custom built control circuits were used to drive the 1460-nm laser module. The pulse duty cycle, repetition rate, output power and working duration are tunable through a LabVIEW (LabVIEW 2015, National Instruments, Austin, TX, USA) based control software. The duty cycle, repetition rate, and working duration is arbitrarily tunable. Limited by the capacities of the commercial laser diode module, the maximum output power and the shortest pulse width are 12 W and 100 μs, respectively. A 650-nm aiming light could be enabled and disabled by the operator from the control software. A foot pedal is used as a trigger for the laser output. Once the laser output is triggered, the 650-nm aiming light will be turned off automatically and the 1460-nm laser will be delivered to the probe according to the preset parameters. However, if the foot pedal is released during the procedure, the infrared laser will be turned off immediately. The photograph of the prototype of the LTC laser system is shown in Fig. 1(b). The whole laser system was assembled into a case with a dimension of 41 × 36 × 15 cm3.
The power, spectral, and temporal characteristics of the 1460-nm laser system were measured. The relationship between output power and driving current is shown in Fig. 8(a) and a linear relationship between them could be found. The lasing threshold current is 6.8 A. The power was measured by a thermal power sensor (S310C, Thorlabs, Newton, NJ, USA). Figure 8(b) shows the laser spectrum with the output power of 1.2 W, 4.1 W and 6.8 W, respectively. The results were measured by an optical spectrum analyzer (AQ6370C, YOKOGAWA, Tokyo, Japan). The wavelength at the peak power is ∼1456 nm. No obvious spectral shift is observed as the output power increases. The spectral bandwidth increases from 3.1 nm at 1.2 W to 3.9 nm at 6.8 W. Compared with the bandwidth of the water absorption peak (∼100 nm), the influence of this spectral broadening on the LTC can be neglected. The temporal features of this laser system were measured by a photodetector (PDA10CF, Thorlabs, Newton, NJ, USA) and an oscilloscope (MSO4104B-L, Tektronix, Beaverton, OR, USA). Figure 8(c) demonstrates the temporal shapes of output pulses with the pulse widths of 100 ms, 200 ms, and 300 ms. Figure 8(d) shows the output pulse trains with the pulse numbers following the Fibonacci sequence, which proves the output pulse number can be precisely controlled by the working duration.
Handheld LTC probe
A cylindrical lens is used to focus the SMA fiber output into a line [Fig. 9(b)]. This design reduces the power loss and keeps the probe compact. The ray-tracing (OpticStudio, Zemax, LLC, Seattle, WA, USA) of the light path is shown in the left of Fig. 9(a). The cylindrical lens (LJ1918L1-C, Thorlabs, Newton, NJ., USA) has a clear aperture of 6 × 4 mm2 and a focal lens of 5.8 mm. The focal plane has dimensions of 10 × 1 mm2 as shown on the right of Fig. 9(a). Figures 9(b) and (c) are the photographs for the probe with the visible aiming light and the 1460-nm laser on, respectively.
Figure 10 shows the 3D model and the photography of the handheld probe with the angled forceps. The probe head and the forceps are connected through an adaptor [Fig. 10(a)]. The line-shaped focal plane is just above the angled platform of the forceps so that the conjunctiva folds grasped by the angled platform can be effectively heated [Fig. 11(b)]. The adaptor also ensures the right angle illumination of the NIR laser light on the conjunctiva fold [Figs 4(b) and 10(a)] and a proper working distance. Figure 10(c) shows the actual probe and angled forceps held in hand. Figure 10(d) shows the configuration for laser delivery to redundant conjunctival tissue. A line-focused laser beam is delivered to the conjunctive fold held by a pair of angled forceps. The laser energy is applied in pulses that confine the peak heating to the conjunctival fold. Multiple pulses are delivered to achieve collagen shrinkage, which can be directly visualized by the surgeon.
Laser parameter selection
A pulsed laser will be used to heat and shrink the conjunctiva. The laser beam will be focused into a line on the tissue and the situation can be modeled as a one-dimensional heat diffusion problem for the calculation of tissue thermal relaxation time21:
where τ is the tissue thermal relaxation time, D is the heat diffusivity of tissue and d is the heat diffusive length. The heat diffusivity of tissue is approximately 1.3 × 10−7 m2 s−122. The thickness of the human conjunctiva is approximately 0.24 mm, and loose conjunctiva folded over when grasped by the surgical forceps should be approximately 0.5 mm thick. For this thickness, the thermal relaxation time is 0.5 s. So the pulse duration needs to be shorter than 0.5 s for the thermal energy to be confined to within this depth.
In addition, instead of using a single pulse, we propose to use the cumulative effect of multiple pulses. This strategy can not only reduce the requirement of peak power for effective shrinkage of the conjunctiva, but also keep the optimal temperature for a suitable period of time. The pulse train is emitted at a fixed repetition rate. This repetition rate can be estimated as following. For a pulse period including the pulse-tissue interaction and laser-off thermal relaxation, a thermal relaxation of 0.5 s indicating the pulse duration should be shorter than 0.5 s for confine the thermal energy in depth, while the laser-off duration should be around or a little bit longer than 0.5 s for effective heat dissipation of this pulse. So the pulse period should be 1 s, corresponding to a repetition rate of 1 Hz. The effectiveness of this estimation has been validated by the temperature variation in Section 3.1. We set a time duration of 4 seconds for the LTC process, which corresponds to 4 pulses at a repetition rate of 1 Hz.
Experiment setup and procedure
Figure 11 shows the setup and procedure of the experiment. Porcine eyes (Animal Technologies, Tyler, TX, USA) were used as samples in this LTC study. Frozen porcine eye balls were thawed in phosphate buffered saline before the experiment. To keep the condition of the tissue consistent, we used each position on the conjunctiva once. Because each porcine eye has a finite area of the conjunctiva and we were using a beam size of 10 × 1 mm2, the laser radiation on a single eye was limited to 6–8 times. Around 30 porcine eyes were used in this study. A 3D-printed marker was used to create a specific, observable region of interest that was denoted by two parallel lines on the conjunctiva surface: blue colors in the microscope images in Fig. 11(c) with a length of 10 mm and a separation of 1 mm. The line-focused light from the LTC probe was targeted to shine between the parallel marker lines on the conjunctiva. The angled forceps which grasp the tissue may introduce false shrinkage due to force applied on the conjunctiva. So, the angle forceps will not be used in the experiments for optimizing the laser parameters. A thermal camera (TiS45, Fluke, Everett, WA, USA) was used to record the temporal temperature change in the region during the LTC. OCT was also used to monitor the tissue structural change during the LTC. Figure 11(a) shows the experimental setup for LTC (bottom) and thermal camera image (top) of the porcine eye during the LTC. The thermal camera allows the record of the images during the LTC. A custom 1310-nm swept-source OCT system with an A-line rate of 50 kHz, a transverse resolution of 15 μm, and an axial resolution of 8.5 μm was used to monitor the LTC process. The detailed information about the OCT system has been described in a previous publication23. Figure 11(b) shows the schematic of the OCT experiment. The OCT B-scan covered a range of 5 mm and repeated B-scan at this same location was recorded. The OCT imaging rate is 50 frames per second. The laser-induced shrinkage was measured [Fig. 11(c)]. The photographs of the targeted area before and after the LTC were taken with a digital microscope (OMAX, Gyeonggi-do, Korea). The shrinkage percentage of the region was obtained by comparing the two images.
References
Meller, D. & Tseng, S. C. G. Conjunctivochalasis: Literature Review and Possible Pathophysiology. Survey of Ophthalmology 43(3), 225–232 (1998).
Liu, D. Conjunctivochalasis: a cause of tearing and its management. Ophthalmic Plastic & Reconstructive Surgery 2(1), 25–28 (1986).
Di Pascuale, M. A., Espana, E. M., Kawakita, T. & Tseng, S. C. G. Clinical characteristics of conjunctivochalasis with or without aqueous tear deficiency. British journal of Ophthalmology 88, 388–392 (2004).
Zhang, X. et al. Bulbar Conjunctival Thickness Measurements With Optical Coherence Tomography in Healthy Chinese Subjects. Investigative Ophthalmology & Visual Science 54, 4705–4709 (2013).
Zhang, Xingru et al. In Vivo Cross-Sectional Observation and Thickness Measurement of Bulbar Conjunctiva Using Optical Coherence Tomography. Investigative Ophthalmology & Visual Science 52, 7787–7791 (2011).
Gumus, K., Crockett, C. H. & Pflugfelder, S. C. Anterior segment optical coherence tomography: a diagnostic instrument for conjunctivochalasis. American Journal of Ophthalmology 150(6), 798–806 (2010).
Nakasato, S., Uemoto, R. & Mizuki, N. Thermocautery for Inferior Conjunctivochalasis. Cornea 31(5), 514–519 (2012).
Ohba, E. et al. Surgical punctal occlusion with a high heat-energy releasing cautery device for severe dry eye with recurrent punctal plug extrusion. American Journal of Ophthalmology 151(3), 483–7.e1 (2011).
Oh, S. J. & Byon, D. Treatment of Conjunctivochalasis Using Bipolar Cautery. J. Korean Ophthalmol Soc. 40(3), 707–711 (1999).
Doss, L. R., Doss, E. L. & Doss, R. P. Paste-pinch-cut conjunctivoplasty: subconjunctival fibrin sealant injection in the repair of conjunctivochalasis. Cornea 31(8), 959 (2012).
Santiago, E. et al. Surgical techniques for the treatment of conjunctivochalasis: paste-pinch-cut conjunctivoplasty versus thermal cautery conjunctivoplasty. Can J Ophthalmol 52, 308–312 (2017).
Nainiwal, S., Kumar, H. & Kumar, A. Laser conjunctivoplasty: a new technique for correction of mild medial ectropion. Orbit. 22(3), 199–201 (2003).
Naval, C. I. N., Cinco, P. R., Yutangco, R. A. & Dinglasan, J. L. Argon Laser Modified. Conjunctivoplasty for Bullous Keratopathy,European journal of Implant and Refractive Surgery 7(5), pp. 279–284 (1995).
Nicoli, S., Ferrari, G. & Quarta, M. Porcine sclera as a model of human sclera for in vitro transport experiments: histology, SEM, and comparative permeability. Mol Vis. 15, 259–266 (2009).
Kuenstnerb, J. Todd & Norris, KarlH. Spectrophotometry of Human Hemoglobin in the near Infrared Region from 1000 to 2500 nm. Journal of Near Infrared Spectroscopy 2(2), 59–65 (1995).
Horecker, L. The absorption spectra of hemoglobin and its derivatives in the visible and near infrared regions. J. Biol. Chem. 148, 173 (1943).
Curcio, JosephA. & Petty, CharlesC. The Near Infrared Absorption Spectrum of Liquid Water. J. Opt. Soc. Am. 41, 302–304 (1951).
Rigby, B. J. & Spikes, J. D. Hydroxyproline and the Shrinkage Temperature of Collagen. Nature 187, 150–151 (1960).
Adrian, R. J. Particle-imaging techniques for experimental fluid mechanics. Annual Review of Fluid Mechanics. 23(1), 261–304 (1991).
Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nature methods 9(7), 676–682 (2012).
Niemz, H. Laser-Tissue interactions: fundamentals and applications, Springer, 1996.
Telenkov, S. A., Youn, J.-I., Goodman, D. M., Welch, A. J. & Milner, T. E. Non-contact measurement of thermal diffusivity in tissue. Phys. Med. Biol. 46, 551 (2001).
Su, J. P. et al. Imaging the anterior eye with dynamic-focus swept-source optical coherence tomography. J. Biomed. Opt. 20(12), 126002 (2015).
Acknowledgements
We would like to thank Wen Chen from TXStar Laser Technology (Shanghai) Co., Ltd and Minzhao Zhang from Shanghai CONSUN Photoelectric Technology Co., Ltd for helping build the laser system and Chao Pang from Heraeus Noblelight for fabricating the handheld probe. National Institutes of Health (R01 EY023285, R01 EY018184); Oregon Health and Science University Foundation; Oregon Clinical & Translational Research Institute (Biomedical Innovation Program).
Author information
Authors and Affiliations
Contributions
G.L., D.H., J.Y., and S.C.P. designed the experiment. J.Y. and G.L. wrote the manuscript. J.Y., R.C., V.G., and T.B. conducted the experiments. G.L., D.H., and S.C.P. supervised the project. All the authors discussed the results and commented on the manuscript.
Corresponding author
Ethics declarations
Competing Interests
Oregon Health & Science University (OHSU), Baylor College of Medicine, David Huang, Gangjun Liu, Jianlong Yang, and Stephen C. Pflugfelder have potential financial interest in pending patent of the device described in the paper. These potential conflicts of interest have been reviewed and managed by OHSU.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Yang, J., Chandwani, R., Gopinatth, V. et al. Near-infrared laser thermal conjunctivoplasty. Sci Rep 8, 3863 (2018). https://doi.org/10.1038/s41598-018-22204-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-22204-0
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.