1.IntroductionDiabetes mellitus is a metabolic disease whereby the human body cannot properly control its levels of glucose.1 Past research has established a clear understanding of how the normal body controls its glucose.1 Currently researchers are focused on methods of monitoring the glucose levels, both invasively and noninvasively. Ultimately, a glucose monitoring system could be used in the design of an artificial control system. The Diabetes Control and Complications Trial Research Group (DCCTRG) conducted a study between 1983 and 1989 to determine if the severity and frequency of complications related to insulin-dependent diabetes mellitus could be decreased with tight monitoring of glucose levels.2 By performing at least four self-monitoring blood glucose treatments a day, the DCCTRG conclusively found a 76 reduction in the development of retinopathy, a 54 reduction in the progression of retinopathy, a 40–50 reduction of nephropathy, and a 60 reduction of clinical neuropathy.2 Therefore, tight glucose monitoring is essential for the reduction of diabetic complications. Since Clark's oxygen electrode was first used in the 1950's to detect glucose,3 many attempts have been made to monitor glucose levels, ranging from other enzymatic techniques4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 to optical techniques.25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 The use of the glucose oxidase enzyme has been well documented and can provide reliable short-term monitoring.3–24 However, as a potential implantable device, the lifetime of the device is too short, and the electrodes are difficult to consistently reproduce. This is illustrated by several recent studies wherein the reliability of several enzymatic devices ranged from only 24 h to 11 days.17 20 21 Cote´ summarized several optical techniques for glucose monitoring in a recent paper.25 Absorption spectroscopy techniques, both in near infrared (NIR) and the middle infrared (MIR), were discussed. There are relative advantages and disadvantages to each frequency range. Water absorbs more in the MIR, but the peaks are less distinct in the NIR. Thus, multivariate analysis is often used with NIR because of peak smearing and absorption interference from other blood constituents.25 Also, the path length must be smaller for MIR because of the higher water absorption. While Heise et al.; first proposed multivariate analysis for MIR, most of the recent research is focused on NIR.30 31 Polarimetry, Raman spectroscopy, and fluorescence techniques for glucose monitoring were also discussed.25 Optical methods are advantageous as there is no dependence on chemical reactions or enzymatic lifetimes.5 34 Also, a fluorescent molecule relating specifically to glucose is not required.25 32 Various methods explored for potential implantation have included reflectance, light scattering, and transmission.23 The approach under investigation in this paper is transmission within the MIR range. The variables that exist with the light transmission approach are: (1) the concentration of the glucose within the sample—of particular interest for practical use, and (2) the path length, or the distance, that the light travels through the sample—a determining factor for the power requirement of the source and for the interference due to the water absorption peak. The concentration of other materials within the sample will affect the level of light transmission at certain frequencies as well. The proper frequencies for the determination of glucose were investigated by researching previous literature and performing new experiments.33 34 36 Zeller, Novak, and Landgraf, Kajiwara et al.;, and Back and Polavarapu have identified several characteristic frequencies in the MIR glucose absorption spectra.33 34 36 Specific frequencies identified in blood by Zeller and co-workers using a demountable liquid cell with a 25 μm path length were 1365, 1152, 1109, 1080, and 1035 cm −1 . Only the 1035 cm −1 peak was limited to glucose in the context of other blood constituents. For example, the 1365 cm −1 peak was identified in albumin, hemoglobin, and other constituents; the 1080 cm −1 peak was identified in albumin and hemoglobin; and the 1109 cm −1 peak was identified in hemoglobin. However, if interstitial fluid is used for glucose detection, thus eliminating the effects of hemoglobin, then the 1109 cm −1 peak may be used. In addition, the presence of large peaks near a characteristic peak may mask its appearance. For example, a strong urea peak at 1160 cm −1 masked the characteristic glucose peak at 1152 cm −1 . While albumin will still cause added absorption at the 1080 cm −1 peak, the frequency is potentially usable because changes in albumin concentration have little effect on the absorption spectra.34 Therefore, more than one peak could determine the glucose concentration. However, Zeller, Novak, and Landgraf used a range of glucose concentrations that were mostly higher than physiological levels, and even then the spectra had to be magnified 20-fold in order to resolve the differences at some of the characteristic peaks. The apparent frequency of the peaks will shift slightly depending on the resolution used and the exact temperature of the sample according to Lambert-Beer's law.34 While there has been discussion of a decrease in water absorption when samples are taken at body temperature (37 °C) versus room temperature (25 °C), the effects are expected to be minimal. There are several sample media that have been explored for the determination of glucose within the human body. For noninvasive external monitoring techniques, blood and tissue have been investigated. Other semi-invasive external techniques draw out interstitial fluid for measurement. The technique described in this paper is designed for use with extracellular fluid, particularly interstitial fluid, whether it is eventually implanted in the peritoneal cavity or used externally. Interstitial fluid closely resembles plasma with similar glucose values.1 The only distinct difference between plasma and interstitial fluid is the decrease in protein in interstitial fluid, about sixfold.1 Plasma and serum are similar as well except that plasma contains fibrinogen, whereas serum does not. Thus, serum was an appropriate substitute for interstitial fluid for experimental purposes. Also, Wientjes and Schoonen have demonstrated that the measurement of glucose in interstitial fluid and blood are similar.24 According to Guyton and Hall, the glucose levels of interstitial fluid are similar to those of plasma.1 However, Cote´ suggested that measurement of glucose in interstitial fluid is several orders of magnitude less than in plasma.25 Guyton and Hall had also noted that intracellular fluid contains a significantly lower level of glucose than the extracellular fluid.1 This may be the reason for the contrary opinion. The primary objective of this study was to investigate the use of two frequencies of light within the middle infrared, one detecting the base line water absorption and the second detecting glucose, to determine the glucose levels in human blood serum samples. While other multivariate approaches have been used, it was hypothesized that a calibration curve using only two frequencies can be constructed. Since other constituents in the interstitial fluid can have an effect on the absorption spectra, a secondary objective was to determine the effects of high and low concentrations of urea and albumin (protein) on glucose absorption. The quantitative relationship between physiological glucose concentrations and the absorption spectra will be demonstrated. 2Materials and MethodsGlucose concentrations ranging from 0 to 450 mg/dl in 50 mg/dl intervals were obtained by the addition of anhydrous glucose (Fisher Scientific) to glucose depleted base human serum (ABT, Sequin, TX). The glucose depleted base human serum assay is documented in Table 1. Table 1
The two serum constituents with the most prominent absorption near the glucose peaks of interest are urea and albumin. Albumin is the main protein in serum.1 Urea and albumin (Sigma-Aldrich) were added to the serum, the samples were mixed thoroughly, and the samples were analyzed (UAB Clinical Outreach Lab) to confirm the amount of urea and albumin in the solution. Samples prepared with urea and albumin are identified in Table 2. The spectra for high and low glucose concentration samples with high and low levels of urea or albumin or both will be compared to control glucose concentration samples. Table 2
Absorption spectra from 1250 to 1000 cm −1 were collected at 25 °C using a Mattson Research Series 1000 Fourier transform infrared (FTIR) spectrometer (Mattson Tech. Inc., Fremont, CA). While a slightly higher absorption difference may be noticed when samples are run at 37 °C (body temperature), the focus of this research was to demonstrate the technique at room temperature. The spectra were collected at a 0.75 cm −1 resolution over 128 sample scans. The spectrometer utilized a liquid nitrogen cooled mercury–cadmium–telluride detector. The transmission cells utilized CaF 2 crystals for windows, as CaF 2 is resistant to aqueous solutions, acids, and alkalis and, more importantly, translucent in the spectral range of interest (Specac, Inc., Smyrna, GA). As the path length of the transmission cell increases, the absorption due to the materials within the cell will also increase. The strongest absorbers of light within this region are water and glucose. Eventually, water absorption peaks will completely mask all resolvable glucose peaks if the path length is continuously increased. A compromise must be reached whereby a path length large enough to allow for resolvable differences between physiological concentrations of glucose and small enough to observe a minimum level of glucose (in this case 5 mg/dl) is reached. There must also be a significant difference between the absorption at 5 mg/dl and the absorption at 450 mg/dl, so that there are enough resolution points for quantitative determination. Only 6 and 12 μm tin spacers were available for this liquid cell (Specac, Inc., Smyrna, GA). Use of up to three 12 μm spacers resulted in a minimum detection limit of 5 mg/dl. The addition of a 6 μm spacer or a 12 μm spacer resulted in a minimum detection limit of 50 and 100 mg/dl, respectively, due to increased water absorption. Thus, three 12 μm spacers were used to construct the liquid cell. Spectra of the empty cell was obtained to validate the path length. By utilizing the fringe patterns resulting from constructive and destructive interference, the nominal path length was calculated according to Eq. (1) where f1 and f2 are the frequencies (cm −1 ) of the two maxima from interference fringes observed in the spectra of the empty cell, n is the number of maxima between f1 and f2, and d is the path length.30The path length of the empty cell was determined to be 39 μm. Due to the tolerance value for the thickness of the spacers, the path length was slightly larger than the expected 36 μm. The nominal path length was used to calculate the K coefficient of absorption according to Eq. (2) where A is the absorption value according to the spectra taken, and d is the calculated path length.Following the determination of the path length, spectra were collected for each of the samples. Three data points were collected for each of the human blood serum samples at the specific concentration intervals. Between each spectra, the cell was purged with distilled water several times in order to remove any possible residue from the previous test. Raw absorption values were recorded at 1035, 1080, 1109, and 1193 cm −1 . The 1193 cm −1 absorption point had no observable change across all spectra of glucose concentrations nor appeared to be related to any other constituent, so it was used as a base line determination of water absorption since the base line of a spectra can vary due to the spectrometer. All raw absorption glucose values were subtracted by the value at 1193 cm −1 . This difference absorption value, which is related to the relative glucose change, was then converted into a K coefficient of absorption value with Eq. (2). Statistical analysis was performed using StatView software (Cary, NC). Linear regression analysis provided the R2 values. P values were calculated for the first, second, third, and fourth order equations. Residual plots were also obtained to determine if the data exhibited randomness about the fitted equations. The estimated error in the regression analysis was determined from the root mean squared residual error values. 3.ResultsFigure 1 presents a comparison between spectra at 5 and 450 mg/dl. All water absorption has been subtracted from each spectra to clearly illustrate absorption due to serum constituents, specifically glucose. A relative increase in absorption due to glucose changes is noticeable at the three frequencies of interest: 1035, 1080, and 1109 cm −1 . The K coefficient of absorption difference, or the absorption due to glucose, is plotted versus the concentration of glucose in Figure 2. Each concentration point is the average of three data points with each error bar representing one standard deviation. The statistical analysis of the data is presented in Table 3. The R2 values suggest potential linearity of the data presented in Figure 2. P values for the coefficients of the first order term were less than 0.05, thus indicating a linear relationship. P values for the coefficients of the second order equation were greater than 0.05 for the 1035 and 1080 cm −1 frequencies indicating that the second order term can be neglected. Similarly, the P values for the third and fourth order equations were greater than 0.4 indicating the higher order terms could also be neglected. The residual plots conducted for the first order fitted lines exhibited a randomness of the data. The estimated error of prediction for glucose concentrations obtained using the 1035 cm −1 frequency was significantly lower than the estimated error for the other two frequencies. Table 3
Urea and albumin changes did not result in notable absorption changes in the high or low glucose spectra. The predicted glucose values that were determined for each of the spectra with varying urea and albumin values were within the estimated error presented in Table 3 for each potential frequency to be used for glucose determination (1035, 1080, and 1109 cm −1 ). The values were random about the clinically determined value and did not exhibit any relationship due to an increase in either urea or albumin. 4.DiscussionA quantitative relationship between the glucose concentrations and the K coefficient of absorption values was determined (Figure 2). The relative change in absorption between high and low glucose values within the physiological region, when water is subtracted, is clearly observed (Figure 1). While there is a larger difference observed at 1109 and 1080 cm −1 , the 1035 cm −1 frequency is more quantitative. A possible cause is that water is a stronger absorber near the 1109 and 1080 cm −1 frequencies and thus causes less discernible results. The 20.6 mg/dl (1.1 mmol/L) error at the 1035 cm −1 frequency is comparable to other research. Cote´ summarized that multivariate NIR research produced error values of 1.67–2.78 mmol/L (30–50 mg/dl).25 Also, high and low concentrations of urea and albumin did not affect the quantitative determination of glucose in serum at 5 or 500 mg/dl concentrations. While Zeller and co-workers demonstrated a qualitative ability to detect glucose using transmission techniques, this study demonstrated the ability to quantitatively determine glucose in serum.33 Most of the concentrations used by Zeller and co-workers were higher than normal levels and revealed the apparent saturation of absorption due to glucose. The results from this study exhibited no saturation within physiological levels and suggested a linear relationship. However, the quantitative ability of the technique is most important, regardless of linearity. The next step will be to transition from the FTIR system to a breadboard design consisting of frequency specific lasers and detectors. The calculated cross section of absorption for the glucose molecule at 1035 cm −1 ( σa=2.2×10−16 cm 2 ) is acceptable for this transition. The required sensitivity of the detector for a successful, breadboard device could be estimated using the values of the path length (d), the glucose concentration, the cross section of absorption coefficient ( σa) , and the excitation power of the required infrared source (∼10 mW). According to this calculation of the sensitivity requirement (hundreds of μW), commercial detectors are available and thus allow a transition to a breadboard design. Currently, a quantum cascade laser will be used in conjunction with a pyroelectric detector. Both of these devices operate at room temperature and could be the beginning of a potentially implantable device. Other spectroscopic techniques besides transmission will also be considered, such as optical cavity ring-down spectroscopy (CRDS).40 CRDS has been reported to be six orders of magnitude more sensitive than FTIR for gas detection.40 Studies are ongoing to determine the applicability of this technique for the detection of glucose. 5.ConclusionsA quantitative relationship between glucose concentrations and corresponding absorption spectra changes within normal physiological ranges was demonstrated. The 1035 cm −1 frequency was the best predictor of glucose with a 20.6 mg/dl error. Additions of urea and albumin at high and low glucose concentrations had no effect on the ability to detect glucose. While this research is a basic spectroscopic study of glucose, it demonstrates the feasibility of using two frequencies of light, one to detect the water background absorption and the other to detect changes, to quantitatively determine the glucose concentrations in serum samples. AcknowledgmentsA Whitaker Foundation Special Opportunity Award provided a partial graduate student stipend for this research. The initial support and guidance of Micah D. Schmidt, Jim L. Davidson, Duco E. Jansen, and Robert L. Galloway (Vanderbilt University) during a senior design project is also acknowledged. REFERENCES
“The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin-dependent diabetes mellitus,”
N. Engl. J. Med. , 329
(14), 977
–986
(1993). Google Scholar
L. C. Clark Jr.
,
“Monitor and control of blood and tissue oxygen tissues,”
Trans. Am. Soc. Artif. Intern. Organs , 41
(2), 41
–57
(1956). Google Scholar
B. D. McKean
and
D. A. Gough
,
“A telemetry-instrumentation system for chronically implanted glucose and oxygen sensors,”
IEEE Trans. Biomed. Eng. , 35
(7), 526
–532
(1988). Google Scholar
P. H. S. Tse
and
D. A. Gough
,
“Time-dependent inactivation of immobilized glucose oxidase and catalase,”
Biotechnol. Bioeng. , 29 705
–13
(1987). Google Scholar
R. C. Mercado
and
F. Moussy
,
“In vitro and in vivo mineralization of nafion membrane used for implantable glucose sensors,”
Biosens. Bioelectron. , 13
(2), 133
–145
(1998). Google Scholar
E. Wilkins
and
P. Atanasov
,
“Integrated implantable device for long-term glucose monitoring,”
Biosens. Bioelectron. , 10
(5), 485
–494
(1995). Google Scholar
J. C. Armour
,
J. Y. Lucisano
,
B. D. McKean
, and
D. A. Gough
,
“Application of chronic intravascular blood glucose sensor in dogs,”
Diabetes , 39 1519
–1526
(1990). Google Scholar
P. H. S. Tse
and
D. A. Gough
,
“Transient response of an enzyme electrode sensor for glucose,”
Anal. Chem. , 59
(19), 2339
–2344
(1987). Google Scholar
M. Shichiri
,
N. Asakawa
,
Y. Tamasaki
,
R. Kawamoni
, and
H. Abe
,
“Telemetry glucose monitoring device with needle-type glucose sensor: A useful tool for blood glucose monitoring in diabetic individuals,”
Diabetes Care , 9
(3), 298
–301
(1986). Google Scholar
F. Moussy
,
D. J. Harrison
, and
R. V. Rajotte
,
“A miniaturized nafion-based glucose sensor,”
Int. J. Artif. Organs , 17
(2), 88
–94
(1994). Google Scholar
T. Miyasaka
,
Y. Taniyama
,
K. Sakai
, and
Y. Yoshimit
,
“Development of a glucose sensor with on/off control of enzyme activity without the effects of protein absorption,”
ASAIO J. , 43 505
–509
(1997). Google Scholar
D. Moatti-Sinat
,
F. Capnon
,
V. Poitout
,
G. Reach
,
D. S. Bindra
,
Y. Zhung
,
G. S. Wilson
, and
D. R. Therenot
,
“Towards continuous glucose monitoring: In vivo evaluation of a miniaturized glucose sensor implanted for several days in rat subcutaneous tissue,”
Diabetologia , 35 224
–230
(1992). Google Scholar
D. A. Baker
and
D. A. Gough
,
“Dynamic delay and maximal dynamic error in continuous biosensors,”
Anal. Chem. , 68
(8), 1292
–1297
(1996). Google Scholar
M. C. Shults
,
R. K. Rodes
,
S. J. Updike
,
B. J. Gilligan
, and
W. N. Reining
,
“A telemetry-instrumentation system for monitoring multiple subcutaneously implanted glucose sensors,”
IEEE Trans. Biomed. Eng. , 41
(10), 937
–942
(1994). Google Scholar
S. Ertefai
and
D. A. Gough
,
“Physiological preparation for studying the response of subcutaneously implanted glucose and oxygen sensors,”
J. Biomed. Eng. , 11 362
–368
(1989). Google Scholar
A. Maran
,
C. Crepaldi
,
A. Tiengo
,
G. Grasso
,
E. Vitali
,
G. Pagano
,
S. Bistoni
,
G. Calabrese
,
F. Santeusanio
,
F. Leonetti
,
M. Ribaudo
,
U. Di Mario
,
G. Annuzzi
,
S. Genovesse
,
G. Ricardi
,
M. Previti
,
D. Cucinotta
,
F. Giorgino
,
A. Bellomo
,
R. Giorgino
,
A. Poscia
, and
M. Varalli
,
“Continuous subcutaneous glucose monitoring in diabetic patients: A multicenter analysis,”
Diabetes Care , 25
(2), 347
–352
(2002). Google Scholar
X. Chen
,
N. Matsumota
,
Y. Hu
, and
G. S. Wilson
,
“Electrochemically mediated electrodeposition/electropolymerization to yield a glucose microbiosensor with improved characteristics,”
Anal. Chem. , 74
(2), 368
–372
(2002). Google Scholar
S. Bharathi
and
M. Nogami
,
“A glucose biosensor based on electrodeposited biocomposites of gold nanoparticles and glucose oxidase enzyme,”
Analyst , 126
(11), 1919
–1922
(2001). Google Scholar
Rm. Rhemrev-Boom
,
R. G. Tiessen
,
A. A. Jonker
,
K. Venema
,
P. Vadgama
, and
J. Korf
,
“A lightweight measuring device for the continuous in vivo monitoring of glucose by means of ultraslow microdialysis in combination with a miniaturized flow-through biosensor,”
Clin. Chim. Acta , 316
(1–2), 1
–10
(2002). Google Scholar
W. K. Ward
,
L. B. Jensen
,
E. Anderson
,
G. Reach
,
J. C. Klein
, and
G. S. Wilson
,
“A new amperometric microsensor: In vitro and short-term in vivo evaluation,”
Biosens. Bioelectron. , 17
(3), 181
–189
(2002). Google Scholar
H. Yang
,
T. D. Chung
,
Y. T. Kim
,
C. A. Choi
,
C. H. Jun
, and
H. C. Kim
,
“Glucose sensor using a microfabricated electrode and electropolymerized bilayer film,”
Biosens. Bioelectron. , 17
(3), 251
–259
(2002). Google Scholar
H. Nawawi
,
B. S. Sazali
,
B. H. Kamaruzaman
,
T. N. Yazid
,
A. A. Jemain
,
F. Ismail
, and
B. A. Khalid
,
“Effect of ambient temperature on analytical and clinical performance of a blood glucose monitoring system: Omnitest sensor glucose meter,”
Ann. Clin. Biochem. , 38
(6), 676
–683
(2001). Google Scholar
K. J. Wientjes
and
A. J. Schoonen
,
“Determination of time delay between blood and interstitial adipose tissue glucose concentration change by microdialysis in healthy volunteers,”
Int. J. Artif. Organs , 24
(12), 884
–889
(2001). Google Scholar
G. Cote´
,
“Noninvasive and minimally invasive optical monitoring technologies,”
J. Nutr. , 131 1596S
–1604S
(2001). Google Scholar
F. M. Ham
,
I. N. Kostanic
,
G. M. Cohen
, and
B. R. Gooch
,
“Determination of glucose concentrations in an aqueous matrix from NIR spectra using optimal time-domain filtering and partial least-squares regression,”
IEEE Trans. Biomed. Eng. , 44
(6), 475
–485
(1997). Google Scholar
M. A. Arnold
and
G. W. Small
,
“Determination of physiological levels of glucose in an aqueous matrix with digitally filtered Fourier transform near-infrared spectra,”
Anal. Chem. , 62 1457
–1464
(1990). Google Scholar
L. A. Marquardt
and
M. A. Aronld
,
“Near-infrared spectroscopic measurement of glucose in a protein matrix,”
Anal. Chem. , 65 3271
–3278
(1993). Google Scholar
M. R. Robinson
,
R. P. Eaton
,
D. M. Haaland
,
G. W. Koepp
,
E. V. Thomas
,
B. R. Stallard
, and
P. L. Robinson
,
“Noninvasive glucose monitoring in diabetic patients: A preliminary evaluation,”
Clin. Chem. , 38
(9), 1618
–1622
(1992). Google Scholar
H. M. Heise
,
R. Marbach
,
T. Koschinsky
, and
F. A. Gries
,
“Noninvasive blood glucose sensors based on near-infrared spectroscopy,”
Artif. Organs , 18
(6), 439
–447
(1994). Google Scholar
H. M. Heise
,
R. Marbach
,
G. Janatsch
, and
J. D. Kruse-Jarres
,
“Multivariate determination of glucose in whole blood by attenuated total reflection infrared spectroscopy,”
Anal. Chem. , 61 2009
–2015
(1989). Google Scholar
J. S. Shultz
,
S. Mansouri
, and
I. J. Goldstein
,
“Affinity sensor: A new technique for developing implantable sensors for glucose and other metabolites,”
Diabetes Care , 245
(5), 245
–253
(1982). Google Scholar
H. Zeller
,
P. Novak
, and
R. Landgraf
,
“Blood glucose measurement by infrared spectroscopy,”
Int. J. Artif. Organs , 12
(2), 129
–135
(1989). Google Scholar
K. Kajiwara
,
H. Fukushima
,
H. Kishikawa
,
K. Wishida
,
Y. Hasiquchi
,
M. Sokakida
,
M. Uehana
, and
M. Shichiri
,
“Spectroscopic quantitative analysis of blood glucose by Fourier transform infrared spectroscopy with an attenuated total reflection prism,”
Med. Prog. Technol. , 18 181
–189
(1992). Google Scholar
K. A. Oberg
and
A. L. Fink
,
“A new attenuated total reflectance Fourier transform spectroscopy method for the study of proteins in solution,”
Anal. Biochem. , 256 92
–106
(1998). Google Scholar
D. M. Back
and
P. L. Polavarapu
,
“Fourier-transform infrared spectroscopy. Sugars, structural changes in aqueous solutions,”
Carbohyd. Res. , 121 308
–311
(1983). Google Scholar
N. Kaiser
,
“Laser absorption spectroscopy with an ATR prism,”
Appl. Phys. , 17 1
–4
(1978). Google Scholar
N. Kaiser
,
“Methods for determining the contents of metabolic products in blood,”
IEEE Trans. Biomed. Eng. , 26 597
–600
(1979). Google Scholar
R. A. Shaw
,
S. Kotowich
,
M. Leroux
and
H.H. Mantsch
,
“Multianalyte serum analysis using mid-infrared spectroscopy,”
Ann. Clin. Biochem. , 35 624
–632
(1998). Google Scholar
|
CITATIONS
Cited by 37 scholarly publications and 2 patents.