The present invention relates to human cancer diagnostic systems and, more particularly, to a photonic structural and chemometric pathology system utilizing unlabeled tissue sections for cancer and pre-cancerous detection, diagnosis, monitoring, and prognosis.
There are a wide range of conventional approaches to cancer and precancerous detection, diagnosis, monitoring, and prognosis. Techniques involving tissue analysis generally utilize tissue sections that are stained or otherwise labeled to increase the contrast of effected areas. Preparing the sections for analysis usually requires time consuming processing including the application stain or other labeling techniques that obscure some aspects of the sample. At present, there are not effective cancer detection and diagnostic techniques available that utilize fresh or frozen tissue sections without prior processing, such as staining or other labeling techniques. This delays and complicates the detection and diagnostic process and makes the process more expensive while also limiting the pathology lab from processing more samples and generating more revenue.
The current pathology practice using conventional light microscopy has several critical limitations, including a) time consuming and costly multiple step tissue processing is required to prepare the sample; b) the lateral resolution is about 0.5 mm or worse and no axial information of the specimen is revealed; and c) the diagnosis based mainly on the manual expert examination of the tissue architectural and morphological changes is subjective and suffers from inter- and intra-observer variations. The PSCP system addresses all the above limitations while achieving a cost effective and objective means of higher accuracy for rapid cancer diagnosis and prognosis.
The issue of overtreatment of prostate cancer has received wide public attention demonstrating a need for better prognostication at time of diagnosis. This challenge results partially from the limited sampling of the cancer and therefore limited routine microscopic clues on the diagnostic needle biopsy. Evolving efforts are being made to determine a gene based score to address this need. With the state-of-the-art risk scoring based on co-registered structural and molecular signatures, PSCP offers better prognosis for determination of need for therapy. Identifying and treating only those cancers which will become clinically significant will considerably reduce the treatment costs and the economic burden of cancer.
There is, therefore, a continuing need for improved cancer and precancerous detection, diagnosis, monitoring, and prognosis techniques. More specifically, there is a continuing need for cancer diagnostic techniques that utilize fresh or frozen tissue sections without requiring time consuming prior staining or other labeling techniques prior to analysis of the tissue sections.
The present invention provides a method for quantifying a microstructure of a biological tissue by utilizing light interaction, which comprises steps of imaging tissue sections with a phase imaging microscope to obtain a Q-DIC image, and analyzing the Q-DIC image to obtain a two-dimensional Q-DIC quantitative map, and thereby the obtained accurate quantitative map can be used for detecting somatic functions of human bodies, including but not limited to detection, diagnosis and prognosis of pathologies and cancers, with an advantage of being quick and accurate.
The present invention also provides a method for quantifying a molecular content of a biological tissue by utilizing the auto-fluorescence thereof, which comprises steps of imaging a tissue section with a fluorescence imaging microscope to obtain a fluorescence microscope image, and computing and factorizing the fluorescence microscope image by the following equation to obtain a fluorescent component data map:
(wherein, I is the fluorescence intensity, w is the weight factor, and C is the content of a fluorescent substance)
With this method, absolute content of fluorescent molecules and a distribution map thereof can be obtained for quick and accurate detection of somatic functions of human bodies, including but not limited to detection, diagnosis and prognosis of pathologies and cancers.
The present further provides a photonic structural and chemometric pathology system for cancer detection, diagnosis, monitoring and/or prognosis that utilizes standard pathology sections without prior staining or other labeling techniques. The same unlabeled tissue section may be imaged, in some cases with the same phase imaging microscope, to obtain phase differential contrast (Q-DIC) images and fluorescence images. The Q-DIC images are analyzed to generate two dimensional Q-DIC quantitative maps, such as morphology, cell mass, and scattering data maps. The fluorescence images are analyzed to generate fluorescence intensity, fluorescence molecular content and metabolism data maps. The combination of Q-DIC and fluorescence intensity digital data maps are then analyzed individually and comparatively to perform cancer diagnosis and prognosis. The system is useful for a wide range of cancers including solid and fibrous tumors and is particularly well suited for providing detection, diagnosis, monitoring and/or prognosis for prostate cancer, lung cancer, breast cancer, etc.
The Q-DIC data maps and fluorescence intensity digital image maps are co-registered, for example (but not necessarily) on a pixel-by-pixel basis, to facilitate image comparison. The co-registered digital image maps may be visually compared side-by-side and in overlay format. Computational analysis may also be performed to data maps individually and comparatively. A prior diagnostic database of Q-DIC data maps and fluorescence intensity maps correlated to known pathology is analyzed to produce a risk assessment scoring methodology, which is applied to patient Q-DIC data maps and fluorescence intensity maps to obtain cancer detection, diagnosis monitoring and/or prognosis.
This approach has the advantage of producing diagnostic and prognostic results within a matter of minutes using fresh or frozen tissue sections that do not have to be stained or otherwise labeled for analysis. Appropriate Q-DIC images and fluorescent images can be obtained from standard pathology sections using a conventional phase and fluorescence imaging microscope. A software application containing the image analysis and risk scoring methodology can be installed and run on any suitable computer, such as a desktop, laptop, server or other computing device. In many cases, the diagnostic system may be cost effectively installed in a medical office, lab, or hospital to provide cancer diagnostic and prognostic results within a matter of minutes during a single patient appointment. This approach also produces a more objective method for cancer detection by removing elements of subjective analysis from interpretations by a human pathologist.
The patent application contains at least one drawing executed in color. Copies of this patent with color drawing(s) will be provided by the Patent and Trademark Office upon request and payment of the necessary fee.
Embodiments of the invention may be realized in a photonic structural and chemometric pathology (PSCP) system for cancer detection, diagnosis, monitoring and/or prognosis that utilizes standard pathology sections from fresh or frozen human tissue without prior staining or other labeling techniques. The appended figures along with the attached slide presentation, which form part of this disclosure, are summarized and referenced below.
The PSCP system is a novel cost-effective, label-free, quantitative, real-time tissue structural and molecular pathology tool for cancer diagnosis and prognosis. Current histopathology, which mainly utilizes the degree of loss of the normal glandular tissue architecture for cancer diagnosis, is subjective and suffers from inter- and intra-observer variations and has difficulty in distinguishing aggressive cancers from indolent ones. PSCP utilizes alternate clues including alterations in tissue structure and cellular metabolism to achieve an objective risk stratification of cancer grade and aggressiveness and may prove a reliable adjunct or an independent system for rapid diagnostic and prognostic assessment.
The PSCP system provides a novel, cost-effective, label-free, quantitative real-time tissue structural and molecular pathology for cancer diagnosis and prognosis. The system integrates quantitative phase imaging and tissue native fluorescence imaging to simultaneously quantify at the cellular and subcellular level the morphology, mass distribution, scattering characteristics, and metabolism of standard, unlabeled pathology tissue sections. Co-registered two dimensional (2D) microscopic maps of tissue morphology, mass distribution, scattering characteristics, and distribution of native fluorescence molecules and metabolism of cancer tissue sections are correlated with histopathology and patient outcomes. Based on these correlations, a risk stratification algorithm based on a state-of-the-art support vector machine (SVM) for identifying the risk score for a patient is utilized to grade cancer and predict patient outcome. The PSCP system integrates spatially-resolved tissue structural and molecular imaging for cancer diagnosis and prognosis for the first time to overcome major challenges faced by conventional pathology including inter- and intra-observer variations related to the subjective diagnosis based on mainly tissue morphology and the difficulty of distinguishing aggressive cancers from indolent ones. Furthermore, PSCP is label-free and can achieve real-time diagnosis by directly analyzing fresh or frozen unstained tissue.
The same tissue slide can be used to create both Q-DIC and fluorescence images. This feature is available in the conventional phase imaging microscope available in many medical offices, labs and hospitals. The Q-DIC and fluorescence images are input into a computer system 20, which may be any suitable general purpose computer configured with software to perform the innovative image processing and diagnostic methodology. More specifically, the computer system 20 includes a QDIC image processor 22 that computes a variety of two dimensional Q-DIC quantitative maps from the base Q-DIC images, such as morphology, cell mass, and scattering characteristic digital data maps. The computer system 20 also includes a fluorescence image processor 24, which computes a content of fluorescent molecules and metabolism digital data map from the fluorescence image. The Q-DIC, fluorescence and their derivative digital data maps are then passed to a diagnostic processor 26 for diagnostic and prognostic analysis.
The Q-DIC and fluorescence digital data maps may be co-registered on a pixel-by-pixel basis to adjust the sizes to provide maps of similar size for comparison. The digital data maps may also be placed side-by-side for visual analysis, placed in an overlay mode for visual analysis, and subjected to individual image and comparative computational analysis. The diagnostic processor 26 also includes a risk assessment scoring methodology and may also include a prior diagnostic database of Q-DIC and fluorescent component data maps for patients with known pathology used to create the risk assessment methodology. Regardless of whether the database itself is included on a particular embodiment, the risk assessment scoring methodology is applied to patient Q-DIC and fluorescent component data maps obtained from a particular tissue sample to produce diagnostic and prognostic results for the patient associated with the tissue sample.
In an illustrative embodiment, the pathology slide preparation station 12 may be any suitable type of slide preparation station and, in many cases, the slides arrive at the analysis location from cold storage or another location. The QDIC microscope may generally be an inverted brightfield, phase contrast, differential interference contrast, and fluorescence microscope with a focus drive.
In block 36, a fluorescence microscope is used to obtain a native fluorescence image of the tissue section. The Q-DIC image and fluorescence image are typically obtained using different settings of the same phase imaging microscope and a fluorescence filter for the fluorescence image.
In block 38, the Q-DIC image is processed to obtain a variety of Q-DIC quantitative maps, such as morphology, cell mass, and scattering characteristic digital data maps. Also in block 38, the fluorescence image is processed to obtain fluorescence intensity, native fluorescence molecular content and metabolism digital data maps. In block 40, a historical database of Q-DIC data maps and fluorescence intensity data maps for subjects with known pathology is analyzed to identify correlations between the map data and pathology diagnosis and prognosis. These correlations are typically reduced to a risk scoring methodology allowing higher risk cancers to be identified based on individual and comparative data map characteristics. The risk scoring methodology includes a range of predictive characteristics, combinations of characteristics, and overall profiles based on analysis of the historical database, which is refined through experience using the system. The effectiveness of the risk scoring methodology is expected to improve as the size and analytical detail of the database increases. It will therefore be beneficial to obtain a large initial data set and continually augment the database as experience with the system develops. The risk scoring methodology may also be adapted to score different types of cancer, different types of tissue, and other variations found to be effective through use and vetting of the system and associated database.
In block 42, the risk scoring methodology is applied to the data obtained from a particular tissue section to obtain specific diagnostic and prognostic results. Since the system utilizes standard pathology sections imaged with a standard phase imaging microscope analyzed with software configured on a general purpose computer, specific diagnostic and prognostic results can be obtained for a patient within minutes utilizing widely available equipment, in many cases within the same appointment where the tissue sample is obtained.
The PSCP system thus integrates quantitative phase imaging and tissue native fluorescence imaging to simultaneously quantify at the cellular and subcellular level the morphology, mass distribution, scattering characteristics, native fluorescence molecular content and metabolism of tissue without any prior tissue staining or other label processing. Co-registered two dimensional (2D) microscopic tissue structural and molecular maps of cancer tissue sections are correlated with histopathology and patient outcomes. Based on these correlations, a risk stratification algorithm for identifying the risk score for a patient is utilized to grade cancer and predict patient outcome. The PSCP system thereby integrates spatially-resolved tissue structural and molecular imaging for cancer diagnosis and prognosis for the first time.
The PSCP approach uses light scattering to sense the sub-wavelength changes in the microenvironment and probes native fluorescent molecules and cellular metabolism based on the auto-fluorescence of biological tissues. Light scattering has been demonstrated to be sensitive to nuclear and cellular alterations as small as a few nanometers and has been successfully applied to detect early carcinoma undetectable by conventional methods. Tissue native fluorescence has also been shown to be an effective probe for tissue physiological state, be it oxidative stress, metabolism, cell growth and cell survival or death.
The PSCP approach extracts objectively the 2D maps of the inherent tissue 3D structural and molecular characteristics directly from light interaction with intact tissue, avoiding subjective human interpretation and post-image pattern analysis (in digital pathology) of microscopic images of stained samples. Compared with auto-fluorescence spectroscopy, the PSCP system provides the 2D microscopic chemometric maps of absolute concentrations of native fluorescent molecules, enabling spatially resolved analysis of molecular activity and leading to significantly improved contrast. The spatially co-resolved tissue structure and molecular activities combining the high sensitivity of light scattering to the former and the high sensitivity of tissue native fluorescence to the latter may lead to more accurate cancer diagnosis and prognosis. This allows the PSCP approach to overcome major challenges faced by conventional pathology including inter- and intra-observer variations inherent in the subjective diagnosis and the difficulty of distinguishing aggressive cancers from indolent ones.
Furthermore, the PSCP approach can be applied to not only standard unstained formalin-fixed paraffin-embedded tissue sections but also fresh or frozen tissue sections without any prior tissue processing. It can achieve real-time diagnosis and be potentially adapted for in-vivo applications. PSCP provides a reliable adjunct or an independent system for rapid diagnostic and prognostic assessment.
It is estimated that more than 1.6 million men and women were diagnosed with cancer and more than a half million men and women died of cancer of all sites in 2012. Prostate cancer (PCa) is the most commonly diagnosed cancer in men. Diagnosis and prognosis of prostate tumors currently relies on the histopathologic evaluation of the biopsy together with clinical findings, such as clinical stage and prostate-specific antigen (PSA) to determine treatment. PCa is among cancers with a well-developed cancer grading system. The Gleason score (GS), the most widely used method of prostate cancer tissue grading today, is the single most important prognostic factor in PCa. It is one determinant of a patient's specific risk of dying due to PCa and strongly influences decisions regarding options for therapy.
The Gleason grading system is based upon the degree of loss of the normal glandular tissue architecture (i.e., shape, size and differentiation of the glands). The subjective microscopic determination of this loss of normal glandular structure caused by the cancer is represented by a Gleason grade (GG), a number ranging from 1 to 5, with 5 being the worst grade possible. The Gleason score (GS) is the sum of GG of the predominant and the next most predominant patterns within prostate cancer tissue. The majority of PCa patients are now diagnosed with tumors of GS 6 and 7. These tumors can take two distinct disease courses—indolent or highly aggressive, leading to death if not treated. Clinicians and patients choose a primary treatment modality based largely on the score, such as surgery, radiation therapy, or cryosurgery with the accompanying morbidity and compromised quality of life; or watchful waiting, risking life with delayed treatment. However, Gleason grading and its score is largely subjective. There exists the possibility of inter-observer as well as intra-observer variations, especially when based on the small thin needle core tissue samples obtained at biopsy. Objective sub-stratification of prostate cancer, in particular, GS 6 and 7 cancers, a prognostically heterogeneous group of tumors, indolent or aggressive groups is much desired. It has been shown that other factors including the amount of Gleason grades 4 and 5, the nuclear structure alteration, and cycle proliferation genes relate to the overall prognosis of the patient. Prostate tumor heterogeneity further complicates use of Gleason grading based on limited samples to evaluate patient prognosis.
Prediction of biologic behavior based on Gleason scores, especially in large volume or multifocal tumor, is limited by the scoring system, which has its own inherent difficulties due to reliance on the volume of tumor. A system that may utilize alternate clues, such as tissue microenvironment and physiological state, may compensate for the great heterogeneity in Gleason scores across individual tumors and prove a more reliable adjunct or even an independent system for diagnostic and prognostic assessment.
The PSCP system addresses these challenges. PSCP integrates quantitative tissue structural and molecular imaging for cancer diagnosis and prognosis for the first time to simultaneously quantify at the cellular/subcellular level important hallmarks of carcinogenesis, including alterations in the tissue microenvironment and cellular metabolism. The tissue microenvironment alterations accompanying carcinogenesis affect the spatial distribution of the refractive index of tissue, resulting in changes in its light scattering characteristics. PSCP uses light scattering to determine the sub-wavelength alterations in the microenvironment. Light scattering has been shown to be highly sensitive to the microarchitecture and composition of tissue and has been used to detect and evaluate sub-wavelength morphologic and biochemical changes in early carcinoma undetectable by other methods.
Quantitative phase imaging has concurrently emerged as a powerful light scattering technique to permit nanometer-scale measurements of structure and motion and to generate a spatially resolved sample scattering characteristics map in a noncontact, non-invasive manner. Another important hallmark of carcinogenesis is alterations in cellular metabolism. Cancer cell metabolism is often shifted from oxidative phosphorylation to aerobic glycolysis as the primary generator of cellular ATP. PSCP probes native fluorescent molecules and cellular metabolism based on auto-fluorescence of biological tissues. Auto-fluorescence has been demonstrated to be an effective probe for tissue physiological state. Monitoring cellular metabolism with auto-fluorescence has been successfully applied to differentiate normal, dysplasia, and cancer in vivo and ex vivo with high sensitivity and specificity.
The tight integration and co-registration of spatially resolved tissue structural and molecular imaging in PSCP offers significant advantages over existing approaches, such as digital pathology or auto-fluorescence spectroscopy. PSCP extracts objectively the inherent tissue sub-wavelength structural and molecular characteristics directly from light interaction with intact tissue, generating co-registered 2D microscopic maps of tissue morphology, mass distribution, scattering characteristics, and metabolism of biological tissue. PSCP also provides 2D microscopic chemometric maps of absolute concentrations of native fluorophores, enabling spatially resolved analysis of cellular metabolism and leading to significantly improved contrast. Test results indicate that tissue structural and molecular imaging complement each other and that the nuclear optical path length (OPL) and scattering characteristics are correlated closely with the aggressiveness of the cancer. The concentration of fluorophores and their ratios, as markers for cellular metabolism, also correlate closely with cancer grade. The combination of the high sensitivity advantage of light scattering to tissue structure and the high sensitivity advantage of co-registered native fluorescence to tissue physiological state (cellular metabolism in particular) in PSCP leads to more accurate cancer diagnosis and prognosis.
Furthermore, PSCP is label-free and achieves cost effective, objective real-time diagnosis and prognosis. In conventional approaches, the standard histopathological examination of the cancer biopsy requires multiple tissue processing steps. It is labor intensive and time consuming with frozen tissue preparation requiring 20-45 min for each tissue sample and routine permanent/fixed tissue sections typically requiring 12-18 hours processing and additional sectioning and staining steps. The time for expert manual assessment by conventional light microscopy may range from 5 minutes to 1 hour depending on the number of slides per specimen part. If there is diagnostic difficulty, immunohistochemistry may be required which can take an additional 24 hours to be performed and add a further 5 minutes to 1 hour for expert manual assessment. The cost for expert manual assessment may be determined based on salaries for a board qualified pathologist or by insurance charges for professional interpretation. In comparison, PSCP can be applied to not only the standard unstained formalin-fixed paraffin-embedded tissue sections but also fresh or frozen tissue specimens and hence eliminates the need for any prior tissue processing. The imaging speed of PSCP is mainly limited by the exposure time needed to record tissue native fluorescence images as the time for data processing utilizing currently available general purpose computer systems is minimal.
A complete PSCP scan can be completed within 2 to 3 minutes over 15 mm×15 mm area of a standard slide (under 20× objective) once the PSCP system is optimized and fully automated. This time can be further reduced using a more sensitive digital camera. The operating cost of PSCP is minimal compared to standard pathology once the system is purchased and setup. One significant advantage of a machine based system is that it will not suffer from day to day variations and will not be subject to human error missing a small focus on a slide. One advantageous use of PSCP is to prescreen the slides automatically and highlight suspicious regions and cases for further expert assessment to reduce the burden on pathologists. Another important advantage of PSCP lies in prognostication.
PSCP integrates quantitative spatially-resolved tissue structural and molecular imaging for the first time for objective cancer diagnosis and prognosis on intact tissue specimens, which simultaneously quantifies at the cellular and subcellular level important hallmarks of carcinogenesis. Varying tissue components differ in refractive index and render different extents of phase delay when light passes through. The quantitative phase map in PSCP may resolve nanometer axial optical path differences. Minute alterations in tissue structure, undetectable by the conventional pathology, are revealed. The tissue structure is quantified in detail beyond the scope of morphology. The cellular mass distribution is given by the phase map after scaling. The tissue scattering characteristics, such as the reduced scattering coefficient and the light scattering anisotropy factor maps computed via the scattering-phase theorem, are utilized to quantify the tissue local microenvironment signature. For example, a small “g” suggests that the size of the scattering tissue is small and the local region is more heterogeneous (such as fragmentation in nuclei) whereas a large “g” means that the size of the scattering tissue is large and the local region tends to be more homogenous. The nuclear OPL is a marker for the cell cycling rate. Both nuclear OPL and g have been found to correlate strongly with cancer aggressiveness. In addition, the co-registered chemometric imaging enables spatially resolved analysis of cellular metabolism, which leads to significantly improved contrast. As the tissue structural and molecular maps are complementary to each other, simultaneous measurement of both achieves higher accuracy in cancer diagnosis and prognosis in comparison to prior techniques.
Chemometric tissue native fluorescence imaging provides the co-registered 2D microscopic maps of absolute concentrations of native fluorescent molecules. The primary oxidation-reduction (redox) reactions in cells to generate energy in the form of ATP are the conversion of NAD+ to its reduced form NADH and the oxidation of flavin adenine dinucleotide (FAD) to FADH2, a process known as oxidative phosphorylation. Cancer cell metabolism is often shifted from oxidative phosphorylation to aerobic glycolysis as the primary generator of cellular ATP. This shift has long been known as the “Warburg effect” and results in a net increase in NADH. The tryptophan and NADH concentrations increase consistently as the cells progress from nonmetastatic/low metastatic to metastatic state. These molecules including tryptophan, NADH and FAD are auto-fluorescent and hence tissue native fluorescence has been extensively investigated for diagnostic purposes. The challenges include: 1) tissue auto-fluorescence is weak and 2) auto-fluorescence is significantly modified by tissue morphology and hemoglobin absorption inside a bulk tissue. In PSCP, a RGB digital camera is used to record auto-fluorescence images of a slide containing a thin (e.g., approximately 4-5 microns thick) slice of tissue excited at UV and blue bands. The absolute concentration map of fluorophores is then computed with non-negative matrix factorization (NMF) of the measured RGB channels and system calibration. By avoiding averaging over the whole interrogated volume and providing spatially resolved concentrations of molecules, the contrast between tissues at different states is significantly enhanced. The shortcomings associated with the conventional tissue native fluorescence spectroscopy are thus addressed by PSCP, which possess the additional advantages a simple, label-free, real-time approach to chemometric mapping and metabolism assessment of the tissue section.
Co-registered quantitative tissue structural and molecular imaging together with a state-of-the-art support vector machine achieves high accuracy in cancer diagnosis and prognosis. Tissue structural and molecular imaging complements each other. Test results show that the nuclear OPL and g maps correlate strongly with biochemical recurrence of prostate cancer, whereas the optical metabolism parameters correlate more closely with cancer. The co-registered structural and molecular maps also open up novel avenues in data analysis to achieve more accurate risk scoring for patients. For example, the computation of fluorescence strength of native fluorophores may be refined by excluding the nonfluorescent nuclear region identified from the phase map. More importantly, the co-registered maps enable the characterization of morphology, mass distribution, scattering characteristics, and metabolism for each individual cell on the slide. The wealth of information on the tissue average characteristics, the single cell level statistics, and the spatial distribution characteristics are used as the input to a state-of-the-art support vector machine (SVM) to develop objective and accurate cancer risk scoring for patients. With these innovations, the performance of PSCP is expected to improve further over the already very promising preliminary results.
An illustrative PSCP microscopic imaging system may be implemented on an inverted epi-fluorescence differential interference contrast (DIC) microscope (Axiovert 40cfl, Zeiss). For quantitative phase imaging, the light source is a Halogen 35 W lamp filtered by narrow-band filters with wavelengths selectable from 400 nm to 700 nm at a step size of 50 nm and measurement at one single wavelength (500 nm) is normally used.
The numerical apertures for the condenser and objective (APlan) are 0.2 and 0.5, respectively. For epi-fluorescence measurement, the light source is a 50 W mercury arc lamp. The microscope is equipped with a focus drive (TOFRA, California). The DIC images are recorded by a low noise CCD camera (INFINITY2-1R, Lumenera) and fluorescence images are recorded by a digital camera (Canon 5D Mark II). A motorized stage may be used to scan the specimen. Custom codes in Matlab (Mathworks) may be utilized to automate the whole system from scanning and focusing the specimen, recording DIC and fluorescence images, to performing data analysis and finally risk-stratification of cancer.
More specifically, in an illustrative embodiment, the Halogen 35 W lamp is filtered by narrow-band filters with wavelengths selectable from 400 nm to 700 nm at a step size of 50 nm and measurement at one single wavelength (500 nm) is used. The numerical apertures for the condenser and objective (×40) is 0.2 and 0.5 respectively. For epi-fluorescence measurement, the light source is a 50 W mercury arc lamp. The microscope is equipped with a focus drive (TOFRA, California). The DIC images are recorded by a low noise CCD camera (INFINITY2-1R, Lumenera) and fluorescence images are recorded by a digital camera (Canon 5D Mark II). A motorized stage is used to scan the specimen. Custom codes in Matlab (Mathworks) automate the overall system from scanning/focusing the specimen, recording DIC and fluorescence images, to performing data analysis and finally risk-stratification of cancer.
A specific example is described below with theoretical development.
Equation (1) enables quantitative phase imaging with a DIC microscope. From Eq. (1), the Laplacian of the phase map can be approximated by Equation (2) shown in
Tissue scattering is characterized by the scattering-phase relationship. A set of equations relates the scattering coefficient μs, the reduced scattering coefficient μ′s, and the anisotropy factor g (normally 0<g<1 for tissue) of light scattering by a bulk turbid medium to the property of a phase map ϕ(ρ) of light transmission. These relationships are termed as the scattering-phase theorem equations (3), (4) and (5) shown in
The scattering-phase theorem is applicable to a slice of homogeneous or inhomogeneous medium. In the latter case, a map of μs, ν′s and g are computed from the phase map using spatial averaging over the local region rather than the whole slice. The squared phase gradient map is computed by the finite difference of the phase map. A proper correction is applied for light diffraction in the microscope in the computation of the phase gradient.
The phase imaging in PSCP may be performed at one single wavelength (e.g., 500 nm) or at multiple wavelengths. When the phase imaging is measured at multiple wavelengths, the wavelength-dependent scattering characteristics can be used to compute the spectral slope of μ′s. The fractal dimension of the scattering structure can be computed from this spectral slope and there is a significant difference in fractal dimension between normal and cancerous tissue utilizing the discriminating power of the spectral slope and the fractal dimension derived from multiple wavelength phase measurements as additional input parameters to the SVM.
Referring to
The first column 802 of images shows the computed OPL maps for the two samples. The scattering property for each individual spheres can be analyzed by applying the scattering-phase theorem to the region in the phase map occupied by the sphere. For example, the region highlighted by white dash lines for the central sphere yields μs=0.234 μm−1, μ′s=0.0202 μm−1, and g=0.91 with an area 61.0 μm2, corresponding to scattering and reduced scattering cross sections of 118 μm2 and 10.2 μm2, respectively. The mean scattering and reduced scattering cross sections for all the spheres contained in the section are 116 μm2 and 9.8 μm2, respectively. These values are in excellent agreement with the theoretical prediction for a polystyrene sphere of the specified size (Csca=125 μm2, C′sca=9.5 μm2 and g=0.92) at the probing wavelength 550 nm predicted by Mie theory.
Referring to table 900 shown in
Specific contributing fluorophores for tissue native fluorescence include: aromatic amino acid residues found in most proteins (λex=200-340 nm, λm=360-370, 455 nm); reduced pyridine nucleotides (NAD(P)H), which are cofactors in cellular metabolism and are found mainly in mitochondria, but are also present in the cytoplasm (λex=360 nm, λem=460 nm); and flavins and flavin nucleotides (riboflavin, flavin mononucleotide, and dinucleotide), which are mostly bound to enzymes as coenzymes of flavoproteins and concentrated in the mitochondria (λex=360,445-470 nm, λem=440,520 nm).
Endogenous porphyrins and age related lipofuscin have also been shown to be sources of cellular auto-fluorescence. Extracellular matrix components collagen and elastin also contribute strongly to tissue native fluorescence.
Referring to
The weights are normalized by Σwij=1 which ensures that the total fluorescence intensity map (summing over all channels from all excitations) produced by the jth component is Cj. The weights wij can be regarded as the spectral signature of the jth component. The weights can be estimated from measurement of pure fluorescent components under the identical experimental condition noting that the weights may differ slightly between pure fluorophores and fluorophores in tissue owing to their different chemical environments. For example, a non-negative matrix factorization (NMF) of the measured RGB channels may be utilized to extract both the weights and the fluorescence map directly from the experimental data. The fluorescence map Cj is then identified to be one of the main species of fluorophores with the spectral signature matching and a priori information on their spatial distribution characteristics. The absolute concentration and distribution of native fluorescent molecules is then given by Cj/Kj where Kj is a calibrated constant representing the total fluorescence intensity produced by one mole of the pure jth species of fluorophore under the identical experimental condition.
This chemometric microscopy quantifies the absolute concentration and distribution of native fluorophores including tryptophan, collagen, NADH, FAD, and porphyrin simultaneously for the first time. The redox ratios such as FAD/NADH and FAD/Tryptophan are obtained by dividing two respective concentration maps. The concentration ratio rather than the fluorescence intensity ratio may be reported as the former does not depend on the measurement condition. Note that strong signals from tryptophan may be observed from both the pure tryptophan sample and the tissue specimens under UV excitation although the illumination light is off its excitation peak (280 nm).
The optical parameters (including tissue average values and the single cell level statistics of at least nuclear OPL, the reduced scattering coefficient, the anisotropy factor, the native fluorophore concentration, and redox ratios) and the clinical parameters (Gleason scores and biochemical recurrence) obtained from both methods are collected and tabulated. The PSCP system correlates these factors with clinical detection of metastasis (usually by radiologic methods—bone scan or MRI/CT scan), and death of disease. The prolonged course of prostate cancer is such that PSA recurrence represents the earliest and most sensitive measure for disease recurrence and is the current standard for most studies. The clinical presence of metastases is an additional strong marker for aggressive disease. However in the clinical setting this is often not sought or documented in the face of ever rising PSA, which is sufficient evidence to warrant treatment. Although death from disease is tracked, these data are often unreliable in the absence of autopsy documentation. Nevertheless the system correlates results with all three markers of disease recurrence 1) PSA recurrence, 2) clinical evidence of metastasis, and 3) death of disease. By comparing the two sets of optical and clinical parameters, the correlation between them is established.
One unique advantage of using TMA and fresh cryopreserved prostate tissue sections from surgical specimens is that these samples are from large cohorts of prostate cancer cases with different Gleason scores and have follow up data covering a mean of 7 years post resection, essential for the development of the PSCP risk score for patient outcome. There is no statistically significant difference between the optical properties of fresh and cryopreserved tissue. Phase and native fluorescence measurement on 191 cases contained on a prostate cancer TMA cohort contains 19 adjacent normal sites, 65 GS6 sites, 82 GS7 sites, 23 GS8 sites, and 8 GS9 sites among which there are a total of 123 cases of PSA non-recurrence (NRec) and 44 cases of PSA recurrence (Rec) cases (either Gleason Score or PSA recurrence status of some cases is unknown).
Table 1300 of
The optical parameters presented in Tables 1200 and 1300 are tissue average characteristics. The single cell level statistics on these parameters may further reveal tissue structural and molecular alterations with cancer grade and aggressiveness and achieve even higher correlations between PSCP and the clinical parameters and patient outcome. This establishes a solid foundation for objective and accurate cancer risk scoring for patients.
A state-of-the-art support vector machine (SVM) technique may be utilized to develop the risk score for prostate cancer where the combination of SVM and a flow peaks algorithm yields high accuracy. The decision function (or discriminant function) from the trained SVM can be utilized as the risk score. When the outcome variable such as tumor grade has more than two categories, a one-against one approach may be applied to all possible two group comparisons.
Based on the correlation study on the cohorts of prostate cancer with the measured parameters, a PSCP risk score for prostate cancer from optically retrieved tissue structural and molecular maps is computed. The risk score may be based on the decision function of the support vector machine with the variables selected by recursive feature elimination. A SVM of the tissue average characteristics and the single cell level statistics on morphology, mass distribution, scattering properties, fluorophore concentrations, and their ratios will be carried out for different prediction problems. The algorithm differentiates prostate cancer from normal prostate including benign prostate hyperplasia, predicts tumor grade, and differentiates GS 6 from 7 tumors and sub-stratifies GS7 tumors. Based on this and information, the algorithm also assess cancer aggressiveness and predict patient outcome.
The system may also calculate receiver operator characteristics (ROC) curves and determine the diagnostic and prognostic performance using accuracy, sensitivity, and specificity for comparisons evaluated. A Gaussian kernel SVM of the tissue average characteristics alone may be computed by using the SVM functions. The accuracy, specificity and sensitivity are computed based on 10 random splits of five-fold cross-validations. The SVM using the Gleason score gives only 58.4±6.3 accuracy, 71.1±10.0 sensitivity, and 28.3±6.1 specificity. The best SVM utilizing a single variable can only achieve 69.6±1.2 accuracy, and all input variables 72.3±2.8 accuracy. Using the best subset of phase features selected by RFE including μ′s, nuclear OPL and g, the SVM achieves 79.5±1.3 accuracy, 79.5±1.5 sensitivity, and 79.3±5.2 specificity. By adding metabolism features, the SVM using the best subset of parameters selected by RFE yields accuracy 81.2±3.2, sensitivity 84.7±3.7, and specificity 73. 3±2.9. This preliminary result illustrates that the SVM utilizing carefully selected variables gives at least a ten-fold increase of accuracy over the best SVM utilizing a single variable and the inclusion of metabolism data is beneficial. The PSCP thus provides improvements in SVM, and more importantly, incorporates single-cell level statistics in addition to the tissue average characteristics into the SVM to produce a final risk scoring algorithm that achieves higher accuracy, sensitivity and specificity that prior approaches.
As another embodiment of the invention, the PSCP system is used for lung cancer tissue matrix array (TMA) (Biomax, USA) imaging. In this arrangement, the channels are UV-R, UV-G and UV-B (under the excitation of ultraviolet light), and are B-R and B-G (under the excitation of blue light). Under both of the above excitations, duration of exposure was fixed at 1.5 s. A set of aqueous suspensions containing a thin slice (˜4 μm thick) of monolayer uniform pure fluorophores (comprising tryptophan, collagen, elastin, NADH, FAD and porphyrin) were prepared and measured individually under the identical experimental condition.
For each specimen, the spatially-resolved fluorescence image based on block principal pivoting NMF algorithm after background removal is the product of W and H. Distribution of a single fluorescence component is included in row H, and fluorescence species are determined by characteristic matrix W.
Primary fluorescence species are determined to be Trp (tryptophan), elastin, NADH and FAD. The fluorescence intensities of the four main components differ significantly in normal lung tissue (n=26), perilesional tissue (n=34) and cancerous tissue (n=58) (in
Typical results in the resolved concentration maps, as shown in
Decrease in tryptophan concentration indicates accelerated decomposition and enhanced cell inducer degradation of tryptophan, similar to what is observed in serum of a patient with lung cancer. Elastin is found to be a main fluorescence in lung tissue. Compared with that in normal tissue, dramatic reduction of disordered fiberous web accompanied with fluorescence is observed in pre-cancerous lesion tissue. Pitzer, et al. found that cells transformed by cancerogenic substances exhibited obvious decrease in NADH and flavin and in dot fluorescence intensity, and that immortalized human bronchial epithelial cells lost spatial positioning of native fluorescence in comparison with human bronchial epithelial cells. The same trend is observed here (see Table 4 and
A potential problem for NMF is that it may lead to a non-unique solution. Such fuzziness is minimized by utilizing spectral signatures of known pure “fluorescence” as initial conjecture for block principal pivoting NMF algorithm of the characteristic matrix. By measuring synthesized data with noise added, stability of this approach is verified. Due to their similar spectral signatures, the reduced NADH and FAD exhibit a difference less than 8% under a noise less than 5%, and the reduced tryptophan and protein exhibit a difference less than 20% under a noise approximating 2%. In our experiments, noise in the measured fluorescence images is 1%.
There are some important benefits for obtaining absolute concentrations of fluorescence components instead of their relative content through individual quantization of fluorescence species. This eliminates potential uncertainty in system response and measurement details, which otherwise might influence the reported data and the observed fluorescence characteristics related to the physiological status of tissue. The objective comparison between different approaches is simplified. Most importantly, the absolute fluorophore concentration is recovered to be relevant biodata to facilitate data interpretation and tissue diagnosis improvement. Diagnosibility through relative fluorescence intensity is much weaker than that through absolute fluorescence component concentration. In fact, absolute concentrations can be used for classifying non-cancers (perilesional normalities, n=66) and tumors (n=52) more accurately than relative fluorescence intensity. The ROC graph shown in
The embodiments illustrated above are for diagnosis of prostate cancer and lung cancer. However, the invention is not limited thereto. The above approaches not only generate good results in diagnosis and prognosis of prostate cancer and lung cancer, but also possess the ability to detect and prognose other cancers and other pathologies. In some embodiments, a method for quantifying a microstructure of a biological tissue by utilizing light interaction obtains a Q-DIC image by imaging a tissue section with a phase imaging microscope; and then analyzes the Q-DIC image, wherein the reduced scattering coefficient μ′s, optical path length OPL and anisotropy factor g obtained by the above equations 1-5 may be used for detection of human bodies and other living bodies, including but not limited to detection, diagnosis and prognosis of pathologies and cancers. In some embodiments, a method for quantifying a molecular content of a biological tissue by utilizing the auto-fluorescence thereof images a tissue section with a fluorescence imaging microscope to obtain a fluorescence microscope image, then compute the fluorescence microscope image by the above equation (6) to obtain an absolute concentration data map of fluorescence components, and then obtains data on tryptophan, collagen, NADH, FAD and porphyrin. The method can also be used for detection of human bodies and other living bodies, including but not limited to detection, diagnosis and prognosis of pathologies and cancers.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2016/090755 | 7/20/2016 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2017/012555 | 1/26/2017 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
20050065440 | Levenson | Mar 2005 | A1 |
20110117025 | Dacosta et al. | May 2011 | A1 |
20110254943 | Ozinsky | Oct 2011 | A1 |
20140050386 | Humayun | Feb 2014 | A1 |
20140119628 | Elad et al. | May 2014 | A1 |
20140205057 | Koehler et al. | Jul 2014 | A1 |
20140268168 | Feldman et al. | Sep 2014 | A1 |
20160266126 | Shipitsin | Sep 2016 | A1 |
Number | Date | Country |
---|---|---|
1317688 | Oct 2001 | CN |
102099671 | Jun 2011 | CN |
102266219 | Dec 2011 | CN |
103918005 | Jul 2014 | CN |
Entry |
---|
Kong et al Phase contrast microscopy of living cells within the whole lens: spatial correlations and morphological dynamics Received Mar. 18, 2012 | Accepted Jul. 28, 2012 | Published Aug. 1, 2012 (Year: 2012). |
Number | Date | Country | |
---|---|---|---|
20200143531 A1 | May 2020 | US |
Number | Date | Country | |
---|---|---|---|
62194462 | Jul 2015 | US |