Combining multiset resolution and segmentation for hyperspectral image analysis of biological tissues

Treatment of recurrent acute tonsillitis-a systematic review and clinical practice recommendations

Background

There is an ongoing debate on the indications for tonsil surgery in both children and adults with recurrent acute tonsillitis. The aim is to provide practical recommendations for diagnostics and treatment for recurrent acute tonsillitis including evidence-based decision making for tonsillectomy.

Methods

A systematic literature search in PubMed, Embase, Web of Science, and ScienceDirect from 2014 until April 2023 resulted in 68 articles. These were the basis for the review and a comprehensive series of consensus statements on the most important diagnostics and indications for both non-surgical and surgical therapy. A consensus paper was circulated among the authors and members of the International Head and Neck Scientific Group until a final agreement was reached for all recommendations.

Results

The differentiation between sore throat and tonsillitis patient episodes is mostly not feasible and hence is not relevant for diagnostic decision making. Diagnostics of a tonsillitis/sore throat episode should always include a classification with a scoring system (Centor, McIssac, FeverPAIN score) to estimate the probability of a bacterial tonsillitis, mainly due to group A streptococcus (GAS). In ambiguous cases, a point-of-care test GAS swab test is helpful. Consecutive counting of the tonsillitis/sore throat episodes is important. In addition, a specific quality of life score (Tonsillectomy Outcome Inventory 14 or Tonsil and Adenoid Health Status Instrument) should be used for each episode. Conservative treatment includes a combination of paracetamol and/or non-steroidal anti-inflammatory drugs. In case of high probability of bacterial tonsillitis, and only in such cases, especially in patients at risk, standard antibiotic treatment is initiated directly or by delayed prescription. Tonsillectomy is indicated and is highly effective if the patient has had ≥7 adequately treated episodes in the preceding year, ≥5 such episodes in each of the preceding 2 years, or ≥3 such episodes in each of the preceding 3 years. An essential part of surgery is standardized pain management because severe postoperative pain can be expected in most patients.

Conclusion

It is necessary to follow a stringent treatment algorithm for an optimal and evidence-based treatment for patients with recurrent acute tonsillitis. This will help decrease worldwide treatment variability, antibiotic overuse, and avoid ineffective tonsillectomy.

Biomolecular alterations detected in multiple sclerosis skin fibroblasts using Fourier transform infrared spectroscopy

Multiple sclerosis (MS) is the leading cause of non-traumatic disability in young adults. New avenues are needed to help predict individuals at risk for developing MS and aid in diagnosis, prognosis, and outcome of therapeutic treatments. Previously, we showed that skin fibroblasts derived from patients with MS have altered signatures of cell stress and bioenergetics, which likely reflects changes in their protein, lipid, and biochemical profiles. Here, we used Fourier transform infrared (FTIR) spectroscopy to determine if the biochemical landscape of MS skin fibroblasts were altered when compared to age- and sex-matched controls (CTRL). More so, we sought to determine if FTIR spectroscopic signatures detected in MS skin fibroblasts are disease specific by comparing them to amyotrophic lateral sclerosis (ALS) skin fibroblasts. Spectral profiling of skin fibroblasts from MS individuals suggests significant alterations in lipid and protein organization and homeostasis, which may be affecting metabolic processes, cellular organization, and oxidation status. Sparse partial least squares-discriminant analysis of spectral profiles show that CTRL skin fibroblasts segregate well from diseased cells and that changes in MS and ALS may be unique. Differential changes in the spectral profile of CTRL, MS, and ALS cells support the development of FTIR spectroscopy to detect biomolecular modifications in patient-derived skin fibroblasts, which may eventually help establish novel peripheral biomarkers.

Unsupervised segmentation of biomedical hyperspectral image data: tackling high dimensionality with convolutional autoencoders

Information about the structure and composition of biopsy specimens can assist in disease monitoring and diagnosis. In principle, this can be acquired from Raman and infrared (IR) hyperspectral images (HSIs) that encode information about how a sample's constituent molecules are arranged in space. Each tissue section/component is defined by a unique combination of spatial and spectral features, but given the high dimensionality of HSI datasets, extracting and utilising them to segment images is non-trivial. Here, we show how networks based on deep convolutional autoencoders (CAEs) can perform this task in an end-to-end fashion by first detecting and compressing relevant features from patches of the HSI into low-dimensional latent vectors, and then performing a clustering step that groups patches containing similar spatio-spectral features together. We showcase the advantages of using this end-to-end spatio-spectral segmentation approach compared to i) the same spatio-spectral technique not trained in an end-to-end manner, and ii) a method that only utilises spectral features (spectral k-means) using simulated HSIs of porcine tissue as test examples. Secondly, we describe the potential advantages/limitations of using three different CAE architectures: a generic 2D CAE, a generic 3D CAE, and a 2D convolutional encoder-decoder architecture inspired by the recently proposed UwU-net that is specialised for extracting features from HSI data. We assess their performance on IR HSIs of real colon samples. We find that all architectures are capable of producing segmentations that show good correspondence with HE stained adjacent tissue slices used as approximate ground truths, indicating the robustness of the CAE-driven spatio-spectral clustering approach for segmenting biomedical HSI data. Additionally, we stress the need for more accurate ground truth information to enable a precise comparison of the advantages offered by each architecture.

High-resolution MRI of the human palatine tonsil and its schematic anatomic 3D reconstruction

The palatine tonsils form an important part of the human immune system. Together with the other lymphoid tonsils of Waldeyer's tonsillar ring, they act as the first line of defense against ingested or inhaled pathogens. Although histologically stained sections of the palatine tonsil are widely available, they represent the tissue only in two dimensions and do not provide reference to three‐dimensional space. Such a representation of a tonsillar specimen based on imaging data as a 3D anatomical reconstruction is lacking both in scientific publications and especially in textbooks. As a first step in this direction, the objective of the present work was to image a resected tonsil specimen with high spatial resolution in a 9.4 T small‐bore pre‐clinical MRI and to combine these data with data from the completely sectioned and H&E stained same palatine tonsil. Based on the information from both image modalities, a 3D anatomical sketch was drawn by a scientific graphic artist. In perspective, such studies could help to overcome the difficulty of capturing the spatial extent and arrangement of anatomical structures from 2D images and to establish a link between three‐dimensional anatomical preparations and two‐dimensional sections or illustrations, as they have been found so far in common textbooks and anatomical atlases.

Use of physiological information based on grayscale images to improve mass spectrometry imaging data analysis from biological tissues

The characterization of cancer tissues by matrix-assisted laser desorption ionization-mass spectrometry images (MALDI-MSI) is of great interest because of the power of MALDI-MS to understand the composition of biological samples and the imaging side that allows for setting spatial boundaries among tissues of different nature based on their compositional differences. In tissue-based cancer research, information on the spatial location of necrotic/tumoral cell populations can be approximately known from grayscale images of the scanned tissue slices. This study proposes as a major novelty the introduction of this physiologically-based information to help in the performance of unmixing methods, oriented to extract the MS signatures and distribution maps of the different tissues present in biological samples. Specifically, the information gathered from grayscale images will be used as a local rank constraint in Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) for the analysis of MALDI-MSI of cancer tissues. The use of this constraint, setting absence of certain kind of tissues only in clear zones of the image, will help to improve the performance of MCR-ALS and to provide a more reliable definition of the chemical MS fingerprint and location of the tissues of interest. The general strategy to address the analysis of MALDI-MSI of cancer tissues will involve the study of the MCR-ALS results and the posterior use of MCR-ALS scores as dimensionality reduction for image segmentation based on K-means clustering. The resolution method will provide the MS signatures and their distribution maps for each tissue in the sample. Then, the resolved distribution maps for each biological component (MCR scores) will be submitted as initial information to K-means clustering for image segmentation to obtain information on the boundaries of the different tissular regions in the samples studied. MCR-ALS prior to K-means not only provides the desired dimensionality reduction, but additionally resolved non-biological signal contributions are not used and the weight given to the different biological components in the segmentation process can be modulated by suitable preprocessing methods.

A novel multivariate curve resolution-alternating least squares (MCR-ALS) methodology for application in hyperspectral Raman imaging analysis

Multivariate curve resolution-alternating least squares (MCR-ALS) applied to hyperspectral Raman imaging is extensively used to spatially and spectrally resolve the individual, pure chemical species within complex, heterogeneous samples. A critical aspect of performing MCR-ALS with hyperspectral Raman imaging is the selection of the number of chemical components within the experimental data. Several methods have previously been proposed to determine the number of chemical components, but it remains a challenging task that if done incorrectly, can lead to the loss of chemical information. In this work, we show that the choice of 'optimal' number of factors in the MCR-ALS model may vary depending on the relative contribution of the targeted species to the overall spectral intensity. In a data set consisting of 27 hyperspectral Raman images of TiO2 polymorphs, it was observed that the more dominant species were best resolved with a parsimonious model. However, species with intensities near the noise level often needed more factors to be resolved than was predicted by standard methods. Based on the observations in this data set, we propose a new method that employs approximate reference spectra for determining optimal model complexity for identifying minor constituents with MCR-ALS.

Potential use of MCR-ALS for the identification of coeliac-related biochemical changes in hyperspectral Raman maps from pediatric intestinal biopsies

Raman hyperspectral imaging is an emerging practice in biological and biomedical research for label free analysis of tissues and cells. Using this method, both spatial distribution and spectral information of analyzed samples can be obtained. The current study reports the first Raman microspectroscopic characterisation of colon tissues from patients with Coeliac Disease (CD). The aim was to assess if Raman imaging coupled with hyperspectral multivariate image analysis is capable of detecting the alterations in the biochemical composition of intestinal tissues associated with CD. The analytical approach was based on a multi-step methodology: duodenal biopsies from healthy and coeliac patients were measured and processed with Multivariate Curve Resolution Alternating Least Squares (MCR-ALS). Based on the distribution maps and the pure spectra of the image constituents obtained from MCR-ALS, interesting biochemical differences between healthy and coeliac patients has been derived. Noticeably, a reduced distribution of complex lipids in the pericryptic space, and a different distribution and abundance of proteins rich in beta-sheet structures was found in CD patients. The output of the MCR-ALS analysis was then used as a starting point for two clustering algorithms (k-means clustering and hierarchical clustering methods). Both methods converged with similar results providing precise segmentation over multiple Raman images of studied tissues.

Discriminating normal regions within cancerous hen ovarian tissue using multivariate hyperspectral image analysis

RATIONALE

Identification of subregions under different pathological conditions on cancerous tissue is of great significance for understanding cancer progression and metastasis. Infrared matrix-assisted laser desorption electrospray ionization mass spectrometry (IR-MALDESI-MS) can be potentially used for diagnostic purposes since it can monitor spatial distribution and abundance of metabolites and lipids in biological tissues. However, the large size and high dimensionality of hyperspectral data make analysis and interpretation challenging. To overcome these barriers, multivariate methods were applied to IR-MALDESI data for the first time, aiming at efficiently resolving mass spectral images, from which these results were then used to identify normal regions within cancerous tissue.

METHODS

Molecular profiles of healthy and cancerous hen ovary tissues were generated by IR-MALDESI-MS. Principal component analysis (PCA) combined with color-coding built a single tissue image which summarizes the high-dimensional data features. Pixels with similar color indicated similar composition. PCA results from healthy tissue were further used to test each pixel in cancerous tissue to determine if it is healthy. Multivariate curve resolution-alternating least squares (MCR-ALS) was used to obtain major spatial features existing in ovary tissues, and group molecules with the same distribution patterns simultaneously.

RESULTS

PCA as the predominating dimensionality reduction approach captured over 90% spectral variances by the first three PCs. The PCA images show the cancerous tissue is more chemically heterogeneous than healthy tissue, where at least four regions with different m/z profiles can be differentiated. PCA modeling assigns top regions of cancerous tissue as healthy-like. MCR-ALS extracted three and four major compounds from healthy and cancerous tissue, respectively. Evaluating similarities of resolved spectra uncovered the chemical components that were distinct in some regions on cancerous tissue, serving as a supplementary way to differentiate healthy and cancerous regions.

CONCLUSIONS

Two unsupervised chemometric methods including PCA and MCR-ALS were applied for resolving and visualizing IR-MALDESI-MS data acquired from hen ovary tissues, improving the interpretation of mass spectrometry imaging results. Then possible normal regions were differentiated from cancerous tissue sections. No prior knowledge is required using either chemometric method, so our approach is readily suitable for unstained tissue samples, which allows one to reveal the molecular events happening during disease progression.

Analyzing Raman spectroscopic data

This chapter is a short introduction into the data analysis pipeline, which is typically utilized to analyze Raman spectra. We empathized in the chapter that this data analysis pipeline must be tailored to the specific application of interest. Nevertheless, the tailored data analysis pipeline consists always of the same general procedures applied sequentially. The utilized procedures correct for artefacts, standardize the measured spectral data and translate the spectroscopic signals into higher level information. These computational procedures can be arranged into separate groups namely data pre-treatment, pre-processing and modeling. Thereby the pre-treatment aims to correct for non-sample-dependent artefacts, like cosmic spikes and contributions of the measurement device. The block of procedures, which needs to be applied next, is called pre-processing. This group consists of smoothing, baseline correction, normalization and dimension reduction. Thereafter, the analysis model is constructed and the performance of the models is evaluated. Every data analysis pipeline should be composed of procedures of these three groups and we describe every group in this chapter. After the description of data pre-treatment, pre-processing and modeling, we summarized trends in the analysis of Raman spectra namely model transfer approaches and data fusion. At the end of the chapter we tried to condense the whole chapter into guidelines for the analysis of Raman spectra.

Understanding the Formation of Heartwood in Larch Using Synchrotron Infrared Imaging Combined With Multivariate Analysis and Atomic Force Microscope Infrared Spectroscopy

Formation of extractive-rich heartwood is a process in live trees that make them and the wood obtained from them more resistant to fungal degradation. Despite the importance of this natural mechanism, little is known about the deposition pathways and cellular level distribution of extractives. Here we follow heartwood formation in Larix gmelinii var. Japonica by use of synchrotron infrared images analyzed by the unmixing method Multivariate Curve Resolution - Alternating Least Squares (MCR-ALS). A subset of the specimens was also analyzed using atomic force microscopy infrared spectroscopy. The main spectral changes observed in the transition zone when going from sapwood to heartwood was a decrease in the intensity of a peak at approximately 1660 cm^-1 and an increase in a peak at approximately 1640 cm^-1. There are several possible interpretations of this observation. One possibility that is supported by the MCR-ALS unmixing is that heartwood formation in larch is a type II or Juglans-type of heartwood formation, where phenolic precursors to extractives accumulate in the sapwood rays. They are then oxidized and/or condensed in the transition zone and spread to the neighboring cells in the heartwood.

Assessment of tissue-specific multifactor effects in environmental -omics studies of heterogeneous biological samples: Combining hyperspectral image information and chemometrics

The use of hyperspectral imaging techniques in biological studies has increased in the recent years. Hyperspectral images (HSI) provide chemical information and preserve the morphology and original structure of heterogeneous biological samples, which can be potentially useful in environmental -omics studies when effects due to several factors, e.g., contaminant exposure, phenotype,…, at a specific tissue level need to be investigated. Yet, no available strategies exist to exploit adequately this kind of information. This work offers a novel chemometric strategy to pass from the raw image information to useful knowledge in terms of statistical assessment of the multifactor effects of interest in -omic studies. To do so, unmixing of the hyperspectral image measurement is carried out to provide tissue-specific information. Afterwards, several specific ANOVA-Simultaneous Component Analysis (ASCA) models are generated to properly assess and interpret the diverse effect of the factors of interest on the spectral fingerprints of the different tissues characterized. The unmixing step is performed by Multivariate Curve Resolution-Alternating Least Squares (MCR-ALS) on multisets of biological images related to each studied condition and provides reliable HSI spectral signatures and related image maps for each specific tissue in the regions imaged. The variability associated with these signatures within a population is obtained through an MCR-based resampling step on representative pixel subsets of the images analyzed. All spectral fingerprints obtained for a particular tissue in the different conditions studied are used to obtain the related ASCA model that will help to assess the significance of the factors studied on the tissue and, if relevant, to describe the associated fingerprint modifications. The potential of the approach is assessed in a real case of study linked to the investigation of the effect of exposure time to chlorpyrifos-oxon (CPO) on ocular tissues of different phenotypes of zebrafish larvae from Raman HSI of eye cryosections. The study allowed the characterization of melanin, crystalline and internal eye tissue and the phenotype, exposure time and the interaction of the two factors were found to be significant in the changes found in all kind of tissues. Factor-related changes in the spectral fingerprint were described and interpreted per each kind of tissue characterized.

Multivariate unmixing approaches on Raman images of plant cell walls: new insights or overinterpretation of results?

BACKGROUND

Plant cell walls are nanocomposites based on cellulose microfibrils embedded in a matrix of polysaccharides and aromatic polymers. They are optimized for different functions (e.g. mechanical stability) by changing cell form, cell wall thickness and composition. To reveal the composition of plant tissues in a non-destructive way on the microscale, Raman imaging has become an important tool. Thousands of Raman spectra are acquired, each one being a spatially resolved molecular fingerprint of the plant cell wall. Nevertheless, due to the multicomponent nature of plant cell walls, many bands are overlapping and classical band integration approaches often not suitable for imaging. Multivariate data analysing approaches have a high potential as the whole wavenumber region of all thousands of spectra is analysed at once.

RESULTS

Three multivariate unmixing algorithms, vertex component analysis, non-negative matrix factorization and multivariate curve resolution-alternating least squares were applied to find the purest components within datasets acquired from micro-sections of spruce wood and Arabidopsis. With all three approaches different cell wall layers (including tiny S1 and S3 with 0.09-0.14 µm thickness) and cell contents were distinguished and endmember spectra with a good signal to noise ratio extracted. Baseline correction influences the results obtained in all methods as well as the way in which algorithm extracts components, i.e. prioritizing the extraction of positive endmembers by sequential orthogonal projections in VCA or performing a simultaneous extraction of non-negative components aiming at explaining the maximum variance in NMF and MCR-ALS. Other constraints applied (e.g. closure in VCA) or a previous principal component analysis filtering step in MCR-ALS also contribute to the differences obtained.

CONCLUSIONS

VCA is recommended as a good preliminary approach, since it is fast, does not require setting many input parameters and the endmember spectra result in good approximations of the raw data. Yet the endmember spectra are more correlated and mixed than those retrieved by NMF and MCR-ALS methods. The latter two give the best model statistics (with lower lack of fit in the models), but care has to be taken about overestimating the rank as it can lead to artificial shapes due to peak splitting or inverted bands.

Tracking hidden organic carbon in rocks using chemometrics and hyperspectral imaging

Finding traces of life or organic components of prebiotic interest in the rock record is an appealing goal for numerous fields in Earth and space sciences. However, this is often hampered by the scarceness and highly heterogeneous distribution of organic compounds within rocks. We assess here an innovative analytical strategy combining Synchrotron radiation-based Fourier-Transform Infrared microspectroscopy (S-FTIR) and multivariate analysis techniques to track and characterize organic compounds at the pore level in complex oceanic rocks. S-FTIR hyperspectral images are analysed individually or as multiple image combinations (multiset analysis) using Principal Component Analyses (PCA) and Multivariate Curve Resolution – Alternating Least Squares (MCR-ALS). This approach allows extracting simultaneously pure organic and mineral spectral signatures and determining their spatial distributions and relationships. MCR-ALS analysis provides resolved S-FTIR signatures of 8 pure mineral and organic components showing the close association at a micrometric scale of organic compounds and secondary clays formed during rock alteration and known to catalyse organic synthesis. These results highlights the potential of the serpentinizing oceanic lithosphere to generate and preserve organic compounds of abiotic origin, in favour of the hydrothermal theory for the origin of life.

Identification of trace additives in polymer materials by attenuated total reflection Fourier transform infrared mapping coupled with multivariate curve resolution

Although multivariate curve resolution (MCR) has been applied to the analysis of Fourier transform infrared (FTIR) imaging, it is still problematic to determine the number of components. The reported methods at present tend to cause the components of low concentration missed. In this paper a new idea was proposed to resolve this problem. First, MCR calculation was repeated by increasing the number of components sequentially, then each retrieved pure spectrum of as-resulted MCR component was directly compared with a real-world pixel spectrum of the local high concentration in the corresponding MCR map. One component was affirmed only if the characteristic bands of the MCR component had been included in its pixel spectrum. This idea was applied to attenuated total reflection (ATR)/FTIR mapping for identifying the trace additives in blind polymer materials and satisfactory results were acquired. The successful demonstration of this novel approach opens up new possibilities for analyzing additives in polymer materials.

Label-Free Molecular Imaging of Biological Cells and Tissues by Linear and Nonlinear Raman Spectroscopic Approaches

Raman spectroscopy is an emerging technique in bioanalysis and imaging of biomaterials owing to its unique capability of generating spectroscopic fingerprints. Imaging cells and tissues by Raman microspectroscopy represents a nondestructive and label-free approach. All components of cells or tissues contribute to the Raman signals, giving rise to complex spectral signatures. Resonance Raman scattering and surface-enhanced Raman scattering can be used to enhance the signals and reduce the spectral complexity. Raman-active labels can be introduced to increase specificity and multimodality. In addition, nonlinear coherent Raman scattering methods offer higher sensitivities, which enable the rapid imaging of larger sampling areas. Finally, fiber-based imaging techniques pave the way towards in vivo applications of Raman spectroscopy. This Review summarizes the basic principles behind medical Raman imaging and its progress since 2012.

Tip in-light on: Advantages, challenges, and applications of combining AFM and Raman microscopy on biological samples

Scanning probe microscopies and spectroscopies, especially AFM and Confocal Raman microscopy are powerful tools to characterize biological materials. They are both non-destructive methods and reveal mechanical and chemical properties on the micro and nano-scale. In the last years the interest for increasing the lateral resolution of optical and spectral images has driven the development of new technologies that overcome the diffraction limit of light. The combination of AFM and Raman reaches resolutions of about 50-150 nm in near-field Raman and 1.7-50 nm in tip enhanced Raman spectroscopy (TERS) and both give a molecular information of the sample and the topography of the scanned surface. In this review, the mentioned approaches are introduced, the main advantages and problems for application on biological samples discussed and some examples for successful experiments given. Finally the potential of colocated AFM and Raman measurements is shown on a case study of cellulose-lignin films: the topography structures revealed by AFM can be related to a certain chemistry by the colocated Raman scan and additionally the mechanical properties be revealed by using the digital pulsed force mode. Microsc. Res. Tech. 80:30-40, 2017. © 2016 Wiley Periodicals, Inc.

Segmentation in dermatological hyperspectral images: dedicated methods

Background

Segmentation of hyperspectral medical images is one of many image segmentation methods which require profiling. This profiling involves either the adjustment of existing, known image segmentation methods or a proposal of new dedicated methods of hyperspectral image segmentation. Taking into consideration the size of analysed data, the time of analysis is of major importance. Therefore, the authors proposed three new dedicated methods of hyperspectral image segmentation with special reference to the time of analysis.

Methods

The segmentation methods presented in this paper were tested and profiled to the images acquired from different hyperspectral cameras including SOC710 Hyperspectral Imaging System, Specim sCMOS-50-V10E. Correct functioning of the method was tested for over 10,000 2D images constituting the sequence of over 700 registrations of the areas of the left and right hand and the forearm.

Results

As a result, three new methods of hyperspectral image segmentation have been proposed: fast analysis of emissivity curves (SKE), 3D segmentation (S3D) and hierarchical segmentation (SH). They have the following features: are fully automatic; allow for implementation of fast segmentation methods; are profiled to hyperspectral image segmentation; use emissivity curves in the model form, can be applied in any type of objects not necessarily biological ones, are faster (SKE—2.3 ms, S3D—1949 ms, SH—844 ms for the computer with Intel^® Core i7 4960X CPU 3.6 GHz) and more accurate (SKE—accuracy 79 %, S3D—90 %, SH—92 %) in comparison with typical methods known from the literature.

Conclusions

Profiling and/or proposing new methods of hyperspectral image segmentation is an indispensable element of developing software. This ensures speed, repeatability and low sensitivity of the algorithm to changing parameters.

Investigation of intervertebral disc degeneration using multivariate FTIR spectroscopic imaging

Traditionally tissue samples are analysed using protein or enzyme specific stains on serial sections to build up a picture of the distribution of components contained within them. In this study we investigated the potential of multivariate curve resolution-alternating least squares (MCR-ALS) to deconvolute 2nd derivative spectra of Fourier transform infrared (FTIR) microscopic images measured in transflectance mode of goat and human paraffin embedded intervertebral disc (IVD) tissue sections, to see if this methodology can provide analogous information to that provided by immunohistochemical stains and bioassays but from a single section. MCR-ALS analysis of non-degenerate and enzymatically in vivo degenerated goat IVDs reveals five matrix components displaying distribution maps matching histological stains for collagen, elastin and proteoglycan (PG), as well as immunohistochemical stains for collagen type I and II. Interestingly, two components exhibiting characteristic spectral and distribution profiles of proteoglycans were found, and relative component/tissue maps of these components (labelled PG1 and PG2) showed distinct distributions in non-degenerate versus mildly degenerate goat samples. MCR-ALS analysis of human IVD sections resulted in comparable spectral profiles to those observed in the goat samples, highlighting the inter species transferability of the presented methodology. Multivariate FTIR image analysis of a set of 43 goat IVD sections allowed the extraction of semi-quantitative information from component/tissue gradients taken across the IVD width of collagen type I, collagen type II, PG1 and PG2. Regional component/tissue parameters were calculated and significant correlations were found between histological grades of degeneration and PG parameters (PG1: p = 0.0003, PG2: p < 0.0001); glycosaminoglycan (GAG) content and PGs (PG1: p = 0.0055, PG2: p = 0.0001); and MRI T2* measurements and PGs (PG1: p = 0.0021, PG2: p < 0.0001). Additionally, component/tissue parameters for collagen type I and II showed significant correlations with total collagen content (p = 0.0204, p = 0.0127). In conclusion, the presented findings illustrate, that the described multivariate FTIR imaging approach affords the necessary chemical specificity to be considered an important tool in the study of IVD degeneration in goat and human IVDs.

Two-hierarchical nonnegative matrix factorization distinguishing the fluorescent targets from autofluorescence for fluorescence imaging

BACKGROUND

Nonnegative matrix factorization (NMF) has been used in blind fluorescence unmixing for multispectral in-vivo fluorescence imaging, which decomposes a mixed source data into a set of constituent fluorescence spectra and corresponding concentrations. However, most classical NMF algorithms have ill convergence problems and they always fail to unmix multiple fluorescent targets from background autofluorescence for the sparse acquisition of multispectral fluorescence imaging, which introduces incomplete measurements and severe discontinuities in multispectral fluorescence emissions across the multiple spectral bands.

METHODS

Observing the spatial distinction between the diffusive autofluorescence and the sparse fluorescent targets, we propose to separate the mixed sparse multispectral data into equality constrained two-hierarchical updating within NMF framework by dividing the concentration matrix of entire endmembers into two hierarchies: the fluorescence targets and the background autofluorescence. Specifically, when updating concentrations of multiple fluorescent targets in the two-hierarchical NMF, we assume that the concentration of autofluorescence is fixed and known, and vice versa. Furthermore, a sparsity constraint is imposed on the concentration matrix components of fluorescence targets only.

RESULTS

Synthetic data sets, in vivo fluorescence imaging data are employed to demonstrate and validate the performance of our approach. The proposed algorithm can achieve more satisfying results of spectral unmixing and autofluorescence removal compared to other state-of-the-art methods, especially for the sparse multispectral fluorescence imaging.

CONCLUSIONS

The proposed algorithm can successfully tackle the sparse acquisition and ill-posed problems in the NMF-based fluorescence unmixing through equality constraint along with partial sparsity constraint during two-hierarchical NMF optimization, at which fixing sparsity constrained target fluorescence can make the update of autofluorescence as accurate as possible and vice versa.