The number of microspheres in Y90 radioembolization directly affects normal tissue radiation exposure.

PURPOSE: In Y90 radioembolization, the number of microspheres infused varies by more than a factor of 20 over the shelf-life of the glass radioembolization device. We investigated the effect of the number of Y90 microspheres on normal liver tissue. METHOD: Healthy pigs received lobar radioembolization with glass Y90 microspheres at 4, 8, 12, and 16 days post-calibration, representing a > 20× range in the number of microspheres deposited per milliliter in tissue. Animals were survived for 1-month post-treatment and the livers were explanted and scanned on a micro CT system to fully characterize the microscopic distribution of individual microspheres. A complete 3D microdosimetric evaluation of each liver was performed with a spatially correlated analysis of histopathologic effect. RESULTS: Through whole-lobe microscopic identification of each microsphere, a consistent number of microspheres per sphere cluster was found at 4, 8, and 12 days postcalibration, despite an 8-fold increase in total microspheres infused from days 4 to 12. The additional microspheres instead resulted in more clusters formed and, therefore, a more homogeneous microscopic absorbed dose. The increased absorbed-dose homogeneity resulted in a greater volume fraction of the liver receiving a potentially toxic absorbed dose based on radiobiologic models. Histopathologic findings in the animals support a possible increase in normal liver toxicity in later treatments with more spheres (i.e., ≥ day 12) compared to early treatments with less spheres (i.e., ≤ day 8). CONCLUSION: The microdosimetric evidence presented supports a recommendation of caution when treating large volumes (e.g., right lobe) using glass 90Y microspheres at more than 8 days post-calibration, i.e., after "2nd week" Monday. The favorable normal tissue microscopic distribution and associated low toxicity of first week therapies may encourage opportunities for dose escalation with glass microspheres and could also be considered for patients with decreased hepatic reserve.

Development of a Customizable Hepatic Arterial Tree and Particle Transport Model for Use in Treatment Planning.

Optimal treatment planning for radioembolization of hepatic cancers produces sufficient dose to tumors for control and dose to normal liver parenchyma that is below the threshold for toxicity. The non-uniform distribution of particles in liver microanatomy complicates the planning process as different functional regions receive different doses. Having realistic and patient-specific models of the arterial tree and microsphere trapping would be useful for developing more optimal treatment plans. We propose a macrocell-based growth method to generate models of the hepatic arterial tree from the proper hepatic artery to the terminal arterioles supplying the capillaries in the parenchyma. We show how these trees can be adapted to match patient values of pressure, flow, and vessel diameters while still conforming to laws controlling vessel bifurcation, changes in pressure, and blood flow. We also introduce a method to model particle transport within the tree that accounts for vessel and particle diameter distributions and show the non-uniform microsphere deposition pattern that results. Potential applications include investigating dose heterogeneity and microsphere deposition patterns.

Verification of a method to detect glass microspheres via micro-CT.

INTRODUCTION: The ability to determine the microscopic distribution of glass microspheres in 90 Y radioembolization has important applications in post-treatment microdosimetry and cluster analysis. Current methods are time-intensive and labor-intensive and thus are typically only applied to small samples. MATERIALS AND METHODS: A high-resolution micro-CT image with a voxel size of 8.74 µm was acquired of phantoms containing ~25 µm-diameter glass microspheres embedded in tissue-equivalent materials that were optically transparent, which allowed true microsphere locations to be determined using transmission light microscopy. A 3-stage algorithm was developed to estimate the number and locations of microspheres in tissue regions. The stages are thresholding the CT image and discarding regions with insufficient voxels, estimating the number of microspheres in each region using the values of the detected and neighboring region voxels and estimating locations for each microsphere using the outputs of the previous two stages. Two different methods for estimating the number of microspheres in each region were derived, as were five methods for localizing microspheres. Metrics for each stage were computed, and the mean absolute error (MAE) between the dose to 72 µm voxels of the true and estimated dose maps created from the microsphere locations was used as the figure of merit for overall algorithm performance. Microsphere locations identified in the optical micrograph were used as the gold standard for the metrics of all stages. The method's utility was then demonstrated using a specimen from a human neuroendocrine tumor (NET) treated with glass 90 Y microspheres. RESULTS: The stage detecting regions containing microspheres found 100% of microspheres inside regions. The number of incorrectly detected regions without microspheres was 1.5% of the total number of regions. In stage 2, with these parameters, nearly 94% of the actual number of spheres in each region was correctly counted, and only 5% of the estimated sphere quantities in each region were false positives. The MAE between the true dose maps and dose maps estimated using the full algorithm with optimal parameter and method choices was 4.2%. A total of 5,713 glass microspheres were identified as being distributed heterogeneously in the NET specimen with a maximum tumor dose of >2500 Gy and 46% of the specimen receiving <20 Gy. CONCLUSIONS: This work developed and evaluated a method to detect and estimate the three-dimensional locations of glass microspheres in whole tissue samples that require less manual effort than traditional methods. This method could be used to gain important insights into the heterogeneity of microsphere distributions that would be useful for improving radioembolization treatment planning.