Open Access Article
Pavel
Matousek
*a and
Nicholas
Stone
*b
aCentral Laser Facility, Research Complex at Harwell, STFC Rutherford Appleton Laboratory, Harwell Oxford, OX11 0QX, UK. E-mail: Pavel.Matousek@stfc.ac.uk
bSchool of Physics, University of Exeter, Exeter EX4 4QL, UK. E-mail: N.Stone@exeter.ac.uk
First published on 12th October 2015
The recently developed array of Raman spectroscopy techniques for deep subsurface analysis of biological tissues unlocks new prospects for medical diagnosis and monitoring of various biological conditions. The central pillars of these methods comprise spatially offset Raman spectroscopy (SORS) and Transmission Raman Spectroscopy facilitating penetration depths into tissue up to two orders of magnitude greater than those achievable with conventional Raman spectroscopy. This article reviews these concepts and discusses their emerging medical applications including non-invasive breast cancer diagnosis, cancer margin evaluation, bone disorder detection and glucose level determination.
Key learning points(1) Principle and context of deep Raman methods(2) The capability and limitations of deep Raman methods (3) Applicability areas of deep Raman methods in medical diagnosis and disease monitoring exemplified on real cases |
The cornerstones of these methods are Spatially Offset Raman Spectroscopy (SORS)1 and Transmission Raman Spectroscopy (TRS).2 The methods utilise the properties of photon diffusion in diffusely scattering (turbid) media in analogy with NIR absorption and fluorescence tomography techniques.3 However with Raman spectroscopy much higher chemical specificity is available and this opens a host of new applications. These advances in Raman approaches stem from earlier research into temporal properties of migrating laser and Raman photons in turbid media using time-gated methods.4 Although the time gated approaches are beneficial in subsurface probing in a number of situations their general applicability is hampered by the requirement for higher instrumental complexity and cost. Their medical use for in vivo applications is further restricted by safety laser intensity limits that are considerably more stringent with short pulsed laser beams than with continuous wave laser beams and penetration depths reached in biological tissues have not been as high as those accessible with SORS. For this reason we limit our review to methods utilising only continuous wave laser beams.
The deployment of SORS enables the retrieval of Raman spectra of sublayers within stratified turbid matrixes. For a two layer sample at least two SORS spectra acquired at different spatial offsets are required to recover the Raman spectra of individual layers. One such spectrum would typically be a zero spatially offset spectrum (equivalent to a conventional backscattering Raman spectrum) and one obtained at a non-zero spatial offset. The numerical recovery of spectra from individual layers is achieved by a scaled subtraction of the two spectra from each other with the multiplication factor adjusted to just cancel the contribution of the undesired layer leaving behind only the contributions from the layer one is recovering.1 For a multilayer sample with n-layers one needs at least n-SORS spectra acquired at different spatial offsets to retrieve individual layer contributions, in analogy with solving linear equations with n-unknowns. Alternatively, an order of magnitude larger number of SORS spectra at different spatial offsets (≥20n) can be acquired and used in conjunction with multivariate methods such as band targeted entropy minimization (BTEM)5 to retrieve the individual layer contributions. This approach could also be applicable to samples with an unknown number of layers. Other decomposition method have also been demonstrated such as an overconstrained extraction algorithm based on fitting with spectral libraries,6 2D correlation analysis.7 There are other multivariate analysis methods that can also be potentially applicable. The reader can find their general description in ref. 8. It should be noted that the Raman signature of individual layers is recovered with no prior knowledge of the composition of any of the layers. In other words the data recovery is performed blind. As such the above processing steps are also amenable to automated data analysis. Apart from analysing effectively layers in tissue the SORS concept can also recover the chemical composition of other zones within tissue of arbitrary shape which are chemically distinct from its surroundings although in its basic form SORS would not provide information on the spatial properties and location of these zones.
With medical samples one would typically use laser excitation wavelengths within the NIR optical window of tissue (e.g. 785 or 830 nm). The optical window is a term often used for a spectral region where tissue is relatively transparent for light and stretches approximately from 650 to 950 nm. Its lower edge is given by the absorption of blood and the higher end by the absorption of water and lipids.9 This range also enables optimum detection systems (silicon based CCD's) to be used and minimises contribution from interfering tissue fluorescence. The Raman signal is typically coupled through imaging optics into an optical fibre bundle and relayed onto a slit of a low f-number spectrograph with fibres re-arranged to match the linear shape of the slit for optimum coupling efficiency. Alternatively a free-space coupling (i.e. without the recourse to optical fibres) can also be utilised. The detection is typically accomplished using a cooled, high performance, deep depletion CCD camera.
A number of SORS variants have been demonstrated in recent years. These include using single point collection and point laser deposition or ring or other patterned illumination and/or collection arrangement. Among these the arrangement with a ring shaped laser illumination and collection within a central zone (‘inverse SORS’) merits a particular note in the context of medical applications.10,11 This is because the laser radiation can be spread over an extended area of the ring enabling to deliver large powers into tissue in situations where laser intensities are restricted, such as in in vivo applications. The ring-shaped laser beam can be formed using a conical lens (axicon). The spatial offset can be varied, for example, by altering the axicon-to-sample distance which changes the radius of the illumination ring.
The applicability of SORS is challenging or impossible in situations where the sample is highly absorbing at the laser or Raman signal wavelengths as this reduces photon propagation distances and as such penetration depths. A high level of fluorescence can also potentially drown weaker Raman signals. Although fluorescence originating from a non-target layer, for example, from samples surface, can be effectively reduced by SORS as outlined above.
The optimum magnitude of spatial offset is either determined empirically by trial and error using actual samples or tissue phantoms or its value can be determined numerically, e.g. using Monte Carlo simulations. A balance needs to be struck between adequate penetration depth, which requires increased spatial offset and the strength of deep layer Raman signals which diminishes with spatial offset; the target being the highest possible value of signal-to-noise ratio of the recovered Raman spectrum from the desired sublayer. This issue was investigated by Maher and Berger12 and Bloomfield et al.13 The latter study also concluded that a longer acquisition time should be spent on collecting spectra from larger spatial offsets for retrieving the optimum signal to noise ratio spectrum from the sublayer in situations where the overall collection times are restricted, e.g. due to patient being able to stay still only over limited time.
000 women.28 To enable early diagnosis, which offers more conservative treatments and better patient outcomes, the UK's National Health Service (NHS) runs a mammographic screening programme involving women over 50 years of age. Since its initiation there has been a 17.5% decrease in deaths from breast cancer for women aged 50–70 years (1996–2005).29 Despite these successes the screening suffers from several limitations; (i) the mammography is only effective with female breasts which are less dense and (ii) the identification of a suspect lesion does not convey information about whether the lesion is benign or malignant. To establish the malignancy status a needle biopsy is typically carried out. Approximately 4.4% (2008/9) of 2.1 million women screened were referred for further tests and out of the referred ones only 18% were found to have malignancies.30 As such 82% of women referred from screening in 2008/9 had investigations including excisional biopsy – with all the associated risks (from increased radiation dose to infection), costs to the NHS and inevitable psychological stress to the patient and their close relatives and friends – when they had no malignancies present.
The provision of a safe, rapid non-invasive method in conjunction with mammography could avoid the biopsy step and associated adverse impact on patient and NHS. The first seeds for this potential were laid by Haka et al.31 who recognised that the chemical content of micro-calcifications associated with benign and malignant calcifications are significantly different and that this difference is reflected also in Raman spectra. The study indicated that microcalcifications can be divided into two types: (i) Type I consisting of calcium oxalate dihydrate, and (ii) Type II consisting of calcium phosphates, mainly calcium hydroxyapatite. Calcium oxalate crystals (Type I) are mainly found in benign ductal cysts and rarely found in carcinoma,32 whereas calcium hydroxyapatite deposits (Type II) are present in proliferative lesions, which can include lesions of both benign and malignant pathology and these can be further subdivided into benign and malignant classes depending on their carbonate contents.31 Baker et al.33 further showed using Fourier transform infrared imaging that for Type II calcifications the chemical composition is not only specific to benign and malignant calcifications but it follows a continuum, from benign through the full range of ductal carcinomas in situ and invasive cancers. This is reflected most profoundly in the reduction of carbonate ion substitutions within calcium hydroxyapatite crystal lattice as the neighbouring tissue is transformed from benign to invasive disease.
The above implies that the observation of Type I calcifications signals the presence of a benign lesion. The detection of Type II calcifications represents an ambiguous situation as these are found both with benign and malignant lesions. However, these can be further differentiated into benign and malignant lesions by their carbonate content. Such chemical characterisation would not be possible reliably by mammography.
The TRS research reached the clinically relevant concentration range with a penetration depth of 20 mm in porcine tissue (see Fig. 2).34 This is about 2.5-times lower that the clinical penetration depth (∼50 mm) for wider screening use.35 It is also worth noting that the depth of 20 mm is approximately two orders of magnitude higher than possible with conventional Raman spectroscopy approaches. Further increase of penetration depth is however achievable by further enhancing the quality of the detected Raman signal by boosting its intensity. This is due the fact that for deep Raman technique the detection sensitivity is ultimately determined by the quality of Raman signal retrieved and this mostly depends on its intensity. The intensity of the detected Raman signal can be boosted by increasing the laser power incident on tissue although potential for photo-damage is possible at high power densities and associated safety illumination levels pose stringent limits here. Nevertheless these can be circumvented to a great degree using extended illumination areas (e.g. possible with inverse SORS and TRS geometries). The signal levels can also be boosted by enhancing Raman collection efficiency of spectrographs (etendue). This can be accomplished by further opening spectrograph slits permitted by higher dispersion gratings (consequently resulting also in narrower spectral ranges) or by purposely degrading spectral resolution in applications where the spectral resolution is not a critical parameter (e.g. when distinguishing calcium oxalate from calcium hydroxyapatite where marked bands are well separated from each other and other main components of tissue). In the case of breast cancer diagnosis these above measures along with the enhanced data processing methods using advanced multivariate analysis tools promise to more than double the existing penetration depths.34,36
![]() | ||
| Fig. 2 The recovered signal achieved for 0.125% (relative volume) of HAP buried within 20 mm of porcine tissue. The 962 cm−1 peak is still clearly recognizable. Acquisition times were 5 s with 10 accumulations. Vertical line, HAP peak.34 (Reprinted with permission from N. Stone and P. Matousek, Cancer Res., 2008, 68, 4424 with permission of the American Association for Cancer Research, Inc.) | ||
As mentioned above in situations where the separation of Type II calcifications into subclasses is required the carbonate content could be evaluated by deriving directly the ratio of phosphate and carbonate Raman bands at ∼960 cm−1 and ∼1070 cm−1, respectively. This can however be challenging in subsurface analysis as the carbonate Raman band at ∼1070 cm−1 is relatively weak and overlapped with Raman collagen bands of tissue. Kerssens et al.37 showed an alternative way of characterising the carbonate content, by direct monitoring of the position and bandwidth of the intense ∼960 cm−1 phosphate Raman band alone. This provides information on the carbonate content as this band's parameters are sensitive to carbonate inclusions in the lattice.
The application of SORS for transcutaneously characterising bones was first demonstrated by Schulmerich et al.43 The researchers attained depths of several millimetres through soft tissue in animal and human cadavers. The technique was then rapidly advanced expanding the penetration depths to above 4 mm.42,44 Okagbare et al. developed a multiple-fibre optic Raman probe allowing the collection of Raman spectra from multiple points around a limb to increase accuracy of recovered bone spectra to enable the monitoring of more subtle changes in composition.45 This was demonstrated with rat tibia phantoms in which the bone had carbonated hydroxyapatite with different carbonate concentrations.
More recently, the research in this area had progressed to include in vivo trials on humans within the operating theatre enabling researchers to perform a direct comparison, for the first time, of in vivo transcutaneous data with those obtained from in vivo exposed bone in surgery on the same patient.46 Preliminary data demonstrate that good correspondence between the exposed data and transcutaneous data can be achieved.
Transcutaneous Raman is also under development for early detection of burn induced heterotopic ossification (HO) of soft tissue,47 a condition associated with major burn injuries and blast traumas. The potential of this approach was demonstrated on mice following a burn injury (see Fig. 3). Raman data showed significantly increased bone mineral signalling in the tenotomy (surgical procedure) compared to control leg at 5 days post injury, with the difference increasing over time. In contrast, micro CT did not demonstrate heterotopic bone until three weeks post injury. Changes in bone mineral and matrix composition of the new bone were also evidenced in the Raman spectra which could facilitate early identification of HO and allow more timely therapy decisions for HO patients.
![]() | ||
| Fig. 3 Raman spectroscopy and cross sectional Micro CT of Achilles tenotomy model on non-injured and on tenotomized leg at 3 months after injury and burn. (A) (top) Raman spectra of tenotomized leg (red) and non-tenotomized control leg (blue) of a burn mouse that had known HO growth. Spectra are normalized to the protein matrix and collagen band at 1600 cm−1 and superimposed to show differences in bone mineral signal at 958 cm−1 (bottom) mineral to matrix ratio from Raman spectra demonstrates increased mineral content in the tenotomized leg (*, pb 0.05). (B) Micro CT confirmation of HO growth in the tenotomized leg seen in representative CT slices (top) and 3D reconstructions (middle). Red arrows indicate HO formation. Gray areas in the reconstructed image indicate ectopic bone. Histologic verification of ectopic bone growth with pentachrome stain (bottom) showing bone as pale yellow. Red arrow indicates HO formation. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)47 (Reprinted from J. R. Peterson, P. I. Okagbare, S. De La Rosa, K. E. Cilwa, J. E. Perosky, O. N. Eboda, A. Donneys, G. L. Su, S. R. Buchman, P. S. Cederna, S. C. Wang, K. M. Kozloff, M. D. Morris, B. Levi, ‘Early detection of burn induced heterotopic ossification using transcutaneous Raman spectroscopy’, Bone, 54, 28–34, Copyright (2013), with permission from Elsevier.) | ||
Recently, Buckley et al applied inverse SORS10,11 to noninvasive diagnosis of osteogenesis imperfecta (‘brittle bone’) condition in vivo.48 The study succeeded in detecting the presence of this condition within a single patient. The obtained spectra were also consistent with Raman spectra obtained ex vivo from the same patient (see Fig. 4). Although the study is only the first step towards delivering a diagnostic method it outlines the potential of SORS in this area. Although it should be noted that specifically in the case of osteogenesis imperfecta, which is genetic condition, effective DNA screening methods also exist although in the UK, for example, these are not yet in routine use.
![]() | ||
| Fig. 4 (a) A Raman spectrum of excised OI bone (black) and excised control bone (red). The inset shows mean ± s.d. (b) A spatially offset Raman spectrum retrieved noninvasively through the skin of the same OI patient (black) and a spatially offset Raman spectrum retrieved noninvasively through the skin of an age-matched control (red). The inset shows mean ± s.d.48 (Reprinted with permission from K. Buckley, J. G. Kerns, P. D. Gikas, H. L. Birch, J. Vinton, R. Keen, A. W. Parker, P. Matousek, A. E. Goodship, IBMS BoneKEy, 2014, 11, 602 with permission of the Nature Publishing Group.) | ||
In the area of the diagnosis of osteoporotic conditions by SORS progress has also been recently made. Buckley et al.49 demonstrated the potential of SORS in this area building on earlier advances by Morris’ group.42 Buckley's study showed that on average, bone fragments from the necks of fractured femora measured ex vivo are more mineralised (by 5–10%) than (cadaveric) non-fractured controls, but the mineralisation distributions of the two cohorts are largely overlapped. SORS in vivo measurements indicated a potential of the presence of similar differences but these were as yet statistically underpowered. The study also identified methodological developments which could be implemented to improve the statistical significance of future experiments that may eventually lead to more sensitive prediction of fragility fractures in vivo.
It should, however, be noted that, in general, this is a difficult area to penetrate due to the widespread of DXA technique, despite DXA existing limitations, and only considerable outperformance of DXA by deep Raman methods would be expected lead to the change of medical practice. Initially, the deep Raman methods are likely to be used as complementary tools to the existing methods rather than replacing them outright.
Although it should be pointed out that this is a very challenging area under intense focus of many groups worldwide utilising a wide range of techniques. As such Raman techniques will face tough competition in the market place and their ultimate success will depend on social economic benefits they would bring as well as the robustness and accuracy of these.
![]() | ||
| Fig. 5 Calibration and validation data sets of SESORS quantitative measurements of BPE functionalized SERS nanotags through bone.53 (Reprinted with permission from B. Sharma, K. Ma, M. R. Glucksberg, R. P. Van Duyne, J. Am. Chem. Soc., 2013, 135, 17290 Copyright (2013) American Chemical Society.) | ||
| This journal is © The Royal Society of Chemistry 2016 |