Towards accelerating model parallelism in distributed deep learning systems

Modern deep neural networks cannot be often trained on a single GPU due to large model size and large data size. Model parallelism splits a model for multiple GPUs, but making it scalable and seamless is challenging due to different information sharing among GPUs with communication overhead. Specifically, we identify two key issues to make the parallelism being inefficient and inaccurate; an efficient pipelining technique is crucial to maximize GPU utilization and normalizations in deep neural networks may affect the performance due to different statistics sharing of mini-batch. In this work, we address these issues by investigating efficient pipelining for model parallelism and effective normalizations in model / data parallelisms when training a model with large mini-batch in multiple GPUs so that the model performance in accuracy can not be compromised. Firstly, we propose a novel method to search for an optimal micro-batch size considering the number of GPUs and memory size for model parallelism. For efficient pipelining, mini-batch is usually divided into smaller batches (called micro-batch). To maximize the utilization of GPU computing resources, training should be performed with the optimal micro-batch size. Our proposed micro-batch size search algorithm achieved increased image throughput by up to 12% and improved trainable mini-batch size by 25% as compared to the conventional model parallelism method. Secondly, we investigate normalizations in distributed deep learning training for different parallelisms. Our experiments using different normalization methods suggested that the performance with batch normalization can be improved by sharing the batch information among GPUs when performing data parallelism. It was also confirmed that group normalization helped minimizing accuracy degradation when performing model parallelism with pipelining and yielded consistent accuracies for diverse mini-batch sizes.

Physics-Guided Deep Scatter Estimation by Weak Supervision for Quantitative SPECT

Accurate scatter estimation is important in quantitative SPECT for improving image contrast and accuracy. With a large number of photon histories, Monte-Carlo (MC) simulation can yield accurate scatter estimation, but is computationally expensive. Recent deep learning-based approaches can yield accurate scatter estimates quickly, yet full MC simulation is still required to generate scatter estimates as ground truth labels for all training data. Here we propose a physics-guided weakly supervised training framework for fast and accurate scatter estimation in quantitative SPECT by using a 100× shorter MC simulation as weak labels and enhancing them with deep neural networks. Our weakly supervised approach also allows quick fine-tuning of the trained network to any new test data for further improved performance with an additional short MC simulation (weak label) for patient-specific scatter modelling. Our method was trained with 18 XCAT phantoms with diverse anatomies / activities and then was evaluated on 6 XCAT phantoms, 4 realistic virtual patient phantoms, 1 torso phantom and 3 clinical scans from 2 patients for 177Lu SPECT with single / dual photopeaks (113, 208 keV). Our proposed weakly supervised method yielded comparable performance to the supervised counterpart in phantom experiments, but with significantly reduced computation in labeling. Our proposed method with patient-specific fine-tuning achieved more accurate scatter estimates than the supervised method in clinical scans. Our method with physics-guided weak supervision enables accurate deep scatter estimation in quantitative SPECT, while requiring much lower computation in labeling, enabling patient-specific fine-tuning capability in testing.

PAIP 2020: Microsatellite instability prediction in colorectal cancer

Microsatellite instability (MSI) refers to alterations in the length of simple repetitive genomic sequences. MSI status serves as a prognostic and predictive factor in colorectal cancer. The MSI-high status is a good prognostic factor in stage II/III cancer, and predicts a lack of benefit to adjuvant fluorouracil chemotherapy in stage II cancer but a good response to immunotherapy in stage IV cancer. Therefore, determining MSI status in patients with colorectal cancer is important for identifying the appropriate treatment protocol. In the Pathology Artificial Intelligence Platform (PAIP) 2020 challenge, artificial intelligence researchers were invited to predict MSI status based on colorectal cancer slide images. Participants were required to perform two tasks. The primary task was to classify a given slide image as belonging to either the MSI-high or the microsatellite-stable group. The second task was tumor area segmentation to avoid ties with the main task. A total of 210 of the 495 participants enrolled in the challenge downloaded the images, and 23 teams submitted their final results. Seven teams from the top 10 participants agreed to disclose their algorithms, most of which were convolutional neural network-based deep learning models, such as EfficientNet and UNet. The top-ranked system achieved the highest F1 score (0.9231). This paper summarizes the various methods used in the PAIP 2020 challenge. This paper supports the effectiveness of digital pathology for identifying the relationship between colorectal cancer and the MSI characteristics.

Development and operation of a digital platform for sharing pathology image data

BACKGROUND

Artificial intelligence (AI) research is highly dependent on the nature of the data available. With the steady increase of AI applications in the medical field, the demand for quality medical data is increasing significantly. We here describe the development of a platform for providing and sharing digital pathology data to AI researchers, and highlight challenges to overcome in operating a sustainable platform in conjunction with pathologists.

METHODS

Over 3000 pathological slides from five organs (liver, colon, prostate, pancreas and biliary tract, and kidney) in histologically confirmed tumor cases by pathology departments at three hospitals were selected for the dataset. After digitalizing the slides, tumor areas were annotated and overlaid onto the images by pathologists as the ground truth for AI training. To reduce the pathologists' workload, AI-assisted annotation was established in collaboration with university AI teams.

RESULTS

A web-based data sharing platform was developed to share massive pathological image data in 2019. This platform includes 3100 images, and 5 pre-processing algorithms for AI researchers to easily load images into their learning models.

DISCUSSION

Due to different regulations among countries for privacy protection, when releasing internationally shared learning platforms, it is considered to be most prudent to obtain consent from patients during data acquisition.

CONCLUSIONS

Despite limitations encountered during platform development and model training, the present medical image sharing platform can steadily fulfill the high demand of AI developers for quality data. This study is expected to help other researchers intending to generate similar platforms that are more effective and accessible in the future.

PAIP 2019: Liver cancer segmentation challenge

Pathology Artificial Intelligence Platform (PAIP) is a free research platform in support of pathological artificial intelligence (AI). The main goal of the platform is to construct a high-quality pathology learning data set that will allow greater accessibility. The PAIP Liver Cancer Segmentation Challenge, organized in conjunction with the Medical Image Computing and Computer Assisted Intervention Society (MICCAI 2019), is the first image analysis challenge to apply PAIP datasets. The goal of the challenge was to evaluate new and existing algorithms for automated detection of liver cancer in whole-slide images (WSIs). Additionally, the PAIP of this year attempted to address potential future problems of AI applicability in clinical settings. In the challenge, participants were asked to use analytical data and statistical metrics to evaluate the performance of automated algorithms in two different tasks. The participants were given the two different tasks: Task 1 involved investigating Liver Cancer Segmentation and Task 2 involved investigating Viable Tumor Burden Estimation. There was a strong correlation between high performance of teams on both tasks, in which teams that performed well on Task 1 also performed well on Task 2. After evaluation, we summarized the top 11 team's algorithms. We then gave pathological implications on the easily predicted images for cancer segmentation and the challenging images for viable tumor burden estimation. Out of the 231 participants of the PAIP challenge datasets, a total of 64 were submitted from 28 team participants. The submitted algorithms predicted the automatic segmentation on the liver cancer with WSIs to an accuracy of a score estimation of 0.78. The PAIP challenge was created in an effort to combat the lack of research that has been done to address Liver cancer using digital pathology. It remains unclear of how the applicability of AI algorithms created during the challenge can affect clinical diagnoses. However, the results of this dataset and evaluation metric provided has the potential to aid the development and benchmarking of cancer diagnosis and segmentation.

Algorithms and Analyses for Joint Spectral Image Reconstruction in Y-90 Bremsstrahlung SPECT

Quantitative yttrium-90 (Y-90) SPECT imaging is challenging due to the nature of Y-90, an almost pure beta emitter that is associated with a continuous spectrum of bremsstrahlung photons that have a relatively low yield. This paper proposes joint spectral reconstruction (JSR), a novel bremsstrahlung SPECT reconstruction method that uses multiple narrow acquisition windows with accurate multi-band forward modeling to cover a wide range of the energy spectrum. Theoretical analyses using Fisher information and Monte-Carlo (MC) simulation with a digital phantom show that the proposed JSR model with multiple acquisition windows has better performance in terms of covariance (precision) than previous methods using multi-band forward modeling with a single acquisition window, or using a single-band forward modeling with a single acquisition window. We also propose an energy-window subset (ES) algorithm for JSR to achieve fast empirical convergence and maximum-likelihood based initialization for all reconstruction methods to improve quantification accuracy in early iterations. For both MC simulation with a digital phantom and experimental study with a physical multi-sphere phantom, our proposed JSR-ES, a fast algorithm for JSR with ES, yielded higher recovery coefficients (RCs) on hot spheres over all iterations and sphere sizes than all the other evaluated methods, due to fast empirical convergence. In experimental study, for the smallest hot sphere (diameter 1.6cm), at the 20th iteration the increase in RCs with JSR-ES was 66 and 31% compared with single wide and narrow band forward models, respectively. JSR-ES also yielded lower residual count error (RCE) on a cold sphere over all iterations than other methods for MC simulation with known scatter, but led to greater RCE compared with single narrow band forward model at higher iterations for experimental study when using estimated scatter.

A Wearable Wrist Band-Type System for Multimodal Biometrics Integrated with Multispectral Skin Photomatrix and Electrocardiogram Sensors

Multimodal biometrics are promising for providing a strong security level for personal authentication, yet the implementation of a multimodal biometric system for practical usage need to meet such criteria that multimodal biometric signals should be easy to acquire but not easily compromised. We developed a wearable wrist band integrated with multispectral skin photomatrix (MSP) and electrocardiogram (ECG) sensors to improve the issues of collectability, performance and circumvention of multimodal biometric authentication. The band was designed to ensure collectability by sensing both MSP and ECG easily and to achieve high authentication performance with low computation, efficient memory usage, and relatively fast response. Acquisition of MSP and ECG using contact-based sensors could also prevent remote access to personal data. Personal authentication with multimodal biometrics using the integrated wearable wrist band was evaluated in 150 subjects and resulted in 0.2% equal error rate ( EER ) and 100% detection probability at 1% FAR (false acceptance rate) ( PD . 1 ), which is comparable to other state-of-the-art multimodal biometrics. An additional investigation with a separate MSP sensor, which enhanced contact with the skin, along with ECG reached 0.1% EER and 100% PD . 1 , showing a great potential of our in-house wearable band for practical applications. The results of this study demonstrate that our newly developed wearable wrist band may provide a reliable and easy-to-use multimodal biometric solution for personal authentication.

Temporary work and depressive symptoms in South Korean workers

Background

In many countries, including South Korea, labour market changes have led to an increase in unstable, temporary jobs. There is evidence that workers in such jobs may experience poorer mental health than those in more stable employment.

Aims

To investigate the association between temporary employment and depressive symptoms in South Korean workers.

Methods

We analysed data from the 2010-2014 Korean Welfare Panel Study (KOWEPS). Employment type was categorized into workers paid per day of labour (day labourers), those on short-term contracts (fixed-term workers) and permanent workers. The association between employment type and depressive symptoms, measured using the Center for Epidemiological Studies Depression scale (CES-D 11), was examined using the generalized estimating equation model.

Results

A total of 3756 workers aged 20-59 were included in the 2010 baseline population. Day labourers had the highest mean CES-D 11 score, followed by fixed-term workers and permanent workers. With the day labourer group as reference, fixed-term workers (β: -1.5027, P < 0.001) and permanent workers (β: -2.1848, P < 0.001) showed statistically significant decreases in depression scores.

Conclusions

Compared with day labourers, fixed-term workers and permanent workers had progressively lower depression scores. The findings of this study suggest that mental health inequalities based on employment type exist in South Korea.

Improved quantitative 90 Y bremsstrahlung SPECT/CT reconstruction with Monte Carlo scatter modeling

PURPOSE

In ⁹⁰ Y microsphere radioembolization (RE), accurate post-therapy imaging-based dosimetry is important for establishing absorbed dose versus outcome relationships for developing future treatment planning strategies. Additionally, accurately assessing microsphere distributions is important because of concerns for unexpected activity deposition outside the liver. Quantitative ⁹⁰ Y imaging by either SPECT or PET is challenging. In ⁹⁰ Y SPECT model based methods are necessary for scatter correction because energy window-based methods are not feasible with the continuous bremsstrahlung energy spectrum. The objective of this work was to implement and evaluate a scatter estimation method for accurate ⁹⁰ Y bremsstrahlung SPECT/CT imaging.

METHODS

Since a fully Monte Carlo (MC) approach to ⁹⁰ Y SPECT reconstruction is computationally very demanding, in the present study the scatter estimate generated by a MC simulator was combined with an analytical projector in the 3D OS-EM reconstruction model. A single window (105 to 195-keV) was used for both the acquisition and the projector modeling. A liver/lung torso phantom with intrahepatic lesions and low-uptake extrahepatic objects was imaged to evaluate SPECT/CT reconstruction without and with scatter correction. Clinical application was demonstrated by applying the reconstruction approach to five patients treated with RE to determine lesion and normal liver activity concentrations using a (liver) relative calibration.

RESULTS

There was convergence of the scatter estimate after just two updates, greatly reducing computational requirements. In the phantom study, compared with reconstruction without scatter correction, with MC scatter modeling there was substantial improvement in activity recovery in intrahepatic lesions (from > 55% to > 86%), normal liver (from 113% to 104%), and lungs (from 227% to 104%) with only a small degradation in noise (13% vs. 17%). Similarly, with scatter modeling contrast improved substantially both visually and in terms of a detectability index, which was especially relevant for the low uptake extrahepatic objects. The trends observed for the phantom were also seen in the patient studies where lesion activity concentrations and lesion-to-liver concentration ratios were lower for SPECT without scatter correction compared with reconstruction with just two MC scatter updates: in eleven lesions the mean uptake was 4.9 vs. 7.1 MBq/mL (P = 0.0547), the mean normal liver uptake was 1.6 vs. 1.5 MBq/mL (P = 0.056) and the mean lesion-to-liver uptake ratio was 2.7 vs. 4.3 (P = 0.0402) for reconstruction without and with scatter correction respectively.

CONCLUSIONS

Quantitative accuracy of ⁹⁰ Y bremsstrahlung imaging can be substantially improved with MC scatter modeling without significant degradation in image noise or intensive computational requirements.

Bounded Self-Weights Estimation Method for Non-Local Means Image Denoising Using Minimax Estimators

A non-local means (NLM) filter is a weighted average of a large number of non-local pixels with various image intensity values. The NLM filters have been shown to have powerful denoising performance, excellent detail preservation by averaging many noisy pixels, and using appropriate values for the weights, respectively. The NLM weights between two different pixels are determined based on the similarities between two patches that surround these pixels and a smoothing parameter. Another important factor that influences the denoising performance is the self-weight values for the same pixel. The recently introduced local James-Stein type center pixel weight estimation method (LJS) outperforms other existing methods when determining the contribution of the center pixels in the NLM filter. However, the LJS method may result in excessively large self-weight estimates since no upper bound is assumed, and the method uses a relatively large local area for estimating the self-weights, which may lead to a strong bias. In this paper, we investigated these issues in the LJS method, and then propose a novel local self-weight estimation methods using direct bounds (LMM-DB) and reparametrization (LMM-RP) based on the Baranchik's minimax estimator. Both the LMM-DB and LMM-RP methods were evaluated using a wide range of natural images and a clinical MRI image together with the various levels of additive Gaussian noise. Our proposed parameter selection methods yielded an improved bias-variance trade-off, a higher peak signal-to-noise (PSNR) ratio, and fewer visual artifacts when compared with the results of the classical NLM and LJS methods. Our proposed methods also provide a heuristic way to select a suitable global smoothing parameters that can yield PSNR values that are close to the optimal values.

The Neuromelanin-related T2* Contrast in Postmortem Human Substantia Nigra with 7T MRI

High field magnetic resonance imaging (MRI)-based delineation of the substantia nigra (SN) and visualization of its inner cellular organization are promising methods for the evaluation of morphological changes associated with neurodegenerative diseases; however, corresponding MR contrasts must be matched and validated with quantitative histological information. Slices from two postmortem SN samples were imaged with a 7 Tesla (7T) MRI with T₁ and T₂* imaging protocols and then stained with Perl’s Prussian blue, Kluver-Barrera, tyrosine hydroxylase, and calbindin immunohistochemistry in a serial manner. The association between T₂* values and quantitative histology was investigated with a co-registration method that accounts for histology slice preparation. The ventral T₂* hypointense layers between the SNr and the crus cerebri extended anteriorly to the posterior part of the crus cerebri, which demonstrates the difficulty with an MRI-based delineation of the SN. We found that the paramagnetic hypointense areas within the dorsolateral SN corresponded to clusters of neuromelanin (NM). These NM-rich zones were distinct from the hypointense ventromedial regions with high iron pigments. Nigral T₂* imaging at 7T can reflect the density of NM-containing neurons as the metal-bound NM macromolecules may decrease T₂* values and cause hypointense signalling in T₂* imaging at 7T.

The Use of Anatomical Information for Molecular Image Reconstruction Algorithms: Attenuation/Scatter Correction, Motion Compensation, and Noise Reduction

PET and SPECT are important tools for providing valuable molecular information about patients to clinicians. Advances in nuclear medicine hardware technologies and statistical image reconstruction algorithms enabled significantly improved image quality. Sequentially or simultaneously acquired anatomical images such as CT and MRI from hybrid scanners are also important ingredients for improving the image quality of PET or SPECT further. High-quality anatomical information has been used and investigated for attenuation and scatter corrections, motion compensation, and noise reduction via post-reconstruction filtering and regularization in inverse problems. In this article, we will review works using anatomical information for molecular image reconstruction algorithms for better image quality by describing mathematical models, discussing sources of anatomical information for different cases, and showing some examples.

Evaluation of dual energy quantitative CT for determining the spatial distributions of red marrow and bone for dosimetry in internal emitter radiation therapy

PURPOSE

To evaluate a three-equation three-unknown dual-energy quantitative CT (DEQCT) technique for determining region specific variations in bone spongiosa composition for improved red marrow dose estimation in radionuclide therapy.

METHODS

The DEQCT method was applied to 80/140 kVp images of patient-simulating lumbar sectional body phantoms of three sizes (small, medium, and large). External calibration rods of bone, red marrow, and fat-simulating materials were placed beneath the body phantoms. Similar internal calibration inserts were placed at vertebral locations within the body phantoms. Six test inserts of known volume fractions of bone, fat, and red marrow were also scanned. External-to-internal calibration correction factors were derived. The effects of body phantom size, radiation dose, spongiosa region segmentation granularity [single (∼17 × 17 mm) region of interest (ROI), 2 × 2, and 3 × 3 segmentation of that single ROI], and calibration method on the accuracy of the calculated volume fractions of red marrow (cellularity) and trabecular bone were evaluated.

RESULTS

For standard low dose DEQCT x-ray technique factors and the internal calibration method, the RMS errors of the estimated volume fractions of red marrow of the test inserts were 1.2-1.3 times greater in the medium body than in the small body phantom and 1.3-1.5 times greater in the large body than in the small body phantom. RMS errors of the calculated volume fractions of red marrow within 2 × 2 segmented subregions of the ROIs were 1.6-1.9 times greater than for no segmentation, and RMS errors for 3 × 3 segmented subregions were 2.3-2.7 times greater than those for no segmentation. Increasing the dose by a factor of 2 reduced the RMS errors of all constituent volume fractions by an average factor of 1.40 ± 0.29 for all segmentation schemes and body phantom sizes; increasing the dose by a factor of 4 reduced those RMS errors by an average factor of 1.71 ± 0.25. Results for external calibrations exhibited much larger RMS errors than size matched internal calibration. Use of an average body size external-to-internal calibration correction factor reduced the errors to closer to those for internal calibration. RMS errors of less than 30% or about 0.01 for the bone and 0.1 for the red marrow volume fractions would likely be satisfactory for human studies. Such accuracies were achieved for 3 × 3 segmentation of 5 mm slice images for: (a) internal calibration with 4 times dose for all size body phantoms, (b) internal calibration with 2 times dose for the small and medium size body phantoms, and (c) corrected external calibration with 4 times dose and all size body phantoms.

CONCLUSIONS

Phantom studies are promising and demonstrate the potential to use dual energy quantitative CT to estimate the spatial distributions of red marrow and bone within the vertebral spongiosa.

Alternating direction method of multiplier for tomography with nonlocal regularizers

The ordered subset expectation maximization (OSEM) algorithm approximates the gradient of a likelihood function using a subset of projections instead of using all projections so that fast image reconstruction is possible for emission and transmission tomography such as SPECT, PET, and CT. However, OSEM does not significantly accelerate reconstruction with computationally expensive regularizers such as patch-based nonlocal (NL) regularizers, because the regularizer gradient is evaluated for every subset. We propose to use variable splitting to separate the likelihood term and the regularizer term for penalized emission tomographic image reconstruction problem and to optimize it using the alternating direction method of multiplier (ADMM). We also propose a fast algorithm to optimize the ADMM parameter based on convergence rate analysis. This new scheme enables more sub-iterations related to the likelihood term. We evaluated our ADMM for 3-D SPECT image reconstruction with a patch-based NL regularizer that uses the Fair potential function. Our proposed ADMM improved the speed of convergence substantially compared to other existing methods such as gradient descent, EM, and OSEM using De Pierro's approach, and the limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm.

Post-reconstruction non-local means filtering methods using CT side information for quantitative SPECT

Quantitative SPECT techniques are important for many applications including internal emitter therapy dosimetry where accurate estimation of total target activity and activity distribution within targets are both potentially important for dose–response evaluations. We investigated non-local means (NLM) post-reconstruction filtering for accurate I-131 SPECT estimation of both total target activity and the 3D activity distribution. We first investigated activity estimation versus number of ordered-subsets expectation–maximization (OSEM) iterations. We performed simulations using the XCAT phantom with tumors containing a uniform and a non-uniform activity distribution, and measured the recovery coefficient (RC) and the root mean squared error (RMSE) to quantify total target activity and activity distribution, respectively. We observed that using more OSEM iterations is essential for accurate estimation of RC, but may or may not improve RMSE. We then investigated various post-reconstruction filtering methods to suppress noise at high iteration while preserving image details so that both RC and RMSE can be improved. Recently, NLM filtering methods have shown promising results for noise reduction. Moreover, NLM methods using high-quality side information can improve image quality further. We investigated several NLM methods with and without CT side information for I-131 SPECT imaging and compared them to conventional Gaussian filtering and to unfiltered methods. We studied four different ways of incorporating CT information in the NLM methods: two known (NLM CT-B and NLM CT-M) and two newly considered (NLM CT-S and NLM CT-H). We also evaluated the robustness of NLM filtering using CT information to erroneous CT. NLM CT-S and NLM CT-H yielded comparable RC values to unfiltered images while substantially reducing RMSE. NLM CT-S achieved −2.7 to 2.6% increase of RC compared to no filtering and NLM CT-H yielded up to 6% decrease in RC while other methods yielded lower RCs than them: Gaussian filtering (up to 11.8% decrease in RC), NLM method without CT (up to 9.5% decrease in RC), and NLM CT-M and NLM CT-B (up to 19.4% decrease in RC). NLM CT-S and NLM CT-H achieved 8.2 to 33.9% and −0.9 to 36% decreased RMSE on tumors compared to no filtering respectively while other methods yielded less reduced or increased RMSE: Gaussian filtering (up to 7.9% increase in RMSE), NLM method without CT (up to 18.3% increase in RMSE), and NLM CT-M and NLM CT-B (up to 31.5% increase in RMSE). NLM CT-S and NLM CT-H also yielded images with tumor shapes that better-matched the true shapes than other methods. All NLM methods using CT information were robust to small misregistration between SPECT and CT, but NLM CT-S and NLM CT-H were more sensitive than NLM CT-M and NLM CT-B to missing CT information.

Bias atlases for segmentation-based PET attenuation correction using PET-CT and MR

This study was to obtain voxel-wise PET accuracy and precision using tissue-segmentation for attenuation correction. We applied multiple thresholds to the CTs of 23 patients to classify tissues. For six of the 23 patients, MR images were also acquired. The MR fat/in-phase ratio images were used for fat segmentation. Segmented tissue classes were used to create attenuation maps, which were used for attenuation correction in PET reconstruction. PET bias images were then computed using the PET reconstructed with the original CT as the reference. We registered the CTs for all the patients and transformed the corresponding bias images accordingly. We then obtained the mean and standard deviation bias atlas using all the registered bias images. Our CT-based study shows that four-class segmentation (air, lungs, fat, other tissues), which is available on most PET-MR scanners, yields 15.1%, 4.1%, 6.6%, and 12.9% RMSE bias in lungs, fat, non-fat soft-tissues, and bones, respectively. An accurate fat identification is achievable using fat/in-phase MR images. Furthermore, we have found that three-class segmentation (air, lungs, other tissues) yields less than 5% standard deviation of bias within the heart, liver, and kidneys. This implies that three-class segmentation can be sufficient to achieve small variation of bias for imaging these three organs. Finally, we have found that inter- and intra-patient lung density variations contribute almost equally to the overall standard deviation of bias within the lungs.

Cardiac motion compensation and resolution modeling in simultaneous PET-MR: a cardiac lesion detection study

Cardiac motion and partial volume effects (PVE) are two of the main causes of image degradation in cardiac PET. Motion generates artifacts and blurring while PVE lead to erroneous myocardial activity measurements. Newly available simultaneous PET-MR scanners offer new possibilities in cardiac imaging as MRI can assess wall contractility while collecting PET perfusion data. In this perspective, we develop a list-mode iterative reconstruction framework incorporating both tagged-MR derived non-rigid myocardial wall motion and position dependent detector point spread function (PSF) directly into the PET system matrix. In this manner, our algorithm performs both motion 'deblurring' and PSF deconvolution while reconstructing images with all available PET counts. The proposed methods are evaluated in a beating non-rigid cardiac phantom whose hot myocardial compartment contains small transmural and non-transmural cold defects. In order to accelerate imaging time, we investigate collecting full and half k-space tagged MR data to obtain tagged volumes that are registered using non-rigid B-spline registration to yield wall motion information. Our experimental results show that tagged-MR based motion correction yielded an improvement in defect/myocardium contrast recovery of 34-206% as compared to motion uncorrected studies. Likewise, lesion detectability improved by respectively 115-136% and 62-235% with MR-based motion compensation as compared to gating and no motion correction and made it possible to distinguish non-transmural from transmural defects, which has clinical significance given the inherent limitations of current single modality imaging in identifying the amount of residual ischemia. The incorporation of PSF modeling within the framework of MR-based motion compensation significantly improved defect/myocardium contrast recovery (5.1-8.5%, p < 0.01) and defect detectability (39-56%, p < 0.01). No statistical difference was found in PET contrast and lesion detectability based on motion fields obtained with half and full k-space tagged data.

Noise properties of motion-compensated tomographic image reconstruction methods

Motion-compensated image reconstruction (MCIR) methods incorporate motion models to improve image quality in the presence of motion. MCIR methods differ in terms of how they use motion information and they have been well studied separately. However, there have been less theoretical comparisions of different MCIR methods. This paper compares the theoretical noise properties of three popular MCIR methods assuming known nonrigid motion. We show the relationship among three MCIR methods-motion-compensated temporal regularization (MTR), the parametric motion model (PMM), and post-reconstruction motion correction (PMC)-for penalized weighted least square cases. These analyses show that PMM and MTR are matrix-weighted sums of all registered image frames, while PMC is a scalar-weighted sum. We further investigate the noise properties of MCIR methods with Poisson models and quadratic regularizers by deriving accurate and fast variance prediction formulas using an "analytical approach." These theoretical noise analyses show that the variances of PMM and MTR are lower than or comparable to the variance of PMC due to the statistical weighting. These analyses also facilitate comparisons of the noise properties of different MCIR methods, including the effects of different quadratic regularizers, the influence of the motion through its Jacobian determinant, and the effect of assuming that total activity is preserved. Two-dimensional positron emission tomography simulations demonstrate the theoretical results.

Correction for collimator-detector response in SPECT using point spread function template

Compensating for the collimator-detector response (CDR) in SPECT is important for accurate quantification. The CDR consists of both a geometric response and a septal penetration and collimator scatter response. The geometric response can be modeled analytically and is often used for modeling the whole CDR if the geometric response dominates. However, for radionuclides that emit medium or high-energy photons such as I-131, the septal penetration and collimator scatter response is significant and its modeling in the CDR correction is important for accurate quantification. There are two main methods for modeling the depth-dependent CDR so as to include both the geometric response and the septal penetration and collimator scatter response. One is to fit a Gaussian plus exponential function that is rotationally invariant to the measured point source response at several source-detector distances. However, a rotationally-invariant exponential function cannot represent the star-shaped septal penetration tails in detail. Another is to perform Monte-Carlo (MC) simulations to generate the depth-dependent point spread functions (PSFs) for all necessary distances. However, MC simulations, which require careful modeling of the SPECT detector components, can be challenging and accurate results may not be available for all of the different SPECT scanners in clinics. In this paper, we propose an alternative approach to CDR modeling. We use a Gaussian function plus a 2-D B-spline PSF template and fit the model to measurements of an I-131 point source at several distances. The proposed PSF-template-based approach is nearly non-parametric, captures the characteristics of the septal penetration tails, and minimizes the difference between the fitted and measured CDR at the distances of interest. The new model is applied to I-131 SPECT reconstructions of experimental phantom measurements, a patient study, and a MC patient simulation study employing the XCAT phantom. The proposed model yields up to a 16.5 and 10.8% higher recovery coefficient compared to the results with the conventional Gaussian model and the Gaussian plus exponential model, respectively.

Spatial resolution properties of motion-compensated tomographic image reconstruction methods

Many motion-compensated image reconstruction (MCIR) methods have been proposed to correct for subject motion in medical imaging. MCIR methods incorporate motion models to improve image quality by reducing motion artifacts and noise. This paper analyzes the spatial resolution properties of MCIR methods and shows that nonrigid local motion can lead to nonuniform and anisotropic spatial resolution for conventional quadratic regularizers. This undesirable property is akin to the known effects of interactions between heteroscedastic log-likelihoods (e.g., Poisson likelihood) and quadratic regularizers. This effect may lead to quantification errors in small or narrow structures (such as small lesions or rings) of reconstructed images. This paper proposes novel spatial regularization design methods for three different MCIR methods that account for known nonrigid motion. We develop MCIR regularization designs that provide approximately uniform and isotropic spatial resolution and that match a user-specified target spatial resolution. Two-dimensional PET simulations demonstrate the performance and benefits of the proposed spatial regularization design methods.

MRI-based nonrigid motion correction in simultaneous PET/MRI

Respiratory and cardiac motion is the most serious limitation to whole-body PET, resulting in spatial resolution close to 1 cm. Furthermore, motion-induced inconsistencies in the attenuation measurements often lead to significant artifacts in the reconstructed images. Gating can remove motion artifacts at the cost of increased noise. This paper presents an approach to respiratory motion correction using simultaneous PET/MRI to demonstrate initial results in phantoms, rabbits, and nonhuman primates and discusses the prospects for clinical application.

METHODS

Studies with a deformable phantom, a free-breathing primate, and rabbits implanted with radioactive beads were performed with simultaneous PET/MRI. Motion fields were estimated from concurrently acquired tagged MR images using 2 B-spline nonrigid image registration methods and incorporated into a PET list-mode ordered-subsets expectation maximization algorithm. Using the measured motion fields to transform both the emission data and the attenuation data, we could use all the coincidence data to reconstruct any phase of the respiratory cycle. We compared the resulting SNR and the channelized Hotelling observer (CHO) detection signal-to-noise ratio (SNR) in the motion-corrected reconstruction with the results obtained from standard gating and uncorrected studies.

RESULTS

Motion correction virtually eliminated motion blur without reducing SNR, yielding images with SNR comparable to those obtained by gating with 5-8 times longer acquisitions in all studies. The CHO study in dynamic phantoms demonstrated a significant improvement (166%-276%) in lesion detection SNR with MRI-based motion correction as compared with gating (P < 0.001). This improvement was 43%-92% for large motion compared with lesion detection without motion correction (P < 0.001). CHO SNR in the rabbit studies confirmed these results.

CONCLUSION

Tagged MRI motion correction in simultaneous PET/MRI significantly improves lesion detection compared with respiratory gating and no motion correction while reducing radiation dose. In vivo primate and rabbit studies confirmed the improvement in PET image quality and provide the rationale for evaluation in simultaneous whole-body PET/MRI clinical studies.

Nonrigid PET motion compensation in the lower abdomen using simultaneous tagged-MRI and PET imaging

PURPOSE

We propose a novel approach for PET respiratory motion correction using tagged-MRI and simultaneous PET-MRI acquisitions.

METHODS

We use a tagged-MRI acquisition followed by motion tracking in the phase domain to estimate the nonrigid deformation of biological tissues during breathing. In order to accurately estimate motion even in the presence of noise and susceptibility artifacts, we regularize the traditional HARP tracking strategy using a quadratic roughness penalty on neighboring displacement vectors (R-HARP). We then incorporate the motion fields estimated with R-HARP in the system matrix of an MLEM PET reconstruction algorithm formulated both for sinogram and list-mode data representations. This approach allows reconstruction of all detected coincidences in a single image while modeling the effect of motion both in the emission and the attenuation maps. At present, tagged-MRI does not allow estimation of motion in the lungs and our approach is therefore limited to motion correction in soft tissues. Since it is difficult to assess the accuracy of motion correction approaches in vivo, we evaluated the proposed approach in numerical simulations of simultaneous PET-MRI acquisitions using the NCAT phantom. We also assessed its practical feasibility in PET-MRI acquisitions of a small deformable phantom that mimics the complex deformation pattern of a lung that we imaged on a combined PET-MRI brain scanner.

RESULTS

Simulations showed that the R-HARP tracking strategy accurately estimated realistic respiratory motion fields for different levels of noise in the tagged-MRI simulation. In simulations of tumors exhibiting increased uptake, contrast estimation was 20% more accurate with motion correction than without. Signal-to-noise ratio (SNR) was more than 100% greater when performing motion-corrected reconstruction which included all counts, compared to when reconstructing only coincidences detected in the first of eight gated frames. These results were confirmed in our proof-of-principle PET-MRI acquisitions, indicating that our motion correction strategy is accurate, practically feasible, and is therefore ready to be tested in vivo.

CONCLUSIONS

This work shows that PET motion correction using motion fields measured with tagged-MRI in simultaneous PET-MRI acquisitions can be made practical for clinical application and that doing so has the potential to remove motion blur in whole-body PET studies of the torso.

A simple regularizer for B-spline nonrigid image registration that encourages local invertibility

Nonrigid image registration is an important task for many medical imaging applications. In particular, for radiation oncology it is desirable to track respiratory motion for thoracic cancer treatment. B-splines are convenient for modeling nonrigid deformations, but ensuring invertibility can be a challenge. This paper describes sufficient conditions for local invertibility of deformations based on B-spline bases. These sufficient conditions can be used with constrained optimization to enforce local invertibility. We also incorporate these conditions into nonrigid image registration methods based on a simple penalty approach that encourages diffeomorphic deformations. Traditional Jacobian penalty methods penalize negative Jacobian determinant values only at grid points. In contrast, our new method enforces a sufficient condition for invertibility directly on the deformation coefficients to encourage invertibility globally over a 3D continuous domain. The proposed penalty approach requires substantially less compute time than Jacobian penalties per iteration.

CDX2 promotes anchorage-independent growth by transcriptional repression of IGFBP-3

CDX2 is a Drosophila caudal-related homeobox transcription factor that is important for the establishment and maintenance of intestinal epithelial cells. We have reported that CDX2 promotes tumorigenicity in a subset of human colorectal cancer cell lines. Here, we present evidence that CDX2 negatively regulates the well-documented growth inhibitor insulin-like growth factor binding protein-3 (IGFBP-3). Specifically, CDX2 binds to the IGFBP-3 gene promoter and can repress IGFBP-3 transcription, protein expression and secretion. Furthermore, inhibition of IGFBP-3 partially rescues the decreased anchorage-independent growth phenotype observed in CDX2 knockout cells. These data demonstrate for the first time that (1) CDX2 can function as a transcriptional repressor, and (2) one mechanism by which CDX2 promotes anchorage-independent growth is by transcriptional repression of IGFBP-3.

CDX2 has tumorigenic potential in the human colon cancer cell lines LOVO and SW48

CDX2 is a Drosophila caudal-related homeobox transcription factor that is important for the establishment and maintenance of intestinal epithelial cells. CDX2 is a marker of colon cancer, with strong staining in up to 90% of colonic adenocarcinomas. CDX2 heterozygous-null mice develop colonic neoplasms, which have suggested that CDX2 is a tumor suppressor. However, CDX2 has not been reported to affect xenograft growth. Furthermore, CDX2 is rarely mutated in colon cancer, which has led to suggestions that it may play only a minor role as a tumor suppressor in colon cancer. To understand the functional contributions of CDX2 to colon cancer, we disrupted CDX2 in LOVO and SW48 human colon cancer cell lines by targeted homologous recombination. Consistent with the literature, disruption of CDX2 enhanced anchorage-dependent cell proliferation. However, homozygous loss of CDX2 led to significant inhibition of anchorage-independent growth in LOVO cells, and cell lethality in SW48 cells. Further analyses revealed that disruption of CDX2 led to anchorage-independent G1 to S growth arrest and anoikis. In vivo xenograft studies confirmed that disruption of CDX2 inhibited LOVO tumor growth. These data demonstrate that CDX2 mediates anchorage-independent growth and survival. Thus, CDX2 has tumorigenic potential in the human colon cancer cell lines LOVO and SW48.

Overexpression of calbindin-D28K induces neurite outgrowth in dopaminergic neuronal cells via activation of p38 MAPK

An MN9D dopaminergic neuronal cell line overexpressing calbindin-D28K (MN9D/Calbindin) was established in order to investigate directly the potential role of calcium-binding protein in neuronal differentiation. Overexpression of calbindin-D28K in MN9D cells resulted in significant increases in the number of neurites, the length of primary neurites, and the total extent of neurites. This robust neurite outgrowth occurred without cessation of cell division. Analysis of immunoblots revealed that this morphological differentiation was accompanied by increased expression of such markers of maturation as the synaptosomal protein SNAP-25. During calbindin-D28K-evoked neurite outgrowth in MN9D cells, phosphorylation of p38 mitogen-activated protein kinase (MAPK) dramatically increased while the levels and extent of phosphorylation of such other MAPKs as c-Jun N-terminal kinase (JNK) or extracellular response kinase (ERK) were not altered. Consequently, calbindin-D28K-induced neurite outgrowth was largely abolished by treatment with a p38 inhibitor, PD 169316, while the level of SNAP-25 in MN9D/Calbindin cells was not altered by this treatment. These data support an idea that calbindin-D28K and its associated p38 signaling pathway play a role in dopaminergic neuronal differentiation.

Stage-dependent regulation of ovarian pituitary adenylate cyclase-activating polypeptide mRNA levels by GnRH in cultured rat granulosa cells

The present study was designed to test whether GnRH regulates pituitary adenylate cyclase-activating polypeptide mRNA levels in a stage-dependent manner during follicle development in the rat ovary. The granulosa cells of preovulatory and immature follicles obtained from PMSG- and estrogen-treated rats, respectively, were cultured in serum-free conditions in the presence of various hormones. GnRH receptor mRNA expression was detected in both preovulatory and immature granulosa cells and was down-regulated by gonadotropins. Treatment of preovulatory granulosa cells with GnRH agonist stimulated pituitary adenylate cyclase-activating polypeptide mRNA levels in a dose-dependent manner. In situ hybridization analysis of cultured preovulatory follicles revealed that GnRH-induced pituitary adenylate cyclase- activating polypeptide signals were detected in granulosa cells, but not thecal cells. In immature granulosa cells, cotreatment with GnRH agonist suppressed FSH-stimulated pituitary adenylate cyclase-activating polypeptide mRNA levels in a dose-dependent manner, whereas treatment with GnRH alone had no effect. Furthermore, treatment with GnRH antagonist inhibited LH-induced pituitary adenylate cyclase-activating polypeptide gene expression in preovulatory granulosa cells, whereas it stimulated FSH-induced pituitary adenylate cyclase-activating polypeptide gene expression in immature granulosa cells. Interestingly, GnRH-stimulated pituitary adenylate cyclase-activating polypeptide mRNA levels in preovulatory granulosa cells was inhibited by arachidonyltri fluoromethyl ketone, an inhibitor of phospholipase A(2), but not by an inhibitor of protein kinase A or C. Lastly, treatment of preovulatory follicles with pituitary adenylate cyclase-activating polypeptide antagonist suppressed GnRH-stimulated progesterone production during 6--9 h of culture. Taken together, these results demonstrate the stage-dependent regulation of pituitary adenylate cyclase-activating polypeptide mRNA levels by GnRH, the stimulatory and inhibitory effect in granulosa cells of preovulatory and immature follicles, respectively.

MPP(+) downregulates mitochondrially encoded gene transcripts and their activities in dopaminergic neuronal cells: protective role of Bcl-2

The effects of neurotoxins on levels of mitochondrially encoded gene transcripts in a dopaminergic neuronal cell line, MN9D, were examined following treatment with 200 microM N-methyl-4-phenylpyridinium (MPP(+)) or 6-hydroxydopamine (6-OHDA). As confirmed by a Northern blot analysis, levels of cytochrome c oxidase subunit 3 (COX III) and ATPase subunit 6 (ATPase 6) transcript were decreased in a time-dependent manner following treatment with MPP(+) but not with 6-OHDA. Accordingly, enzymatic activity of cytochrome c oxidase (COX) and the intracellular ATP content were also decreased in MPP(+)-treated cells while these remained unaltered in 6-OHDA-treated cells. In the cell death paradigm induced by MPP(+), overexpression of Bcl-2 in MN9D cells (MN9D/Bcl-2) significantly blocked MPP(+)-induced downregulation of COX III and ATPase 6 transcripts. In MN9D/Bcl-2 cells, MPP(+)-induced downregulation of COX activity and the intracellular level of ATP was also blocked. Treatment with a pan-caspase inhibitor, however, neither prevented MPP(+)-induced downregulation of COX activity nor affected intracellular level of ATP in MN9D cells. Taken together, our present data suggest that Bcl-2 may play a regulatory role in energy metabolism by preventing downregulation of mitochondrially encoded gene(s) at a point distinct from its known anticaspase activity in MPP(+)-induced dopaminergic neuronal death.

Gonadotropin regulation of NGFI-B messenger ribonucleic acid expression during ovarian follicle development in the rat

NGFI-B is an immediate-early gene that encodes an orphan nuclear receptor. The present study was designed to examine the localization and gonadotropin regulation of NGFI-B expression in the rat ovary. Northern blot analysis of ovaries obtained from prepubertal rats revealed the increased expression of NGFI-B during prepubertal development. Treatment of immature rats with PMSG, however, decreased ovarian NGFI-B expression. The major cell types expressing NGFI-B messenger RNA were thecal cells of follicles in different sizes. In contrast, treatment of PMSG-primed rats with human (h) CG resulted in the rapid and transient stimulation of ovarian NGFI-B messenger RNA, reaching a peak within 1 h. In situ hybridization analysis revealed that hCG treatment induced the expression of NGFI-B in granulosa cells of preovulatory follicles. Treatment of cultured preovulatory follicles in vitro with LH further confirmed the time- and dose-dependent stimulation of NGFI-B messenger RNA and protein. LH-stimulated NGFI-B expression in preovulatory follicles was abolished by alpha-amanitin, but was superinduced by cycloheximide. Furthermore, treatment of adult cycling rats with pentobarbital abolished NGFI-B expression on proestrus, and exogenous administration of hCG restored it, indicating the role of the preovulatory surge of LH in the stimulation of NGFI-B expression. These results demonstrate the cell type-specific expression and gonadotropin induction of NGFI-B in granulosa cells of preovulatory follicles and suggest a role for NGFI-B in the ovulatory process.

Expression of messenger ribonucleic acid for the antiapoptosis gene P11 in the rat ovary: gonadotropin stimulation in granulosa cells of preovulatory follicles

P11, a member of the S100 family of calcium-binding proteins, has been shown to interact with BAD (Bcl-xL/Bcl-2-associated death promoter) in the yeast two-hybrid protein-protein interaction assay. Because overexpression of P11 dampens the proapoptotic activity of BAD in transfected cells, we tested the possibility that the expression of this antiapoptotic protein may be regulated by gonadotropins and other survival factors in the ovary. Northern blot analysis of ovaries obtained from prepubertal rats revealed an increased expression of P11 messenger RNA (mRNA) during prepubertal development in the theca cells of preantral and early antral follicles. Treatment of immature rats with PMSG did not affect P11 expression, whereas treatment of PMSG-primed rats with an ovulatory dose of human (h)CG stimulated ovarian P11 mRNA within 6-9 h in the granulosa cells of preovulatory follicles. Treatment of cultured preovulatory follicles in vitro with LH further confirmed the time-dependent stimulation of P11 by gonadotropins. In addition, treatment of cultured preovulatory follicles with MDL-12,330A, an inhibitor of adenylate cyclase, inhibited LH-stimulated P11 mRNA, whereas treatment with forskolin, an adenylate cyclase activator, but not the protein kinase C activator, 2-O-tetradecanol-phorbal-13-acetate, mimicked the LH action, suggesting the role of adenylate cyclase activation in P11 expression. Treatment with other follicle survival factors, including the epidermal growth factor, the basic fibroblast growth factor, and interleukin-1beta, could also stimulate P11 expression in cultured preovulatory follicles. These results demonstrate the expression of P11 mRNA in theca cells of different-sized follicles and in granulosa cells of preovulatory follicles following gonadotropin stimulation, and suggest that P11 may mediate, at least partially, the survival action of gonadotropins during the ovulatory process.

Three distinct types of GnRH receptor characterized in the bullfrog

It has been proposed recently that two types of GnRH receptors (GnRHR) exist in a particular species. Here we present data demonstrating that at least three types of GnRHR are expressed in a single diploid species, the bullfrog. Three different cDNAs, encoding distinct types of bullfrog GnRHR (bfGnRHR-1, bfGnRHR-2, and bfGnRHR-3), were isolated from pituitary and hindbrain of the bullfrog. BfGnRHR-1 mRNA was expressed predominantly in pituitary, whereas bfGnRHR-2 and -3 mRNAs were expressed in brain. The bfGnRHR-1, bfGnRHR-2, and bfGnRHR-3 proteins have an amino acid identity of approximately 30% to approximately 35% with mammalian GnRHRs and approximately 40% to approximately 50% with nonmammalian GnRHRs. Interestingly, bfGnRHR-2 has an 85% amino acid homology with Xenopus GnRHR. Less than 53% amino acid identity was observed among the three bfGnRHRs. All isolated cDNAs encode functional receptors because their transient expression in COS-7 cells resulted in a ligand-dependent increase in inositol phosphate production. Notably, all three receptors exhibited a differential ligand selectivity. For all receptors, cGnRH-II has a higher potency than mGnRH. In addition, salmon GnRH also has a strikingly high potency to stimulate all three receptors. In conclusion, we demonstrated the presence of three GnRHRs in the bullfrog. Their expression in pituitary and brain suggests that bfGnRHRs play an important role in the regulation of reproductive functions in the bullfrog.

Identification and characterization of a 19q12 amplicon in esophageal adenocarcinomas reveals cyclin E as the best candidate gene for this amplicon

Genomic DNA amplification in tumors is frequently associated with an increased gene copy number of oncogenes or other cancer-related genes. We have used a two-dimensional whole-genome scanning technique to identify gene amplification events in esophageal adenocarcinomas. A multicopy genomic fragment from a tumor two-dimensional gel was cloned, and genomic amplification encompassing this fragment was confirmed by Southern blot analysis. The corresponding DNA sequence was matched by BLAST to a BAC contig, which allowed the use of electronic-PCR to localize this amplicon to 19q12. Sequence tagged site-amplification mapping, an approach recently implemented in our laboratory (Lin, L. et al., Cancer Res., 60: 1341-1347,2000), was used to characterize the amplicon. Genomic DNA from 65 esophageal and 11 gastric cardia adenocarcinomas were investigated for 19q12 amplification using quantitative PCR at 11 sequence tagged site markers neighboring the cloned fragment. The amplicon was narrowed from >8 cM to a minimal critical region spanning <0.8 cM, between D19S919 and D19S882. This region includes the cyclin E gene. Fourteen expressed sequence tags (ESTs) covering the minimal region were then assayed for potential gene overexpression using quantitative reverse transcription-PCR. Seven of the selected ESTs were found to be both amplified and overexpressed. Among these seven ESTs, cyclin E showed the highest frequency of gene amplification and overexpression in the tumors examined, which allowed us to finalize the core-amplified region to <300 kb. These results indicate that cyclin E is the likely target gene selected by the amplification event at 19q12. The fact that cyclin E overexpression was found only in the amplified tumors examined indicates that gene amplification underlies the cyclin E gene overexpression. Our study represents the first extensive analysis of the 19q12 amplicon, and is the first to physically map the core-amplified domain to a region of <300 kb that includes cyclin E. Amplification of 19q12 was found neither in the 28 esophageal squamous cancers nor in the 39 lung adenocarcinomas examined but was observed in 13.8% of esophageal and 9.1% of gastric cardia adenocarcinomas.

The significance of CD44 variants expression in colorectal cancer and its regional lymph nodes

CD44 is a cell adhesion molecule with numerous isoforms created by mRNA alternative splicing. Expression of CD44 variants has been suggested to play a potential role in tumor progression and metastasis. We designed primers CD44V, CD44V6/7, CD44R1 and CD44V6-10 to analyze and compare the roles of each CD44 variants. Expressions of CD44 variants were investigated in normal colonic mucosa, the lymph nodes which was histopathologically free of cancer cell, and cancer tissues of 44 human colorectal cancer patients by RT-PCR method. The expression of CD44V was observed in 28 out of 39 (71.8%) tumors and 7 out of 11 (63.6%) N1 normal regional lymph nodes, and CD44V6/7 was observed in 28 out of 39 (71.8%) tumors and 9 out of 11 (81.8%) N1 normal regional lymph nodes. The expressions of CD44V and CD44V6/7 were most frequently observed compared with any other CD44 variants. In normal colonic mucosa, the expression of CD44 variants are low but in cancer tissue and its regional lymph node, the expression of CD44V and CD44V6/7 were significantly higher and more frequent than any other CD44 variants (p<0.05). These results suggest that CD44V and CD44V6/7 can be a molecular marker for colorectal cancer and its micrometastasis to the regional normal lymph node.

Progesterone increases mRNA levels of pituitary adenylate cyclase-activating polypeptide (PACAP) and type I PACAP receptor (PAC(1)) in the rat hypothalamus

Pituitary adenylate cyclase-activating polypeptide (PACAP) regulates pituitary hormone biosynthesis and secretion through its cognate receptors. PACAP also plays an important role in the regulation of ovarian steroid biosynthesis. If so, there might be a feedback regulation of hypothalamic PACAP synthesis by the pituitary and by ovarian steroids. In the present study, we used RNase protection assays to determine changes in mRNA levels of PACAP and type I PACAP receptor (PAC(1)) under the conditions of ovariectomy and replacement with ovarian steroids. Progesterone (P) alone or in combination with estradiol (E) induced significant increases in PACAP mRNA level in the medial basal hypothalamus (MBH) and PAC(1) mRNA levels in MBH and the preoptic area (POA). This finding suggests that feedback regulation takes place between the ovary and hypothalamic PACAP neurons. P is known to be a major regulatory feedback factor for hypothalamic luteinizing hormone-releasing hormone (LHRH) neurons, but P receptor is not present in these neurons. Therefore, we examined a possible involvement of PACAP in the feedback regulatory pathway of P to LHRH neurons. After an antisense PAC(1) oligodeoxynucleotide (ODN) was i.c.v.-injected into the third ventricle of E and P-treated rats, LHRH mRNA levels were determined. The ODN markedly decreased the P-induced increase in the LHRH mRNA level. Taken together, the present data suggest that PACAP may play a role as a mediator in the regulation of LHRH synthetic machinery by stimulatory feedback of P.

Expression of pituitary adenylate cyclase activating polypeptide (PACAP) and PACAP type I A receptor mRNAs in granulosa cells of preovulatory follicles of the rat ovary

Pituitary adenylate cyclase activating polypeptide (PACAP) was isolated from ovine hypothalamus and known to stimulate the production of cAMP in anterior pituitary cells. In the recent report, the expression of PACAP was detected in preovulatory follicles, and treatment with PACAP stimulated the production of progesterone and prostaglandin E(2) through the action of AC and PLC pathways in the ovary. PACAP binds to three type receptors. Type I A receptor is coupled to adenylate cyclase (AC) and phospholipase C (PLC) pathways, while type I B and type II receptors are only coupled to AC. Thus, the present study aimed to evaluate the temporal expression of PACAP and its type I A receptor mRNAs in the rat ovary after treatment with pregnant mare's serum gonadotropin and human chorionic gonadotropin (hCG). Northern blot analysis showed that PACAP transcripts were transiently expressed from 3-9 hr after hCG treatment, reaching a maximum at 6 hr. During these time points, PACAP mRNAs were specifically and strongly expressed in granulosa cells and cumulus cells of large preovulatory follicles and interstitial glandular cells. Type I A receptor mRNAs were also transiently expressed in granulosa cells of large preovulatory follicles from 3-9 hr after hCG treatment. PACAP and its type I A receptor mRNAs were expressed in the same preovulatory follicles. These results demonstrate that PACAP acts as an autoregulator or pararegulator through type I A receptor in granulosa cells and cumulus cells of large preovulatory follicles. Thus, we suggest that PACAP may have a critical role in granulosa cells of preovulatory follicles for the preparation of ovulation.

Involvement of progesterone in gonadotrophin-induced pituitary adenylate cyclase-activating polypeptide gene expression in pre-ovulatory follicles of rat ovary

The present study was designed to determine whether progesterone might have a role in gonadotrophin-induced pituitary adenylate cyclase-activating polypeptide (Pacap) gene expression in rat ovary. Northern blot analysis revealed that treatment of pregnant mare's serum gonadotrophin (PMSG)-primed immature rats with the progestin antagonist RU486 or an inhibitor of 3beta-hydroxysteroid dehydrogenase epostane, 1 h before HCG, resulted in a dose-dependent inhibition of the HCG-induced Pacap gene expression. In-situ hybridization demonstrated that the number of pre-ovulatory follicles expressing Pacap mRNA in their granulosa cells was greatly reduced in ovaries treated with RU486. Moreover, the suppressive effect of RU486 or epostane on the LH-induced Pacap gene expression in cultured pre-ovulatory follicles was reversed by co-treatment with the synthetic progestin R5020. We further cloned the 5'-flanking region of the rat Pacap gene and identified the presence of a consensus progesterone receptor element. When luciferase fusion genes containing Pacap gene promoter were transiently transfected into granulosa cells of pre-ovulatory follicles, luciferase activity was markedly stimulated by LH. Treatment with RU486 or epostane resulted in partial suppression of LH-stimulated PACAP promoter activity. Taken together, these results indicate that progesterone, acting through progesterone receptors, plays a role in gonadotrophin induction of Pacap gene expression in granulosa cells of pre-ovulatory follicles, and thereby may be involved in the process of ovulation.

Stage-specific expression of pituitary adenylate cyclase-activating polypeptide type I receptor messenger ribonucleic acid during ovarian follicle development in the rat

Expression of pituitary adenylate cyclase-activating polypeptide (PACAP), a neuropeptide with considerable homology to vasoactive intestinal peptide, has been shown to be stimulated by gonadotropins in the ovary. The present studies further evaluated the cell-type specific expression and gonadotropin regulation of PACAP type I receptor (PACAPR) messenger RNA in immature rat ovaries and in cultured preovulatory follicles. Northern blot analysis of ovaries obtained from prepubertal rats revealed the increased expression of PACAPR during prepubertal development. The major cell types expressing PACAPR messenger RNA were granulosa cells of large preantral follicles. Treatment of immature rats with PMSG caused a decrease in ovarian PACAPR expression. In contrast, treatment with human (h) CG at 2 days after PMSG treatment stimulated ovarian PACAPR messenger RNA within 3-6 h in granulosa cells of preovulatory follicles. Treatment of cultured preovulatory follicles in vitro with LH further confirmed the time- and dose-dependent stimulation of PACAPR by gonadotropins in granulosa cells of preovulatory follicles. Moreover, RNase protection assay revealed that the short variant of ovarian PACAPR was the predominant form stimulated during prepubertal development and by gonadotropins. These results demonstrate the expression of PACAPR messenger RNA in granulosa cells of growing follicles and of preovulatory follicles stimulated by gonadotropins, and suggest that PACAP may play a role in the growth of developing follicles and in ovulation as an autocrine/paracrine factor.

Characterization of the antiapoptotic Bcl-2 family member myeloid cell leukemia-1 (Mcl-1) and the stimulation of its message by gonadotropins in the rat ovary

The majority of ovarian follicles undergo atresia mediated by apoptosis. Bcl-2-related proteins act as regulators of apoptosis via the formation of dimers with proteins inside and outside the Bcl-2 family. Previous studies have identified BAD as a proapoptotic Bcl-2 family member expressed in the ovary. It is known that BAD phosphorylation induced by survival factors leads to its preferential binding to 14-3-3 and suppression of the death-inducing function of BAD. To identify ovarian binding partners for hypophosphorylated BAD, we performed a yeast two-hybrid screening of a rat ovary complementary DNA library using as bait a mutant BAD incapable of binding to 14-3-3. Screening of yeast transformants yielded positive clones encoding the rat ortholog of Mcl-1 (myeloid cell leukemia-1), an antiapoptotic Bcl-2 protein. Amino acid sequence analysis revealed that rat and human Mcl-1 showed a complete conservation of the Bcl-2 homology domains BH1, BH2, and BH3. In the yeast two-hybrid system, Mcl-1 binds to the hypophosphorylated mutant of BAD and interacts preferentially with different proapoptotic (Bax, Bak, Bok, Bik, and BOD) compared with antiapoptotic Bcl-2 family members (Bcl-2, Bcl-xL, Bcl-w, Bfl-1, CED-9, and BHRF-1). Northern blot hybridization demonstrated expression of Mcl-1 transcripts of 2.3 and 3.7 kb in the ovary and diverse other rat tissues. In immature rats, PMSG treatment led to a transient increase in the 2.3-kb Mcl-1 transcript, peaking at 6 h after injection and returning to baseline levels after 24 h. Moreover, the same transcript was induced in the PMSG-primed preovulatory rat ovary 6 h after the administration of ovulatory doses of either hCG or FSH. In situ hybridization studies revealed that the gonadotropin stimulation of ovarian Mcl-1 message occurs in both granulosa and thecal cells. In conclusion, rat Mcl-1 was identified as an ovarian BAD-interacting protein and the message for the antiapoptotic Mcl-1 protein was induced after treatment with gonadotropins in granulosa and thecal cells of growing follicles.

Type 4 cyclic adenosine monophosphate-specific phosphodiesterases are expressed in discrete subcellular compartments during rat spermiogenesis

The type 4 cAMP-specific phosphodiesterases (PDE4) are a family of closely related enzymes with similar catalytic domains and divergent amino- and carboxyl-terminus domains. Multiple PDE proteins with heterogeneous amino termini are derived from each gene. To understand the significance of this heterogeneity, the expression and localization of variants derived from PDE4A and PDE4D genes was investigated during spermatogenesis in the rat. RNase protection analysis with mRNA for testes at different ages of development showed that two transcripts (PDE4D1 and PDE4D2) are expressed at day 10 and 15 of age and become undetectable thereafter. An additional PDE4D transcript appears at day 30 and increased during testid maturation. This latter transcript codes for a long variant of the PDE4D gene and is expressed in germ cells as demonstrated by RNase protection with RNA from isolated pachytene spermatocytes and round spermatids. The presence of a corresponding PDE4D protein with a molecular mass of 98 kDa was established by immunoprecipitation and Western blot analysis with antibodies specific for PDE4D and by immunoaffinity chromatography purification of the 98 kDa variant from isolated germ cells. PDE4A transcripts were also expressed in pachytene spermatocytes and round spermatids. Two polypeptides encoded by these PDE4A transcripts were expressed in pachytene spermatocytes, reached a maximum in round spermatids, and declined thereafter. Immunofluorescence analysis demonstrated a localization of the PDE4D protein in the manchette and in a periacrosomal region of the developing spermatid, a localization confirmed by immunogold electron microscopy. Conversely, the PDE4A was mostly soluble in the cytoplasm of round spermatids. These data demonstrate that PDE4D and PDE4A variants are expressed at different stages and localized in distinct subcellular structures of developing spermatids. Different properties of the mRNAs derived from the two genes and localization signals are responsible for the temporal and spatial expression of the different PDE4 isoenzymes.

Keratinocyte growth factor promotes the survival, growth, and differentiation of preantral ovarian follicles

OBJECTIVE

To determine the effect of treatment with keratinocyte growth factor (KGF) on the survival of cells in cultured preantral follicles and on the growth and differentiation of preantral follicles.

DESIGN

Preantral follicles (140-150 microm) were dissected mechanically from the ovaries of 14-day-old rats and cultured for 24 hours with and without KGF. Genomic DNA was extracted, labeled with [32P]-dideoxyadenosine triphosphate, and fractionated through agarose gels. For growth studies, the follicles were cultured individually in 96-well dishes. After 72 hours, the follicles were collected and their protein or DNA content was evaluated and their inhibin-alpha content was determined.

RESULT(S)

Keratinocyte growth factor suppressed apoptosis in cultured preantral follicles by 60%. Treatment with KGF or FSH increased follicle diameter by 8% and 16%, respectively, and combined treatment with KGF and FSH increased follicle diameter by 26%. Western blot analysis demonstrated increased expression of inhibin-alpha content after treatment with KGF (2-fold), treatment with FSH (4-fold), and combined treatment with FSH and KGF (12-fold), demonstrating the effect of KGF on preantral follicle differentiation.

CONCLUSION(S)

Treatment with KGF promotes the survival, growth, and differentiation of cultured preantral follicles. Keratinocyte growth factor produced by theca cells may play a role in the progression of early follicle development.

Gonadotropin stimulation of pituitary adenylate cyclase-activating polypeptide (PACAP) messenger ribonucleic acid in the rat ovary and the role of PACAP as a follicle survival factor

Pituitary adenylate cyclase-activating polypeptide (PACAP), a novel neuropeptide with considerable homology to vasoactive intestinal peptide and GH-releasing hormone, exists in two biologically active forms, PACAP-38 and -27. The presence of PACAP in the ovary has been demonstrated, where it stimulates steroidogenesis and cAMP accumulation in cultured granulosa cells. In the present study, gonadotropin regulation of PACAP gene expression was examined in PMSG/human (h)CG-treated immature rat ovaries and cultured preovulatory follicles. Northern blot analysis of ovaries obtained from PMSG/hCG-treated immature animals revealed the transient induction of PACAP transcripts by hCG, reaching a maximum at 6 h. The major cell types expressing PACAP messenger RNA were granulosa cells of preovulatory follicles and some theca/interstitial cells. In preovulatory follicles cultured in serum-free medium, PACAP transcripts were transiently induced by LH and FSH, reaching a maximum 6-9 h after stimulation in granulosa cells but not in theca cells. Treatment with cycloheximide or alpha-amanitin abolished LH-induced PACAP transcripts, indicating that new protein synthesis and transcription are necessary. Treatment with MDL-12,330A, an inhibitor of adenylate cyclase, inhibited LH-induced PACAP messenger RNA, and forskolin mimicked the LH action, implying the role of adenylate cyclase activation. In contrast, treatment with chelerythrine, an inhibitor of protein kinase C, and 2-O-tetradecanol-phorbol-13-acetate had no effect. We further tested the role of PACAP in follicle apoptosis using apoptotic DNA fragmentation analysis. Treatment with PACAP-38 suppressed follicle apoptosis in a dose-dependent manner. Moreover, the LH suppression of follicle apoptosis was partially blocked by cotreatment with PACAP-38 antagonist, indicating mediation by endogenous PACAP-38. These results suggest that PACAP, transiently induced by the gonadotropin surge, could be a local regulator of a number of events and may act as a follicle survival factor during the periovulatory period.

Restricted expression of WT1 messenger ribonucleic acid in immature ovarian follicles: uniformity in mammalian and avian species and maintenance during reproductive senescence

WT1 is a zinc finger protein with transcriptional repressor activity on several growth factor and growth factor receptor genes. In the ovary, a potential role for WT1 in the suppression of the development of immature follicles has been demonstrated. Here, gel retardation assays further showed that recombinant WT1 protein interacted with consensus DNA sequences in the inhibin-alpha gene promoter. We investigated the pattern of WT1 expression in a wide variety of species and also over the reproductive life span in rats. In chicken ovaries, Northern blot analysis revealed the presence of WT1 transcript in small healthy white follicles (1-5 mm in diameter) and its absence in small yellow (6-12 mm in diameter) or larger follicles (F1-F5). In pig and monkey ovaries, WT1 expression was limited to granulosa cells of preantral follicles, as shown by in situ hybridization analysis. In rats, Northern blot analyses demonstrated the presence of WT1 transcript in the ovaries of young (3-mo-old) and middle-aged (9-mo-old) rats on the proestrous day, with a decrease in old (12-mo-old) rats in persistent estrus. In situ hybridization analysis further suggested that the decrease in WT1 expression in aging ovaries was associated with fewer immature follicles. Thus, WT1 expression is restricted to immature follicles in diverse avian and mammalian species and over the reproductive life span in rats. These data demonstrated that WT1 is a marker for immature follicles and suggested a potential role of this transcriptional repressor in the slow growth of early follicles.

Paracrine mechanisms of ovarian follicle apoptosis

It has become evident that the physiological removal of cells through apoptosis during embryonic and postnatal development of multicellular organisms is a mandatory process to maintain a homeostatic state of the individual. In the ovary, massive cell death occurs during neonatal and postnatal life as an integral part of the normal ovarian development. At birth, mammalian ovaries are endowed with a fixed number of non-growing follicles that will be gradually recruited into a growing pool during reproductive life. Once follicles start growth they are either selected for ovulation or, for the majority of them, removed by apoptosis. Thus, removal of excess ovarian cells by apoptosis is necessary for normal development of the ovary. Despite the important role of follicle atresia in the maintenance of normal follicle development, studies on the hormonal control of follicle cell demise during follicle growth have not been possible until the recent development of apoptosis detection methods. Recent biochemical analysis has revealed the occurrence of internucleosomal DNA fragmentation, a hallmark of apoptosis, in atretic follicles and has facilitated the investigation into the intra-ovarian hormonal regulation of follicle atresia. This review summarizes the recent advances in the intra-ovarian hormonal mechanisms that control follicle apoptosis.

Transport of light ions in matter

A recent set of light ion experiments are analyzed using the Green's function method of solving the Boltzmann equation for ions of high charge and energy (the GRNTRN transport code) and the NUCFRG2 fragmentation database generator code. Although the NUCFRG2 code reasonably represents the fragmentation of heavy ions, the effects of light ion fragmentation requires a more detailed nuclear model including shell structure and short range correlations appearing as tightly bound clusters in the light ion nucleus. The most recent NUCFRG2 code is augmented with a quasielastic alpha knockout model and semiempirical adjustments (up to 30 percent in charge removal) in the fragmentation process allowing reasonable agreement with the experiments to be obtained. A final resolution of the appropriate cross sections must await the full development of a coupled channel reaction model in which shell structure and clustering can be accurately evaluated.

Preantral ovarian follicles in serum-free culture: suppression of apoptosis after activation of the cyclic guanosine 3',5'-monophosphate pathway and stimulation of growth and differentiation by follicle-stimulating hormone

Progression of preantral follicle development is essential to further follicle maturation and ovulation, but there are few models for studying the regulation of preantral follicle survival and growth. We have evaluated preantral follicle survival in vivo and in vitro, and have developed a serum-free rat follicle culture system that can be used to characterize the regulation of preantral follicle growth and differentiation. Analysis of ovarian cell DNA fragmentation during the first wave of follicle growth in the infantile rat indicated negligible apoptosis up to day 16 of age. However, a major increase in apoptosis was found by day 18, a time point associated with the appearance of large antral follicles. In situ analysis confirmed that apoptotic DNA fragments were limited to antral follicles. Culture of individual preantral follicles mechanically dissected from ovaries of 12- or 14-day-old rats in serum-free conditions led to major increases in follicle cell apoptosis, similar to that seen in cultures of antral and preovulatory follicles. In contrast to antral and preovulatory follicles, treatment of preantral follicles with gonadotropins or cAMP analogs did not prevent apoptosis. However, treatment with 8-bromo-cGMP or 10% serum suppressed apoptosis by 75% in cultured preantral follicles. In situ analysis identified granulosa cells as the cell type susceptible to apoptosis regulation. Taking advantage of the ability of the cGMP analog to suppress apoptosis, we evaluated the potential of FSH as a growth factor. In the absence of serum, FSH treatment for 48 h did not affect follicle size compared to controls; however, treatment with the cGMP analog together with FSH increased follicle diameter (13%; P < 0.01) and viable cells (2.4-fold; P < 0.01) compared to control values. Immunoblot analysis further indicated that the inhibin-alpha content of the cultured follicles was increased by treatment with the combination of FSH and 8-bromo-cGMP, demonstrating the induction of follicle cell differentiation during culture. Therefore, we demonstrated that activation of the cGMP pathway promotes the survival of cultured preantral follicles and that in the presence of alpha cGMP analog, FSH is a growth and differentiation factor for preantral follicles. The present serum-free follicle culture model system will be useful in further evaluation of the regulation of growth and differentiation of preantral follicles.

Early onset of deoxyribonucleic acid fragmentation during atresia of preovulatory ovarian follicles in rats

Atretic demise rather than ovulation is the ultimate fate of the vast majority of ovarian follicles in mammals, affecting 70-99.9% of the follicles in various species. Recent studies have established that atretic degeneration of follicles is an apoptotic process, heralded by endonuclease degradation of DNA at internucleosomal sites, which generates DNA fragments in size multiples of 185-200 bp that are seen as distinct ladder bands after agarose gel electrophoresis. Using the well-characterized model of inducing atresia of preovulatory follicles in vivo by hypophysectomy and analyzing DNA fragmentation by autoradiography of size-fractionated DNA labeled at the 3' ends by [12P] dideoxy-ATP, we have examined the timing of atretic changes. DNA degradation was related to morphological signs of atresia, ovulability, and changes in follicular steroidogenesis. Rats were hypophysectomized on the morning of the day of proestrus, after which largest follicles were collected at various times. DNA fragmentation was analyzed in groups of five follicles. The increase in DNA fragments of low molecular weight up to 4 h after hypophysectomy was negligible (101 +/- 10%; 0 h time = 100%) but progressed 8, 12, 24, 48, and 72 h after hypophysectomy (143 +/- 20%, 168 +/- 27%, 235 +/- 29%, 3299 +/- 1075%, and 2249 +/- 805%, respectively; p < 0.03, n = 5). At 48 and 72 h, the extent of DNA degradation was higher than that observed in follicles cultured in a serum-free medium for 24 h. Likewise, staining of DNA by 4',6-diamido-2-phenylindole hydrochloride revealed apoptotic nuclei at 8 h after hypophysectomy (p < 0.01), and the percentage of such nuclei progressively increased afterwards. Thus, the increase in DNA fragmentation appeared concomitantly with atretic changes observed in previous studies (a decrease in ovulability at 6 h, and a spontaneous increase in progesterone accumulation and decrease in androgen and estrogen in follicles explanted 6 h after hypophysectomy) and preceded atresia detectable by morphological changes at 24 h. Detection of internucleosomal DNA degradation in preovulatory follicles early in the atretic cascade underscores the central role of apoptosis in ovarian follicle atresia.

Tumor necrosis factor-alpha and its second messenger, ceramide, stimulate apoptosis in cultured ovarian follicles

In the mammalian ovary, only a small fraction of follicles fully mature and ovulate, while most of them die via apoptosis. Multiple factors promoting follicle survival have been identified, but intraovarian mediators of apoptosis are poorly known. Tumor necrosis factor-alpha (TNF alpha) is a cytokine capable of inducing apoptosis in diverse cell types, and the apoptotic effect of TNF alpha is, partially, coupled to the sphingomyelin signaling pathway with ceramide as a second messenger. Because TNF alpha has been localized in the rat ovary, and TNF alpha treatment increases granulosa cell ceramide production, we studied the effect of treatment with TNF alpha and ceramide on follicle apoptosis. Immature rats were implanted with diethylstilbestrol to stimulate the development of early antral follicles. Follicles were isolated and cultured in a serum-free medium for 24 h with or without hormone treatments. During culture, spontaneous follicle apoptosis occurred (10-fold increase in DNA fragmentation), which was partially blocked by 100 ng/ml FSH (60% suppression). The effect of FSH was counteracted by TNF alpha in a dose-dependent manner, with the maximal effect at 100 ng/ml TNF alpha (90% reversal of FSH action). In situ analysis indicated that the granulosa cell is the follicle cell type undergoing DNA fragmentation. A membrane-permeable ceramide analog, C2-ceramide N-acetyl sphingosine, mimicked the effect of TNF alpha and was able to completely abolish the action of FSH at 50 microM. In contrast, another ceramide analog, C2-dihydroceramide N-acetyl dihydrosphingosine, did not alter the effect of FSH, verifying the specificity of ceramide action. To study the mechanism of TNF alpha and ceramide action, the effect of sodium aurathiomalate (ATM), an inhibitor of interleukin-1 beta-converting enzyme/ced-3-related cystine proteases known to be essential in the execution of mammalian cell apoptosis, was studied. Treatment with ATM (1 mM) prevented the apoptosis-inducing effect of both TNF alpha and ceramide, suggesting a role for cysteine proteases in mediating follicle apoptosis. Treatment with either TNF alpha or ceramide increased both basal and FSH-stimulated progesterone production by cultured follicles. Concomitant treatment by ATM did not alter the stimulatory effect of TNF alpha or ceramide on progesterone production, ruling out nonspecific toxic effect of the inhibitor and indicating that the apoptotic and steroidogenic pathways are independent. In summary, treatment with TNF alpha or its second messenger, ceramide, stimulates apoptosis of early antral follicles in culture, suggesting a potential role for TNF alpha as an intraovarian regulator of follicle atresia by acting through the ceramide signaling pathway.

Oocyte maturation involves compartmentalization and opposing changes of cAMP levels in follicular somatic and germ cells: studies using selective phosphodiesterase inhibitors

The second messenger cAMP has been implicated in the regulation of mammalian and amphibian oocyte maturation. Although a decrease in intraoocyte levels of cAMP precedes germinal vesicle breakdown (GVBD), the gonadotropin induction of ovulation and oocyte maturation is associated with major increases of cAMP in ovarian follicles. In the mammalian system, isolated oocytes undergo spontaneous maturation in vitro but this process is blocked by treatment with a phosphodiesterase (PDE) inhibitor, IBMX, which increases intraoocyte cAMP levels. In contrast, the same inhibitor, when added to cultured follicles for a brief time, increases follicle cAMP levels, followed by the induction of GVBD. To resolve the paradoxical actions of this PDE inhibitor on the maturation of isolated and follicle-enclosed oocytes, we hypothesized that meiotic maturation requires opposing fluctuations of cAMP levels in the somatic granulosa and germ cells. Such opposing fluctuations may result from selective expression and regulation of PDEs in the somatic and germ cell compartments of the follicle. To test this hypothesis, PDE activity was manipulated in different follicular cells using type-specific inhibitors. The impact of the ensuing changes in cAMP levels in the two compartments was monitored by the induction of GVBD. In isolated oocytes, spontaneous GVBD was blocked by two inhibitors of type 3 PDE (cGMP-inhibited: CGI-PDE), milrinone and cilostamide. In contrast, treatment with an inhibitor for type 4 PDE (cAMP-specific), rolipram, was ineffective. These findings suggest that the oocyte expresses type 3 but not type 4 PDE and that increases in intraoocyte cAMP suppress GVBD. This hypothesis was confirmed by in situ hybridization studies with PDE3 and PDE4 probes. PDE3B mRNA was concentrated in oocytes while PDE4D was mainly expressed in granulosa cells. In cultured follicles, LH treatment induced oocyte maturation but the gonadotropin action was blocked by inhibitors of type 3 but not the type 4 PDE inhibitors. Furthermore, treatment with the type 4, but not the type 3, PDE inhibitor mimics the action of LH and induces oocyte maturation, presumably by increasing cAMP levels in granulosa cells. Our findings indicate that PDE subtypes 4 and 3 are located in follicle somatic and germ cells, respectively. Preferential inhibition of PDE 3 in the oocyte may lead to a delay in oocyte maturation without affecting the cAMP-induced ovulatory process in the somatic cells. Conversely, selective suppression of granulosa cell cAMP-PDE may enhance the gonadotropin induction of ovulation and oocyte maturation. Thus, in addition to the well-recognized differential expression and regulation of adenylate cyclase in the somatic and germ cell compartments of the follicle, we suggest that selective regulation and expression of PDEs may be involved in the regulation of cAMP levels and control of oocyte maturation in the preovulatory mammalian follicle.