Dynamic contrast-enhanced MRI-based early detection of acute renal transplant rejection

Kidney Segmentation from Dynamic Contrast-Enhanced Magnetic Resonance Imaging Integrating Deep Convolutional Neural Networks and Level Set Methods

The dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) technique has taken on a significant and increasing role in diagnostic procedures and treatments for patients who suffer from chronic kidney disease. Careful segmentation of kidneys from DCE-MRI scans is an essential early step towards the evaluation of kidney function. Recently, deep convolutional neural networks have increased in popularity in medical image segmentation. To this end, in this paper, we propose a new and fully automated two-phase approach that integrates convolutional neural networks and level set methods to delimit kidneys in DCE-MRI scans. We first develop two convolutional neural networks that rely on the U-Net structure (UNT) to predict a kidney probability map for DCE-MRI scans. Then, to leverage the segmentation performance, the pixel-wise kidney probability map predicted from the deep model is exploited with the shape prior information in a level set method to guide the contour evolution towards the target kidney. Real DCE-MRI datasets of 45 subjects are used for training, validating, and testing the proposed approach. The valuation results demonstrate the high performance of the two-phase approach, achieving a Dice similarity coefficient of 0.95 ± 0.02 and intersection over union of 0.91 ± 0.03, and 1.54 ± 1.6 considering a 95% Hausdorff distance. Our intensive experiments confirm the potential and effectiveness of that approach over both UNT models and numerous recent level set-based methods.

Variational Approach for Joint Kidney Segmentation and Registration from DCE-MRI Using Fuzzy Clustering with Shape Priors

The dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) technique has great potential in the diagnosis, therapy, and follow-up of patients with chronic kidney disease (CKD). Towards that end, precise kidney segmentation from DCE-MRI data becomes a prerequisite processing step. Exploiting the useful information about the kidney's shape in this step mandates a registration operation beforehand to relate the shape model coordinates to those of the image to be segmented. Imprecise alignment of the shape model induces errors in the segmentation results. In this paper, we propose a new variational formulation to jointly segment and register DCE-MRI kidney images based on fuzzy c-means clustering embedded within a level-set (LSet) method. The image pixels' fuzzy memberships and the spatial registration parameters are simultaneously updated in each evolution step to direct the LSet contour toward the target kidney. Results on real medical datasets of 45 subjects demonstrate the superior performance of the proposed approach, reporting a Dice similarity coefficient of 0.94 ± 0.03, Intersection-over-Union of 0.89 ± 0.05, and 2.2 ± 2.3 in 95-percentile of Hausdorff distance. Extensive experiments show that our approach outperforms several state-of-the-art LSet-based methods as well as two UNet-based deep neural models trained for the same task in terms of accuracy and consistency.

Level-Set-Based Kidney Segmentation from DCE-MRI Using Fuzzy Clustering with Population-Based and Subject-Specific Shape Statistics

The segmentation of dynamic contrast-enhanced magnetic resonance images (DCE-MRI) of the kidney is a fundamental step in the early and noninvasive detection of acute renal allograft rejection. In this paper, a new and accurate DCE‑MRI kidney segmentation method is proposed. In this method, fuzzy c-means (FCM) clustering is embedded into a level set method, with the fuzzy memberships being iteratively updated during the level set contour evolution. Moreover, population‑based shape (PB-shape) and subject-specific shape (SS-shape) statistics are both exploited. The PB-shape model is trained offline from ground-truth kidney segmentations of various subjects, whereas the SS-shape model is trained on the fly using the segmentation results that are obtained for a specific subject. The proposed method was evaluated on the real medical datasets of 45 subjects and reports a Dice similarity coefficient (DSC) of 0.953 ± 0.018, an intersection-over-union (IoU) of 0.91 ± 0.033, and 1.10 ± 1.4 in the 95-percentile of Hausdorff distance (HD95). Extensive experiments confirm the superiority of the proposed method over several state-of-the-art level set methods, with an average improvement of 0.7 in terms of HD95. It also offers an HD95 improvement of 9.5 and 3.8 over two deep neural networks based on the U-Net architecture. The accuracy improvements have been experimentally found to be more prominent on low-contrast and noisy images.

Kidney segmentation from DCE-MRI converging level set methods, fuzzy clustering and Markov random field modeling

Early diagnosis of transplanted kidney function requires precise Kidney segmentation from Dynamic Contrast-Enhanced Magnetic Resonance Imaging images as a preliminary step. In this regard, this paper aims to propose an automated and accurate DCE-MRI kidney segmentation method integrating fuzzy c-means (FCM) clustering and Markov random field modeling into a level set formulation. The fuzzy memberships, kidney's shape prior model, and spatial interactions modeled using a second-order MRF guide the LS contour evolution towards the target kidney. Several experiments on real medical data of 45 subjects have shown that the proposed method can achieve high and consistent segmentation accuracy regardless of where the LS contour was initialized. It achieves an accuracy of 0.956 ± 0.019 in Dice similarity coefficient (DSC) and 1.15 ± 1.46 in 95% percentile of Hausdorff distance (HD95). Our quantitative comparisons confirm the superiority of the proposed method over several LS methods with an average improvement of more than 0.63 in terms of HD95. It also offers HD95 improvements of 9.62 and 3.94 over two deep neural networks based on the U-Net model. The accuracy improvements are experimentally found to be more profound on low-contrast images as well as DCE-MRI images with high noise levels.

On the potential use of dynamic contrast-enhanced (DCE) MRI parameters as radiomic features of cervical cancer

PURPOSE

To evaluate whether the analysis of high-temporal resolution DCE-MRI by various tracer kinetic models could yield useful radiomic features in discriminating cervix carcinoma and normal cervix tissue.

METHODS

Forty-three patients (median age 51 yr; range 26-78 yr) diagnosed with cervical cancer based on postoperative pathology were enrolled in this study with informed consent. DCE-MRI data with temporal resolution of 2 s were acquired and analyzed using the Tofts (TOFTS), extended Tofts (EXTOFTS), conventional two-compartment (CC), adiabatic tissue homogeneity (ATH), and distributed parameter (DP) models. Ability of all kinetic parameters in distinguishing tumor from normal tissue was assessed using Mann-Whitney U test and receiver operating characteristic (ROC) curves. Repeatability of parameter estimates due to sampling of arterial input functions (AIFs) was also studied using intraclass correlation (ICC) analysis.

RESULTS

Fractional extravascular, extracellular volume (Ve) of all models were significantly smaller in cervix carcinoma than normal cervix tissue, and were associated with large values of area under ROC curve (AUC 0.884-0.961). Capillary permeability PS derived from the ATH, CC, and DP models also yielded large AUC values (0.730, 0.860, and 0.797). Transfer constant Ktrans derived from TOFTS and EXTOFTS models yielded smaller AUC (0.587 and 0.701). Repeatability of parameters derived from all models was robust to AIF sampling, with ICC coefficients typically larger than 0.80.

CONCLUSIONS

With the use of high-temporal resolution DCE-MRI, all tracer kinetic models could reflect pathophysiological differences between cervix carcinoma and normal tissue (with significant differences in Ve and PS) and potentially yield radiomic features with diagnostic value.

A Novel CNN-Based CAD System for Early Assessment of Transplanted Kidney Dysfunction

This paper introduces a deep-learning based computer-aided diagnostic (CAD) system for the early detection of acute renal transplant rejection. For noninvasive detection of kidney rejection at an early stage, the proposed CAD system is based on the fusion of both imaging markers and clinical biomarkers. The former are derived from diffusion-weighted magnetic resonance imaging (DW-MRI) by estimating the apparent diffusion coefficients (ADC) representing the perfusion of the blood and the diffusion of the water inside the transplanted kidney. The clinical biomarkers, namely: creatinine clearance (CrCl) and serum plasma creatinine (SPCr), are integrated into the proposed CAD system as kidney functionality indexes to enhance its diagnostic performance. The ADC maps are estimated for a user-defined region of interest (ROI) that encompasses the whole kidney. The estimated ADCs are fused with the clinical biomarkers and the fused data is then used as an input to train and test a convolutional neural network (CNN) based classifier. The CAD system is tested on DW-MRI scans collected from 56 subjects from geographically diverse populations and different scanner types/image collection protocols. The overall accuracy of the proposed system is 92.9% with 93.3% sensitivity and 92.3% specificity in distinguishing non-rejected kidney transplants from rejected ones. These results demonstrate the potential of the proposed system for a reliable non-invasive diagnosis of renal transplant status for any DW-MRI scans, regardless of the geographical differences and/or imaging protocol.

Computer-Aided Diagnostic System for Early Detection of Acute Renal Transplant Rejection Using Diffusion-Weighted MRI

OBJECTIVE

Early diagnosis of acute renal transplant rejection (ARTR) is critical for accurate treatment. Although the current gold standard, diagnostic technique is renal biopsy, it is not preferred due to its invasiveness, long recovery time (1-2 weeks), and potential for complications, e.g., bleeding and/or infection.

METHODS

This paper presents a computer-aided diagnostic (CAD) system for early ARTR detection using (3D + b-value) diffusion-weighted (DW) magnetic resonance imaging (MRI) data. The CAD process starts from kidney tissue segmentation with an evolving geometric (level-set-based) deformable model. The evolution is guided by a voxel-wise stochastic speed function, which follows from a joint kidney-background Markov-Gibbs random field model accounting for an adaptive kidney shape prior and on-going kidney-background visual appearances. A B-spline-based three-dimensional data alignment is employed to handle local deviations due to breathing and heart beating. Then, empirical cumulative distribution functions of apparent diffusion coefficients of the segmented DW-MRI at different b-values are collected as discriminatory transplant status features. Finally, a deep-learning-based classifier with stacked nonnegative constrained autoencoders is employed to distinguish between rejected and nonrejected renal transplants.

RESULTS

In our initial "leave-one-subject-out" experiment on 100 subjects, [Formula: see text] of the subjects were correctly classified. The subsequent four-fold and ten-fold cross-validations gave the average accuracy of [Formula: see text] and [Formula: see text], respectively.

CONCLUSION

These results demonstrate the promise of this new CAD system to reliably diagnose renal transplant rejection.

SIGNIFICANCE

The technology presented here can significantly impact the quality of care of renal transplant patients since it has the potential to replace the gold standard in kidney diagnosis, biopsy.

3D kidney segmentation from abdominal diffusion MRI using an appearance-guided deformable boundary

A new technique for more accurate automatic segmentation of the kidney from its surrounding abdominal structures in diffusion-weighted magnetic resonance imaging (DW-MRI) is presented. This approach combines a new 3D probabilistic shape model of the kidney with a first-order appearance model and fourth-order spatial model of the diffusion-weighted signal intensity to guide the evolution of a 3D geometric deformable model. The probabilistic shape model was built from labeled training datasets to produce a spatially variant, independent random field of region labels. A Markov-Gibbs random field spatial model with up to fourth-order interactions was adequate to capture the inhomogeneity of renal tissues in the DW-MRI signal. A new analytical approach estimated the Gibbs potentials directly from the DW-MRI data to be segmented, in order that the segmentation procedure would be fully automatic. Finally, to better distinguish the kidney object from the surrounding tissues, marginal gray level distributions inside and outside of the deformable boundary were modeled with adaptive linear combinations of discrete Gaussians (first-order appearance model). The approach was tested on a cohort of 64 DW-MRI datasets with b-values ranging from 50 to 1000 s/mm2. The performance of the presented approach was evaluated using leave-one-subject-out cross validation and compared against three other well-known segmentation methods applied to the same DW-MRI data using the following evaluation metrics: 1) the Dice similarity coefficient (DSC); 2) the 95-percentile modified Hausdorff distance (MHD); and 3) the percentage kidney volume difference (PKVD). High performance of the new approach was confirmed by the high DSC (0.95±0.01), low MHD (3.9±0.76) mm, and low PKVD (9.5±2.2)% relative to manual segmentation by an MR expert (a board certified radiologist).

Automatic blood vessels segmentation based on different retinal maps from OCTA scans

The retinal vascular network reflects the health of the retina, which is a useful diagnostic indicator of systemic vascular. Therefore, the segmentation of retinal blood vessels is a powerful method for diagnosing vascular diseases. This paper presents an automatic segmentation system for retinal blood vessels from Optical Coherence Tomography Angiography (OCTA) images. The system segments blood vessels from the superficial and deep retinal maps for normal and diabetic cases. Initially, we reduced the noise and improved the contrast of the OCTA images by using the Generalized Gauss-Markov random field (GGMRF) model. Secondly, we proposed a joint Markov-Gibbs random field (MGRF) model to segment the retinal blood vessels from other background tissues. It integrates both appearance and spatial models in addition to the prior probability model of OCTA images. The higher order MGRF (HO-MGRF) model in addition to the 1^st-order intensity model are used to consider the spatial information in order to overcome the low contrast between vessels and other tissues. Finally, we refined the segmentation by extracting connected regions using a 2D connectivity filter. The proposed segmentation system was trained and tested on 47 data sets, which are 23 normal data sets and 24 data sets for diabetic patients. To evaluate the accuracy and robustness of the proposed segmentation framework, we used three different metrics, which are Dice similarity coefficient (DSC), absolute vessels volume difference (VVD), and area under the curve (AUC). The results on OCTA data sets (DSC=95.04±3.75%, VVD=8.51±1.49%, and AUC=95.20±1.52%) show the promise of the proposed segmentation approach.

3D Kidney Segmentation from Abdominal Images Using Spatial-Appearance Models

Kidney segmentation is an essential step in developing any noninvasive computer-assisted diagnostic system for renal function assessment. This paper introduces an automated framework for 3D kidney segmentation from dynamic computed tomography (CT) images that integrates discriminative features from the current and prior CT appearances into a random forest classification approach. To account for CT images' inhomogeneities, we employ discriminate features that are extracted from a higher-order spatial model and an adaptive shape model in addition to the first-order CT appearance. To model the interactions between CT data voxels, we employed a higher-order spatial model, which adds the triple and quad clique families to the traditional pairwise clique family. The kidney shape prior model is built using a set of training CT data and is updated during segmentation using not only region labels but also voxels' appearances in neighboring spatial voxel locations. Our framework performance has been evaluated on in vivo dynamic CT data collected from 20 subjects and comprises multiple 3D scans acquired before and after contrast medium administration. Quantitative evaluation between manually and automatically segmented kidney contours using Dice similarity, percentage volume differences, and 95th-percentile bidirectional Hausdorff distances confirms the high accuracy of our approach.

Non-invasive approaches in the diagnosis of acute rejection in kidney transplant recipients, part II: omics analyses of urine and blood samples

Kidney transplantation (KTx) represents the best available treatment for patients with end-stage renal disease. Still, the full benefits of KTx are undermined by acute rejection (AR). The diagnosis of AR ultimately relies on transplant needle biopsy. However, such an invasive procedure is associated with a significant risk of complications and is limited by sampling error and interobserver variability. In the present review, we summarize the current literature about non-invasive approaches for the diagnosis of AR in kidney transplant recipients (KTRs), including in vivo imaging, gene-expression profiling and omics analyses of blood and urine samples. Most imaging techniques, such as contrast-enhanced ultrasound and magnetic resonance, exploit the fact that blood flow is significantly lowered in case of AR-induced inflammation. In addition, AR-associated recruitment of activated leucocytes may be detectable by 18F-fluorodeoxyglucose positron emission tomography. In parallel, urine biomarkers, including CXCL9/CXCL10 or a three-gene signature of CD3ε, CXCL10 and 18S RNA levels, have been identified. None of these approaches has yet been adopted in the clinical follow-up of KTRs, but standardization of analysis procedures may help assess reproducibility and comparative diagnostic yield in large, prospective, multicentre trials.

Non-invasive approaches in the diagnosis of acute rejection in kidney transplant recipients. Part I. In vivo imaging methods

Kidney transplantation (KTx) represents the best available treatment for patients with end-stage renal disease. Still, full benefits of KTx are undermined by acute rejection (AR). The diagnosis of AR ultimately relies on transplant needle biopsy. However, such an invasive procedure is associated with a significant risk of complications and is limited by sampling error and interobserver variability. In the present review, we summarize the current literature about non-invasive approaches for the diagnosis of AR in kidney transplant recipients (KTRs), including in vivo imaging, gene expression profiling and omics analyses of blood and urine samples. Most imaging techniques, like contrast-enhanced ultrasound and magnetic resonance, exploit the fact that blood flow is significantly lowered in case of AR-induced inflammation. In addition, AR-associated recruitment of activated leukocytes may be detectable by ¹⁸F-fluoro-deoxy-glucose positron emission tomography. In parallel, urine biomarkers, including CXCL9/CXCL10 or a three-gene signature of CD3ε, IP-10 and 18S RNA levels, have been identified. None of these approaches has been adopted yet in the clinical follow-up of KTRs, but standardization of procedures may help assess reproducibility and compare diagnostic yields in large prospective multicentric trials.

Models and methods for analyzing DCE-MRI: a review

PURPOSE

To present a review of most commonly used techniques to analyze dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), discusses their strengths and weaknesses, and outlines recent clinical applications of findings from these approaches.

METHODS

DCE-MRI allows for noninvasive quantitative analysis of contrast agent (CA) transient in soft tissues. Thus, it is an important and well-established tool to reveal microvasculature and perfusion in various clinical applications. In the last three decades, a host of nonparametric and parametric models and methods have been developed in order to quantify the CA's perfusion into tissue and estimate perfusion-related parameters (indexes) from signal- or concentration-time curves. These indexes are widely used in various clinical applications for the detection, characterization, and therapy monitoring of different diseases.

RESULTS

Promising theoretical findings and experimental results for the reviewed models and techniques in a variety of clinical applications suggest that DCE-MRI is a clinically relevant imaging modality, which can be used for early diagnosis of different diseases, such as breast and prostate cancer, renal rejection, and liver tumors.

CONCLUSIONS

Both nonparametric and parametric approaches for DCE-MRI analysis possess the ability to quantify tissue perfusion.

Image registration for triggered and non-triggered DTI of the human kidney: reduced variability of diffusion parameter estimation

BACKGROUND

To investigate if non-rigid image-registration reduces motion artifacts in triggered and non-triggered diffusion tensor imaging (DTI) of native kidneys. A secondary aim was to determine, if improvements through registration allow for omitting respiratory-triggering.

METHODS

Twenty volunteers underwent coronal DTI of the kidneys with nine b-values (10-700 s/mm(2) ) at 3 Tesla. Image-registration was performed using a multimodal nonrigid registration algorithm. Data processing yielded the apparent diffusion coefficient (ADC), the contribution of perfusion (FP ), and the fractional anisotropy (FA). For comparison of the data stability, the root mean square error (RMSE) of the fitting and the standard deviations within the regions of interest (SDROI ) were evaluated.

RESULTS

RMSEs decreased significantly after registration for triggered and also for non-triggered scans (P < 0.05). SDROI for ADC, FA, and FP were significantly lower after registration in both medulla and cortex of triggered scans (P < 0.01). Similarly the SDROI of FA and FP decreased significantly in non-triggered scans after registration (P < 0.05). RMSEs were significantly lower in triggered than in non-triggered scans, both with and without registration (P < 0.05).

CONCLUSION

Respiratory motion correction by registration of individual echo-planar images leads to clearly reduced signal variations in renal DTI for both triggered and particularly non-triggered scans. Secondarily, the results suggest that respiratory-triggering still seems advantageous.

Individually wide range of renal motion evaluated by four-dimensional computed tomography

Objectives

Assessment of physiologic renal motion in order to optimize abdominal intensity-modulated radiation therapy and stereotactic body radiation therapy.

Methods and materials

Twenty patients with a median age of 47 years underwent computed tomography simulation and four-dimensional computed tomography acquisition. Thirty-nine kidneys were contoured during ten phases of respiration to estimate renal motion.

Results

Kidney motion was not related to age (p = 0.42), sex (p = 0.28), height (p = 0.75), or body weight (p = 0.63). The average +/- standard deviation (SD) of movement of the center of gravity for all subjects was 11.1 +/- 4.8 mm in the cranio-caudal (CC) direction (range, 2.5-20.5 mm), 3.6 +/- 2.1 mm in the anterior-posterior (AP) direction (range, 0.6-8.0 mm), and 1.7 +/- 1.4 mm in the right-left (RL) direction (range, 0.4-5.9 mm). Renal motion strongly correlated with the respiratory phases (r > 0.97 and p < 0.01 in all three directions).

Conclusions

Renal motion was independent of age, sex, height, or body weight. Renal motion in all directions was strongly respiration dependent, but motion in the cranio-caudal direction showed wide individual variation. In a clinical setting, it will be necessary to evaluate renal respiratory motion separately in each individual.