Imaging the dynamics of murine uterine contractions in early pregnancy†.

Uterine muscle contractility is essential for reproductive processes including sperm and embryo transport, and during the uterine cycle to remove menstrual effluent. Even still, uterine contractions have primarily been studied in the context of preterm labor. This is partly due to a lack of methods for studying the uterine muscle contractility in the intact organ. Here, we describe an imaging-based method to evaluate mouse uterine contractility of both the longitudinal and circular muscles in the cycling stages and in early pregnancy. By transforming the image-based data into three-dimensional spatiotemporal contractility maps, we calculate waveform characteristics of muscle contractions, including amplitude, frequency, wavelength, and velocity. We report that the native organ is highly contractile during the progesterone-dominant diestrus stage of the cycle when compared to the estrogen-dominant proestrus and estrus stages. We also observed that during the first phase of uterine embryo movement when clustered embryos move toward the middle of the uterine horn, contractions are dynamic and non-uniform between different segments of the uterine horn. In the second phase of embryo movement, contractions are more uniform and rhythmic throughout the uterine horn. Finally, in Lpar3-/- uteri, which display faster embryo movement, we observe global and regional increases in contractility. Our method provides a means to understand the wave characteristics of uterine smooth muscle in response to modulators and in genetic mutants. Better understanding uterine contractility in the early pregnancy stages is critical for the advancement of artificial reproductive technologies and a possibility of modulating embryo movement during clinical embryo transfers.

Clinical-Genomic Risk Group Classification of Suspicious Lesions on Prostate Multiparametric-MRI.

The utilization of multi-parametric MRI (mpMRI) in clinical decisions regarding prostate cancer patients' management has recently increased. After biopsy, clinicians can assess risk using National Comprehensive Cancer Network (NCCN) risk stratification schema and commercially available genomic classifiers, such as Decipher. We built radiomics-based models to predict lesions/patients at low risk prior to biopsy based on an established three-tier clinical-genomic classification system. Radiomic features were extracted from regions of positive biopsies and Normally Appearing Tissues (NAT) on T2-weighted and Diffusion-weighted Imaging. Using only clinical information available prior to biopsy, five models for predicting low-risk lesions/patients were evaluated, based on: 1: Clinical variables; 2: Lesion-based radiomic features; 3: Lesion and NAT radiomics; 4: Clinical and lesion-based radiomics; and 5: Clinical, lesion and NAT radiomic features. Eighty-three mpMRI exams from 78 men were analyzed. Models 1 and 2 performed similarly (Area under the receiver operating characteristic curve were 0.835 and 0.838, respectively), but radiomics significantly improved the lesion-based performance of the model in a subset analysis of patients with a negative Digital Rectal Exam (DRE). Adding normal tissue radiomics significantly improved the performance in all cases. Similar patterns were observed on patient-level models. To the best of our knowledge, this is the first study to demonstrate that machine learning radiomics-based models can predict patients' risk using combined clinical-genomic classification.

Deep Learning for Per-Fraction Automatic Segmentation of Gross Tumor Volume (GTV) and Organs at Risk (OARs) in Adaptive Radiotherapy of Cervical Cancer.

Background/Hypothesis: MRI-guided online adaptive radiotherapy (MRI-g-OART) improves target coverage and organs-at-risk (OARs) sparing in radiation therapy (RT). For patients with locally advanced cervical cancer (LACC) undergoing RT, changes in bladder and rectal filling contribute to large inter-fraction target volume motion. We hypothesized that deep learning (DL) convolutional neural networks (CNN) can be trained to accurately segment gross tumor volume (GTV) and OARs both in planning and daily fractions' MRI scans. Materials/Methods: We utilized planning and daily treatment fraction setup (RT-Fr) MRIs from LACC patients, treated with stereotactic body RT to a dose of 45-54 Gy in 25 fractions. Nine structures were manually contoured. MASK R-CNN network was trained and tested under three scenarios: (i) Leave-one-out (LOO), using the planning images of N- 1 patients for training; (ii) the same network, tested on the RT-Fr MRIs of the "left-out" patient, (iii) including the planning MRI of the "left-out" patient as an additional training sample, and tested on RT-Fr MRIs. The network performance was evaluated using the Dice Similarity Coefficient (DSC) and Hausdorff distances. The association between the structures' volume and corresponding DSCs was investigated using Pearson's Correlation Coefficient, r. Results: MRIs from fifteen LACC patients were analyzed. In the LOO scenario the DSC for Rectum, Femur, and Bladder was >0.8, followed by the GTV, Uterus, Mesorectum and Parametrium (0.6-0.7). The results for Vagina and Sigmoid were suboptimal. The performance of the network was similar for most organs when tested on RT-Fr MRI. Including the planning MRI in the training did not improve the segmentation of the RT-Fr MRI. There was a significant correlation between the average organ volume and the corresponding DSC (r = 0.759, p = 0.018). Conclusion: We have established a robust workflow for training MASK R-CNN to automatically segment GTV and OARs in MRI-g-OART of LACC. Albeit the small number of patients in this pilot project, the network was trained to successfully identify several structures while challenges remain, especially in relatively small organs. With the increase of the LACC cases, the performance of the network will improve. A robust auto-contouring tool would improve workflow efficiency and patient tolerance of the OART process.

A Quality Control System for Automated Prostate Segmentation on T2-Weighted MRI.

Computer-aided detection and diagnosis (CAD) systems have the potential to improve robustness and efficiency compared to traditional radiological reading of magnetic resonance imaging (MRI). Fully automated segmentation of the prostate is a crucial step of CAD for prostate cancer, but visual inspection is still required to detect poorly segmented cases. The aim of this work was therefore to establish a fully automated quality control (QC) system for prostate segmentation based on T2-weighted MRI. Four different deep learning-based segmentation methods were used to segment the prostate for 585 patients. First order, shape and textural radiomics features were extracted from the segmented prostate masks. A reference quality score (QS) was calculated for each automated segmentation in comparison to a manual segmentation. A least absolute shrinkage and selection operator (LASSO) was trained and optimized on a randomly assigned training dataset (N = 1756, 439 cases from each segmentation method) to build a generalizable linear regression model based on the radiomics features that best estimated the reference QS. Subsequently, the model was used to estimate the QSs for an independent testing dataset (N = 584, 146 cases from each segmentation method). The mean ± standard deviation absolute error between the estimated and reference QSs was 5.47 ± 6.33 on a scale from 0 to 100. In addition, we found a strong correlation between the estimated and reference QSs (rho = 0.70). In conclusion, we developed an automated QC system that may be helpful for evaluating the quality of automated prostate segmentations.

Segmentation of prostate and prostate zones using deep learning : A multi-MRI vendor analysis.

PURPOSE: Develop a deep-learning-based segmentation algorithm for prostate and its peripheral zone (PZ) that is reliable across multiple MRI vendors. METHODS: This is a retrospective study. The dataset consisted of 550 MRIs (Siemens-330, General Electric[GE]-220). A multistream 3D convolutional neural network is used for automatic segmentation of the prostate and its PZ using T2-weighted (T2-w) MRI. Prostate and PZ were manually contoured on axial T2­w. The network uses axial, coronal, and sagittal T2­w series as input. The preprocessing of the input data includes bias correction, resampling, and image normalization. A dataset from two MRI vendors (Siemens and GE) is used to test the proposed network. Six different models were trained, three for the prostate and three for the PZ. Of the three, two were trained on data from each vendor separately, and a third (Combined) on the aggregate of the datasets. The Dice coefficient (DSC) is used to compare the manual and predicted segmentation. RESULTS: For prostate segmentation, the Combined model obtained DSCs of 0.893 ± 0.036 and 0.825 ± 0.112 (mean ± standard deviation) on Siemens and GE, respectively. For PZ, the best DSCs were from the Combined model: 0.811 ± 0.079 and 0.788 ± 0.093. While the Siemens model underperformed on the GE dataset and vice versa, the Combined model achieved robust performance on both datasets. CONCLUSION: The proposed network has a performance comparable to the interexpert variability for segmenting the prostate and its PZ. Combining images from different MRI vendors on the training of the network is of paramount importance for building a universal model for prostate and PZ segmentation.

Classification of suspicious lesions on prostate multiparametric MRI using machine learning.

We present a radiomics-based approach developed for the SPIE-AAPM-NCI PROSTATEx challenge. The task was to classify clinically significant prostate cancer in multiparametric (mp) MRI. Data consisted of a "training dataset" (330 suspected lesions from 204 patients) and a "test dataset" (208 lesions/140 patients). All studies included T2-weighted (T2-W), proton density-weighted, dynamic contrast enhanced, and diffusion-weighted imaging. Analysis of the images was performed using the MIM imaging platform (MIM Software, Cleveland, Ohio). Prostate and peripheral zone contours were manually outlined on the T2-W images. A workflow for rigid fusion of the aforementioned images to T2-W was created in MIM. The suspicious lesion was outlined using the high b-value image. Intensity and texture features were extracted on four imaging modalities and characterized using nine histogram descriptors: 10%, 25%, 50%, 75%, 90%, mean, standard deviation, kurtosis, and skewness (216 features). Three classification methods were used: classification and regression trees (CART), random forests, and adaptive least absolute shrinkage and selection operator (LASSO). In the held out by the organizers test dataset, the areas under the curve (AUCs) were: 0.82 (random forests), 0.76 (CART), and 0.76 (adaptive LASSO). AUC of 0.82 was the fourth-highest score of 71 entries (32 teams) and the highest for feature-based methods.

A new method for the analysis of soft tissues with data acquired under field conditions.

Analyzing soft-tissue structures is particularly challenging due to the lack of homologous landmarks that can be reliably identified across time and specimens. This is particularly true when data are to be collected under field conditions. Here, we present a method that combines photogrammetric techniques and geometric morphometrics methods (GMM) to quantify soft tissues for their subsequent volumetric analysis. We combine previously developed methods for landmark data acquisition and processing with a custom program for volumetric computations. Photogrammetric methods are a particularly powerful tool for field studies as they allow for image acquisition with minimal equipment requirements and for the acquisition of the spatial coordinates of points (anatomical landmarks or others) from these images. For our method, a limited number of homologous landmarks, i.e., points that can be found on any specimen independent of space and time, and further distinctive points, which may vary over time, space and subject, are identified on two-dimensional photographs and their three-dimensional coordinates estimated using photogrammetric methods. The three-dimensional configurations are oriented by the spatial principal components (PCs) of the homologous points. Crucially, this last step orients the configuration such that x and y-information (PC1 and PC2 coordinates) constitute an anatomically-defined plane with the z-values (PC3 coordinate) in the direction of interest for volume computation. The z-coordinates are then used to estimate the volume of the tissue. We validate our method using a physical, geometric model of known dimensions and physical (wax) models designed to approximate perineal swellings in female macaques. To demonstrate the usefulness and potential of our method, we use it to estimate the volumes of Barbary macaque sexual swellings recorded in the field with video images. By analyzing both the artificial data and real monkey swellings, we validate our method's accuracy and illustrate its potential for application in important areas of biological research.