A Multicompartmental Diffusion Model for Improved Assessment of Whole-Body Diffusion-weighted Imaging Data and Evaluation of Prostate Cancer Bone Metastases

Purpose To develop a multicompartmental signal model for whole-body diffusion-weighted imaging (DWI) and apply it to study the diffusion properties of normal tissue and metastatic prostate cancer bone lesions in vivo. Materials and Methods This prospective study (ClinicalTrials.gov: NCT03440554) included 139 men with prostate cancer (mean age, 70 years ± 9 [SD]). Multicompartmental models with two to four tissue compartments were fit to DWI data from whole-body scans to determine optimal compartmental diffusion coefficients. Bayesian information criterion (BIC) and model-fitting residuals were calculated to quantify model complexity and goodness of fit. Diffusion coefficients for the optimal model (having lowest BIC) were used to compute compartmental signal-contribution maps. The signal intensity ratio (SIR) of bone lesions to normal-appearing bone was measured on these signal-contribution maps and on conventional DWI scans and compared using paired t tests (α = .05). Two-sample t tests (α = .05) were used to compare compartmental signal fractions between lesions and normal-appearing bone. Results Lowest BIC was observed from the four-compartment model, with optimal compartmental diffusion coefficients of 0, 1.1 × 10^-3, 2.8 × 10^-3, and >3.0 ×10^-2 mm²/sec. Fitting residuals from this model were significantly lower than from conventional apparent diffusion coefficient mapping (P < .001). Bone lesion SIR was significantly higher on signal-contribution maps of model compartments 1 and 2 than on conventional DWI scans (P < .008). The fraction of signal from compartments 2, 3, and 4 was also significantly different between metastatic bone lesions and normal-appearing bone tissue (P ≤ .02). Conclusion The four-compartment model best described whole-body diffusion properties. Compartmental signal contributions from this model can be used to examine prostate cancer bone involvement. Keywords: Whole-Body MRI, Diffusion-weighted Imaging, Restriction Spectrum Imaging, Diffusion Signal Model, Bone Metastases, Prostate Cancer Clinical trial registration no. NCT03440554 Supplemental material is available for this article. © RSNA, 2023 See also commentary by Margolis in this issue.

An unsupervised learning approach to identify novel signatures of health and disease from multimodal data

BACKGROUND

Modern medicine is rapidly moving towards a data-driven paradigm based on comprehensive multimodal health assessments. Integrated analysis of data from different modalities has the potential of uncovering novel biomarkers and disease signatures.

METHODS

We collected 1385 data features from diverse modalities, including metabolome, microbiome, genetics, and advanced imaging, from 1253 individuals and from a longitudinal validation cohort of 1083 individuals. We utilized a combination of unsupervised machine learning methods to identify multimodal biomarker signatures of health and disease risk.

RESULTS

Our method identified a set of cardiometabolic biomarkers that goes beyond standard clinical biomarkers. Stratification of individuals based on the signatures of these biomarkers identified distinct subsets of individuals with similar health statuses. Subset membership was a better predictor for diabetes than established clinical biomarkers such as glucose, insulin resistance, and body mass index. The novel biomarkers in the diabetes signature included 1-stearoyl-2-dihomo-linolenoyl-GPC and 1-(1-enyl-palmitoyl)-2-oleoyl-GPC. Another metabolite, cinnamoylglycine, was identified as a potential biomarker for both gut microbiome health and lean mass percentage. We identified potential early signatures for hypertension and a poor metabolic health outcome. Additionally, we found novel associations between a uremic toxin, p-cresol sulfate, and the abundance of the microbiome genera Intestinimonas and an unclassified genus in the Erysipelotrichaceae family.

CONCLUSIONS

Our methodology and results demonstrate the potential of multimodal data integration, from the identification of novel biomarker signatures to a data-driven stratification of individuals into disease subtypes and stages-an essential step towards personalized, preventative health risk assessment.

Relationship between kurtosis and bi-exponential characterization of high b-value diffusion-weighted imaging: application to prostate cancer

BACKGROUND

High b-value diffusion-weighted imaging has application in the detection of cancerous tissue across multiple body sites. Diffusional kurtosis and bi-exponential modeling are two popular model-based techniques, whose performance in relation to each other has yet to be fully explored.

PURPOSE

To determine the relationship between excess kurtosis and signal fractions derived from bi-exponential modeling in the detection of suspicious prostate lesions.

MATERIAL AND METHODS

This retrospective study analyzed patients with normal prostate tissue (n = 12) or suspicious lesions (n = 13, one lesion per patient), as determined by a radiologist whose clinical care included a high b-value diffusion series. The observed signal intensity was modeled using a bi-exponential decay, from which the signal fraction of the slow-moving component was derived ( SFs). In addition, the excess kurtosis was calculated using the signal fractions and ADCs of the two exponentials ( KCOMP). As a comparison, the kurtosis was also calculated using the cumulant expansion for the diffusion signal ( KCE).

RESULTS

Both K and KCE were found to increase with SFs within the range of SFs commonly found within the prostate. Voxel-wise receiver operating characteristic performance of SFs, KCE, and KCOMP in discriminating between suspicious lesions and normal prostate tissue was 0.86 (95% confidence interval [CI] = 0.85 - 0.87), 0.69 (95% CI = 0.68-0.70), and 0.86 (95% CI = 0.86-0.87), respectively.

CONCLUSION

In a two-component diffusion environment, KCOMP is a scaled value of SFs and is thus able to discriminate suspicious lesions with equal precision . KCE provides a computationally inexpensive approximation of kurtosis but does not provide the same discriminatory abilities as SFs and KCOMP.

Voxel Level Radiologic-Pathologic Validation of Restriction Spectrum Imaging Cellularity Index with Gleason Grade in Prostate Cancer

PURPOSE

Restriction spectrum imaging (RSI-MRI), an advanced diffusion imaging technique, can potentially circumvent current limitations in tumor conspicuity, in vivo characterization, and location demonstrated by multiparametric magnetic resonance imaging (MP-MRI) techniques in prostate cancer detection. Prior reports show that the quantitative signal derived from RSI-MRI, the cellularity index, is associated with aggressive prostate cancer as measured by Gleason grade (GG). We evaluated the reliability of RSI-MRI to predict variance with GG at the voxel-level within clinically demarcated prostate cancer regions.

EXPERIMENTAL DESIGN

Ten cases were processed using whole mount sectioning after radical prostatectomy. Regions of tumor were identified by an uropathologist. Stained prostate sections were scanned at high resolution (75 μm/pixel). A grid of tiles corresponding to voxel dimensions was graded using the GG system. RSI-MRI cellularity index was calculated from presurgical prostate MR scans and presented as normalized z-score maps. In total, 2,795 tiles were analyzed and compared with RSI-MRI cellularity.

RESULTS

RSI-MRI cellularity index was found to distinguish between prostate cancer and benign tumor (t = 25.48, P < 0.00001). Significant differences were also found between benign tissue and prostate cancer classified as low-grade (GG = 3; t = 11.56, P < 0.001) or high-grade (GG ≥ 4; t = 24.03, P < 0.001). Furthermore, RSI-MRI differentiated between low and high-grade prostate cancer (t = 3.23; P = 0.003).

CONCLUSIONS

Building on our previous findings of correlation between GG and the RSI-MRI among whole tumors, our current study reveals a similar correlation at voxel resolution within tumors. Because it can detect variations in tumor grade with voxel-level precision, RSI-MRI may become an option for planning targeted procedures where identifying the area with the most aggressive disease is important. Clin Cancer Res; 22(11); 2668-74. ©2016 AACR.

Polygenic hazard score to guide screening for aggressive prostate cancer: development and validation in large scale cohorts

OBJECTIVES

To develop and validate a genetic tool to predict age of onset of aggressive prostate cancer (PCa) and to guide decisions of who to screen and at what age.

DESIGN

Analysis of genotype, PCa status, and age to select single nucleotide polymorphisms (SNPs) associated with diagnosis. These polymorphisms were incorporated into a survival analysis to estimate their effects on age at diagnosis of aggressive PCa (that is, not eligible for surveillance according to National Comprehensive Cancer Network guidelines; any of Gleason score ≥7, stage T3-T4, PSA (prostate specific antigen) concentration ≥10 ng/L, nodal metastasis, distant metastasis). The resulting polygenic hazard score is an assessment of individual genetic risk. The final model was applied to an independent dataset containing genotype and PSA screening data. The hazard score was calculated for these men to test prediction of survival free from PCa.

SETTING

Multiple institutions that were members of international PRACTICAL consortium.

PARTICIPANTS

All consortium participants of European ancestry with known age, PCa status, and quality assured custom (iCOGS) array genotype data. The development dataset comprised 31 747 men; the validation dataset comprised 6411 men.

MAIN OUTCOME MEASURES

Prediction with hazard score of age of onset of aggressive cancer in validation set.

RESULTS

In the independent validation set, the hazard score calculated from 54 single nucleotide polymorphisms was a highly significant predictor of age at diagnosis of aggressive cancer (z=11.2, P<10-16). When men in the validation set with high scores (>98th centile) were compared with those with average scores (30th-70th centile), the hazard ratio for aggressive cancer was 2.9 (95% confidence interval 2.4 to 3.4). Inclusion of family history in a combined model did not improve prediction of onset of aggressive PCa (P=0.59), and polygenic hazard score performance remained high when family history was accounted for. Additionally, the positive predictive value of PSA screening for aggressive PCa was increased with increasing polygenic hazard score.

CONCLUSIONS

Polygenic hazard scores can be used for personalised genetic risk estimates that can predict for age at onset of aggressive PCa.

Restriction spectrum imaging: An evolving imaging biomarker in prostate MRI

Restriction spectrum imaging (RSI) is a novel diffusion-weighted MRI technique that uses the mathematically distinct behavior of water diffusion in separable microscopic tissue compartments to highlight key aspects of the tissue microarchitecture with high conspicuity. RSI can be acquired in less than 5 min on modern scanners using a surface coil. Multiple field gradients and high b-values in combination with postprocessing techniques allow the simultaneous resolution of length-scale and geometric information, as well as compartmental and nuclear volume fraction filtering. RSI also uses a distortion correction technique and can thus be fused to high resolution T2-weighted images for detailed localization, which improves delineation of disease extension into critical anatomic structures. In this review, we discuss the acquisition, postprocessing, and interpretation of RSI for prostate MRI. We also summarize existing data demonstrating the applicability of RSI for prostate cancer detection, in vivo characterization, localization, and targeting.

LEVEL OF EVIDENCE

5 J. Magn. Reson. Imaging 2017;45:323-336.

Restriction spectrum imaging improves MRI-based prostate cancer detection

PURPOSE

To compare the diagnostic performance of restriction spectrum imaging (RSI), with that of conventional multi-parametric (MP) magnetic resonance imaging (MRI) for prostate cancer (PCa) detection in a blinded reader-based format.

METHODS

Three readers independently evaluated 100 patients (67 with proven PCa) who underwent MP-MRI and RSI within 6 months of systematic biopsy (N = 67; 23 with targeting performed) or prostatectomy (N = 33). Imaging was performed at 3 Tesla using a phased-array coil. Readers used a five-point scale estimating the likelihood of PCa present in each prostate sextant. Evaluation was performed in two separate sessions, first using conventional MP-MRI alone then immediately with MP-MRI and RSI in the same session. Four weeks later, another scoring session used RSI and T2-weighted imaging (T2WI) without conventional diffusion-weighted or dynamic contrast-enhanced imaging. Reader interpretations were then compared to prostatectomy data or biopsy results. Receiver operating characteristic curves were performed, with area under the curve (AUC) used to compare across groups.

RESULTS

MP-MRI with RSI achieved higher AUCs compared to MP-MRI alone for identifying high-grade (Gleason score greater than or equal to 4 + 3=7) PCa (0.78 vs. 0.70 at the sextant level; P < 0.001 and 0.85 vs. 0.79 at the hemigland level; P = 0.04). RSI and T2WI alone achieved AUCs similar to MP-MRI for high-grade PCa (0.71 vs. 0.70 at the sextant level). With hemigland analysis, high-grade disease results were similar when comparing RSI + T2WI with MP-MRI, although with greater AUCs compared to the sextant analysis (0.80 vs. 0.79).

CONCLUSION

Including RSI with MP-MRI improves PCa detection compared to MP-MRI alone, and RSI with T2WI achieves similar PCa detection as MP-MRI.

Prostate diffusion imaging with distortion correction

PURPOSE

Diffusion imaging in the prostate is susceptible to distortion from B0 inhomogeneity. Distortion correction in prostate imaging is not routinely performed, resulting in diffusion images without accurate localization of tumors. We performed and evaluated distortion correction for diffusion imaging in the prostate.

MATERIALS AND METHODS

28 patients underwent pre-operative MRI (T2, Gadolinium perfusion, diffusion at b=800 s/mm(2)). The restriction spectrum protocol parameters included b-values of 0, 800, 1500, and 4000 s/mm(2) in 30 directions for each nonzero b-value. To correct for distortion, forward and reverse trajectories were collected at b=0 s/mm(2). Distortion maps were generated to reflect the offset of the collected data versus the corrected data. Whole-mount histology was available for correlation.

RESULTS

Across the 27 patients evaluated (excluding one patient due to data collection error), the average root mean square distortion distance of the prostate was 3.1 mm (standard deviation, 2.2mm; and maximum distortion, 12 mm).

CONCLUSION

Improved localization of prostate cancer by MRI will allow better surgical planning, targeted biopsies and image-guided treatment therapies. Distortion distances of up to 12 mm due to standard diffusion imaging may grossly misdirect treatment decisions. Distortion correction for diffusion imaging in the prostate improves tumor localization.

Polygenic Overlap Between C-Reactive Protein, Plasma Lipids, and Alzheimer Disease

BACKGROUND

Epidemiological findings suggest a relationship between Alzheimer disease (AD), inflammation, and dyslipidemia, although the nature of this relationship is not well understood. We investigated whether this phenotypic association arises from a shared genetic basis.

METHODS AND RESULTS

Using summary statistics (P values and odds ratios) from genome-wide association studies of >200 000 individuals, we investigated overlap in single-nucleotide polymorphisms associated with clinically diagnosed AD and C-reactive protein (CRP), triglycerides, and high- and low-density lipoprotein levels. We found up to 50-fold enrichment of AD single-nucleotide polymorphisms for different levels of association with C-reactive protein, low-density lipoprotein, high-density lipoprotein, and triglyceride single-nucleotide polymorphisms using a false discovery rate threshold <0.05. By conditioning on polymorphisms associated with the 4 phenotypes, we identified 55 loci associated with increased AD risk. We then conducted a meta-analysis of these 55 variants across 4 independent AD cohorts (total: n=29 054 AD cases and 114 824 healthy controls) and discovered 2 genome-wide significant variants on chromosome 4 (rs13113697; closest gene, HS3ST1; odds ratio=1.07; 95% confidence interval=1.05-1.11; P=2.86×10(-8)) and chromosome 10 (rs7920721; closest gene, ECHDC3; odds ratio=1.07; 95% confidence interval=1.04-1.11; P=3.38×10(-8)). We also found that gene expression of HS3ST1 and ECHDC3 was altered in AD brains compared with control brains.

CONCLUSIONS

We demonstrate genetic overlap between AD, C-reactive protein, and plasma lipids. By conditioning on the genetic association with the cardiovascular phenotypes, we identify novel AD susceptibility loci, including 2 genome-wide significant variants conferring increased risk for AD.

MRI-Derived Restriction Spectrum Imaging Cellularity Index is Associated with High Grade Prostate Cancer on Radical Prostatectomy Specimens

Purpose: We evaluate a novel magnetic resonance imaging (MRI) technique to improve detection of aggressive prostate cancer (PCa).

Materials and Methods: We performed a retrospective analysis of pre-surgical prostate MRI scans using an advanced diffusion-weighted imaging technique called restriction spectrum imaging (RSI), which can be presented as a normalized z-score statistic. Scans were acquired prior to radical prostatectomy. Prostatectomy specimens were processed using whole-mount sectioning and regions of interest (ROIs) were drawn around individual PCa tumors. Corresponding ROIs were drawn on the MRI imaging and paired with ROIs in regions with no pathology. RSI z-score and conventional apparent diffusion coefficient (ADC) values were recorded for each ROI. Paired t-test, ANOVA, and logistic regression analyses were performed.

Results: We evaluated 28 patients with 64 ROIs (28 benign and 36 PCa). The mean difference in RSI z-score (PCa ROI–Benign ROI) was 2.17 (SE = 0.11; p < 0.001) and in ADC was 551 mm²/s (SE = 80 mm²/s; paired t-test, p < 0.001). The differences in the means among all groups (benign, primary Gleason 3, and primary Gleason 4) was significant for both RSI z-score (F_3,64 = 97.7, p < 0.001) and ADC (F_3,64 = 13.9, p < 0.001). A t-test was performed on only PCa tumor ROIs (n = 36) to determine PCa aggressiveness (Gleason 3 vs. Gleason 4) revealing that RSI z-score was still significant (p = 0.03), whereas, ADC values were no longer significant (p = 0.08). In multivariable analysis adjusting for age and race, RSI z-score was associated with PCa aggressiveness (OR 10.3, 95% CI: 1.4–78.0, p = 0.02) while ADC trended to significance (p = 0.07).

Conclusion: The RSI-derived normalized cellularity index is associated with aggressive PCa as determined by pathologic Gleason scores. Further utilization of RSI techniques may serve to enhance standardized reporting systems for PCa in the future.

Relative capability of MR imaging and FDG PET to depict changes associated with prodromal and early Alzheimer disease

PURPOSE

To quantify the effect sizes of regional metabolic and morphometric measures in patients with preclinical and mild Alzheimer disease (AD) to aid in the identification of noninvasive biomarkers for the early detection of AD.

MATERIALS AND METHODS

The study was conducted with institutional review board approval and in compliance with HIPAA regulations. Written informed consent was obtained from each participant or participant's legal guardian. Fluorine 18 fluorodeoxyglucose (FDG) positron emission tomography (PET) and magnetic resonance (MR) imaging data were analyzed from 80 healthy control (HC) subjects, 68 individuals with AD, and 156 with amnestic mild cognitive impairment (MCI), 69 of whom had single-domain amnestic MCI. Regions of interest (ROIs) were derived after coregistering FDG PET and MR images by using high-throughput, subject-specific procedures. The Cohen d effect sizes were calculated for 42 predefined ROIs across the brain. Statistical comparison of the largest overall effect sizes for MR imaging and PET was performed. Metabolic effect sizes were determined with and without accounting for regional atrophy. Discriminative accuracy of ROIs showing the largest effect sizes were compared by calculating receiver operating characteristic curves.

RESULTS

For all disease groups, the hippocampus showed the largest morphometric effect size and the entorhinal cortex showed the largest metabolic effect size. In mild AD, the Cohen d effect size for hippocampal volume (1.92) was significantly larger (P < .05) than that for entorhinal metabolism (1.43). Regression of regional atrophy substantially reduced most metabolic effects. For all group comparisons, the areas under the receiver operating characteristic curves were significantly larger for hippocampal volume than for entorhinal metabolism.

CONCLUSION

The current results show no evidence that FDG PET is more sensitive than MR imaging to the degeneration occurring in preclinical and mild AD, suggesting that an MR imaging finding may be a more practical clinical biomarker for early detection of AD.

Structural MRI biomarkers for preclinical and mild Alzheimer's disease

Noninvasive MRI biomarkers for Alzheimer's disease (AD) may enable earlier clinical diagnosis and the monitoring of therapeutic effectiveness. To assess potential neuroimaging biomarkers, the Alzheimer's Disease Neuroimaging Initiative is following normal controls (NC) and individuals with mild cognitive impairment (MCI) or AD. We applied high-throughput image analyses procedures to these data to demonstrate the feasibility of detecting subtle structural changes in prodromal AD. Raw DICOM scans (139 NC, 175 MCI, and 84 AD) were downloaded for analysis. Volumetric segmentation and cortical surface reconstruction produced continuous cortical surface maps and region-of-interest (ROI) measures. The MCI cohort was subdivided into single- (SMCI) and multiple-domain MCI (MMCI) based on neuropsychological performance. Repeated measures analyses of covariance were used to examine group and hemispheric effects while controlling for age, sex, and, for volumetric measures, intracranial vault. ROI analyses showed group differences for ventricular, temporal, posterior and rostral anterior cingulate, posterior parietal, and frontal regions. SMCI and NC differed within temporal, rostral posterior cingulate, inferior parietal, precuneus, and caudal midfrontal regions. With MMCI and AD, greater differences were evident in these regions and additional frontal and retrosplenial cortices; evidence for non-AD pathology in MMCI also was suggested. Mesial temporal right-dominant asymmetries were evident and did not interact with diagnosis. Our findings demonstrate that high-throughput methods provide numerous measures to detect subtle effects of prodromal AD, suggesting early and later stages of the preclinical state in this cross-sectional sample. These methods will enable a more complete longitudinal characterization and allow us to identify changes that are predictive of conversion to AD.

Neurons detect increases and decreases in oxygen levels using distinct guanylate cyclases

Homeostatic sensory systems detect small deviations in temperature, water balance, pH, and energy needs to regulate adaptive behavior and physiology. In C. elegans, a homeostatic preference for intermediate oxygen (O2) levels requires cGMP signaling through soluble guanylate cyclases (sGCs), proteins that bind gases through an associated heme group. Here we use behavioral analysis, functional imaging, and genetics to show that reciprocal changes in O2 levels are encoded by sensory neurons that express alternative sets of sGCs. URX sensory neurons are activated by increases in O2 levels, and require the sGCs gcy-35 and gcy-36. BAG sensory neurons are activated by decreases in O2 levels, and require the sGCs gcy-31 and gcy-33. The sGCs are instructive O2 sensors, as forced expression of URX sGC genes causes BAG neurons to detect O2 increases. Both sGC expression and cell-intrinsic dynamics contribute to the differential roles of URX and BAG in O2-dependent behaviors.

Alzheimer disease: quantitative structural neuroimaging for detection and prediction of clinical and structural changes in mild cognitive impairment

PURPOSE

To use structural magnetic resonance (MR) images to identify a pattern of regional atrophy characteristic of mild Alzheimer disease (AD) and to investigate whether presence of this pattern prospectively can aid prediction of 1-year clinical decline and increased structural loss in mild cognitive impairment (MCI).

MATERIALS AND METHODS

The study was conducted with institutional review board approval and compliance with HIPAA regulations. Written informed consent was obtained from each participant. High-throughput volumetric segmentation and cortical surface reconstruction methods were applied to MR images from 84 subjects with mild AD, 175 with MCI, and 139 healthy control (HC) subjects. Stepwise linear discriminant analysis was used to identify regions that best can aid discrimination of HC subjects from subjects with AD. A classifier trained on data from HC subjects and those with AD was applied to data from subjects with MCI to determine whether presence of phenotypic AD atrophy at baseline was predictive of clinical decline and structural loss.

RESULTS

Atrophy in mesial and lateral temporal, isthmus cingulate, and orbitofrontal areas aided discrimination of HC subjects from subjects with AD, with fully cross-validated sensitivity of 83% and specificity of 93%. Subjects with MCI who had phenotypic AD atrophy showed significantly greater 1-year clinical decline and structural loss than those who did not and were more likely to have progression to probable AD (annual progression rate of 29% for subjects with MCI who had AD atrophy vs 8% for those who did not).

CONCLUSION

Semiautomated, individually specific quantitative MR imaging methods can be used to identify a pattern of regional atrophy in MCI that is predictive of clinical decline. Such information may aid in prediction of patient prognosis and increase the efficiency of clinical trials.

Multi-modal imaging predicts memory performance in normal aging and cognitive decline

This study (n=161) related morphometric MR imaging, FDG-PET and APOE genotype to memory scores in normal controls (NC), mild cognitive impairment (MCI) and Alzheimer's disease (AD). Stepwise regression analyses focused on morphometric and metabolic characteristics of the episodic memory network: hippocampus, entorhinal, parahippocampal, retrosplenial, posterior cingulate, precuneus, inferior parietal, and lateral orbitofrontal cortices. In NC, hippocampal metabolism predicted learning; entorhinal metabolism predicted recognition; and hippocampal metabolism predicted recall. In MCI, thickness of the entorhinal and precuneus cortices predicted learning, while parahippocampal metabolism predicted recognition. In AD, posterior cingulate cortical thickness predicted learning, while APOE genotype predicted recognition. In the total sample, hippocampal volume and metabolism, cortical thickness of the precuneus, and inferior parietal metabolism predicted learning; hippocampal volume and metabolism, parahippocampal thickness and APOE genotype predicted recognition. Imaging methods appear complementary and differentially sensitive to memory in health and disease. Medial temporal and parietal metabolism and morphometry best explained memory variance. Medial temporal characteristics were related to learning, recall and recognition, while parietal structures only predicted learning.

A distributed chemosensory circuit for oxygen preference in C. elegans

The nematode Caenorhabditis elegans has complex, naturally variable behavioral responses to environmental oxygen, food, and other animals. C. elegans detects oxygen through soluble guanylate cyclase homologs (sGCs) and responds to it differently depending on the activity of the neuropeptide receptor NPR-1: npr-1(lf) and naturally isolated npr-1(215F) animals avoid high oxygen and aggregate in the presence of food; npr-1(215V) animals do not. We show here that hyperoxia avoidance integrates food with npr-1 activity through neuromodulation of a distributed oxygen-sensing network. Hyperoxia avoidance is stimulated by sGC-expressing oxygen-sensing neurons, nociceptive neurons, and ADF sensory neurons. In npr-1(215V) animals, the switch from weak aerotaxis on food to strong aerotaxis in its absence requires close regulation of the neurotransmitter serotonin in the ADF neurons; high levels of ADF serotonin promote hyperoxia avoidance. In npr-1(lf) animals, food regulation is masked by increased activity of the oxygen-sensing neurons. Hyperoxia avoidance is also regulated by the neuronal TGF-beta homolog DAF-7, a secreted mediator of crowding and stress responses. DAF-7 inhibits serotonin synthesis in ADF, suggesting that ADF serotonin is a convergence point for regulation of hyperoxia avoidance. Coalitions of neurons that promote and repress hyperoxia avoidance generate a subtle and flexible response to environmental oxygen.

Nitric oxide binding to prokaryotic homologs of the soluble guanylate cyclase beta1 H-NOX domain

The heme cofactor in soluble guanylate cyclase (sGC) is a selective receptor for NO, an important signaling molecule in eukaryotes. The sGC heme domain has been localized to the N-terminal 194 amino acids of the beta1 subunit of sGC and is a member of a family of conserved hemoproteins, called the H-NOX family (Heme-Nitric Oxide and/or OXygen-binding domain). Three new members of this family have now been cloned and characterized, two proteins from Legionella pneumophila (L1 H-NOX and L2 H-NOX) and one from Nostoc punctiforme (Np H-NOX). Like sGC, L1 H-NOX forms a 5-coordinate Fe(II)-NO complex. However, both L2 H-NOX and Np H-NOX form temperature-dependent mixtures of 5- and 6-coordinate Fe(II)-NO complexes; at low temperature, they are primarily 6-coordinate, and at high temperature, the equilibrium is shifted toward a 5-coordinate geometry. This equilibrium is fully reversible with temperature in the absence of free NO. This process is analyzed in terms of a thermally labile proximal Fe(II)-His bond and suggests that in both the 5- and 6-coordinate Fe(II)-NO complexes of L2 H-NOX and Np H-NOX, NO is bound in the distal heme pocket of the H-NOX fold. NO dissociation kinetics for L1 H-NOX and L2 H-NOX have been determined and support a model in which NO dissociates from the distal side of the heme in both 5- and 6-coordinate complexes.

Characterization of functional heme domains from soluble guanylate cyclase

Soluble guanylate cyclase (sGC) is a heterodimeric, nitric oxide (NO)-sensing hemoprotein composed of two subunits, alpha1 and beta1. NO binds to the heme cofactor in the beta1 subunit, forming a five-coordinate NO complex that activates the enzyme several hundred-fold. In this paper, the heme domain has been localized to the N-terminal 194 residues of the beta1 subunit. This fragment represents the smallest construct of the beta1 subunit that retains the ligand-binding characteristics of the native enzyme, namely, tight affinity for NO and no observable binding of O(2). A functional heme domain from the rat beta2 subunit has been localized to the first 217 amino acids beta2(1-217). These proteins are approximately 40% identical to the rat beta1 heme domain and form five-coordinate, low-spin NO complexes and six-coordinate, low-spin CO complexes. Similar to sGC, these constructs have a weak Fe-His stretch [208 and 207 cm(-)(1) for beta1(1-194) and beta2(1-217), respectively]. beta2(1-217) forms a CO complex that is very similar to sGC and has a high nu(CO) stretching frequency at 1994 cm(-)(1). The autoxidation rate of beta1(1-194) was 0.073/min, while the beta2(1-217) was substantially more stable in the ferrous form with an autoxidation rate of 0.003/min at 37 degrees C. This paper has identified and characterized the minimum functional ligand-binding heme domain derived from sGC, providing key details toward a comprehensive characterization.

Crystal structure of an oxygen-binding heme domain related to soluble guanylate cyclases

Soluble guanylate cyclases are nitric oxide-responsive signaling proteins in which the nitric oxide sensor is a heme-binding domain of unknown structure that we have termed the heme-NO and oxygen binding (H-NOX) domain. H-NOX domains are also found in bacteria, either as isolated domains, or are fused through a membrane-spanning region to methyl-accepting chemotaxis proteins. We have determined the crystal structure of an oxygen-binding H-NOX domain of one such signaling protein from the obligate anaerobe Thermoanaerobacter tengcongensis at 1.77-angstroms resolution, revealing a protein fold unrelated to known structures. Particularly striking is the structure of the protoporphyrin IX group, which is distorted from planarity to an extent not seen before in protein-bound heme groups. Comparison of the structure of the H-NOX domain in two different crystal forms suggests a mechanism whereby alteration in the degree of distortion of the heme group is coupled to changes on the molecular surface of the H-NOX domain and potentially to changes in intermolecular interactions.

Spectroscopic characterization of the soluble guanylate cyclase-like heme domains from Vibrio cholerae and Thermoanaerobacter tengcongensis

Soluble guanylate cyclase (sGC) is a nitric oxide- (NO-) sensing hemoprotein that has been found in eukaryotes from Drosophila to humans. Prokaryotic proteins with significant homology to the heme domain of sGC have recently been identified through genomic analysis. Characterization of two of these proteins is reported here. The first is a 181 amino acid protein cloned from Vibrio cholerae (VCA0720) that is encoded in a histidine kinase-containing operon. The ferrous unligated form of VCA0720 is 5-coordinate, high-spin. The CO complex is low-spin, 6-coordinate, and the NO complex is high-spin and 5-coordinate. These ligand-binding properties are very similar to those of sGC. The second protein is the N-terminal 188 amino acids of Tar4 (TtTar4H), a predicted methyl-accepting chemotaxis protein (MCP) from the strict anaerobe Thermoanaerobacter tengcongensis. TtTar4H forms a low-spin, 6-coordinate ferrous-oxy complex, the first of this sGC-related family that binds O2. TtTar4H has ligand-binding properties similar to those of the heme-containing O2 sensors such as AxPDEA1. sGC does not bind O2 despite having a porphyrin with a histidyl ligand like the globins. The results reported here, with sequence-related proteins from prokaryotes but in the same family as the sGC heme domain, show that these proteins have evolved to discriminate between ligands such as NO and O2; hence, we term this family H-NOX domains (heme-nitric oxide/oxygen).

Oxygen sensation and social feeding mediated by a C. elegans guanylate cyclase homologue

Specialized oxygen-sensing cells in the nervous system generate rapid behavioural responses to oxygen. We show here that the nematode Caenorhabditis elegans exhibits a strong behavioural preference for 5-12% oxygen, avoiding higher and lower oxygen levels. 3',5'-cyclic guanosine monophosphate (cGMP) is a common second messenger in sensory transduction and is implicated in oxygen sensation. Avoidance of high oxygen levels by C. elegans requires the sensory cGMP-gated channel tax-2/tax-4 and a specific soluble guanylate cyclase homologue, gcy-35. The GCY-35 haem domain binds molecular oxygen, unlike the haem domains of classical nitric-oxide-regulated guanylate cyclases. GCY-35 and TAX-4 mediate oxygen sensation in four sensory neurons that control a naturally polymorphic social feeding behaviour in C. elegans. Social feeding and related behaviours occur only when oxygen exceeds C. elegans' preferred level, and require gcy-35 activity. Our results suggest that GCY-35 is regulated by molecular oxygen, and that social feeding can be a behavioural strategy for responding to hyperoxic environments.

RNA polymerase II-dependent positional effects on mRNA 3' end processing in the adenovirus major late transcription unit

During the early phase of adenovirus infection, the promoter-proximal L1 poly(A) site in the major late transcription unit is used preferentially despite the fact that the distal L3 poly(A) site is stronger (i.e. it competes better for processing factors and is cleaved at a faster rate, in vitro). Previous work had established that this was due at least in part to the stable binding of the processing factor, cleavage and polyadenylation specificity factor, to the L1 poly(A) site as mediated by specific regulatory sequences. It is now demonstrated that in addition, the L1 poly(A) site has a positional advantage because of its 5' location in the transcription unit. We also show that preferential processing of a particular poly(A) site in a complex transcription unit is dependent on RNA polymerase II. Our results are consistent with recent reports demonstrating that the processing factors cleavage and polyadenylation specificity factor and cleavage stimulatory factor are associated with the RNA polymerase II holoenzyme; thus, processing at a weak poly(A) site like L1 can be enhanced by virtue of its being the first site to be transcribed.