Molecular and anatomical signatures of sleep deprivation in the mouse brain.

Sleep deprivation (SD) leads to a suite of cognitive and behavioral impairments, and yet the molecular consequences of SD in the brain are poorly understood. Using a systematic immediate-early gene (IEG) mapping to detect neuronal activation, the consequences of SD were mapped primarily to forebrain regions. SD was found to both induce and suppress IEG expression (and thus neuronal activity) in subregions of neocortex, striatum, and other brain regions. Laser microdissection and cDNA microarrays were used to identify the molecular consequences of SD in seven brain regions. In situ hybridization (ISH) for 222 genes selected from the microarray data and other sources confirmed that robust molecular changes were largely restricted to the forebrain. Analysis of the ISH data for 222 genes (publicly accessible at http://sleep.alleninstitute.org) provided a molecular and anatomic signature of the effects of SD on the brain. The suprachiasmatic nucleus (SCN) and the neocortex exhibited differential regulation of the same genes, such that in the SCN genes exhibited time-of-day effects while in the neocortex, genes exhibited only SD and waking (W) effects. In the neocortex, SD activated gene expression in areal-, layer-, and cell type-specific manner. In the forebrain, SD preferentially activated excitatory neurons, as demonstrated by double-labeling, except for striatum which consists primarily of inhibitory neurons. These data provide a characterization of the anatomical and cell type-specific signatures of SD on neuronal activity and gene expression that may account for the associated cognitive and behavioral effects.

Divergent and nonuniform gene expression patterns in mouse brain.

Considerable progress has been made in understanding variations in gene sequence and expression level associated with phenotype, yet how genetic diversity translates into complex phenotypic differences remains poorly understood. Here, we examine the relationship between genetic background and spatial patterns of gene expression across seven strains of mice, providing the most extensive cellular-resolution comparative analysis of gene expression in the mammalian brain to date. Using comprehensive brainwide anatomic coverage (more than 200 brain regions), we applied in situ hybridization to analyze the spatial expression patterns of 49 genes encoding well-known pharmaceutical drug targets. Remarkably, over 50% of the genes examined showed interstrain expression variation. In addition, the variability was nonuniformly distributed across strain and neuroanatomic region, suggesting certain organizing principles. First, the degree of expression variance among strains mirrors genealogic relationships. Second, expression pattern differences were concentrated in higher-order brain regions such as the cortex and hippocampus. Divergence in gene expression patterns across the brain could contribute significantly to variations in behavior and responses to neuroactive drugs in laboratory mouse strains and may help to explain individual differences in human responsiveness to neuroactive drugs.

Clustering of spatial gene expression patterns in the mouse brain and comparison with classical neuroanatomy.

Spatial gene expression profiles provide a novel means of exploring the structural organization of the brain. Computational analysis of these patterns is made possible by genome-scale mapping of the C57BL/6J mouse brain in the Allen Brain Atlas. Here we describe methodology used to explore the spatial structure of gene expression patterns across a set of 3041 genes chosen on the basis of consistency across experimental observations (N=2). The analysis was performed on smoothed, co-registered 3D expression volumes for each gene obtained by aggregating cellular resolution image data. Following dimensionality and noise reduction, voxels were clustered according to similarity of expression across the gene set. We illustrate the resulting parcellations of the mouse brain for different numbers of clusters (K) and quantitatively compare these parcellations with a classically-defined anatomical reference atlas at different levels of granularity, revealing a high degree of correspondence. These observations suggest that spatial localization of gene expression offers substantial promise in connecting knowledge at the molecular level with higher-level information about brain organization.

Genomic anatomy of the hippocampus.

Availability of genome-scale in situ hybridization data allows systematic analysis of genetic neuroanatomical architecture. Within the hippocampus, electrophysiology and lesion and imaging studies demonstrate functional heterogeneity along the septotemporal axis, although precise underlying circuitry and molecular substrates remain uncharacterized. Application of unbiased statistical component analyses to genome-scale hippocampal gene expression data revealed robust septotemporal molecular heterogeneity, leading to the identification of a large cohort of genes with robust regionalized hippocampal expression. Manual mapping of heterogeneous CA3 pyramidal neuron expression patterns demonstrates an unexpectedly complex molecular parcellation into a relatively coherent set of nine expression domains in the septal/temporal and proximal/distal axes with reciprocal, nonoverlapping boundaries. Unique combinatorial profiles of adhesion molecules within these domains suggest corresponding differential connectivity, which is demonstrated for CA3 projections to the lateral septum using retrograde labeling. This complex, discrete molecular architecture provides a novel paradigm for predicting functional differentiation across the full septotemporal extent of the hippocampus.

Seed-based transrectal ultrasound-fluoroscopy registration method for intraoperative dosimetry analysis of prostate brachytherapy.

Prostate brachytherapy is an effective treatment option for early-stage prostate cancer. During a prostate brachytherapy procedure, transrectal ultrasound (TRUS) and fluoroscopy imaging modalities complement each other by providing good visualization of soft tissue and implanted seeds, respectively. Therefore, the registration of these two imaging modalities, which are readily available in the operating room, could facilitate intraoperative dosimetry, thus enabling physicians to implant additional seeds into the underdosed portions of the prostate while the patient is still on the operating table. It is desirable to register TRUS and fluoroscopy images by using the seeds as fiducial markers. Although the locations of all the implanted seeds can be reconstructed from three fluoroscopy images, only a fraction of these seeds can be located in TRUS images. It is challenging to register the TRUS and fluoroscopy images by using the identified seeds, since the correspondence between them is unknown. Furthermore, misdetection of nonseed structures as seeds can lead to the inclusion of spurious points in the data set. We developed a new method called iterative optimal assignment (IOA) to overcome these challenges in TRUS-fluoroscopy registration. By using the Hungarian method in an optimization framework, IOA computes a set of transformation parameters that yield the one-to-one correspondence with minimum cost. We have evaluated our registration method at varying noise levels, seed detection rates, and number of spurious points using data collected from 25 patients. We have found that IOA can perform registration with an average root mean square error of about 0.2 cm even when the seed detection rate is only 10%. We believe that IOA can offer a robust solution to seed-based TRUS-fluoroscopy registration, thus making intraoperative dosimetry possible.

Semiautomatic 3-D prostate segmentation from TRUS images using spherical harmonics.

Prostate brachytherapy quality assessment procedure should be performed while the patient is still on the operating table since this would enable physicians to implant additional seeds immediately into the prostate if necessary thus reducing the costs and increasing patient outcome. Seed placement procedure is readily performed under fluoroscopy and ultrasound guidance. Therefore, it has been proposed that seed locations be reconstructed from fluoroscopic images and prostate boundaries be identified in ultrasound images to perform dosimetry in the operating room. However, there is a key hurdle that needs to be overcome to perform the ultrasound and fluoroscopy-based dosimetry: it is highly time-consuming for physicians to outline prostate boundaries in ultrasound images manually, and there is no method that enables physicians to identify three-dimensional (3-D) prostate boundaries in postimplant ultrasound images in a fast and robust fashion. In this paper, we propose a new method where the segmentation is defined in an optimization framework as fitting the best surface to the underlying images under shape constraints. To derive these constraints, we modeled the shape of the prostate using spherical harmonics of degree eight and performed statistical analysis on the shape parameters. After user initialization, our algorithm identifies the prostate boundaries on the average in 2 min. For algorithm validation, we collected 30 postimplant prostate volume sets, each consisting of axial transrectal ultrasound images acquired at 1-mm increments. For each volume set, three experts outlined the prostate boundaries first manually and then using our algorithm. By treating the average of manual boundaries as the ground truth, we computed the segmentation error. The overall mean absolute distance error was 1.26 +/- 0.41 mm while the percent volume overlap was 83.5 +/- 4.2. We found the segmentation error to be slightly less than the clinically-observed interobserver variability.

Parametric shape modeling using deformable superellipses for prostate segmentation.

Automatic prostate segmentation in ultrasound images is a challenging task due to speckle noise, missing boundary segments, and complex prostate anatomy. One popular approach has been the use of deformable models. For such techniques, prior knowledge of the prostate shape plays an important role in automating model initialization and constraining model evolution. In this paper, we have modeled the prostate shape using deformable superellipses. This model was fitted to 594 manual prostate contours outlined by five experts. We found that the superellipse with simple parametric deformations can efficiently model the prostate shape with the Hausdorff distance error (model versus manual outline) of 1.32 +/- 0.62 mm and mean absolute distance error of 0.54 +/- 0.20 mm. The variability between the manual outlinings and their corresponding fitted deformable superellipses was significantly less than the variability between human experts with p-value being less than 0.0001. Based on this deformable superellipse model, we have developed an efficient and robust Bayesian segmentation algorithm. This algorithm was applied to 125 prostate ultrasound images collected from 16 patients. The mean error between the computer-generated boundaries and the manual outlinings was 1.36 +/- 0.58 mm, which is significantly less than the manual interobserver distances. The algorithm was also shown to be fairly insensitive to the choice of the initial curve.

Quantitative image analysis: software systems in drug development trials.

Multi-dimensional image analysis is being used increasingly to arrive at surrogate end-points for drug development trials. Various imaging modalities such as computed tomography (CT), magnetic resonance imaging (MRI), positron emission tomography (PET) and ultrasound are used to analyze treatments for diseases such as cancer, multiple sclerosis, osteoarthritis, and Alzheimer's disease. However, extracting information from images can be tedious and is prone to high user variability. The medical image analysis community is moving towards advanced software systems specifically designed for drug development trials. These systems can automatically identify the anatomy of interest in medical images (segmentation methods), can compare the anatomy over time or between patients (registration methods) and allow the quantitative extraction of anatomical features and the integration of the data and results into a database management system, automatically tracking the changes made to the data (audit trail generation). In this article, we present a case study using a prototype system that is used for quantifying multiple sclerosis lesions from multivariate MRI.

Effects of video digitization in pubic arch interference assessment for prostate brachytherapy.

Recently, it has been shown that prior to surgery a transrectal ultrasound (TRUS) study of the prostate and pubic arch can effectively determine pubic arch interference (PAI), a major stumbling block for the prostate brachytherapy (radioactive seed implantation) procedure. This PAI determination is currently being done with digital images taken directly from an ultrasound (US) machine. However, 70-75% of US machines used in prostate brachytherapy do not have a method to save or transfer digital image data for external use. To allow PAI assessment regardless of US platform and to keep costs to a minimum, we need to digitize the images from the US video output when there is no direct digital transfer capability. D/A and A/D conversions can introduce quantization error and other noises in these digitized images. The purpose of this work is to assess the image degradation caused by digitization and quantitatively evaluate whether after digitization it is still possible to accurately assess PAI. We used a PAI assessment algorithm (developed in previous research by our group) to predict the location of the pubic arch on both digital images and those captured after digitization. These predicted arch locations were compared to the "true" position of the pubic arch as established during surgery. Despite apparent image degradation due to the D/A and A/D conversions, we found no statistically significant difference between the accuracy of the predicted arch locations from the digitized images and those from the digital images. By demonstrating equally accurate determination of pubic arch locations using digital and digitized images, we conclude that TRUS-based PAI assessment can be easily and inexpensively performed in clinics where it is needed.

Tomosynthesis-based localization of radioactive seeds in prostate brachytherapy.

Accurately assessing the quality of prostate brachytherapy intraoperatively would be valuable for improved clinical outcome by ensuring the delivery of a prescribed tumoricidal radiation dose to the entire prostate gland. One necessary step towards this goal is the robust and rapid localization of implanted seeds. Several methods have been developed to locate seeds from x-ray projection images, but they fail to detect completely-overlapping seeds, thus necessitating manual intervention. To overcome this limitation, we have developed a new method where (1) a three-dimensional volume is reconstructed from x-ray projection images using a brachytherapy-specific tomosynthesis reconstruction algorithm with built-in blur compensation and (2) the seeds are located in this reconstructed volume. In contrast to other projection-based methods, our method can detect completely overlapping seeds. Our simulation results indicate that we can locate all implanted seeds in the prostate using a tomosynthesis angle of 30 degrees and seven projection images. The mean localization error is 1.27 mm for a case with 100 seeds. We have also tested our method using a prostate phantom with 61 implanted seeds and succeeded in locating all seeds automatically. We believe this new method can be useful for the intraoperative quality assessment of prostate brachytherapy in the future.

Ultrasonography and fluoroscopic fusion for prostate brachytherapy dosimetry.

PURPOSE: To investigate the feasibility of performing postimplant and intraoperative dosimetry for prostate brachytherapy by fusing transrectal ultrasound (TRUS) and fluoroscopic data. METHODS AND MATERIALS: Registration of ultrasound (prostate boundary) and fluoroscopic (seed) data requires spatial markers that are detectable by both imaging modalities. In this study, the needle tips were considered as such fiducials. Prostate phantoms were implanted with the seeds, and four localization needles were inserted. In the TRUS frame of reference, the longitudinal coordinate of the needle tip was determined by advancing the needle until the echo from its tip just registered at a known probe depth. The tip's transverse coordinates were determined from the associated TRUS slice. The three-dimensional needle tip positions were also calculated in the fluoroscopic coordinate system using a seed reconstruction method. The transformation between the TRUS and fluoroscopy coordinate systems was established by the least-squares solution using the singular value decomposition. RESULTS: With three of four needle tips as fiducials and the one remaining needle as a test target, the mean fiducial registration error was 0.8 mm and the test target registration error was 2.5 mm. When all four points were used for registration, the errors decreased to 1.1 mm. A comparison between the proposed method and CT-based dosimetry yielded a percentage of prostate volume receiving 100% and 150% of the prescribed minimal peripheral dose and minimal dose received by 90% of the prostate gland that agreed within 0.4%, 2.7%, and 4.2%, respectively. CONCLUSION: The combination of TRUS and fluoroscopy is a feasible alternative to the currently used CT-based postimplant dosimetry. Furthermore, because of online imaging capability, the method lends itself to real-time intraoperative applications.