A deformable atlas of the laboratory mouse

Learning a stable approximation of an existing but unknown inverse mapping: application to the half-time circular Radon transform

Supervised deep learning-based methods have inspired a new wave of image reconstruction methods that implicitly learn effective regularization strategies from a set of training data. While they hold potential for improving image quality, they have also raised concerns regarding their robustness. Instabilities can manifest when learned methods are applied to find approximate solutions to ill-posed image reconstruction problems for which a unique and stable inverse mapping does not exist, which is a typical use case. In this study, we investigate the performance of supervised deep learning-based image reconstruction in an alternate use case in which a stable inverse mapping is known to exist but is not yet analytically available in closed form. For such problems, a deep learning-based method can learn a stable approximation of the unknown inverse mapping that generalizes well to data that differ significantly from the training set. The learned approximation of the inverse mapping eliminates the need to employ an implicit (optimization-based) reconstruction method and can potentially yield insights into the unknown analytic inverse formula. The specific problem addressed is image reconstruction from a particular case of radially truncated circular Radon transform (CRT) data, referred to as 'half-time' measurement data. For the half-time image reconstruction problem, we develop and investigate a learned filtered backprojection method that employs a convolutional neural network to approximate the unknown filtering operation. We demonstrate that this method behaves stably and readily generalizes to data that differ significantly from training data. The developed method may find application to wave-based imaging modalities that include photoacoustic computed tomography.

Bioluminescence tomography reconstruction in conjunction with an organ probability map as an anatomical reference

To alleviate the ill-posedness of bioluminescence tomography (BLT) reconstruction, anatomical information from computed tomography (CT) or magnetic resonance imaging (MRI) is usually adopted to improve the reconstruction quality. With the anatomical information, different organs could be segmented and assigned with appropriate optical parameters, and the reconstruction could be confined into certain organs. However, image segmentation is a time-consuming and challenging work, especially for the low-contrast organs. In this paper, we present a BLT reconstruction method in conjunction with an organ probability map to effectively incorporate the anatomical information. Instead of using a segmentation with a fixed organ map, an organ probability map is established by registering the CT image of the mouse to the statistical mouse atlas with the constraints of the mouse surface and high-contrast organs (bone and lung). Then the organ probability map of the low-contrast organs, such as the liver and kidney, is determined automatically. After discretization of the mouse torso, a heterogeneous model is established as the input for reconstruction, in which the optical parameter of each node is calculated according to the organ probability map. To take the advantage of the sparse Bayesian Learning (SBL) method in recovering block sparse signals in inverse problems, which is common in BLT applications where the target distribution has the characteristic of sparsity and block structure, a two-step method in conjunction with the organ probability map is presented. In the first step, a fast sparse algorithm, L1-LS, is used to reveal the source distribution on the organ level. In the second step, the bioluminescent source is reconstructed on the pixel level based on the SBL method. Both simulation and in vivo experiments are conducted, and the results demonstrate that the organ probability map in conjunction with the proposed two-step BLT reconstruction method is feasible to accurately reconstruct the localization of the bioluminescent light source.

Deep learning-based segmentation of the thorax in mouse micro-CT scans

For image-guided small animal irradiations, the whole workflow of imaging, organ contouring, irradiation planning, and delivery is typically performed in a single session requiring continuous administration of anaesthetic agents. Automating contouring leads to a faster workflow, which limits exposure to anaesthesia and thereby, reducing its impact on experimental results and on animal wellbeing. Here, we trained the 2D and 3D U-Net architectures of no-new-Net (nnU-Net) for autocontouring of the thorax in mouse micro-CT images. We trained the models only on native CTs and evaluated their performance using an independent testing dataset (i.e., native CTs not included in the training and validation). Unlike previous studies, we also tested the model performance on an external dataset (i.e., contrast-enhanced CTs) to see how well they predict on CTs completely different from what they were trained on. We also assessed the interobserver variability using the generalized conformity index ([Formula: see text]) among three observers, providing a stronger human baseline for evaluating automated contours than previous studies. Lastly, we showed the benefit on the contouring time compared to manual contouring. The results show that 3D models of nnU-Net achieve superior segmentation accuracy and are more robust to unseen data than 2D models. For all target organs, the mean surface distance (MSD) and the Hausdorff distance (95p HD) of the best performing model for this task (nnU-Net 3d_fullres) are within 0.16 mm and 0.60 mm, respectively. These values are below the minimum required contouring accuracy of 1 mm for small animal irradiations, and improve significantly upon state-of-the-art 2D U-Net-based AIMOS method. Moreover, the conformity indices of the 3d_fullres model also compare favourably to the interobserver variability for all target organs, whereas the 2D models perform poorly in this regard. Importantly, the 3d_fullres model offers 98% reduction in contouring time.

A preclinical micro-computed tomography database including 3D whole body organ segmentations

The gold-standard of preclinical micro-computed tomography (μCT) data processing is still manual delineation of complete organs or regions by specialists. However, this method is time-consuming, error-prone, has limited reproducibility, and therefore is not suitable for large-scale data analysis. Unfortunately, robust and accurate automated whole body segmentation algorithms are still missing. In this publication, we introduce a database containing 225 murine 3D whole body μCT scans along with manual organ segmentation of most important organs including heart, liver, lung, trachea, spleen, kidneys, stomach, intestine, bladder, thigh muscle, bone, as well as subcutaneous tumors. The database includes native and contrast-enhanced, regarding spleen and liver, μCT data. All scans along with organ segmentation are freely accessible at the online repository Figshare. We encourage researchers to reuse the provided data to evaluate and improve methods and algorithms for accurate automated organ segmentation which may reduce manual segmentation effort, increase reproducibility, and even reduce the number of required laboratory animals by reducing a source of variability and having access to a reliable reference group.

Automated quantification of bioluminescence images

We developed a computer-aided analysis tool for quantitatively determining bioluminescent reporter distributions inside small animals. The core innovations are a body-fitting animal shuttle and a statistical mouse atlas, both of which are spatially aligned and scaled according to the animal’s weight, and hence provide data congruency across animals of varying size and pose. In conjunction with a multispectral bioluminescence tomography technique capitalizing on the spatial framework of the shuttle, the in vivo biodistribution of luminescent reporters can rapidly be calculated and, thus, enables operator-independent and computer-driven data analysis. We demonstrate its functionality by quantitatively monitoring a bacterial infection, where the bacterial organ burden was determined and validated with the established serial-plating method. In addition, the statistical mouse atlas was validated and compared to existing techniques providing an anatomical reference. The proposed data analysis tool promises to increase data throughput and data reproducibility and accelerate human disease modeling in mice.

Analysis of bioluminescence images of bacterial distributions in living animals is mostly manual and semiquantitative. Here, the authors present an analysis platform featuring an animal mold, a probabilistic organ atlas, and a mirror gantry to perform automatic in vivo bioluminescence quantification.

Bioluminescence tomography with structural information estimated via statistical mouse atlas registration

Due to an ill-posed and underestimated characteristic of bioluminescence tomography (BLT) reconstruction, a priori anatomical information obtained from computed tomography (CT) or magnetic resonance imaging (MRI), is usually incorporated to improve the reconstruction accuracy. The organs need to be segmented, which is time-consuming and challenging, especially for the low-contrast CT images. In this paper, we present a BLT reconstruction method based on a statistical mouse atlas to improve the efficiency of heterogeneous model generation and the accuracy of target localization. The low-contrast CT image of the mouse was first registered to the statistical mouse atlas model with the constraints of mouse surface and high-contrast organs (bone and lung). Then the other organs, such as the liver and kidney, were determined automatically through the statistical mouse atlas model. The estimated organs were then discretized into tetrahedral meshes for BLT reconstruction. The linearized Bregman method was used to solve the sparse inverse problem of BLT by minimizing the regularization function (L1 norm plus L2 norm with smooth factor). Both numerical simulations and in vivo experiments were conducted, and the results demonstrate that even though the localization of the estimated organs may not be exactly accurate, the proposed method is feasible to reconstruct the bioluminescent source effectively and accurately with the estimated organs. This method would greatly benefit the bioluminescent light source localization for hybrid BLT/CT systems.

Deformable torso phantoms of Chinese adults for personalized anatomy modelling

In recent years, there has been increasing demand for personalized anatomy modelling for medical and industrial applications, such as ergonomics device development, clinical radiological exposure simulation, biomechanics analysis, and 3D animation character design. In this study, we constructed deformable torso phantoms that can be deformed to match the personal anatomy of Chinese male and female adults. The phantoms were created based on a training set of 79 trunk computed tomography (CT) images (41 males and 38 females) from normal Chinese subjects. Major torso organs were segmented from the CT images, and the statistical shape model (SSM) approach was used to learn the inter-subject anatomical variations. To match the personal anatomy, the phantoms were registered to individual body surface scans or medical images using the active shape model method. The constructed SSM demonstrated anatomical variations in body height, fat quantity, respiratory status, organ geometry, male muscle size, and female breast size. The masses of the deformed phantom organs were consistent with Chinese population organ mass ranges. To validate the performance of personal anatomy modelling, the phantoms were registered to the body surface scan and CT images. The registration accuracy measured from 22 test CT images showed a median Dice coefficient over 0.85, a median volume recovery coefficient (RCvlm ) between 0.85 and 1.1, and a median averaged surface distance (ASD) < 1.5 mm. We hope these phantoms can serve as computational tools for personalized anatomy modelling for the research community.

Metabolic profiling of follistatin overexpression: a novel therapeutic strategy for metabolic diseases

BACKGROUND

Follistatin (Fst) promotes brown adipocyte characteristics in adipose tissues.

METHODS

Abdominal fat volume (CT scan), glucose clearance (GTT test), and metabolomics analysis (mass spectrometry) of adipose tissues from Fst transgenic (Fst-Tg) and wild type (WT) control mice were analyzed. Oxygen consumption (Seahorse Analyzer) and lipidomics (gas chromatography) was analyzed in 3T3-L1 cells.

RESULTS

Fst-Tg mice show significant decrease in abdominal fat content, increased glucose clearance, improved plasma lipid profiles and significant changes in several conventional metabolites compared to the WT mice. Furthermore, overexpression of Fst in 3T3-L1 cells resulted in up regulation of key brown/beige markers and changes in lipidomics profiles.

CONCLUSION

Fst modulates key factors involved in promoting metabolic syndrome and could be used for therapeutic intervention.

Estimation of subject coregistration errors during multimodal preclinical imaging using separate instruments: origins and avoidance of artifacts

We use high-resolution [Formula: see text] data in multiple experiments to estimate the sources of error during coregistration of images acquired on separate preclinical instruments. In combination with experiments with phantoms, we completed in vivo imaging on mice, aimed at identifying the possible sources of registration errors, caused either by transport of the animal, movement of the animal itself, or methods of coregistration. The same imaging cell was used as a holder for phantoms and animals. For all procedures, rigid coregistration was carried out using a common landmark coregistration system, placed inside the imaging cell. We used the fiducial registration error and the target registration error to analyze the coregistration accuracy. We found that moving an imaging cell between two preclinical devices during a multimodal procedure gives an error of about [Formula: see text] at most. Therefore, it could not be considered a source of coregistration errors. Errors linked to spontaneous movements of the animal increased with time, to nearly 1 mm at most, excepted for body parts that were properly restrained. This work highlights the importance of animal intrinsic movements during a multiacquisition procedure and demonstrates a simple method to identify and quantify the sources of error during coregistration.

PET Image Reconstruction and Deformable Motion Correction Using Unorganized Point Clouds

Quantitative positron emission tomography imaging often requires correcting the image data for deformable motion. With cyclic motion, this is traditionally achieved by separating the coincidence data into a relatively small number of gates, and incorporating the inter-gate image transformation matrices into the reconstruction algorithm. In the presence of non-cyclic deformable motion, this approach may be impractical due to a large number of required gates. In this paper, we propose an alternative approach to iterative image reconstruction with correction for deformable motion, wherein unorganized point clouds are used to model the imaged objects in the image space, and motion is corrected for explicitly by introducing a time-dependence into the point coordinates. The image function is represented using constant basis functions with finite support determined by the boundaries of the Voronoi cells in the point cloud. We validate the quantitative accuracy and stability of the proposed approach by reconstructing noise-free and noisy projection data from digital and physical phantoms. The point-cloud-based maximum likelihood expectation maximization (MLEM) and one-pass list-mode ordered-subset expectation maximization (OSEM) algorithms are validated. The results demonstrate that images reconstructed using the proposed method are quantitatively stable, with noise and convergence properties comparable to image reconstruction based on the use of rectangular and radially-symmetric basis functions.

Noninvasive Assessment of Elimination and Retention using CT-FMT and Kinetic Whole-body Modeling

Fluorescence-mediated tomography (FMT) is a quantitative three-dimensional imaging technique for preclinical research applications. The combination with micro-computed tomography (µCT) enables improved reconstruction and analysis. The aim of this study is to assess the potential of µCT-FMT and kinetic modeling to determine elimination and retention of typical model drugs and drug delivery systems.

We selected four fluorescent probes with different but well-known biodistribution and elimination routes: Indocyanine green (ICG), hydroxyapatite-binding OsteoSense (OS), biodegradable nanogels (NG) and microbubbles (MB). µCT-FMT scans were performed in twenty BALB/c nude mice (5 per group) at 0.25, 2, 4, 8, 24, 48 and 72 h after intravenous injection. Longitudinal organ curves were determined using interactive organ segmentation software and a pharmacokinetic whole-body model was implemented and applied to compute physiological parameters describing elimination and retention.

ICG demonstrated high initial hepatic uptake which decreased rapidly while intestinal accumulation appeared for around 8 hours which is in line with the known direct uptake by hepatocytes followed by hepatobiliary elimination. Complete clearance from the body was observed at 48 h. NG showed similar but slower hepatobiliary elimination because these nanoparticles require degradation before elimination can take place. OS was strongly located in the bones in addition to high signal in the bladder at 0.25 h indicating fast renal excretion. MB showed longest retention in liver and spleen and low signal in the kidneys likely caused by renal elimination or retention of fragments. Furthermore, probe retention was found in liver (MB, NG and OS), spleen (MB) and kidneys (MB and NG) at 72 h which was confirmed by ex vivo data. The kinetic model enabled robust extraction of physiological parameters from the organ curves.

In summary, µCT-FMT and kinetic modeling enable differentiation of hepatobiliary and renal elimination routes and allow for the noninvasive assessment of retention sites in relevant organs including liver, kidney, bone and spleen.

A Novel Mouse Segmentation Method Based on Dynamic Contrast Enhanced Micro-CT Images

With the development of hybrid imaging scanners, micro-CT is widely used in locating abnormalities, studying drug metabolism, and providing structural priors to aid image reconstruction in functional imaging. Due to the low contrast of soft tissues, segmentation of soft tissue organs from mouse micro-CT images is a challenging problem. In this paper, we propose a mouse segmentation scheme based on dynamic contrast enhanced micro-CT images. With a homemade fast scanning micro-CT scanner, dynamic contrast enhanced images were acquired before and after injection of non-ionic iodinated contrast agents (iohexol). Then the feature vector of each voxel was extracted from the signal intensities at different time points. Based on these features, the heart, liver, spleen, lung, and kidney could be classified into different categories and extracted from separate categories by morphological processing. The bone structure was segmented using a thresholding method. Our method was validated on seven BALB/c mice using two different classifiers: a support vector machine classifier with a radial basis function kernel and a random forest classifier. The results were compared to manual segmentation, and the performance was assessed using the Dice similarity coefficient, false positive ratio, and false negative ratio. The results showed high accuracy with the Dice similarity coefficient ranging from 0.709 ± 0.078 for the spleen to 0.929 ± 0.006 for the kidney.

Development of computational small animal models and their applications in preclinical imaging and therapy research

The development of multimodality preclinical imaging techniques and the rapid growth of realistic computer simulation tools have promoted the construction and application of computational laboratory animal models in preclinical research. Since the early 1990s, over 120 realistic computational animal models have been reported in the literature and used as surrogates to characterize the anatomy of actual animals for the simulation of preclinical studies involving the use of bioluminescence tomography, fluorescence molecular tomography, positron emission tomography, single-photon emission computed tomography, microcomputed tomography, magnetic resonance imaging, and optical imaging. Other applications include electromagnetic field simulation, ionizing and nonionizing radiation dosimetry, and the development and evaluation of new methodologies for multimodality image coregistration, segmentation, and reconstruction of small animal images. This paper provides a comprehensive review of the history and fundamental technologies used for the development of computational small animal models with a particular focus on their application in preclinical imaging as well as nonionizing and ionizing radiation dosimetry calculations. An overview of the overall process involved in the design of these models, including the fundamental elements used for the construction of different types of computational models, the identification of original anatomical data, the simulation tools used for solving various computational problems, and the applications of computational animal models in preclinical research. The authors also analyze the characteristics of categories of computational models (stylized, voxel-based, and boundary representation) and discuss the technical challenges faced at the present time as well as research needs in the future.