Near-infrared photons: a non-invasive probe for studying bone blood flow regulation in humans
© Binzoni and Spinelli; licensee BioMed Central Ltd. 2015
Received: 1 October 2014
Accepted: 22 June 2015
Published: 25 July 2015
The study of bone blood flow regulation in humans has always represented a difficult task for the clinician and the researcher. Classical measurement techniques imply the presence of ionizing radiation or contrast agents, or they are slow or cannot be repeated too often in time. In the present review, we would like to give a perspective on how the optical approach might overcome some of these problems and give unique solutions to the study of bone blood flow regulation. We hope that the present contribution will encourage the scientific community to put a greater attention on this approach.
KeywordsBone blood flow regulation Human Bone neurovascular system Blood volume regulation Near infrared light Biophotonics
The study of bone blood flow regulation in humans has always been a challenging topic. In fact, known techniques allowing the monitoring of blood flow appear to be invasive (e.g., utilization of ionizing radiation or contrast agents) or do not permit fast measurements (e.g., to follow blood flow pulsations) over repeated or long periods of time. Moreover, these techniques are often very expensive and cannot be used at the patient bedside [1–4]. In this frame, optical technologies, based on near infrared light, might represent a unique solution allowing to surmount some of the these difficulties and open the possibility to approach from a new angle the study of blood flow regulation in bone.
Light in the 650- to 950-nm range has the singular property to deeply penetrate biological tissues such as brain, muscle, fat, or bone. Moreover, due to the phenomenon of diffusion, light propagates in all directions, and then, it is re-emitted also from the same surface where it is injected. Photons composing this light interact in various ways with living tissues, and when detected, they carry information about blood flow and other parameters related to it. In the following sections, we will explain, by reporting the related scientific literature, how non-invasive and contrast-agent free optical techniques can be exploited to investigate the regulatory mechanisms of the bone neurovascular system in humans.
The optical instrumentation
Optical instruments utilized to investigate bone blood flow in humans
Possible source types
Possible number of light wavelengths
Light source intensity
Detected parameters for light (possible detector types)
Measured physiological parameters
Intensity (PMT, APD)
LED, CW laser, white light
1, 2, 3, …, to continuous
Intensity (PMT, APD)
Hb, HbO2, ΔHb, ΔHbO2, (a.u.) and %SO2
Temporal correlation (SPD)
BF (cm2 s −1)
2, 3, …
Intensity and intensity phase shifts (PMT, APD)
Hb, HbO2, (μM) and %SO2
2, 3, …
Photon DTOF (SPD)
Hb, HbO2, (μM) and %SO2
Doppler frequency shift (PMT, APD, SPD)
BF, BS, #rbc (a.u.)
The aim of the present contribution is to show to a non-technical reader the potentialities of the different optical instrumentation when applied to the study of human bone blood flow. For this reason, we will not discuss technical or theoretical advantages/problems related to the different approaches, because this is out of scope. The reader interested on a more technical side of this presentation can refer to specialized literature [5, 6].
Investigations on bone blood flow by the optical approach
For this section, an organization per instrument has been conceived, emphasizing what the different optical techniques can actually measure, in the framework of human bone blood flow investigations. As a matter of fact, the specific parameters of interest for the investigator directly depend on the type of optical instrument. We think that this approach is closer to the real needs of the experimenter or the clinicians that first of all wants to know what it is possible or not to measure with the optical approach, without being lost in technical details or often confusing chronological history.
Near infrared spectroscopy
Near infrared spectroscopy (NIRS) allows the assessment of oxy-(HbO2) and deoxy-hemoglobin (Hb) concentration changes (ΔHbO2 and ΔHb, respectively). This is made possible by the fact that light is absorbed differently by HbO2 compared to Hb. Some clever mathematical algorithm, exploiting this phenomenon, allows one to obtain ΔHbO2 and ΔHb as instrument output.
An early attempt to monitor non-invasively bone blood perfusion, by optical techniques, has been performed by NIRS in human tibia . The aim was to define an indirect perfusion index allowing future investigations on blood flow regulation in special environmental conditions such as in microgravity. The rates of deoxy-hemoglobin concentration changes were utilized as an index allowing to evaluate the reperfusion speed after a short ischemic period. It was shown (13 males) that the return to preischemic conditions is faster for the tibialis muscle than for the tibia (from ∼1.8 times at 72 years old to ∼3.5 times at 25 years old) and that this kinetic becomes slower with aging in both tissues (∼2 times for tibia from 25 to 72 year old). In a subsequent NIRS study , these results have been confirmed (13 males) and it was also shown that, in tibia bone, resting blood oxygen saturation (%SO2) remains constant (∼85 %) with aging, while in tibialis muscle, it decreases (−0.35 % per year with an %SO2 of ∼80 % at 25 years old). This gives us some inside on the age dependence of bone vascular system and related oxydative metabolism, information that cannot be easily obtained by other classical measurement techniques. Always using NIRS, it has been shown that, as expected, tibia oxygen consumption does not change during submaximal (80 %) isometric contractions (foot dorsiflexion, 10 males and 5 females) . On the other side, despite the rigid bone structure, NIRS has allowed to demonstrate that orthostatic stress (tilt bed) may induce blood volume changes in human tibia bone/bone marrow  (6 males). This phenomenon is probably related to a concurrent extravascular fluid displacement in the tibia and a typical half time of 1 min (estimated from Figure two in ) is necessary to reach a ΔHb + ΔHbO2 steady state level after tilting. This value is independent of the magnitude of the rapid increase in the orthostatic pressure. In the same study, a sudden change (at ∼28 mmHg) in the slopes of ΔHb and ΔHbO2, as function of the orthostatic pressure, was also revealed. This was interpreted as the consequence of the possible presence of a direct and/or indirect venous-arteriolar reflex. Since then, NIRS instruments specifically dedicated to bone investigation have been developed . This has allowed to better define the region of interest monitored by the NIRS instruments and, as a consequence, by any other optical instrument presented here (thanks to the similarity in the optical geometries; see, e.g., Fig. 1).
Attempts to apply NIRS in clinical settings have been made on spinal cord injured (SCI, five males and five females) patients . Even if it does not represent an extensive study, it has already been possible to observe, by analyzing the ΔHb and kinetics, that normal subjects (five males and five females) have faster post-ischemic hyperhemic response in tibia (baseline-to-peak Hb in 10.1 s in normal compared to 16.6 s in SCI) and faster return to baseline compared to SCI . The result may be explained by an impaired blood flow in SCI.
Time resolved spectroscopy
The NIRS instrumentation considered in the previous section does not give access to the absolute values of Hb and HbO2 concentrations but only to their changes, often expressed in arbitrary units. This prevents the clinician to easily compare different patients. For this reason, a dedicated time resolved spectroscopy (TRS) system giving Hb, HbO2 in absolute values has been developed . In TRS, by exploiting very short laser pulses (durations of tens of picoseconds), it is possible to measure the distribution of time-of-flight of photons arrived to the detector, after having traveled through the tissue. From these data, by means of a dedicated algorithm, one can obtain Hb and HbO2. By using this system, it has been shown in human calcaneous (range 25–72 years) that %SO2 remains constant (∼70 %) for increasing age. This is in accordance with previous NIRS findings (see above) in human tibia  but at a lower %SO2 value. The lower %SO2 may be explained by the different vascular structure of the two bones determining different blood flows (%SO2 may depend, e.g., on blood speed) and different Hb and HbO2 fractions. In this frame, values for Hb+HbO2 were found to be constant (∼23 μM) in the same age range, implying that, in calcaneous, the vascular volume is independent from age.
It would be interesting to know if in general the vascular system of aging bone in healthy subjects undergoes blood flow regulatory changes. Considering the possible link existing for example between impaired bone blood flow and osteoporosis , this knowledge may reveal important for the comprehension of the mechanisms of initiation of this pathology. Optical techniques might be in this case a non-invasive and attractive early screening tool.
Photoplethysmography (PPG) has been introduced with the aim to have a more direct mean to assess blood flow, and the fundamental improvement has been the capability to obtain mean blood flow (BF) changes (ΔBF) in arbitrary units (a.u.). PPG technique is based on the fact that total hemoglobin concentration variations also changes light tissue absorption properties. By modeling this phenomenon, it is possible to obtain ΔBF. Compared to NIRS, it seems that PPG is also sensitive to blood flow speed, and for this reason, PPG should generate a signal that better represents blood flow changes. However, this potential mechanism has not yet been completely clarified and care must be taken on this point. Thanks to PPG it has been possible to monitor fast BF changes induced by heart pulsations in human patella  (12 males and 8 females). In a subsequent study, the same authors were able to demonstrate that the amplitude of the pulsatile BF (ΔBF, peak-to-peak) in patella is markedly reduced (∼[−26 %) in patient (9 males and 13 females) with patellofemoral pain syndrome when the knee is being flexed (from 20° to 90°), while in normal subjects (18 males and 15 females), ΔBF show no distinctive patterns . This corroborate the hypothesis of a potential ischemic mechanism involved in this pathology.
Other PPG experimentation have demonstrated that in healthy humans (five males and four females), local cooling (5-min ice pack on the tibia) decreases the amplitude of the tibial ΔBF pulsations (∼]−33 %; from figure 3 in ). This might be a possible evidence of the presence of a sympathetic reflex in bone. It has also been observed that increasing local external pressure (from −50 to 30 mmHg) increases ΔBF pulsations amplitudes (from ∼[−3.2 to ∼1.4 a.u.) in the tibia of healthy subjects (six males and six females), with a sudden ΔBF decrease (∼]−0.5 a.u.) for pressure values exceeding 30 mmHg (from figure 3 in ). The exact reason of this behavior has not yet been clarified. To finish, it was observed always in healthy subjects that post-exercise increase in ΔBF pulsations in patellar bone (“hyperhemia”) depends on the intensity of the quadriceps femoris contraction  (28 males and 14 females). This resulted in an increase of 22 and 61 % in ΔBF for an isotonic exercise intensity of 60 and 75 % of one repetition maximum, respectively (duration of high was half of that of the low intensity exercise).
Intensity-modulated spectroscopy (IMS; also called frequency domain spectrometer) is based on the modulation of the amplitude of a CW laser light source. The analysis of the changes in the modulation amplitude at the output of the observed tissue region of interest allows one to estimate Hb, HbO2, and %SO2. Blood flow pulsations in patella has been indirectly monitored by IMS by observing Hb, HbO2, and %SO2  variations. In this case, the authors have observed a periodic increase in Hb+HbO2 and %SO2 at each heartbeat of ∼3.1and ∼0.7 %, respectively (8 males). This result may be explained, e.g., by an increase in tissue blood volume mainly determined by the entry of the arterial blood in the tissue.
Diffuse correlation spectroscopy
All the instrumentation cited in the previous paragraphs try to derive information on blood flow through indirect measurements or through BF changes. For this reason, diffuse correlation spectroscopy (DCS) and laser-Doppler flowmetry (LDF) have been developed to obtain measurements directly proportional to BF. Note that DCS and LDF do not allow to obtain Hb, HbO2, or %SO2. DCS is based on the fact that photons traveling through the investigated tissue undergo a frequency Doppler-shift when interacting with moving red blood cells. By analyzing the photo-electric current produced by these photons at the output of the photodetector (through an autocorrelation operation), one obtains BF. BF obtained by DCS is expressed in special cm2 s −1 units. It has been demonstrated that the measured parameter is proportional to BF. In principle, this allows to obtain comparable values between different subjects. In fact, DCS has permitted to define a “normal” BF reference value (5.0×10−9 cm2 s −1; range 4.2×10−9 cm2 s −1 to 7.4×10−9 cm2 s −1) for the human (15 males and 17 females) manubrium ; the aim of the authors is to compare in the future this BF value with that of patients with cancers.
LDF directly assess BF in arbitrary units, but the supplementary advantage is that it allows to monitor two important parameters determining BF, i.e., mean blood speed (BS, in a.u.) and the number moving red blood cells (#rbc, in a.u.) . As DCS, also LDF is based on the fact that photons traveling through the investigated tissue undergo a frequency Doppler-shift when interacting with moving red blood cells. By analyzing the photo-electric current produced by these photons at the output of the photodetector (through the power density spectrum assessment) one obtains BF, BS, and #rbc. Thanks to this particularity, it has been shown using ECG-gated LDF that heart-generated BF pulsations in patella, clavicle, tibial malleolus, and tibial diaphysis are determined by BS changes at constant #rbc  (5 males). This signifies that fast BF pulsations do not mechanically induce in this case fast blood volumes changes. ECG-gated LDF has allowed also to observe a secondary peak, in the BF pulsation, probably due to the closure of the cardiac valves. To improve the sensitivity of the LDF technique and to reach deeper tissue regions of interest, a new instrument based on single-photon counting has been developed . This low-cost LDF allows to follow in real time BF kinetics in human tibia, during an ischemia-reperfusion protocol, with a time resolution of ∼1.6 s.
Perspectives and conclusions
Compared to sophisticated imaging techniques such as MRI, PET, CT, etc, the use of NIRS, TRS, PPG, IMS, DCS, and LDF may seem very simplified approaches to the study of bone blood flow. Actually, this is a misleading point of view. Optical techniques have unique characteristics that might allow original investigations on bone blood flow regulation in humans. In fact, we must not forget that these techniques are portable, have imaging capabilities, and can be used in special environments such as during water immersion , in flight simulators , in hyperbaric chambers [27, 28], during high altitude hiking , in microgravity, or during parabolic flights [29, 30]. Moreover, the non-invasiveness of the approach allows us to study the bone neurovascular system without the unwanted influence of, e.g., the mental stress that might modulate the activity of the autonomous nervous system and thus potentially bias the experimental results . The large spectrum of sampling rates (from fraction of second to seconds, depending on the specific implementation of a given hardware), acquisition times (from a fraction of second to hours), and the repeatability of the measurements, without incurring in potential risks determined by ionizing radiation and/or contrast agents allow exploiting optical techniques for both functional studies and monitoring purposes, enabling the conception of new investigations considered impossible until now for ethical reasons. The follow-up of acute or long term effects induced by pharmaceutical products modulating bone blood flow might also represent an interesting application.
It is difficult at this stage to decide which is the better optical technique. This question has not yet been completely answered by the scientific community. Probably, the best choice still depends on the precise experimental goal. For example, if very fast acquisition rates are needed, and only changes in BF values suffice, probably the PPG technique is superior compared to LDF or DCS. If absolute (a.u.) BF values are required, then LDF or DCS may represent the best choice. However, depending on the future technical improvements in the hardware components composing these instruments, this situation may rapidly change.
To shortly highlight some pros and cons of the proposed modalities in terms of their capability to assess bone blood flow, one can say for example that LDF or DCS are the only techniques allowing to directly obtain this parameter. PPG has the advantage of the high acquisition speed, but it can detect only changes in blood flow, and the mechanism of detection has not yet been fully clarified (living an open question on the linearity of the instrument response). NIRS, TRS, and IMS give only indirect information on blood flow, but they have the advantage to give other important physiological parameters (linked to blood flow) such as %SO2. However, in terms of explicit bone blood flow monitoring, NIRS, TRS, and IMS are probably the less suitable approaches.
Thus, besides the obvious potential utilization of the optical techniques for clinical purposes, the aim of the present review is to encourage orthopedists and scientist to exploit this new possibility to deepen the knowledge on bone blood flow regulation mechanisms in humans. A topic that certainly deserves the attention of the scientific and medical community.
In conclusion, the large spectrum of applications in the domain of human bone blood flow regulation, made possible by the advent of the optical techniques, potentially opens a new domain of study in physiological anthropology. In fact, nutrition, fitness, aging, sex, and growing are all parameters that may influence human bone blood flow. The access to low-cost, portable, and non-invasive optical techniques may represent in this sense a new powerful tool allowing large scale screenings and a better comprehension of long/short term bone blood flow regulation and its adaptability to human modern environment.
- Brookes M, Revell WJ. Blood supply of bone. Scientific Aspects. London: Springer; 1998.View ArticleGoogle Scholar
- Dyke JP, Aaron RK. Noninvasive methods of measuring bone blood perfusion. Ann New York Acad Sci. 2010; 1192:95–102.View ArticleGoogle Scholar
- Schoutens A, Arlet J, Gardeniers JWM, Hughes SPF. Bone circulation and vascularization in normal and pathological conditions. New York: Plenum Press; 1993.View ArticleGoogle Scholar
- Lan SM, Wu YN, Wu PC, Sun CK, Shieh DB, Lin RM. Advances in noninvasive functional imaging of bone. Acad Radiol. 2014; 21:281–301.PubMedView ArticleGoogle Scholar
- Boas DA, Pitris C, Ramanujam N. Handbook of biomedical optics. Boca Raton, London, New York: CRC Press, Taylor and Francis Group; 2011.View ArticleGoogle Scholar
- Delpy DT, Cope M. Quantification in tissue near-infrared spectroscopy. Philos Trans R Soc Lond B Biol Sci. 1997; 352:649–59.PubMed CentralView ArticleGoogle Scholar
- Binzoni T, Blanchi S, Fasel JH, Bounameaux H, Hiltbrand E, Delpy D. Human tibia bone marrow blood perfusion by non-invasive near infrared spectroscopy: a new tool for studies on microgravity. Philos Trans R Soc Lond B Biol Sci. 2002; 9:183–4.Google Scholar
- Binzoni T, Leung T, Hollis V, Bianchi S, Fasel JH, Bounameaux H, et al. Human tibia bone marrow: defining a model for the study of haemodynamics as a function of age by near infrared spectroscopy. J Physiol Anthropol Appl Human Sci. 2003; 22:211–8.PubMedView ArticleGoogle Scholar
- Klasing M, Zange J. In vivo quantitative near- infrared spectroscopy in skeletal muscle and bone during rest and isometric exercise. Proc SPIE. 2003; 5138:318–22.View ArticleGoogle Scholar
- Binzoni T, Leung TS, Courvoisier C, Giust R, Tribillon G, Gharbi T, et al. Blood volume and haemoglobin oxygen content changes in human bone marrow during orthostatic stress. J Physiol Anthropol. 2006; 25:1–6.PubMedView ArticleGoogle Scholar
- Aziz SM, Khambatta F, Vaithianathan T, Thomas JC, Clark JM, Marshall R. A near infrared instrument to monitor relative hemoglobin concentrations of human bone tissue in vitro and in vivo. Rev Sci Instrum. 2010; 81:043111.PubMedView ArticleGoogle Scholar
- McCarthy I. Application of near infrared spectroscopy in the assessment of bone perfusion. J Bone Joint Surg Br. 2012; 94-B(Supp VIII):36.Google Scholar
- Khakha RS, Bloomer Z, Bain D, Nicholson G, Gall A, Ferguson-Pell M. Differences in the reactive hyperhemia response of bone in able bodied and spinal injured persons using near infrared spectroscopy. J Bone Joint Surg Br. 2006; 88-B(Supp II):305.Google Scholar
- Pifferi A, Torricelli A, Taroni P, Bassi A, Chikoidze E, Giambattistelli E, et al. Optical biopsy of bone tissue: a step toward the diagnosis of bone pathologies. J Biomed Opt. 2004; 9:474–80.PubMedView ArticleGoogle Scholar
- Näslund J, Pettersson J, Lundeberg T, Linnarsson D, Lindberg LG. Non-invasive continuous estimation of blood flow changes in human patellar bone. Med Biol Eng Comput. 2006; 44:501–9.PubMedView ArticleGoogle Scholar
- Näslund J, Waldén M, Lindberg LG. Decreased pulsatile blood flow in the patella in patellofemoral pain syndrome. Am J Sports Med. 2007; 35:1668–73.PubMedView ArticleGoogle Scholar
- Mateus J, Hargens AR. Photoplethysmography for non-invasive in vivo measurement of bone hemodynamics. Physiol Meas. 2012; 33:1027–42.PubMedView ArticleGoogle Scholar
- Mateus J, Hargens AR. Bone hemodynamic responses to changes in external pressure. Bone. 2013; 52:604–10.PubMedView ArticleGoogle Scholar
- Näslund JE, Näslund S, Lundeberg E, Lindberg LG, Lund I. Bone blood flow is influenced by muscle contractions. J Biomed Sci Eng. 2011; 4:490–6.View ArticleGoogle Scholar
- Farzam P, Zirak P, Binzoni T, Durduran T. Pulsatile and steady-state hemodynamics of the human patella bone by diffuse optical spectroscopy. Physiol Meas. 2013; 34:839–57.PubMedView ArticleGoogle Scholar
- Farzam P, Lindner C, Weigel U, Suarez M, Urbano-Ispizua A, Durduran T. Noninvasive characterization of the healthy human manubrium using diffuse optical spectroscopies. Physiol Meas. 2014; 35:1469–91.PubMedView ArticleGoogle Scholar
- Binzoni T, Boggett D, Van De Ville D. Laser-doppler flowmetry at large interoptode spacing in human tibia diaphysis: Monte carlo simulations and preliminary experimental results. Physiol Meas. 2011; 32:33–53.View ArticleGoogle Scholar
- Binzoni T, Tchernin D, Hyacinthe JN, Van De Ville D, Richiardi J. Pulsatile blood flow in human bone assessed by laser-doppler flowmetry and the interpretation of photoplethysmographic signals. Physiol Meas. 2013; 34:25–40.View ArticleGoogle Scholar
- Binzoni T, Van De Ville D, Sanguineti B. Time domain algorithm for single-photon laser-doppler flowmetry at large interoptode spacing in human bone. Appl Opt. 2014; 53:7017–24.PubMedView ArticleGoogle Scholar
- Binzoni T, Tchernin D, Richiardi J, Van De Ville D, Hyacinthe JN. Haemodynamic responses to temperature changes of human skeletal muscle studied by laser-doppler flowmetry. Physiol Meas. 2012; 33:1181–97.PubMedView ArticleGoogle Scholar
- Harrison M, Neary JP, Albert WJ, Veillette DW, McKenzie NP, Croll JC. Trapezius muscle metabolism measured with nirs in helicopter pilots flying a simulator. Aviat Space Environ Med. 2007; 78:110–6.PubMedView ArticleGoogle Scholar
- Kot J. Medical equipment for multiplace hyperbaric chambers. part i: Devices for monitoring and cardiac support. Eur J Underwater Hyperb Med. 2005; 6:115–20.Google Scholar
- Litscher G, Schwarz G, Ratzenhofer-Komenda B, Kovac H, Gabor S, Smolle-Jüttner FM. Transcranial cerebral oximetry in the hyperbaric environment. Biomed Tech (Berl). 1997; 42:38–41.View ArticleGoogle Scholar
- Zhang Q, Yan X, Strangman GE. Development of motion resistant instrumentation for ambulatory near-infrared spectroscopy. J Biomed Opt. 2011; 16:087008.PubMed CentralPubMedView ArticleGoogle Scholar
- Schneider S, Abeln V, Askew CD, Vogt T, Hoffmann U, Denise P, et al. Changes in cerebral oxygenation during parabolic flight. Eur J Appl Physiol. 2013; 113:1617–23.PubMedView ArticleGoogle Scholar
- Binzoni T, Van De Ville D. Noninvasive probing of the neurovascular system in human bone/bone marrow using near-infrared light. Innovative Opt Health Sci. 2011; 4:183–9.View ArticleGoogle Scholar
- Allen J. Photoplethysmography and its application in clinical physiological measurement. Physiol Meas. 2007; 28:1–39.View ArticleGoogle Scholar
- Scholkmann F, Kleiser S, Metz AJ, Zimmermann R, Mata Pavia J, Wolf U, Wolf M. A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology. Physiol Meas. 2014; 85(Pt 1):6–27.Google Scholar
- Durduran T, Yodh AG. Diffuse correlation spectroscopy for non-invasive, micro-vascular cerebral blood flow measurement. Neuroimage. 2014; 85(Pt 1):51–63.PubMed CentralPubMedView ArticleGoogle Scholar
- Durduran T, Choe R, Baker W, Yodh AG. Diffuse optics for tissue monitoring and tomography. Rep Prog Phys. 2010; 73:076701–44.PubMed CentralPubMedView ArticleGoogle Scholar
- Torricelli A, Contini D, Pifferi A, Caffini M, Re R, Zucchelli L, Spinelli L. Time domain functional nirs imaging for human brain mapping. Neuroimage. 2014; 85(Pt 1):28–50.PubMedView ArticleGoogle Scholar
- Humeau A, Steenbergen W, Nilsson H, Strömberg T. Laser doppler perfusion monitoring and imaging: novel approaches. Med Biol Eng Comput. 2007; 45:421–35.PubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License(http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.