Abstract

Bioimaging using endogenous cell fluorescence, without any external biomarkers makes it possible to explore cells and tissues in their original native state, also in vivo. In order to be informative, this label-free method requires careful multispectral or hyperspectral recording of autofluorescence images followed by unsupervised extraction (unmixing) of biochemical signatures. The unmixing is difficult due to the scarcity of biochemically pure regions in cells and also because autofluorescence is weak compared with signals from labelled cells, typically leading to low signal to noise ratio. Here, we solve the problem of unsupervised hyperspectral unmixing of cellular autofluorescence by introducing the Robust Dependent Component Analysis (RoDECA). This approach provides sophisticated and statistically robust quantitative biochemical analysis of cellular autofluorescence images. We validate our method on artificial images, where the addition of varying known level of noise has allowed us to quantify the accuracy of our RoDECA analysis in a way that can be applied to real biological datasets. The same unsupervised statistical minimisation is then applied to imaging of mouse retinal photoreceptor cells where we establish the identity of key endogenous fluorophores (free NADH, FAD and lipofuscin) and derive the corresponding molecular abundance maps. The pre-processing methodology of image datasets is also presented, which is essential for the spectral unmixing analysis, but mostly overlooked in the previous studies.

Highlights

  • IntroductionThe hyperspectral dataset of image pixel i, (i = 1, ..., Nde)2s3c.rTibheisdmbaytarixmcaatnrixbeyc=ons[iydkei r]e, dwahsearesuymki is of the the noiseless signal matrix x and the matrix of image noise n, where y = x + n

  • Background of hyperspectral unmixingIn hyperspectral imaging, two-dimensional images (N pixels each; here N ≈ 106) of the same sample are captured in L different spectral channels (L = 18 in our study)

  • For the unmixing to be accurate, the number of components needs to be smaller than the number of channels (p < L) Algebraically, in the linear mixing model (LMM), the noiseless signals xare expressed with the aid of an endmember matrix M = [Mkj], is weighted with abundance fractions specified in the abundance matrix s = [sji], (j = 1, ..., p, i = 1, ..., N) according to p x = [xki] = ∑Mkjsji = Ms

Read more

Summary

Introduction

The hyperspectral dataset of image pixel i, (i = 1, ..., Nde)2s3c.rTibheisdmbaytarixmcaatnrixbeyc=ons[iydkei r]e, dwahsearesuymki is of the the noiseless signal matrix x and the matrix of image noise n, where y = x + n. For the unmixing to be accurate, the number of components needs to be smaller than the number of channels (p < L) Algebraically, in the LMM, the noiseless signals xare expressed with the aid of an endmember matrix M = [Mkj] , (where, k = 1, ..., L, j = 1, ..., p) is weighted with abundance fractions specified in the abundance matrix s = [sji], (j = 1, ..., p, i = 1, ..., N) according to p x = [xki] = ∑Mkjsji = Ms.

Objectives
Methods
Results
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call