1
|
Bobholz SA, Lowman AK, Connelly JM, Duenweg SR, Winiarz A, Nath B, Kyereme F, Brehler M, Bukowy J, Coss D, Lupo JM, Phillips JJ, Ellingson BM, Krucoff MO, Mueller WM, Banerjee A, LaViolette PS. Noninvasive Autopsy-Validated Tumor Probability Maps Identify Glioma Invasion Beyond Contrast Enhancement. Neurosurgery 2024:00006123-990000000-01091. [PMID: 38501824 DOI: 10.1227/neu.0000000000002898] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2023] [Accepted: 01/09/2024] [Indexed: 03/20/2024] Open
Abstract
BACKGROUND AND OBJECTIVES This study identified a clinically significant subset of patients with glioma with tumor outside of contrast enhancement present at autopsy and subsequently developed a method for detecting nonenhancing tumor using radio-pathomic mapping. We tested the hypothesis that autopsy-based radio-pathomic tumor probability maps would be able to noninvasively identify areas of infiltrative tumor beyond traditional imaging signatures. METHODS A total of 159 tissue samples from 65 subjects were aligned to MRI acquired nearest to death for this retrospective study. Demographic and survival characteristics for patients with and without tumor beyond the contrast-enhancing margin were computed. An ensemble algorithm was used to predict pixelwise tumor presence from pathological annotations using segmented cellularity (Cell), extracellular fluid, and cytoplasm density as input (6 train/3 test subjects). A second level of ensemble algorithms was used to predict voxelwise Cell, extracellular fluid, and cytoplasm on the full data set (43 train/22 test subjects) using 5-by-5 voxel tiles from T1, T1 + C, fluid-attenuated inversion recovery, and apparent diffusion coefficient as input. The models were then combined to generate noninvasive whole brain maps of tumor probability. RESULTS Tumor outside of contrast was identified in 41.5% of patients, who showed worse survival outcomes (hazard ratio = 3.90, P < .001). Tumor probability maps reliably tracked nonenhancing tumor on a range of local and external unseen data, identifying tumor outside of contrast in 69% of presurgical cases that also showed reduced survival outcomes (hazard ratio = 1.67, P = .027). CONCLUSION This study developed a multistage model for mapping gliomas using autopsy tissue samples as ground truth, which was able to identify regions of tumor beyond traditional imaging signatures.
Collapse
Affiliation(s)
- Samuel A Bobholz
- Department of Radiology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Allison K Lowman
- Department of Radiology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Jennifer M Connelly
- Department of Neurology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Savannah R Duenweg
- Department of Biophysics, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Aleksandra Winiarz
- Department of Biophysics, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Biprojit Nath
- Department of Biophysics, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Fitzgerald Kyereme
- Department of Radiology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Michael Brehler
- Department of Radiology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - John Bukowy
- Department of Electrical Engineering and Computer Science, Milwaukee School of Engineering, Milwaukee , Wisconsin , USA
| | - Dylan Coss
- Department of Pathology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Janine M Lupo
- Department of Radiology and Biomedical Imaging, University of California, San Francisco , California , USA
- UCSF/UC Berkeley Graduate Program in Bioengineering, University of California, San Francisco and Berkeley , California , USA
| | - Joanna J Phillips
- Department of Neurological Surgery, University of California, San Francisco , California , USA
- Department of Pathology, University of California, San Francisco , California , USA
| | - Benjamin M Ellingson
- UCLA Brain Tumor Imaging Laboratory, Department of Radiological Sciences, David Geffen School of Medicine, University of California Los Angeles, Los Angeles , California , USA
| | - Max O Krucoff
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Wade M Mueller
- Department of Neurosurgery, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Anjishnu Banerjee
- Department of Biostatistics, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| | - Peter S LaViolette
- Department of Radiology, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
- Department of Biophysics, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
- Department of Biomedical Engineering, Medical College of Wisconsin, Milwaukee , Wisconsin , USA
| |
Collapse
|
2
|
Duenweg SR, Brehler M, Lowman AK, Bobholz SA, Kyereme F, Winiarz A, Nath B, Iczkowski KA, Jacobsohn KM, LaViolette PS. Quantitative Histomorphometric Features of Prostate Cancer Predict Patients Who Biochemically Recur Following Prostatectomy. J Transl Med 2023; 103:100269. [PMID: 37898290 PMCID: PMC10872376 DOI: 10.1016/j.labinv.2023.100269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Revised: 10/10/2023] [Accepted: 10/19/2023] [Indexed: 10/30/2023] Open
Abstract
Prostate cancer is the most commonly diagnosed cancer in men, accounting for 27% of the new male cancer diagnoses in 2022. If organ-confined, removal of the prostate through radical prostatectomy is considered curative; however, distant metastases may occur, resulting in a poor patient prognosis. This study sought to determine whether quantitative pathomic features of prostate cancer differ in patients who biochemically experience biological recurrence after surgery. Whole-mount prostate histology from 78 patients was analyzed for this study. In total, 614 slides were hematoxylin and eosin stained and digitized to produce whole slide images (WSI). Regions of differing Gleason patterns were digitally annotated by a genitourinary fellowship-trained pathologist, and high-resolution tiles were extracted from each annotated region of interest for further analysis. Individual glands within the prostate were identified using automated image processing algorithms, and histomorphometric features were calculated on a per-tile basis and across WSI and averaged by patients. Tiles were organized into cancer and benign tissues. Logistic regression models were fit to assess the predictive value of the calculated pathomic features across tile groups and WSI; additionally, models using clinical information were used for comparisons. Logistic regression classified each pathomic feature model at accuracies >80% with areas under the curve of 0.82, 0.76, 0.75, and 0.72 for all tiles, cancer only, noncancer only, and across WSI. This was comparable with standard clinical information, Gleason Grade Groups, and CAPRA score, which achieved similar accuracies but areas under the curve of 0.80, 0.77, and 0.70, respectively. This study demonstrates that the use of quantitative pathomic features calculated from digital histology of prostate cancer may provide clinicians with additional information beyond the traditional qualitative pathologist assessment. Further research is warranted to determine possible inclusion in treatment guidance.
Collapse
Affiliation(s)
- Savannah R Duenweg
- Departments of Biophysics, Medical College of Wisconsin, Milwaukee, Wisconsin
| | - Michael Brehler
- Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin
| | | | | | | | - Aleksandra Winiarz
- Departments of Biophysics, Medical College of Wisconsin, Milwaukee, Wisconsin
| | - Biprojit Nath
- Departments of Biophysics, Medical College of Wisconsin, Milwaukee, Wisconsin
| | | | | | - Peter S LaViolette
- Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin; Biomedical Engineering, Medical College of Wisconsin, Milwaukee, Wisconsin.
| |
Collapse
|
3
|
Kylies D, Zimmermann M, Haas F, Schwerk M, Kuehl M, Brehler M, Czogalla J, Hernandez LC, Konczalla L, Okabayashi Y, Menzel J, Edenhofer I, Mezher S, Aypek H, Dumoulin B, Wu H, Hofmann S, Kretz O, Wanner N, Tomas NM, Krasemann S, Glatzel M, Kuppe C, Kramann R, Banjanin B, Schneider RK, Urbschat C, Arck P, Gagliani N, van Zandvoort M, Wiech T, Grahammer F, Sáez PJ, Wong MN, Bonn S, Huber TB, Puelles VG. Expansion-enhanced super-resolution radial fluctuations enable nanoscale molecular profiling of pathology specimens. Nat Nanotechnol 2023; 18:336-342. [PMID: 37037895 PMCID: PMC10115634 DOI: 10.1038/s41565-023-01328-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/12/2022] [Accepted: 01/13/2023] [Indexed: 06/19/2023]
Abstract
Expansion microscopy physically enlarges biological specimens to achieve nanoscale resolution using diffraction-limited microscopy systems1. However, optimal performance is usually reached using laser-based systems (for example, confocal microscopy), restricting its broad applicability in clinical pathology, as most centres have access only to light-emitting diode (LED)-based widefield systems. As a possible alternative, a computational method for image resolution enhancement, namely, super-resolution radial fluctuations (SRRF)2,3, has recently been developed. However, this method has not been explored in pathology specimens to date, because on its own, it does not achieve sufficient resolution for routine clinical use. Here, we report expansion-enhanced super-resolution radial fluctuations (ExSRRF), a simple, robust, scalable and accessible workflow that provides a resolution of up to 25 nm using LED-based widefield microscopy. ExSRRF enables molecular profiling of subcellular structures from archival formalin-fixed paraffin-embedded tissues in complex clinical and experimental specimens, including ischaemic, degenerative, neoplastic, genetic and immune-mediated disorders. Furthermore, as examples of its potential application to experimental and clinical pathology, we show that ExSRRF can be used to identify and quantify classical features of endoplasmic reticulum stress in the murine ischaemic kidney and diagnostic ultrastructural features in human kidney biopsies.
Collapse
Affiliation(s)
- Dominik Kylies
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Research Center On Rare Kidney Diseases (RECORD), University Hospital Erlangen, Erlangen, Germany
| | - Marina Zimmermann
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Institute of Medical Systems Biology, Center for Biomedical AI (bAIome), Center for Molecular Neurobiology Hamburg (ZMNH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Fabian Haas
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Maria Schwerk
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Malte Kuehl
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Institute of Medical Systems Biology, Center for Biomedical AI (bAIome), Center for Molecular Neurobiology Hamburg (ZMNH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Michael Brehler
- Institute of Medical Systems Biology, Center for Biomedical AI (bAIome), Center for Molecular Neurobiology Hamburg (ZMNH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Jan Czogalla
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Lola C Hernandez
- Cell Communication and Migration Laboratory, Department of Biochemistry and Molecular Cell Biology (IBMZ), Center for Experimental Medicine, Hamburg, Germany
| | - Leonie Konczalla
- Department of General, Visceral and Thoracic Surgery, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Mildred Scheel Cancer Career Center HaTriCS4, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Yusuke Okabayashi
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Division of Nephrology and Hypertension, Department of Internal Medicine, The Jikei University School of Medicine, Tokyo, Japan
| | | | - Ilka Edenhofer
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Sam Mezher
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Hande Aypek
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Bernhard Dumoulin
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Hui Wu
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Smilla Hofmann
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Oliver Kretz
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Nicola Wanner
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Nicola M Tomas
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Susanne Krasemann
- Institute of Neuropathology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Markus Glatzel
- Institute of Neuropathology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Christoph Kuppe
- Institute of Experimental Medicine and Systems Biology and Division of Nephrology and Clinical Immunology, RWTH Aachen University Medical Faculty, Aachen, Germany
| | - Rafael Kramann
- Institute of Experimental Medicine and Systems Biology and Division of Nephrology and Clinical Immunology, RWTH Aachen University Medical Faculty, Aachen, Germany
| | - Bella Banjanin
- Department of Developmental Biology, Erasmus Medical Center, Rotterdam, The Netherlands
- Oncode Institute, Erasmus Medical Center Cancer Institute, Rotterdam, The Netherlands
| | - Rebekka K Schneider
- Department of Developmental Biology, Erasmus Medical Center, Rotterdam, The Netherlands
- Oncode Institute, Erasmus Medical Center Cancer Institute, Rotterdam, The Netherlands
- Institute for Cell and Tumor Biology, RWTH Aachen University, Aachen, Germany
| | - Christopher Urbschat
- Department of Obstetrics and Fetal Medicine, Division of Experimental Feto-Maternal Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Petra Arck
- Department of Obstetrics and Fetal Medicine, Division of Experimental Feto-Maternal Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Nicola Gagliani
- Department of General, Visceral and Thoracic Surgery, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- I. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Marc van Zandvoort
- Department of Genetics and Cell Biology, Maastricht University, School for Oncology and Reproduction GROW, School for Mental Health and Neuroscience MHeNS, and School for Cardiovascular Diseases CARIM, Maastricht University, Maastricht, The Netherlands
- Institute for Molecular Cardiovascular Research (IMCAR), RWTH Aachen University, Aachen, Germany
| | - Thorsten Wiech
- Institute of Pathology, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Florian Grahammer
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Pablo J Sáez
- Cell Communication and Migration Laboratory, Department of Biochemistry and Molecular Cell Biology (IBMZ), Center for Experimental Medicine, Hamburg, Germany
| | - Milagros N Wong
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
- Department of Pathology, Aarhus University Hospital, Aarhus, Denmark
| | - Stefan Bonn
- Institute of Medical Systems Biology, Center for Biomedical AI (bAIome), Center for Molecular Neurobiology Hamburg (ZMNH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Tobias B Huber
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany
| | - Victor G Puelles
- III. Department of Medicine, University Medical Center Hamburg-Eppendorf, Hamburg, Germany.
- Hamburg Center for Kidney Health (HCKH), University Medical Center Hamburg-Eppendorf, Hamburg, Germany.
- Department of Clinical Medicine, Aarhus University, Aarhus, Denmark.
- Department of Pathology, Aarhus University Hospital, Aarhus, Denmark.
| |
Collapse
|
4
|
Duenweg SR, Brehler M, Bobholz SA, Lowman AK, Winiarz A, Kyereme F, Nencka A, Iczkowski KA, LaViolette PS. Comparison of a machine and deep learning model for automated tumor annotation on digitized whole slide prostate cancer histology. PLoS One 2023; 18:e0278084. [PMID: 36928230 PMCID: PMC10019669 DOI: 10.1371/journal.pone.0278084] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 03/04/2023] [Indexed: 03/18/2023] Open
Abstract
One in eight men will be affected by prostate cancer (PCa) in their lives. While the current clinical standard prognostic marker for PCa is the Gleason score, it is subject to inter-reviewer variability. This study compares two machine learning methods for discriminating between cancerous regions on digitized histology from 47 PCa patients. Whole-slide images were annotated by a GU fellowship-trained pathologist for each Gleason pattern. High-resolution tiles were extracted from annotated and unlabeled tissue. Patients were separated into a training set of 31 patients (Cohort A, n = 9345 tiles) and a testing cohort of 16 patients (Cohort B, n = 4375 tiles). Tiles from Cohort A were used to train a ResNet model, and glands from these tiles were segmented to calculate pathomic features to train a bagged ensemble model to discriminate tumors as (1) cancer and noncancer, (2) high- and low-grade cancer from noncancer, and (3) all Gleason patterns. The outputs of these models were compared to ground-truth pathologist annotations. The ensemble and ResNet models had overall accuracies of 89% and 88%, respectively, at predicting cancer from noncancer. The ResNet model was additionally able to differentiate Gleason patterns on data from Cohort B while the ensemble model was not. Our results suggest that quantitative pathomic features calculated from PCa histology can distinguish regions of cancer; however, texture features captured by deep learning frameworks better differentiate unique Gleason patterns.
Collapse
Affiliation(s)
- Savannah R Duenweg
- Department of Biophysics, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Michael Brehler
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Samuel A Bobholz
- Department of Biophysics, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Allison K Lowman
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Aleksandra Winiarz
- Department of Biophysics, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Fitzgerald Kyereme
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Andrew Nencka
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Kenneth A Iczkowski
- Department of Pathology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| | - Peter S LaViolette
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
- Department of Biomedical Engineering, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America
| |
Collapse
|
5
|
Lowman A, Bobholz S, Brehler M, Duenweg S, Kyereme F, Cochran E, Coss D, Connelly J, Mueller W, Agarwal M, Banerjee A, LaViolette P. NIMG-18. A GROUND TRUTH COMPARISON OF PATHOLOGICALLY CONFIRMED GLIOBLASTOMA MARGINS TO CONTRAST ENHANCEMENT AT AUTOPSY. Neuro Oncol 2022. [DOI: 10.1093/neuonc/noac209.636] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022] Open
Abstract
Abstract
PURPOSE
Glioblastoma is one of the most common and deadly adult brain tumors. Current standard treatment is surgical resection followed by radiation and concomitant chemotherapy (chemoRT). Glioblastoma progression is monitored using MRI, primarily relying on post-contrast T1-weighted imaging (T1C). Unfortunately, tumor invasion is known to extend beyond traditional contrast enhancement. T1-subtraction (T1S) maps have been introduced as a better tumor volume estimate. In this study we compare T1S map derived tumor annotations to a radiologist for identifying histologically confirmed tumor in patients with differing treatment histories at autopsy.
METHODS
Ten patients with autopsy confirmed glioblastoma and MRI within 1 month of death were recruited for this study. Seven patients received chemoRT combined (chemoRT+) and three patients received no treatment beyond surgery (chemoRT-). Patient’s brains were sliced axially in the same orientation as their final MRI using a patient-specific slicing jig. Large tissue samples were taken, processed, embedded in paraffin, stained for hematoxylin and eosin, and digitized at 40x resolution. Digital images were annotated for infiltrative tumor, pseudopalisading necrosis, and necrosis without palisading cells. T1S and radiologist annotations were created for each patient using their final MRI (mean 18 days prior to death). The annotated histology images were aligned and resampled into MRI space using custom software and the overlap of pathologically confirmed tumor and MRI derived annotations was compared.
RESULTS
T1S maps alone were significantly better at identifying areas of histologically confirmed tumor in chemoRT+ patients compared to chemoRT- patients (p=0.043). T1S derived annotations overlapped with 52% of histologically confirmed tumor in the chemoRT- patients and 78% in the chemoRT+. The radiologist drawn tumor masks were more accurate in chemoRT+ patients, identifying 61% confirmed tumor (trending, p=0.097, chemoRT+=61%).
CONCLUSION
These results demonstrate the difficulty of identifying pathologically confirmed tumor outside contrast enhancement in glioblastoma patients, even in the untreated state.
Collapse
Affiliation(s)
| | | | | | | | | | | | - Dylan Coss
- Medical College of Wisconsin , Milwaukee, WI , USA
| | | | - Wade Mueller
- Medical College of Wisconsins , Milwaukee, WI , USA
| | | | | | | |
Collapse
|
6
|
McGarry SD, Brehler M, Bukowy JD, Lowman AK, Bobholz SA, Duenweg SR, Banerjee A, Hurrell SL, Malyarenko D, Chenevert TL, Cao Y, Li Y, You D, Fedorov A, Bell LC, Quarles CC, Prah MA, Schmainda KM, Taouli B, LoCastro E, Mazaheri Y, Shukla‐Dave A, Yankeelov TE, Hormuth DA, Madhuranthakam AJ, Hulsey K, Li K, Huang W, Huang W, Muzi M, Jacobs MA, Solaiyappan M, Hectors S, Antic T, Paner GP, Palangmonthip W, Jacobsohn K, Hohenwalter M, Duvnjak P, Griffin M, See W, Nevalainen MT, Iczkowski KA, LaViolette PS. Multi-Site Concordance of Diffusion-Weighted Imaging Quantification for Assessing Prostate Cancer Aggressiveness. J Magn Reson Imaging 2022; 55:1745-1758. [PMID: 34767682 PMCID: PMC9095769 DOI: 10.1002/jmri.27983] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2021] [Revised: 10/21/2021] [Accepted: 10/22/2021] [Indexed: 01/01/2023] Open
Abstract
BACKGROUND Diffusion-weighted imaging (DWI) is commonly used to detect prostate cancer, and a major clinical challenge is differentiating aggressive from indolent disease. PURPOSE To compare 14 site-specific parametric fitting implementations applied to the same dataset of whole-mount pathologically validated DWI to test the hypothesis that cancer differentiation varies with different fitting algorithms. STUDY TYPE Prospective. POPULATION Thirty-three patients prospectively imaged prior to prostatectomy. FIELD STRENGTH/SEQUENCE 3 T, field-of-view optimized and constrained undistorted single-shot DWI sequence. ASSESSMENT Datasets, including a noise-free digital reference object (DRO), were distributed to the 14 teams, where locally implemented DWI parameter maps were calculated, including mono-exponential apparent diffusion coefficient (MEADC), kurtosis (K), diffusion kurtosis (DK), bi-exponential diffusion (BID), pseudo-diffusion (BID*), and perfusion fraction (F). The resulting parametric maps were centrally analyzed, where differentiation of benign from cancerous tissue was compared between DWI parameters and the fitting algorithms with a receiver operating characteristic area under the curve (ROC AUC). STATISTICAL TEST Levene's test, P < 0.05 corrected for multiple comparisons was considered statistically significant. RESULTS The DRO results indicated minimal discordance between sites. Comparison across sites indicated that K, DK, and MEADC had significantly higher prostate cancer detection capability (AUC range = 0.72-0.76, 0.76-0.81, and 0.76-0.80 respectively) as compared to bi-exponential parameters (BID, BID*, F) which had lower AUC and greater between site variation (AUC range = 0.53-0.80, 0.51-0.81, and 0.52-0.80 respectively). Post-processing parameters also affected the resulting AUC, moving from, for example, 0.75 to 0.87 for MEADC varying cluster size. DATA CONCLUSION We found that conventional diffusion models had consistent performance at differentiating prostate cancer from benign tissue. Our results also indicated that post-processing decisions on DWI data can affect sensitivity and specificity when applied to radiological-pathological studies in prostate cancer. LEVEL OF EVIDENCE 1 TECHNICAL EFFICACY: Stage 3.
Collapse
Affiliation(s)
- Sean D. McGarry
- Department of BiophysicsMedical College of WisconsinMilwaukeeWisconsinUSA
| | - Michael Brehler
- Department of RadiologyMedical College of WisconsinMilwaukeeWIUSA
| | - John D. Bukowy
- Department of Electrical Engineering and Computer ScienceMilwaukee School of EngineeringMilwaukeeWIUSA
| | | | - Samuel A. Bobholz
- Department of BiophysicsMedical College of WisconsinMilwaukeeWisconsinUSA
| | | | - Anjishnu Banerjee
- Division of BiostatisticsMedical College of WisconsinMilwaukeeWisconsinUSA
| | - Sarah L. Hurrell
- Department of RadiologyMedical College of WisconsinMilwaukeeWIUSA
| | | | | | - Yue Cao
- Department of RadiologyUniversity of MichiganAnn ArborMichiganUSA,Department of Radiation OncologyUniversity of MichiganAnn ArborMichiganUSA
| | - Yuan Li
- Department of Radiation OncologyUniversity of MichiganAnn ArborMichiganUSA
| | - Daekeun You
- Department of Radiation OncologyUniversity of MichiganAnn ArborMichiganUSA
| | - Andrey Fedorov
- Department of RadiologyBrigham and Women's HospitalBostonMassachusettsUSA
| | - Laura C. Bell
- Division of Neuroimaging ResearchBarrow Neurological InstitutePhoenixArizonaUSA
| | - C. Chad Quarles
- Division of Neuroimaging ResearchBarrow Neurological InstitutePhoenixArizonaUSA
| | - Melissa A. Prah
- Department of BiophysicsMedical College of WisconsinMilwaukeeWisconsinUSA
| | | | - Bachir Taouli
- Department of RadiologyIcahn School of Medicine at Mount SinaiNew YorkNew YorkUSA
| | - Eve LoCastro
- Department of Medical PhysicsMemorial Sloan Kettering Cancer CenterNew YorkNew YorkUSA
| | - Yousef Mazaheri
- Department of Medical PhysicsMemorial Sloan Kettering Cancer CenterNew YorkNew YorkUSA,Department of RadiologyMemorial Sloan Kettering Cancer CenterNew YorkNew YorkUSA
| | - Amita Shukla‐Dave
- Department of Medical PhysicsMemorial Sloan Kettering Cancer CenterNew YorkNew YorkUSA,Department of RadiologyMemorial Sloan Kettering Cancer CenterNew YorkNew YorkUSA
| | - Thomas E. Yankeelov
- Department of Biomedical Engineering, Diagnostic Medicine, Oncology, Oden Institute for Computational Engineering and Sciences, Livestrong Cancer InstitutesThe University of TexasAustinTexasUSA
| | - David A. Hormuth
- Department of Biomedical Engineering, Diagnostic Medicine, Oncology, Oden Institute for Computational Engineering and Sciences, Livestrong Cancer InstitutesThe University of TexasAustinTexasUSA
| | | | - Keith Hulsey
- Department of RadiologyThe University of Texas Southwestern Medical CenterDallasTexasUSA
| | - Kurt Li
- International School of BeavertonAlohaOregonUSA
| | - Wei Huang
- Advanced Imaging Research CenterOregon Health Sciences UniversityPortlandOregonUSA
| | - Wei Huang
- Department of PathologyOregon Health and Science UniversityMadisonWisconsinUSA
| | - Mark Muzi
- Department of Radiology, Neurology, and Radiation OncologyUniversity of WashingtonSeattleWashingtonUSA
| | - Michael A. Jacobs
- The Russell H. Morgan Department of Radiology and Radiological Science and Sidney Kimmel Comprehensive Cancer CenterJohns Hopkins University School of MedicineBaltimoreMarylandUSA
| | - Meiyappan Solaiyappan
- The Russell H. Morgan Department of Radiology and Radiological Science and Sidney Kimmel Comprehensive Cancer CenterJohns Hopkins University School of MedicineBaltimoreMarylandUSA
| | - Stefanie Hectors
- Department of biomedical engineering and imaging instituteWeill Cornell Medical CollegeNew York CityNew YorkUSA
| | - Tatjana Antic
- Department of PathologyUniversity of ChicagoChicagoIllinoisUSA
| | | | - Watchareepohn Palangmonthip
- Department of PathologyMedical College of WisconsinMilwaukeeWisconsinUSA,Department of PathologyChiang Mai UniversityChiang MaiThailand
| | - Kenneth Jacobsohn
- Department of UrologyMedical College of WisconsinMilwaukeeWisconsinUSA
| | - Mark Hohenwalter
- Department of RadiologyMedical College of WisconsinMilwaukeeWIUSA
| | - Petar Duvnjak
- Department of RadiologyMedical College of WisconsinMilwaukeeWIUSA
| | - Michael Griffin
- Department of RadiologyMedical College of WisconsinMilwaukeeWIUSA
| | - William See
- Department of UrologyMedical College of WisconsinMilwaukeeWisconsinUSA
| | | | | | - Peter S. LaViolette
- Department of RadiologyMedical College of WisconsinMilwaukeeWIUSA,Department of Biomedical EngineeringMedical College of WisconsinMilwaukeeWisconsinUSA
| |
Collapse
|
7
|
Bobholz SA, Lowman AK, Brehler M, Kyereme F, Duenweg SR, Sherman J, McGarry SD, Cochran EJ, Connelly J, Mueller WM, Agarwal M, Banerjee A, LaViolette PS. Radio-Pathomic Maps of Cell Density Identify Brain Tumor Invasion beyond Traditional MRI-Defined Margins. AJNR Am J Neuroradiol 2022; 43:682-688. [PMID: 35422419 PMCID: PMC9089258 DOI: 10.3174/ajnr.a7477] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Accepted: 02/07/2022] [Indexed: 11/07/2022]
Abstract
BACKGROUND AND PURPOSE Currently, contrast-enhancing margins on T1WI are used to guide treatment of gliomas, yet tumor invasion beyond the contrast-enhancing region is a known confounding factor. Therefore, this study used postmortem tissue samples aligned with clinically acquired MRIs to quantify the relationship between intensity values and cellularity as well as to develop a radio-pathomic model to predict cellularity using MR imaging data. MATERIALS AND METHODS This single-institution study used 93 samples collected at postmortem examination from 44 patients with brain cancer. Tissue samples were processed, stained with H&E, and digitized for nuclei segmentation and cell density calculation. Pre- and postgadolinium contrast T1WI, T2 FLAIR, and ADC images were collected from each patient's final acquisition before death. In-house software was used to align tissue samples to the FLAIR image via manually defined control points. Mixed-effects models were used to assess the relationship between single-image intensity and cellularity for each image. An ensemble learner was trained to predict cellularity using 5 × 5 voxel tiles from each image, with a two-thirds to one-third train-test split for validation. RESULTS Single-image analyses found subtle associations between image intensity and cellularity, with a less pronounced relationship in patients with glioblastoma. The radio-pathomic model accurately predicted cellularity in the test set (root mean squared error = 1015 cells/mm2) and identified regions of hypercellularity beyond the contrast-enhancing region. CONCLUSIONS A radio-pathomic model for cellularity trained with tissue samples acquired at postmortem examination is able to identify regions of hypercellular tumor beyond traditional imaging signatures.
Collapse
Affiliation(s)
- S A Bobholz
- From the Departments of Biophysics (S.A.B., S.R.D., J.S., S.D.M.)
| | | | - M Brehler
- Radiology (A.L., M.B., M.A., P.S.L.)
| | | | - S R Duenweg
- From the Departments of Biophysics (S.A.B., S.R.D., J.S., S.D.M.)
| | - J Sherman
- From the Departments of Biophysics (S.A.B., S.R.D., J.S., S.D.M.)
| | - S D McGarry
- From the Departments of Biophysics (S.A.B., S.R.D., J.S., S.D.M.)
| | | | | | | | - M Agarwal
- Radiology (A.L., M.B., M.A., P.S.L.)
| | | | - P S LaViolette
- Radiology (A.L., M.B., M.A., P.S.L.)
- Biomedical Engineering (P.S.L.), Medical College of Wisconsin, Milwaukee, Wisconsin
| |
Collapse
|
8
|
Bobholz S, Lowman A, Brehler M, Sherman J, Duenweg S, Kyereme F, Connelly J, Cochran E, Mueller W, Banerjee A, LaViolette P. INNV-18. EFFECT OF DURATION OF TUMOR TREATING FIELDS THERAPY ON CELLULARITY DISTRIBUTIONS BEYOND T1W MRI CONTRAST ENHANCING MARGIN AT AUTOPSY IN GLIOMA PATIENTS: PRELIMINARY RESULTS. Neuro Oncol 2021. [DOI: 10.1093/neuonc/noab196.429] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Abstract
Tumor treating fields (TTFields) are thought to disrupt the cell division process in glioma. This study uses autopsy tissue samples from patients recruited for brain donation to determine whether TTFields treatment duration and usage affects cellularity distributions beyond the gadolinium contrast enhancing region on T1-weighted (T1w) magnetic resonance imaging (MRI). At autopsy, 43 tissue samples from 18 patients who had undergone TTFields treatment (inclusion criteria: >50% compliance and 25-day duration) were collected from brain slices sectioned in-line with slices from the patients’ final MRI prior to death. Nuclei were segmented on the digitized hematoxylin and eosin (HE) stained tissue samples, which were then used to compute cell count across the slides. Tissue was then aligned to the patient’s final MRI scan prior to death using a custom in-house software, where manually defined control points were used to compute a nonlinear transform to warp the tissue to match the MRI. Histogram features including mean, median, 90th percentile, variance, and skewness, were computed from the cellularity distributions within regions outside the traditional tumor margin defined by T1w contrast enhancement for each subject. General linear models were fit to assess the effects of TTFields usage (percent of day) and duration (in days) on each histogram feature, controlling for age, overall survival, and time between TTFields treatment and death, as well as time between MRI acquisition and death. Longer TTFields treatment duration and overall survival was associated with significantly decreased skewness in the non-enhancing cellularity distributions (p< 0.05) while associations with other metrics remained statistically insignificant. These preliminary results in a small patient cohort suggest that TTFields duration may reduce the presence of high cellularity tails in the non-contrast enhancing distributions. Additional research in larger patient populations is warranted to better understand these findings given the number of confounding factors.
Collapse
Affiliation(s)
| | | | | | - John Sherman
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | | | | | | - Wade Mueller
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | |
Collapse
|
9
|
Lowman A, Bobholz S, Brehler M, Duenweg S, Sherman J, Kyereme F, Connelly J, Cochran E, Mueller W, Banerjee A, LaViolette P. NIMG-30. PATHOLOGICAL VALIDATION OF CONTRAST ENHANCEMENT AT AUTOPSY BETWEEN PATIENTS TREATED AND UNTREATED WITH CHEMOTHERAPY AND RADIATION. Neuro Oncol 2021. [DOI: 10.1093/neuonc/noab196.530] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Abstract
Tumor heterogeneity in glioblastoma complicates delineation of active tumor using standard MR imaging. Pseudo-progression following treatment with chemotherapy and radiation (chemoRT) further complicates how tumors appear. T1-weighted subtraction maps (T1S) have been used to better identify subtly enhancing regions containing infiltrative tumor. This study examines the differences in tumor appearance between patients treated with chemoRT compared to a cohort opting out at autopsy, to understand how chemoRT changes contrast enhancement on MRI. Ten patients diagnosed with glioblastoma were recruited for whole brain donation for this study. Three patients received no treatment (chemoRT-), and seven received a combination of chemoRT and additional treatments (chemoRT+), including but not limited to bevacizumab (Bev) and tumor treating fields (TTF). Large tissue samples were taken at autopsy from whole brain samples sliced axially to align with the last clinical MRI using patient-specific 3D-printed slicing jigs. All tissue samples were hematoxylin and eosin (HE) stained and digitized at 40X resolution (27 total samples). The whole slide images (WSI) were annotated to outline regions containing necrosis without pseudopalisading cells, tumor with pseudopalisading necrosis, and infiltrative tumor. T1S were created for each patient by subtracting intensity normalized T1-weighted images from T1 post-contrast images (T1C). The annotated WSIs were aligned and resampled into MRI space using a custom software. A mixed effect model was used to compare the T1S intensity between chemoRT +/- cohorts within each pathological annotation, incorporating a random effect for subject. Mean T1S intensity was greater in untreated subjects when compared to chemoRT-subjects within each pathological annotation, including necrosis without pseudopalisading cells, tumor with pseudopalisading necrosis, and infiltrative tumor (p< 0.001). We show that chemoRT reduces the contrast enhancement in all aspects of pathologically validated tumor compartments, including infiltrative tumor. Future research is needed to examine if other patterns are evident in additional MR sequences.
Collapse
Affiliation(s)
| | | | | | | | - John Sherman
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | | | | - Wade Mueller
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | |
Collapse
|
10
|
Bobholz S, Lowman A, Brehler M, Sherman J, Duenweg S, Kyereme F, Connelly J, Cochran E, Mueller W, Banerjee A, LaViolette P. NIMG-42. MP-MRI-BASED TUMOR PROBABILITY MAPS TRAINED USING AUTOPSY TISSUE SAMPLES AS GROUND TRUTH NON-INVASIVELY PREDICT INFILTRATIVE TUMOR BEYOND THE CONTRAST ENHANCING REGION. Neuro Oncol 2021. [DOI: 10.1093/neuonc/noab196.541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Abstract
Infiltrative glioma beyond contrast enhancement on MRI is often difficult to identify with conventional imaging. In this study, we use large-format autopsy samples aligned to multi-parametric MRI to test the hypothesis that radio-pathomic machine learning models are able to accurately identify areas of infiltrative tumor beyond the contrast enhancing region. At autopsy, 140 tissue samples from 62 brain cancer patients were collected from brain slices sectioned to align with the patients’ last clinical MRI prior to death. Cell, extra-cellular fluid (ECF), and cytoplasm densities were computed from digitized, hematoxylin and eosin-stained samples, and a subset of 20 slides from 9 patients were annotated for tumor presence by a pathologist-trained technician. In-house custom software was used to align the tissue samples to the patients’ last clinical imaging, which included pre- and post-contrast T1, FLAIR, and ADC images. Bagging random forest models were then trained to predict cellularity, ECF, and cytoplasm density using 5-by-5 voxel tiles from each MRI as input. A 2/3-1/3 train-test split was used to validate model generalizability. A naïve Bayes classifier was trained to predict tumor class using cellularity, ECF, and cytoplasm segmentations within the annotation data set, again using a 2/3-1/3 train-test split to validate performance. The random forest models each accurately predicted cellularity, ECF, and cytoplasm density within the test data set, with root-mean-squared error values for each falling within one standard deviation of the ground truth. The histology-based tumor prediction model accurately predicted tumor, with a test set ROC AUC of 0.86. When using whole brain cellularity, ECF, and cytoplasm predictions from the random forest models as inputs for the naïve Bayes classifier, tumor probability maps identified regions of infiltrative tumor beyond contrast enhancement. Our results suggest that radio-pathomic maps of tumor probability accurately identify regions of infiltrative tumor beyond currently accepted MRI signatures.
Collapse
Affiliation(s)
| | | | | | - John Sherman
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | | | | | | - Wade Mueller
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | |
Collapse
|
11
|
Bobholz SA, Lowman AK, Barrington A, Brehler M, McGarry S, Cochran EJ, Connelly J, Mueller WM, Agarwal M, O'Neill D, Nencka AS, Banerjee A, LaViolette PS. Radiomic Features of Multiparametric MRI Present Stable Associations With Analogous Histological Features in Patients With Brain Cancer. ACTA ACUST UNITED AC 2021; 6:160-169. [PMID: 32548292 PMCID: PMC7289245 DOI: 10.18383/j.tom.2019.00029] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Magnetic resonance (MR)-derived radiomic features have shown substantial predictive utility in modeling different prognostic factors of glioblastoma and other brain cancers. However, the biological relationship underpinning these predictive models has been largely unstudied, and the generalizability of these models had been called into question. Here, we examine the localized relationship between MR-derived radiomic features and histology-derived “histomic” features using a data set of 16 patients with brain cancer. Tile-based radiomic features were collected on T1, post-contrast T1, FLAIR, and diffusion-weighted imaging (DWI)-derived apparent diffusion coefficient (ADC) images acquired before patient death, with analogous histomic features collected for autopsy samples coregistered to the magnetic resonance imaging. Features were collected for each original image, as well as a 3D wavelet decomposition of each image, resulting in 837 features per MR and histology image. Correlative analyses were used to assess the degree of association between radiomic–histomic pairs for each magnetic resonance imaging. The influence of several confounds was also assessed using linear mixed-effect models for the normalized radiomic–histomic distance, testing for main effects of different acquisition field strengths. Results as a whole were largely heterogeneous, but several features showed substantial associations with their histomic analogs, particularly those derived from the FLAIR and postcontrast T1W images. These features with the strongest association typically presented as stable across field strengths as well. These data suggest that a subset of radiomic features can consistently capture texture information on underlying tissue histology.
Collapse
|
12
|
Lowman A, Bobholz S, Connelly JM, Cochran E, Mueller W, Brehler M, Kyereme F, Sherman J, Duenweg S, LaViolette PS. Measurement of treatment dependent glioblastoma cell density in T1- weighted contrast enhancement at autopsy. J Clin Oncol 2021. [DOI: 10.1200/jco.2021.39.15_suppl.e14035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
e14035 Background: With an average overall survival of 12-18 months, glioblastoma has a particularly grim diagnosis. Standard treatment of glioblastoma, following detection on MRI, is surgical resection followed by radiation therapy and chemotherapy and is monitored through MR imaging. Glioblastoma has a unique heterogenous nature that complicates visualization of subtly enhancing tumor. This study used autopsy tissue samples taken from glioblastoma patients with varying treatment, to examine the effects of treatment on cell density within regions of contrast enhancement, using T1-weighted subtraction maps (T1S) from the last MR images to death. Methods: Eight patients diagnosed with glioblastoma at autopsy were recruited for this study. Two patients had no treatment and six received a combination of chemo-radiation and other treatments, including but not limited to bevacizumab (Bev) and TTFields therapy. At autopsy, whole brain samples were sliced axially aligned to the patient’s final MRI to death. Time between last MRI and death ranged from 4-27 days (mean 16 days). Overall survival (OS) ranged from 4-538 days (mean 307 days). Large tissue samples were taken from regions of suspected tumor or treatment effect, for a total of 18 tissue samples. Tissue samples were processed, H&E stained, and digitized at 40X resolution (Huron Slide Scanner). Cell density (cells/mm2) was calculated using digital histology. T1S were created for each patient by subtracting intensity normalized T1 weighted images from T1 post contrast images (T1C). Digital histology was aligned and resampled into MRI space using manual control point registration. Mixed effect models were used to compare differences in cell density across contrast enhancement (T1SE vs. NE) as well as across treatment groups (treatment vs. no treatment). Results: Cellularity was compared across regions of T1S enhancement (T1SE) and non-enhancement (NE) within manually selected regions of interest. Cell density compared between regions of T1SE and NE was not different (p=0.219). Total cell density was increased in patients who had received treatment compared to no treatment in both regions of T1SE and NE (p=0.014). Conclusions: Overall, cell density was increased in patients who had received treatment after diagnosis of glioblastoma. Additional research is needed to examine the extent of treatment’s effect on cellularity of glioblastoma. This work begins to characterize the use of T1S in evaluating glioblastoma tumor burden in patients with varying treatment histories.[Table: see text]
Collapse
|
13
|
Ibrahim ESH, Baruah D, Croisille P, Stojanovska J, Rubenstein JC, Frei A, Schlaak RA, Lin CY, Pipke JL, Lemke A, Xu Z, Klaas A, Brehler M, Flister MJ, Laviolette PS, Gore EM, Bergom C. Cardiac Magnetic Resonance for Early Detection of Radiation Therapy-Induced Cardiotoxicity in a Small Animal Model. JACC CardioOncol 2021; 3:113-130. [PMID: 33912843 PMCID: PMC8078846 DOI: 10.1016/j.jaccao.2020.12.006] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Background Over half of all cancer patients receive radiation therapy (RT). However, radiation exposure to the heart can cause cardiotoxicity. Nevertheless, there is a paucity of data on RT-induced cardiac damage, with limited understanding of safe regional RT doses, early detection, prevention and management. A common initial feature of cardiotoxicity is asymptomatic dysfunction, which if left untreated may progress to heart failure. The current paradigm for cardiotoxicity detection and management relies primarily upon assessment of ejection fraction (EF). However, cardiac injury can occur without a clear change in EF. Objectives To identify magnetic resonance imaging (MRI) markers of early RT-induced cardiac dysfunction. Methods We investigated the effect of RT on global and regional cardiac function and myocardial T1/T2 values at two timepoints post-RT using cardiac MRI in a rat model of localized cardiac RT. Rats who received image-guided whole-heart radiation of 24Gy were compared to sham-treated rats. Results The rats maintained normal global cardiac function post-RT. However, a deterioration in strain was particularly notable at 10-weeks post RT, and changes in circumferential strain were larger than changes in radial or longitudinal strain. Compared to sham, circumferential strain changes occurred at the basal, mid-ventricular and apical levels (p<0.05 for all at both 8-weeks and 10-weeks post-RT), most of the radial strain changes occurred at the mid-ventricular (p=0.044 at 8-weeks post-RT) and basal (p=0.018 at 10-weeks post-RT) levels, and most of the longitudinal strain changes occurred at the apical (p=0.002 at 8-weeks post-RT) and basal (p=0.035 at 10-weeks post-RT) levels. Regionally, lateral myocardial segments showed the greatest worsening in strain measurements, and histologic changes supported these findings. Despite worsened myocardial strain post-RT, myocardial tissue displacement measures were maintained, or even increased. T1/T2 measurements showed small non-significant changes post-RT compared to values in non-irradiated rats. Conclusions Our findings suggest MRI regional myocardial strain is a sensitive imaging biomarker for detecting RT-induced subclinical cardiac dysfunction prior to compromise of global cardiac function.
Collapse
Affiliation(s)
- El-Sayed H Ibrahim
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Cardiovascular Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Cancer Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Biomedical Engineering, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Dhiraj Baruah
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Pierre Croisille
- Jean-Monnet University, 10 Rue Trefilerie, 42100 Saint-Etienne, France
| | | | - Jason C Rubenstein
- Department of Medicine, Division of Cardiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Anne Frei
- Department of Radiation Oncology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Rachel A Schlaak
- Department of Pharmacology and Toxicology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Chieh-Yu Lin
- Department of Pathology & Immunology, Washington University School of Medicine, St. Louis, Missouri, USA
| | - Jamie L Pipke
- Department of Radiation Oncology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Angela Lemke
- Department of Physiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Zhiqiang Xu
- Department of Radiation Oncology, Washington University School of Medicine, St. Louis, Missouri, USA
| | - Amanda Klaas
- Department of Radiation Oncology, Washington University School of Medicine, St. Louis, Missouri, USA
| | - Michael Brehler
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Michael J Flister
- Cardiovascular Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Cancer Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Physiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Genomic Sciences and Precision Medicine Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Peter S Laviolette
- Department of Radiology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Biomedical Engineering, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Elizabeth M Gore
- Department of Radiation Oncology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA
| | - Carmen Bergom
- Cardiovascular Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Cancer Center, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Radiation Oncology, Medical College of Wisconsin, Milwaukee, Wisconsin, USA.,Department of Radiation Oncology, Washington University School of Medicine, St. Louis, Missouri, USA
| |
Collapse
|
14
|
Bobholz S, Lowman A, Connelly J, Cochran E, Mueller W, McGarry S, Brehler M, Gliszinski C, Wilczynski A, LaViolette P. NIMG-23. COMPARISON OF A RADIO-PATHOMIC MODEL VERSUS A RADIOLOGY-ONLY TUMOR SEGMENTATION MODEL FOR THE DETECTION OF INFILTRATIVE TUMOR IN GLIOMA PATIENTS. Neuro Oncol 2020. [DOI: 10.1093/neuonc/noaa215.636] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Abstract
This study used large format autopsy tissue samples to compare radio-pathomic maps of brain cancer to a current tumor segmentation algorithm. We hypothesized that an MRI-based machine learning model trained with actual histology rather than radiologist annotations cellularity would 1) improve delineation between tumor and treatment effect, and 2) detect abnormal pathology beyond the contrast-enhancing tumor region. Seventeen patients with pathologically confirmed glioma were included in this study. At autopsy, 43 tissue samples were collected from 17 subjects from whole brain slices sectioned to align with the last axial MRI prior to death. Cellularity was calculated using a color deconvolution on 40X digitized H&E stained slides from the tissue samples. In-house custom software was used to align tissue samples and cellularity information to the FLAIR image using manually defined control points. The DeepMedic algorithm was trained to segment tumors using the BraTs 2017 dataset, and then applied to our patients in order to create automated tumor probability maps. An MRI-based ensemble algorithm using a 5x5 voxel searchlight (input: T1, T1C, FLAIR, ADC) was used to predict cellularity at each voxel, using tissue samples from 14 subjects as ground truth. Both models were applied to 3 withheld test subjects in order to compare tumor probability and cellularity predictions to the pathological ground truth. The mutual information between tumor probability and actual cellularity was 1X10-15, relatively low compared to the rad-path predicted cellularity (=0.16), despite the tumor prediction model accurately highlighting regions of contrast enhancement. Additionally, the radio-pathomic ensemble model correctly identify regions of hypercellularity beyond the tumor segmentation model as well as regions of within the segmented tumor area. This study demonstrates the utility of training machine learning models with pathological ground truth rather than radiologist annotations for predicting localized tumor information, particularly in the post-treatment stage.
Collapse
Affiliation(s)
| | | | | | | | - Wade Mueller
- Medical College of Wisconsin, Milwaukee, WI, USA
| | - Sean McGarry
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | | | | | |
Collapse
|
15
|
Bobholz S, Lowman A, Connelly J, Cochran E, Mueller W, McGarry S, Brehler M, LaViolette P. INNV-10. EFFECTS OF TTFIELDS USAGE AND DURATION OF USAGE ON CELLULARITY AND Ki67 DISTRIBUTIONS AT AUTOPSY. Neuro Oncol 2020. [DOI: 10.1093/neuonc/noaa215.493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Abstract
This study used autopsy tissue samples taken from patients with glioblastoma who had undergone Tumor-Treating Fields (TTFields) therapy to examine the effects of treatment usage and duration on cellular and mitotic activity distributions. We hypothesized that treatment duration and percent compliance would be associated with a change in cellularity. Three tissue samples were collected from 10 patients with glioblastoma at autopsy. These samples were specifically collected in regions that demonstrated clear contrast enhancement on the patients’ last clinical imaging, in order to capture regions of suspected active tumor. Samples were stained with both hematoxylin and eosin (HE) and Ki67 immuno-histochemistry (IHC) and digitized at 40X resolution using a digital microscope. A color deconvolution algorithm was used to separate stains, which allowed nuclei segmentation. Number of cells and percent of Ki67 positive cells were computed across 100x100 superpixels. Spearman’s correlations were used to examine the association between TTFields usage (measured as a percent) and duration of use (in days) and signatures of active tumor proliferation, measured as median cellularity and proportion of Ki67 positive cells across each patients’ three samples. Trending negative associations with median cellularity were observed for both TTFields usage (R=-0.60, p=0.065) and duration (R=-0.61,p=0.063). The magnitude of the treatment duration effect increased when controlling for time between diagnosis and death, despite a p-value drop due to the lost statistical degree of freedom (R=-0.62, p=0.073). No significant effects were observed for Ki67 positive staining. These results generally suggest a possible effect of TTFields therapy duration, where longer treatments and increased compliance lead to lower tumor proliferation. However, these results should be interpreted with caution, as larger, more statistically powerful studies will be better suited to assess the generalizability of the preliminary trends seen in these RESULTS:
Collapse
Affiliation(s)
| | | | | | | | - Wade Mueller
- Medical College of Wisconsin, Milwaukee, WI, USA
| | - Sean McGarry
- Medical College of Wisconsin, Milwaukee, WI, USA
| | | | | |
Collapse
|
16
|
McGarry SD, Bukowy JD, Iczkowski KA, Lowman AK, Brehler M, Bobholz S, Nencka A, Barrington A, Jacobsohn K, Unteriner J, Duvnjak P, Griffin M, Hohenwalter M, Keuter T, Huang W, Antic T, Paner G, Palangmonthip W, Banerjee A, LaViolette PS. Radio-pathomic mapping model generated using annotations from five pathologists reliably distinguishes high-grade prostate cancer. J Med Imaging (Bellingham) 2020; 7:054501. [PMID: 32923510 PMCID: PMC7479263 DOI: 10.1117/1.jmi.7.5.054501] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2020] [Accepted: 08/20/2020] [Indexed: 12/19/2022] Open
Abstract
Purpose: Our study predictively maps epithelium density in magnetic resonance imaging (MRI) space while varying the ground truth labels provided by five pathologists to quantify the downstream effects of interobserver variability. Approach: Clinical imaging and postsurgical tissue from 48 recruited prospective patients were used in our study. Tissue was sliced to match the MRI orientation and whole-mount slides were stained and digitized. Data from 28 patients ( n = 33 slides) were sent to five pathologists to be annotated. Slides from the remaining 20 patients ( n = 123 slides) were annotated by one of the five pathologists. Interpathologist variability was measured using Krippendorff's alpha. Pathologist-specific radiopathomic mapping models were trained using a partial least-squares regression using MRI values to predict epithelium density, a known marker for disease severity. An analysis of variance characterized intermodel means difference in epithelium density. A consensus model was created and evaluated using a receiver operator characteristic classifying high grade versus low grade and benign, and was statistically compared to apparent diffusion coefficient (ADC). Results: Interobserver variability ranged from low to acceptable agreement (0.31 to 0.69). There was a statistically significant difference in mean predicted epithelium density values ( p < 0.001 ) between the five models. The consensus model outperformed ADC (areas under the curve = 0.80 and 0.71, respectively, p < 0.05 ). Conclusion: We demonstrate that radiopathomic maps of epithelium density are sensitive to the pathologist annotating the dataset; however, it is unclear if these differences are clinically significant. The consensus model produced the best maps, matched the performance of the best individual model, and outperformed ADC.
Collapse
Affiliation(s)
- Sean D McGarry
- Medical College of Wisconsin, Department of Biophysics, Milwaukee, Wisconsin, United States
| | - John D Bukowy
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Kenneth A Iczkowski
- Medical College of Wisconsin, Department of Pathology, Milwaukee, Wisconsin, United States
| | - Allison K Lowman
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Michael Brehler
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Samuel Bobholz
- Medical College of Wisconsin, Department of Biophysics, Milwaukee, Wisconsin, United States
| | - Andrew Nencka
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Alex Barrington
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Kenneth Jacobsohn
- Medical College of Wisconsin, Department of Urological Surgery, Milwaukee, Wisconsin, United States
| | - Jackson Unteriner
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Petar Duvnjak
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Michael Griffin
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Mark Hohenwalter
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States
| | - Tucker Keuter
- Medical College of Wisconsin, Department of Biostatistics, Milwaukee, Wisconsin, United States
| | - Wei Huang
- University of Wisconsin-Madison, Department of Pathology, Madison, Wisconsin, United States
| | - Tatjana Antic
- University of Chicago, Department of Pathology, Chicago, Illinois, United States
| | - Gladell Paner
- University of Chicago, Department of Pathology, Chicago, Illinois, United States
| | - Watchareepohn Palangmonthip
- Medical College of Wisconsin, Department of Pathology, Milwaukee, Wisconsin, United States.,Chiang Mai University, Department of Pathology, Faculty of Medicine, Chiang Mai, Thailand
| | - Anjishnu Banerjee
- Medical College of Wisconsin, Department of Biostatistics, Milwaukee, Wisconsin, United States
| | - Peter S LaViolette
- Medical College of Wisconsin, Department of Radiology, Milwaukee, Wisconsin, United States.,Medical College of Wisconsin, Department of Biomedical Engineering, Milwaukee, Wisconsin, United States
| |
Collapse
|
17
|
Kaplan JT, Ramsay JW, Cameron SE, Seymore KD, Brehler M, Thawait GK, Zbijewski WB, Siewerdsen JH, Brown TN. Association Between Knee Anatomic Metrics and Biomechanics for Male Soldiers Landing With Load. Am J Sports Med 2020; 48:1389-1397. [PMID: 32255657 DOI: 10.1177/0363546520911608] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
BACKGROUND Anterior cruciate ligament (ACL) injury is a military occupational hazard that may be attributed to an individual's knee biomechanics and joint anatomy. This study sought to determine if greater flexion when landing with load resulted in knee biomechanics thought to decrease ACL injury risk and whether knee biomechanics during landing relate to knee anatomic metrics. HYPOTHESIS Anatomic metrics regarding the slope and concavity of the tibial plateau will exhibit a significant relation to the increased anterior shear force on the knee and decreased knee flexion posture during landing with body-borne load. STUDY DESIGN Descriptive laboratory study. METHODS Twenty male military personnel completed a drop landing task with 3 load conditions: light (~6 kg), medium (15% body weight), and heavy (30% body weight). Participants were divided into groups based on knee flexion exhibited when landing with the heavy load (high- and low-Δflexion). Tibial slopes and depth were measured on weightbearing volumetric images of the knee obtained with a prototype cone beam computed tomography system. Knee biomechanics were submitted to a linear mixed model to evaluate the effect of landing group and load, with the anatomic metrics considered covariates. RESULTS Load increased peak proximal anterior tibial shear force (P = .034), knee flexion angle (P = .024), and moment (P = .001) during landing. Only the high flexion group increased knee flexion (P < .001) during weighted landings with medium and heavy loads. The low flexion group used greater knee abduction angle (P = .030) and peak proximal anterior tibial shear force (P = .034) when landing with load. Anatomic metrics did not differ between groups, but ratio of medial-to-lateral tibial slope and medial tibial depth predicted peak proximal anterior tibial shear force (P = .009) and knee flexion (P = .034) during landing, respectively. CONCLUSION Increasing knee flexion is an attainable strategy to mitigate risk of ACL injury, but certain individuals may be predisposed to knee forces and biomechanics that load the ACL during weighted landings. CLINICAL RELEVANCE The ability to screen individuals for anatomic metrics that predict knee flexion may identify soldiers and athletes who require additional training to mitigate the risk of lower extremity injury.
Collapse
Affiliation(s)
- Jonathan T Kaplan
- Combat Capabilities Development Command Soldier Center, Natick, Massachusetts, USA
| | - John W Ramsay
- Combat Capabilities Development Command Soldier Center, Natick, Massachusetts, USA
| | | | - Kayla D Seymore
- Department of Kinesiology, Boise State University, Boise, Idaho, USA
| | - Michael Brehler
- Russel H. Morgan Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA
| | - Gaurav K Thawait
- Russel H. Morgan Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA
| | - Wojciech B Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
| | - Jeffrey H Siewerdsen
- Russel H. Morgan Department of Radiology, Johns Hopkins University, Baltimore, Maryland, USA.,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, Maryland, USA
| | - Tyler N Brown
- Department of Kinesiology, Boise State University, Boise, Idaho, USA
| |
Collapse
|
18
|
Subramanian S, Brehler M, Cao Q, Quevedo Gonzalez FJ, Breighner RE, Carrino JA, Wright T, Yorkston J, Siewerdsen JH, Zbijewski W. Quantitative Evaluation of Bone Microstructure using High-Resolution Extremity Cone-Beam CT with a CMOS Detector. Proc SPIE Int Soc Opt Eng 2019; 10953. [PMID: 31814656 DOI: 10.1117/12.2515504] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
Purpose A high-resolution cone-beam CT (CBCT) system for extremity imaging has been developed using a custom complementary metal-oxide-semiconductor (CMOS) x-ray detector. The system has spatial resolution capability beyond that of recently introduced clinical orthopedic CBCT. We evaluate performance of this new scanner in quantifying trabecular microstructure in subchondral bone of the knee. Methods The high-resolution scanner uses the same mechanical platform as the commercially available Carestream OnSight 3D extremity CBCT, but replaces the conventional amorphous silicon flat-panel detector (a-Si:H FPD with 0.137 mm pixels and a ~0.7 mm thick scintillator) with a Dalsa Xineos3030 CMOS detector (0.1 mm pixels and a custom 0.4 mm scintillator). The CMOS system demonstrates ~40% improved spatial resolution (FWHM of a ~0.1 mm tungsten wire) and ~4× faster scan time than FPD-based extremity CBCT (FPD-CBCT). To investigate potential benefits of this enhanced spatial resolution in quantitative assessment of bone microstructure, 26 trabecular core samples were obtained from four cadaveric tibias and imaged using FPD-CBCT (75 μm voxels), CMOS-CBCT (75 μm voxels), and reference micro-CT (μCT, 15 μm voxels). CBCT bone segmentations were obtained using local Bernsen's thresholding combined with global histogram-based pre-thresholding; μCT segmentation involved Otsu's method. Measurements of trabecular thickness (Tb.Th), spacing (Tb.Sp), number (Tb.N) and bone volume (BV/TV) were performed in registered regions of interest in the segmented CBCT and μCT reconstructions. Results CMOS-CBCT achieved noticeably improved delineation of trabecular detail compared to FPD-CBCT. Correlations with reference μCT for metrics of bone microstructure were better for CMOS-CBCT than FPD-CBCT, in particular for Tb.Th (increase in Pearson correlation from 0.84 with FPD-CBCT to 0.96 with CMOS-CBCT) and Tb.Sp (increase from 0.80 to 0.85). This improved quantitative performance of CMOS-CBCT is accompanied by a reduction in scan time, from ~60 sec for a clinical high resolution protocol on FPD-CBCT to ~17 sec for CMOS-CBCT. Conclusion The CMOS-based extremity CBCT prototype achieves improved performance in quantification of bone microstructure, while retaining other diagnostic capabilities of its FPD-based precursor, including weight-bearing imaging. The new system offers a promising platform for quantitative imaging of skeletal health in osteoporosis and osteoarthritis.
Collapse
Affiliation(s)
- S Subramanian
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - M Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - Q Cao
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | | | - R E Breighner
- Department of Radiology and Imaging, Hospital for Special Surgery, New York, NY USA
| | - J A Carrino
- Department of Radiology and Imaging, Hospital for Special Surgery, New York, NY USA
| | - T Wright
- Biomechanics Laboratory, Hospital for Special Surgery, New York, NY USA
| | | | - J H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA.,Russell H. Morgan Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - W Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| |
Collapse
|
19
|
Brehler M, Islam A, Vogelsang L, Yang D, Sehnert W, Shakoor D, Demehri S, Siewerdsen JH, Zbijewski W. Coupled Active Shape Models for Automated Segmentation and Landmark Localization in High-Resolution CT of the Foot and Ankle. Proc SPIE Int Soc Opt Eng 2019; 10953. [PMID: 31337927 DOI: 10.1117/12.2515022] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Purpose We develop an Active Shape Model (ASM) framework for automated bone segmentation and anatomical landmark localization in weight-bearing Cone-Beam CT (CBCT). To achieve a robust shape model fit in narrow joint spaces of the foot (0.5 - 1 mm), a new approach for incorporating proximity constraints in ASM (coupled ASM, cASM) is proposed. Methods In cASM, shape models of multiple adjacent foot bones are jointly fit to the CBCT volume. This coupling enables checking for proximity between the evolving shapes to avoid situations where a conventional single-bone ASM might erroneously fit to articular surfaces of neighbouring bones. We used 21 extremity CBCT scans of the weight-bearing foot to compare segmentation and landmark localization accuracy of ASM and cASM in leave-one-out validation. Each scan was used as a test image once; shape models of calcaneus, talus, navicular, and cuboid were built from manual surface segmentations of the remaining 20 scans. The models were augmented with seven anatomical landmarks used for common measurements of foot alignment. The landmarks were identified in the original CBCT volumes and mapped onto mean bone shape surfaces. ASM and cASM were run for 100 iterations, and the number of principal shape components was increased every 10 iterations. Automated landmark localization was achieved by applying known point correspondences between landmark vertices on the mean shape and vertices of the final active shape segmentation of the test image. Results Root Mean Squared (RMS) error of bone surface segmentation improved from 3.6 mm with conventional ASM to 2.7 mm with cASM. Furthermore, cASM achieved convergence (no change in RMS error with iteration) after ~40 iterations of shape fitting, compared to ~60 iterations for ASM. Distance error in landmark localization was 25% to 55% lower (depending on the landmark) with cASM than with ASM. The importance of using a coupled model is underscored by the finding that cASM detected and corrected collisions between evolving shapes in 50% to 80% (depending on the bone) of shape model fits. Conclusion The proposed cASM framework improves accuracy of shape model fits, especially in complexes of tightly interlocking, articulated joints. The approach enables automated anatomical analysis in volumetric imaging of the foot and ankle, where narrow joint spaces challenge conventional shape models.
Collapse
Affiliation(s)
- M Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - A Islam
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | | | - D Yang
- Carestream Health, Rochester, NY USA
| | - W Sehnert
- Carestream Health, Rochester, NY USA
| | - D Shakoor
- Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - S Demehri
- Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - J H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA.,Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - W Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| |
Collapse
|
20
|
Shakoor D, Osgood GM, Brehler M, Zbijewski WB, de Cesar Netto C, Shafiq B, Orapin J, Thawait GK, Shon LC, Demehri S. Cone-beam CT measurements of distal tibio-fibular syndesmosis in asymptomatic uninjured ankles: does weight-bearing matter? Skeletal Radiol 2019; 48:583-594. [PMID: 30242446 DOI: 10.1007/s00256-018-3074-6] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Revised: 08/31/2018] [Accepted: 09/09/2018] [Indexed: 02/02/2023]
Abstract
OBJECTIVE To evaluate the influence of weight-bearing (WB) load in standard axial ankle syndesmotic measurements using cone beam CT (CBCT) examination of asymptomatic uninjured ankles. MATERIALS AND METHODS In this IRB approved, prospective study, patients with previous unilateral ankle fractures were recruited. We simultaneously scanned the injured ankles and asymptomatic contralateral ankles of 27 patients in both WB and NWB modes. For this study, only asymptomatic contralateral ankles with normal plain radiographs were included. Twelve standardized syndesmosis measurements at two axial planes (10 mm above the tibial plafond and 5 mm below the talar dome) were obtained by two expert readers using a custom CBCT viewer with the capability for geometric measurements between user-identified anatomical landmarks. Inter-reader reliability between two readers was obtained using the intra-class correlation coefficient (ICC). We compared the WB and NWB measurements using paired t test. RESULTS Significant agreement was observed between two readers for both WB and NWB measurements (p <0.05). ICC values for WB and NWB measurements had a range of 50-95 and 31-71 respectively. Mean values of the medial clear space on WB images (1.75, 95% confidence interval [95% CI]: 1.6, 1.9) were significantly lower than on NWB images (2.05, 95% CI: 1.8, 2.2) measurements (p <0.001). There was no significant difference between the remaining WB and NWB measurements. CONCLUSION Measurements obtained from WB images are reliable. Except for the medial clear space, no significant difference in syndesmotic measurements were observed during the WB mode of CBCT acquisition, implying that the tibio-fibular relationship remains unchanged when the physiological axial weight-bearing load is applied.
Collapse
Affiliation(s)
- Delaram Shakoor
- Russell H. Morgan Department of Radiology and Radiological Sciences, Johns Hopkins University, 601 North Caroline Street, Baltimore, MD, USA.
| | - Greg M Osgood
- Department of Orthopedic Surgery, Johns Hopkins University, Baltimore, MD, USA
| | - Michael Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Wojciech B Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Cesar de Cesar Netto
- Department of Orthopedic Surgery, MedStar Union Memorial Hospital, Baltimore, MD, USA
| | - Babar Shafiq
- Department of Orthopedic Surgery, Johns Hopkins University, Baltimore, MD, USA
| | - Jakrapong Orapin
- Department of Orthopedic Surgery, MedStar Union Memorial Hospital, Baltimore, MD, USA
| | - Gaurav K Thawait
- Russell H. Morgan Department of Radiology and Radiological Sciences, Johns Hopkins University, 601 North Caroline Street, Baltimore, MD, USA
| | - Lew C Shon
- Department of Orthopedic Surgery, MedStar Union Memorial Hospital, Baltimore, MD, USA
| | - Shadpour Demehri
- Russell H. Morgan Department of Radiology and Radiological Sciences, Johns Hopkins University, 601 North Caroline Street, Baltimore, MD, USA
| |
Collapse
|
21
|
Brehler M, Thawait G, Kaplan J, Ramsay J, Tanaka MJ, Demehri S, Siewerdsen JH, Zbijewski W. Atlas-based algorithm for automatic anatomical measurements in the knee. J Med Imaging (Bellingham) 2019; 6:026002. [PMID: 31259202 PMCID: PMC6582228 DOI: 10.1117/1.jmi.6.2.026002] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 06/04/2019] [Indexed: 11/14/2022] Open
Abstract
We present an algorithm for automatic anatomical measurements in tomographic datasets of the knee. The algorithm uses a set of atlases, each consisting of a knee image, surface segmentations of the bones, and locations of landmarks required by the anatomical metrics. A multistage volume-to-volume and surface-to-volume registration is performed to transfer the landmarks from the atlases to the target volume. Manual segmentation of the target volume is not required in this approach. Metrics were computed from the transferred landmarks of a best-matching atlas member (different for each bone), identified based on a mutual information criterion. Leave-one-out validation of the algorithm was performed on 24 scans of the knee obtained using extremity cone-beam computed tomography. Intraclass correlation (ICC) between the algorithm and the expert who generated atlas landmarks was above 0.95 for all metrics. This compares favorably to inter-reader ICC, which varied from 0.19 to 0.95, depending on the metric. Absolute agreement with the expert was also good, with median errors below 0.25 deg for measurements of tibial slope and static alignment, and below 0.2 mm for tibial tuberosity-trochlear groove distance and medial tibial depth. The automatic approach is anticipated to improve measurement workflow and mitigate the effects of operator experience and training on reliability of the metrics.
Collapse
Affiliation(s)
- Michael Brehler
- Johns Hopkins University, Department of Biomedical Engineering, Baltimore, Maryland, United States
| | - Gaurav Thawait
- Johns Hopkins University, Russell H. Morgan Department of Radiology, Baltimore, Maryland, United States
| | - Jonathan Kaplan
- U.S. Army Natick Soldier Systems Center, Natick, Massachusetts, United States
| | - John Ramsay
- U.S. Army Natick Soldier Systems Center, Natick, Massachusetts, United States
| | - Miho J. Tanaka
- Johns Hopkins University, Department of Orthopaedic Surgery, Baltimore, Maryland, United States
| | - Shadpour Demehri
- Johns Hopkins University, Russell H. Morgan Department of Radiology, Baltimore, Maryland, United States
| | - Jeffrey H. Siewerdsen
- Johns Hopkins University, Department of Biomedical Engineering, Baltimore, Maryland, United States
- Johns Hopkins University, Russell H. Morgan Department of Radiology, Baltimore, Maryland, United States
| | - Wojciech Zbijewski
- Johns Hopkins University, Department of Biomedical Engineering, Baltimore, Maryland, United States
| |
Collapse
|
22
|
Tilley S, Jacobson M, Cao Q, Brehler M, Sisniega A, Zbijewski W, Stayman JW. Penalized-Likelihood Reconstruction With High-Fidelity Measurement Models for High-Resolution Cone-Beam Imaging. IEEE Trans Med Imaging 2018; 37:988-999. [PMID: 29621002 PMCID: PMC5889122 DOI: 10.1109/tmi.2017.2779406] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/19/2023]
Abstract
We present a novel reconstruction algorithm based on a general cone-beam CT forward model, which is capable of incorporating the blur and noise correlations that are exhibited in flat-panel CBCT measurement data. Specifically, the proposed model may include scintillator blur, focal-spot blur, and noise correlations due to light spread in the scintillator. The proposed algorithm (GPL-BC) uses a Gaussian Penalized-Likelihood objective function, which incorporates models of blur and correlated noise. In a simulation study, GPL-BC was able to achieve lower bias as compared with deblurring followed by FDK as well as a model-based reconstruction method without integration of measurement blur. In the same study, GPL-BC was able to achieve better line-pair reconstructions (in terms of segmented-image accuracy) as compared with deblurring followed by FDK, a model-based method without blur, and a model-based method with blur but not noise correlations. A prototype extremities quantitative cone-beam CT test-bench was used to image a physical sample of human trabecular bone. These data were used to compare reconstructions using the proposed method and model-based methods without blur and/or correlation to a registered CT image of the same bone sample. The GPL-BC reconstructions resulted in more accurate trabecular bone segmentation. Multiple trabecular bone metrics, including trabecular thickness (Tb.Th.) were computed for each reconstruction approach as well as the CT volume. The GPL-BC reconstruction provided the most accurate Tb.Th. measurement, 0.255 mm, as compared with the CT derived value of 0.193 mm, followed by the GPL-B reconstruction, the GPL-I reconstruction, and then the FDK reconstruction (0.271 mm, 0.309 mm, and 0.335 mm, respectively).
Collapse
|
23
|
Cao Q, Brehler M, Sisniega A, Tilley S, Shiraz Bhruwani MM, Stayman JW, Yorkston J, Siewerdsen JH, Zbijewski W. High-Resolution Extremity Cone-Beam CT with a CMOS Detector: Evaluation of a Clinical Prototype in Quantitative Assessment of Bone Microarchitecture. Proc SPIE Int Soc Opt Eng 2018; 10573:105730R. [PMID: 31346302 PMCID: PMC6657686 DOI: 10.1117/12.2293810] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
PURPOSE A prototype high-resolution extremity cone-beam CT (CBCT) system based on a CMOS detector was developed to support quantitative in vivo assessment of bone microarchitecture. We compare the performance of CMOS CBCT to an amorphous silicon (a-Si:H) FPD extremity CBCT in imaging of trabecular bone. METHODS The prototype CMOS-based CBCT involves a DALSA Xineos3030 detector (99 μm pixels) with 400 μm-thick CsI scintillator and a compact 0.3 FS rotating anode x-ray source. We compare the performance of CMOS CBCT to an a-Si:H FPD scanner built on a similar gantry, but using a Varian PaxScan2530 detector with 0.137 mm pixels and a 0.5 FS stationary anode x-ray source. Experimental studies include measurements of Modulation Transfer Function (MTF) for the detectors and in 3D image reconstructions. Image quality in clinical scenarios is evaluated in scans of a cadaver ankle. Metrics of trabecular microarchitecture (BV/TV, Bone Volume/Total Volume, TbSp, Trabecular Spacing, and TbTh, trabecular thickness) are obtained in a human ulna using CMOS CBCT and a-Si:H FPD CBCT and compared to gold standard μCT. RESULTS The CMOS detector achieves ~40% increase in the f20 value (frequency at which MTF reduces to 0.20) compared to the a-Si:H FPD. In the reconstruction domain, the FWHM of a 127 μm tungsten wire is also improved by ~40%. Reconstructions of a cadaveric ankle reveal enhanced modulation of trabecular structures with the CMOS detector and soft-tissue visibility that is similar to that of the a-Si:H FPD system. Correlations of the metrics of bone microarchitecture with gold-standard μCT are improved with CMOS CBCT: from 0.93 to 0.98 for BV/TV, from 0.49 to 0.74 for TbTh, and from 0.9 to 0.96 for TbSp. CONCLUSION Adoption of a CMOS detector in extremity CBCT improved spatial resolution and enhanced performance in metrics of bone microarchitecture compared to a conventional a-Si:H FPD. The results support development of clinical applications of CMOS CBCT in quantitative imaging of bone health.
Collapse
Affiliation(s)
- Q Cao
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - M Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - A Sisniega
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - S Tilley
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - M M Shiraz Bhruwani
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - J W Stayman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | | | - J H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
- Russell H. Morgan Department of Radiology, Johns Hopkins University, Baltimore, MD USA 21287
| | - W Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| |
Collapse
|
24
|
Brehler M, Cao Q, Moseley KF, Osgood G, Morris C, Demehri S, Yorkston J, Siewerdsen JH, Zbijewski W. Robust Quantitative Assessment of Trabecular Microarchitecture in Extremity Cone-Beam CT Using Optimized Segmentation Algorithms. Proc SPIE Int Soc Opt Eng 2018; 10578. [PMID: 31337926 DOI: 10.1117/12.2293346] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Purpose In-vivo evaluation of bone microarchitecture remains challenging because of limited resolution of conventional orthopaedic imaging modalities. We investigate the performance of flat-panel detector extremity Cone-Beam CT (CBCT) in quantitative analysis of trabecular bone. To enable accurate morphometry of fine trabecular bone architecture, advanced CBCT pre-processing and segmentation algorithms are developed. Methods The study involved 35 transilliac bone biopsy samples imaged on extremity CBCT (voxel size 75 μm, imaging dose ~13 mGy) and gold standard μCT (voxel size 7.67 μm). CBCT image segmentation was performed using (i) global Otsu's thresholding, (ii) Bernsen's local thresholding, (iii) Bernsen's local thresholding with additional histogram-based global pre-thresholding, and (iv) the same as (iii) but combined with contrast enhancement using a Laplacian Pyramid. Correlations between extremity CBCT with the different segmentation algorithms and gold standard μCT were investigated for measurements of Bone Volume over Total Volume (BV/TV), Trabecular Thickness (Tb.Th), Trabecular Spacing (Tb.Sp), and Trabecular Number (Tb.N). Results The combination of local thresholding with global pre-thresholding and Laplacian contrast enhancement outperformed other CBCT segmentation methods. Using this optimal segmentation scheme, strong correlation between extremity CBCT and μCT was achieved, with Pearson coefficients of 0.93 for BV/TV, 0.89 for Tb.Th, 0.91 for Tb.Sp, and 0.88 for Tb.N (all results statistically significant). Compared to a simple global CBCT segmentation using Otsu's algorithm, the advanced segmentation method achieved ~20% improvement in the correlation coefficient for Tb.Th and ~50% improvement for Tb.Sp. Conclusions Extremity CBCT combined with advanced image pre-processing and segmentation achieves high correlation with gold standard μCT in measurements of trabecular microstructure. This motivates ongoing development of clinical applications of extremity CBCT in in-vivo evaluation of bone health e.g. in early osteoarthritis and osteoporosis.
Collapse
Affiliation(s)
- M Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - Q Cao
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - K F Moseley
- Division of Endocrinology, Diabetes and Metabolism, Johns Hopkins University, Baltimore, MD USA
| | - G Osgood
- Department of Orthopedics, Johns Hopkins University, Baltimore, MD USA
| | - C Morris
- Department of Orthopedics, Johns Hopkins University, Baltimore, MD USA
| | - S Demehri
- Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | | | - J H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA.,Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - W Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| |
Collapse
|
25
|
Cao Q, Sisniega A, Brehler M, Stayman JW, Yorkston J, Siewerdsen JH, Zbijewski W. Modeling and evaluation of a high-resolution CMOS detector for cone-beam CT of the extremities. Med Phys 2017; 45:114-130. [PMID: 29095489 DOI: 10.1002/mp.12654] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2017] [Revised: 10/19/2017] [Accepted: 10/23/2017] [Indexed: 12/12/2022] Open
Abstract
PURPOSE Quantitative assessment of trabecular bone microarchitecture in extremity cone-beam CT (CBCT) would benefit from the high spatial resolution, low electronic noise, and fast scan time provided by complementary metal-oxide semiconductor (CMOS) x-ray detectors. We investigate the performance of CMOS sensors in extremity CBCT, in particular with respect to potential advantages of thin (<0.7 mm) scintillators offering higher spatial resolution. METHODS A cascaded systems model of a CMOS x-ray detector incorporating the effects of CsI:Tl scintillator thickness was developed. Simulation studies were performed using nominal extremity CBCT acquisition protocols (90 kVp, 0.126 mAs/projection). A range of scintillator thickness (0.35-0.75 mm), pixel size (0.05-0.4 mm), focal spot size (0.05-0.7 mm), magnification (1.1-2.1), and dose (15-40 mGy) was considered. The detectability index was evaluated for both CMOS and a-Si:H flat-panel detector (FPD) configurations for a range of imaging tasks emphasizing spatial frequencies associated with feature size aobj. Experimental validation was performed on a CBCT test bench in the geometry of a compact orthopedic CBCT system (SAD = 43.1 cm, SDD = 56.0 cm, matching that of the Carestream OnSight 3D system). The test-bench studies involved a 0.3 mm focal spot x-ray source and two CMOS detectors (Dalsa Xineos-3030HR, 0.099 mm pixel pitch) - one with the standard CsI:Tl thickness of 0.7 mm (C700) and one with a custom 0.4 mm thick scintillator (C400). Measurements of modulation transfer function (MTF), detective quantum efficiency (DQE), and CBCT scans of a cadaveric knee (15 mGy) were obtained for each detector. RESULTS Optimal detectability for high-frequency tasks (feature size of ~0.06 mm, consistent with the size of trabeculae) was ~4× for the C700 CMOS detector compared to the a-Si:H FPD at nominal system geometry of extremity CBCT. This is due to ~5× lower electronic noise of a CMOS sensor, which enables input quantum-limited imaging at smaller pixel size. Optimal pixel size for high-frequency tasks was <0.1 mm for a CMOS, compared to ~0.14 mm for an a-Si:H FPD. For this fine pixel pitch, detectability of fine features could be improved by using a thinner scintillator to reduce light spread blur. A 22% increase in detectability of 0.06 mm features was found for the C400 configuration compared to C700. An improvement in the frequency at 50% modulation (f50 ) of MTF was measured, increasing from 1.8 lp/mm for C700 to 2.5 lp/mm for C400. The C400 configuration also achieved equivalent or better DQE as C700 for frequencies above ~2 mm-1 . Images of cadaver specimens confirmed improved visualization of trabeculae with the C400 sensor. CONCLUSIONS The small pixel size of CMOS detectors yields improved performance in high-resolution extremity CBCT compared to a-Si:H FPDs, particularly when coupled with a custom 0.4 mm thick scintillator. The results indicate that adoption of a CMOS detector in extremity CBCT can benefit applications in quantitative imaging of trabecular microstructure in humans.
Collapse
Affiliation(s)
- Qian Cao
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21205, USA
| | - Alejandro Sisniega
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21205, USA
| | - Michael Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21205, USA
| | - J Webster Stayman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21205, USA
| | | | - Jeffrey H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21205, USA.,Russell H Morgan Department of Radiology, Johns Hopkins University, Baltimore, 21205, USA
| | - Wojciech Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21205, USA
| |
Collapse
|
26
|
Brehler M, Thawait G, Shyr W, Ramsay J, Siewerdsen JH, Zbijewski W. Atlas-based automatic measurements of the morphology of the tibiofemoral joint. Proc SPIE Int Soc Opt Eng 2017. [PMID: 28638170 DOI: 10.1117/12.2255566] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
PURPOSE Anatomical metrics of the tibiofemoral joint support assessment of joint stability and surgical planning. We propose an automated, atlas-based algorithm to streamline the measurements in 3D images of the joint and reduce user-dependence of the metrics arising from manual identification of the anatomical landmarks. METHODS The method is initialized with coarse registrations of a set of atlas images to the fixed input image. The initial registrations are then refined separately for the tibia and femur and the best matching atlas is selected. Finally, the anatomical landmarks of the best matching atlas are transformed onto the input image by deforming a surface model of the atlas to fit the shape of the tibial plateau in the input image (a mesh-to-volume registration). We apply the method to weight-bearing volumetric images of the knee obtained from 23 subjects using an extremity cone-beam CT system. Results of the automated algorithm were compared to an expert radiologist for measurements of Static Alignment (SA), Medial Tibial Slope (MTS) and Lateral Tibial Slope (LTS). RESULTS Intra-reader variability as high as ~10% for LTS and 7% for MTS (ratio of standard deviation to the mean in repeated measurements) was found for expert radiologist, illustrating the potential benefits of an automated approach in improving the precision of the metrics. The proposed method achieved excellent registration of the atlas mesh to the input volumes. The resulting automated measurements yielded high correlations with expert radiologist, as indicated by correlation coefficients of 0.72 for MTS, 0.8 for LTS, and 0.89 for SA. CONCLUSIONS The automated method for measurement of anatomical metrics of the tibiofemoral joint achieves high correlation with expert radiologist without the need for time consuming and error prone manual selection of landmarks.
Collapse
Affiliation(s)
- M Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - G Thawait
- Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - W Shyr
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| | - J Ramsay
- Natick Soldier Research, Development and Engineering Center (NSRDEC), Natick, MA USA
| | - J H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA.,Department of Radiology, Johns Hopkins University, Baltimore, MD USA
| | - W Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA
| |
Collapse
|
27
|
Cao Q, Brehler M, Sisniega A, Stayman JW, Yorkston J, Siewerdsen JH, Zbijewski W. High-resolution extremity cone-beam CT with a CMOS detector: Task-based optimization of scintillator thickness. Proc SPIE Int Soc Opt Eng 2017; 10132:1013210. [PMID: 28989220 PMCID: PMC5630149 DOI: 10.1117/12.2255695] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
PURPOSE CMOS x-ray detectors offer small pixel sizes and low electronic noise that may support the development of novel high-resolution imaging applications of cone-beam CT (CBCT). We investigate the effects of CsI scintillator thickness on the performance of CMOS detectors in high resolution imaging tasks, in particular in quantitative imaging of bone microstructure in extremity CBCT. METHODS A scintillator thickness-dependent cascaded systems model of CMOS x-ray detectors was developed. Detectability in low-, high- and ultra-high resolution imaging tasks (Gaussian with FWHM of ~250 μm, ~80 μm and ~40 μm, respectively) was studied as a function of scintillator thickness using the theoretical model. Experimental studies were performed on a CBCT test bench equipped with DALSA Xineos3030 CMOS detectors (99 μm pixels) with CsI scintillator thicknesses of 400 μm and 700 μm, and a 0.3 FS compact rotating anode x-ray source. The evaluation involved a radiographic resolution gauge (0.6-5.0 lp/mm), a 127 μm tungsten wire for assessment of 3D resolution, a contrast phantom with tissue-mimicking inserts, and an excised fragment of human tibia for visual assessment of fine trabecular detail. RESULTS Experimental studies show ~35% improvement in the frequency of 50% MTF modulation when using the 400 μm scintillator compared to the standard nominal CsI thickness of 700 μm. Even though the high-frequency DQE of the two detectors is comparable, theoretical studies show a 14% to 28% increase in detectability index (d'2) of high- and ultrahigh resolution tasks, respectively, for the detector with 400 μm CsI compared to 700 μm CsI. Experiments confirm the theoretical findings, showing improvements with the adoption of 400 μm panel in the visibility of the radiographic pattern (2× improvement in peak-to-through distance at 4.6 lp/mm) and a 12.5% decrease in the FWHM of the tungsten wire. Reconstructions of the tibial plateau reveal enhanced visibility of trabecular structures with the CMOS detector with 400 μm scinitllator. CONCLUSION Applications on CMOS detectors in high resolution CBCT imaging of trabecular bone will benefit from using a thinner scintillator than the current standard in general radiography. The results support the translation of the CMOS sensor with 400 μm CsI onto the clinical prototype of CMOS-based extremity CBCT.
Collapse
Affiliation(s)
- Q Cao
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - M Brehler
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - A Sisniega
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | - J W Stayman
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| | | | - J H Siewerdsen
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
- Russell H. Morgan Department of Radiology, Johns Hopkins University, Baltimore, MD USA 21287
| | - W Zbijewski
- Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD USA 21205
| |
Collapse
|
28
|
Cao Q, Brehler M, Sisniega A, Marinetto E, Zyazin A, Peters I, Stayman J, Yorkston J, Siewerdsen J, Zbijewski W. WE-AB-207A-01: BEST IN PHYSICS (IMAGING): High-Resolution Cone-Beam CT of the Extremities and Cancellous Bone Architecture with a CMOS Detector. Med Phys 2016. [DOI: 10.1118/1.4957754] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
|
29
|
Görres J, Brehler M, Franke J, Vetter SY, Grützner PA, Meinzer HP, Wolf I. Articular surface segmentation using active shape models for intraoperative implant assessment. Int J Comput Assist Radiol Surg 2016; 11:1661-72. [DOI: 10.1007/s11548-015-1316-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2015] [Accepted: 10/13/2015] [Indexed: 10/21/2022]
|
30
|
Brehler M, Görres J, Franke J, Barth K, Vetter SY, Grützner PA, Meinzer HP, Wolf I, Nabers D. Intra-operative adjustment of standard planes in C-arm CT image data. Int J Comput Assist Radiol Surg 2015; 11:495-504. [PMID: 26316065 DOI: 10.1007/s11548-015-1281-3] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2015] [Accepted: 08/13/2015] [Indexed: 11/26/2022]
Abstract
PURPOSE With the help of an intra-operative mobile C-arm CT, medical interventions can be verified and corrected, avoiding the need for a post-operative CT and a second intervention. An exact adjustment of standard plane positions is necessary for the best possible assessment of the anatomical regions of interest but the mobility of the C-arm causes the need for a time-consuming manual adjustment. In this article, we present an automatic plane adjustment at the example of calcaneal fractures. METHODS We developed two feature detection methods (2D and pseudo-3D) based on SURF key points and also transferred the SURF approach to 3D. Combined with an atlas-based registration, our algorithm adjusts the standard planes of the calcaneal C-arm images automatically. The robustness of the algorithms is evaluated using a clinical data set. Additionally, we tested the algorithm's performance for two registration approaches, two resolutions of C-arm images and two methods for metal artifact reduction. RESULTS For the feature extraction, the novel 3D-SURF approach performs best. As expected, a higher resolution ([Formula: see text] voxel) leads also to more robust feature points and is therefore slightly better than the [Formula: see text] voxel images (standard setting of device). Our comparison of two different artifact reduction methods and the complete removal of metal in the images shows that our approach is highly robust against artifacts and the number and position of metal implants. CONCLUSIONS By introducing our fast algorithmic processing pipeline, we developed the first steps for a fully automatic assistance system for the assessment of C-arm CT images.
Collapse
Affiliation(s)
- Michael Brehler
- Division of Medical and Biological Informatics (E130), German Cancer Research Center, Im Neuenheimer Feld 280, 69120, Heidelberg, Germany.
| | - Joseph Görres
- Division of Medical and Biological Informatics (E130), German Cancer Research Center, Im Neuenheimer Feld 280, 69120, Heidelberg, Germany
| | - Jochen Franke
- BG Trauma Center, Ludwig-Guttmann-Straße 13, 67071, Ludwigshafen am Rhein, Germany.
| | - Karl Barth
- Siemens AG, Healthcare Sector, Henkestr. 127, 91052, Erlangen, Germany.
| | - Sven Y Vetter
- BG Trauma Center, Ludwig-Guttmann-Straße 13, 67071, Ludwigshafen am Rhein, Germany
| | - Paul A Grützner
- BG Trauma Center, Ludwig-Guttmann-Straße 13, 67071, Ludwigshafen am Rhein, Germany
| | - Hans-Peter Meinzer
- Division of Medical and Biological Informatics (E130), German Cancer Research Center, Im Neuenheimer Feld 280, 69120, Heidelberg, Germany
| | - Ivo Wolf
- Mannheim University of Applied Sciences, Paul-Wittsack-Str. 10, 68163, Mannheim, Germany.
| | - Diana Nabers
- Division of Medical and Biological Informatics (E130), German Cancer Research Center, Im Neuenheimer Feld 280, 69120, Heidelberg, Germany
| |
Collapse
|