3D visualization software to analyze topological outcomes of topoisomerase reactions.

The action of various DNA topoisomerases frequently results in characteristic changes in DNA topology. Important information for understanding mechanistic details of action of these topoisomerases can be provided by investigating the knot types resulting from topoisomerase action on circular DNA forming a particular knot type. Depending on the topological bias of a given topoisomerase reaction, one observes different subsets of knotted products. To establish the character of topological bias, one needs to be aware of all possible topological outcomes of intersegmental passages occurring within a given knot type. However, it is not trivial to systematically enumerate topological outcomes of strand passage from a given knot type. We present here a 3D visualization software (TopoICE-X in KnotPlot) that incorporates topological analysis methods in order to visualize, for example, knots that can be obtained from a given knot by one intersegmental passage. The software has several other options for the topological analysis of mechanisms of action of various topoisomerases.

X-ray microtomography of Larger Foraminifera

X-ray microtomography has become a new tool in earth sciences to obtain non-destructive 3D-image data from geological objects in which variations in mineralogy, chemical composition and/or porosity create sufficient x-ray density contrasts.We present here first, preliminary results of an application to the external and internal morphology of Permian to Recent Larger Foraminifera. We use a SkyScan-1072 high-resolution desk-top micro-CT system. The system has a conical x-ray source with a spot size of about 5µm that runs at 20-100kV, 0-250µA, resulting in a maximal resolution of 5µm. X-ray transmission images are captured by a scintillator coupled via fibre optics to a 1024x1024 pixel 12-bit CCD. The object is placed between the x-ray source and the scintillator on a stub that rotates 360°around its vertical axis in steps as small as 0.24 degrees. Sample size is limited to 2 cm due to the absorption of geologic material for x-rays. The transmission images are back projected using a Feldkamp algorithm into a vertical stack of up to 1000 1Kx1K images that represent horizontal cuts of the object. This calculation takes 2 to several hours on a Double-Processor 2.4GHz PC. The stack of images (.bmp) can be visualized with any 3D-imaging software, used to produce cuts of Larger Foraminifera. Among other applications, the 3D-imaging software furnished by SkyScan can produce 3D-models by defining a threshold density value to distinguish "solid" from "void. Several models with variable threshold values and colors can be imbricated, rotated and cut together. The best results were obtained with microfossils devoid of chamber-filling cements (Permian, Eocene, Recent). However, even slight differences in cement mineralogy/composition can result in surprisingly good x-ray density contrasts.X-ray microtomography may develop into a powerful tool for larger microfossils with a complex internal structure, because it is non-destructive, requires no preparation of the specimens, and produces a true 3D-image data set. We will use these data sets in the future to produce cuts in any direction to compare them with arbitrary cuts of complex microfossils in thin sections. Many groups of benthic and planktonic foraminifera may become more easily determinable in thin section by this way.

Rapid prototyping of compliant human aortic roots for assessment of valved stents.

Adequate in-vitro training in valved stents deployment as well as testing of the latter devices requires compliant real-size models of the human aortic root. The casting methods utilized up to now are multi-step, time consuming and complicated. We pursued a goal of building a flexible 3D model in a single-step procedure. We created a precise 3D CAD model of a human aortic root using previously published anatomical and geometrical data and printed it using a novel rapid prototyping system developed by the Fab@Home project. As a material for 3D fabrication we used common house-hold silicone and afterwards dip-coated several models with dispersion silicone one or two times. To assess the production precision we compared the size of the final product with the CAD model. Compliance of the models was measured and compared with native porcine aortic root. Total fabrication time was 3 h and 20 min. Dip-coating one or two times with dispersion silicone if applied took one or two extra days, respectively. The error in dimensions of non-coated aortic root model compared to the CAD design was <3.0% along X, Y-axes and 4.1% along Z-axis. Compliance of a non-coated model as judged by the changes of radius values in the radial direction by 16.39% is significantly different (P<0.001) from native aortic tissue--23.54% at the pressure of 80-100 mmHg. Rapid prototyping of compliant, life-size anatomical models with the Fab@Home 3D printer is feasible--it is very quick compared to previous casting methods.

Coltop3D: A new software for structural analysis with high resolution 3D point clouds and DEM

Coltop3D is a software that performs structural analysis by using digital elevation model (DEM) and 3D point clouds acquired with terrestrial laser scanners. A color representation merging slope aspect and slope angle is used in order to obtain a unique code of color for each orientation of a local slope. Thus a continuous planar structure appears in a unique color. Several tools are included to create stereonets, to draw traces of discontinuities, or to compute automatically density stereonet. Examples are shown to demonstrate the efficiency of the method.

JULIDE: a software tool for 3D reconstruction and statistical analysis of autoradiographic mouse brain sections.

In this article we introduce JULIDE, a software toolkit developed to perform the 3D reconstruction, intensity normalization, volume standardization by 3D image registration and voxel-wise statistical analysis of autoradiographs of mouse brain sections. This software tool has been developed in the open-source ITK software framework and is freely available under a GPL license. The article presents the complete image processing chain from raw data acquisition to 3D statistical group analysis. Results of the group comparison in the context of a study on spatial learning are shown as an illustration of the data that can be obtained with this tool.

Helical tomotherapy (HT) for the treatment of anal canal cancer: preliminary clinical results, and dosimetric comparison between HT and intensity-modulated or 3D conformal radiotherapy

Background: To report a single-center experience in 19 patients (pts) with anal canal cancer treated with helical tomotherapy (HT) and concurrent chemotherapy, and compare the dosimetric results with fixed-field intensitymodulated radiotherapy (IMRT) and 3D conformal radiotherapy (3D RT). Materials and Methods: Between 2007 and 2008, 19 consecutive pts were treated with HT and concurrent CT for anal canal cancer. Median age was 59 years (range, 38−83), and female/male ratio was 14/5. The majority of the pts had T2 or T3 tumours (68.4%), and 52.6% had positive lymph nodes. In all 19 pts, pelvic and inguinal nodes, and tumour irradiation was given using HT upto a median dose of 36 Gy (1.8 Gy/fr) followed by a 1-week gap. A boost dose of 23.4 Gy (1.8 Gy/fr) was delivered to the tumour and involved nodes using 3DRT (n = 12), HT (n = 6), or IMRT (n = 1). Simultaneous integrated boost was used in none of the pts. All but one patient with a T1N0 tumour received concomitant mitomycin/5- fluorouracil (n = 12) or mitomycin/capecitabin (n = 7) CT. Toxicity was scored according to the Common Terminology Criteria for Adverse Events (NCICTCAE v3.0). HT plans and treatments were generated using Tomotherapy, Inc., software and hardware; and 3D or IMRT boost plans with the CMS treatment planning system (TPS), using 6−18 MV photons from a Siemens Primus accelerator. For dosimetric comparison, computed tomography data sets of 10 pts were imported into the TPS, and 3D and 5-field step-andshoot IMRT plans were generated for each case. Plans were optimized with the aim of assessing organs at risk (OAR) and healthy-tissue sparing while enforcing highly conformal target coverage, and evaluated by dose-volume histograms (DVH) of planning target volumes (PTV) and OAR. Results: With a median follow-up of 13 months (range, 3−18), all pts are alive and well; except one patient developing local recurrence at 12 months. No patient developed grade 3 or more acute toxicity. No unplanned treatment interruption was necessary because of toxicity. With 360-degree-of-freedom beam projection, HT showed an advantage over 3D or IMRT plans in terms of dose conformity around the PTV, and dose gradients were steeper outside the PTV, resulting in reduced doses to OARs. Using HT, acute toxicity was acceptable, and seemed to be better than historical standards. Conclusion: We conclude that HT combined with concurrent chemotherapy for anal canal cancer is effective and tolerable. Compared to 3DRT or 5-field IMRT, there is better conformity around the PTV, and OAR sparing.

Variability of ascending aorta diameter measurements as assessed with electrocardiography-gated multidetector computerized tomography and computer assisted diagnosis software.

Recently, morphometric measurements of the ascending aorta have been done with ECG-gated multidector computerized tomography (MDCT) to help the development of future novel transcatheter therapies (TCT); nevertheless, the variability of such measurements remains unknown. Thirty patients referred for ECG-gated CT thoracic angiography were evaluated. Continuous reformations of the ascending aorta, perpendicular to the centerline, were obtained automatically with a commercially available computer aided diagnosis (CAD). Then measurements of the maximal diameter were done with the CAD and manually by two observers (separately). Measurements were repeated one month later. The Bland-Altman method, Spearman coefficients, and a Wilcoxon signed-rank test were used to evaluate the variability, the correlation, and the differences between observers. The interobserver variability for maximal diameter between the two observers was up to 1.2 mm with limits of agreement [-1.5, +0.9] mm; whereas the intraobserver limits were [-1.2, +1.0] mm for the first observer and [-0.8, +0.8] mm for the second observer. The intraobserver CAD variability was 0.8 mm. The correlation was good between observers and the CAD (0.980-0.986); however, significant differences do exist (P<0.001). The maximum variability observed was 1.2 mm and should be considered in reports of measurements of the ascending aorta. The CAD is as reproducible as an experienced reader.

Evolution of a debris flow channel monitored using a 3D terrestrial laser scanner

Like numerous torrents in mountainous regions, the Illgraben creek (canton of Wallis, SW Switzerland) produces almost every year several debris flows. The total area of the active catchment is only 4.7 km², but large events ranging from 50'000 to 400'000 m³ are common (Zimmermann 2000). Consequently, the pathway of the main channel often changes suddenly. One single event can for instance fill the whole river bed and dig new several-meters-deep channels somewhere else (Bardou et al. 2003). The quantification of both, the rhythm and the magnitude of these changes, is very important to assess the variability of the bed's cross section and long profile. These parameters are indispensable for numerical modelling, as they should be considered as initial conditions. To monitor the channel evolution an Optech ILRIS 3D terrestrial laser scanner (LIDAR) was used. LIDAR permits to make a complete high precision 3D model of the channel and its surroundings by scanning it from different view points. The 3D data are treated and interpreted with the software Polyworks from Innovmetric Software Inc. Sequential 3D models allow for the determination of the variation in the bed's cross section and long profile. These data will afterwards be used to quantify the erosion and the deposition in the torrent reaches. To complete the chronological evolution of the landforms, precise digital terrain models, obtained by high resolution photogrammetry based on old aerial photographs, will be used. A 500 m long section of the Illgraben channel was scanned on 18th of August 2005 and on 7th of April 2006. These two data sets permit identifying the changes of the channel that occurred during the winter season. An upcoming scanning campaign in September 2006 will allow for the determination of the changes during this summer. Preliminary results show huge variations in the pathway of the Illgraben channel, as well as important vertical and lateral erosion of the river bed. Here we present the results of a river bank on the left (north-western) flank of the channel (Figure 1). For the August 2005 model the scans from 3 viewpoints were superposed, whereas the April 2006 3D image was obtained by combining 5 separate scans. The bank was eroded. The bank got eroded essentially on its left part (up to 6.3 m), where it is hit by the river and the debris flows (Figures 2 and 3). A debris cone has also formed (Figure 3), which suggests that a part of the bank erosion is due to shallow landslides. They probably occur when the river erosion creates an undercut slope. These geometrical data allow for the monitoring of the alluvial dynamics (i.e. aggradation and degradation) on different time scales and the influence of debris flows occurrence on these changes. Finally, the resistance against erosion of the bed's cross section and long profile will be analysed to assess the variability of these two key parameters. This information may then be used in debris flow simulation.

A 2D/3D image analysis system to track fluorescently labeled structures in rod-shaped cells: application to measure spindle pole asymmetry during mitosis.

BACKGROUND: The yeast Schizosaccharomyces pombe is frequently used as a model for studying the cell cycle. The cells are rod-shaped and divide by medial fission. The process of cell division, or cytokinesis, is controlled by a network of signaling proteins called the Septation Initiation Network (SIN); SIN proteins associate with the SPBs during nuclear division (mitosis). Some SIN proteins associate with both SPBs early in mitosis, and then display strongly asymmetric signal intensity at the SPBs in late mitosis, just before cytokinesis. This asymmetry is thought to be important for correct regulation of SIN signaling, and coordination of cytokinesis and mitosis. In order to study the dynamics of organelles or large protein complexes such as the spindle pole body (SPB), which have been labeled with a fluorescent protein tag in living cells, a number of the image analysis problems must be solved; the cell outline must be detected automatically, and the position and signal intensity associated with the structures of interest within the cell must be determined. RESULTS: We present a new 2D and 3D image analysis system that permits versatile and robust analysis of motile, fluorescently labeled structures in rod-shaped cells. We have designed an image analysis system that we have implemented as a user-friendly software package allowing the fast and robust image-analysis of large numbers of rod-shaped cells. We have developed new robust algorithms, which we combined with existing methodologies to facilitate fast and accurate analysis. Our software permits the detection and segmentation of rod-shaped cells in either static or dynamic (i.e. time lapse) multi-channel images. It enables tracking of two structures (for example SPBs) in two different image channels. For 2D or 3D static images, the locations of the structures are identified, and then intensity values are extracted together with several quantitative parameters, such as length, width, cell orientation, background fluorescence and the distance between the structures of interest. Furthermore, two kinds of kymographs of the tracked structures can be established, one representing the migration with respect to their relative position, the other representing their individual trajectories inside the cell. This software package, called "RodCellJ", allowed us to analyze a large number of S. pombe cells to understand the rules that govern SIN protein asymmetry. CONCLUSIONS: "RodCell" is freely available to the community as a package of several ImageJ plugins to simultaneously analyze the behavior of a large number of rod-shaped cells in an extensive manner. The integration of different image-processing techniques in a single package, as well as the development of novel algorithms does not only allow to speed up the analysis with respect to the usage of existing tools, but also accounts for higher accuracy. Its utility was demonstrated on both 2D and 3D static and dynamic images to study the septation initiation network of the yeast Schizosaccharomyces pombe. More generally, it can be used in any kind of biological context where fluorescent-protein labeled structures need to be analyzed in rod-shaped cells. AVAILABILITY: RodCellJ is freely available under http://bigwww.epfl.ch/algorithms.html, (after acceptance of the publication).

Development of a system for 3D high-resolution seismic reflection profiling on lakes

A high-resolution three-dimensional (3D) seismic reflection system for small-scale targets in lacustrine settings has been developed. Its main characteristics include navigation and shot-triggering software that fires the seismic source at regular distance intervals (max. error of 0.25 m) with real-time control on navigation using differential GPS (Global Positioning System). Receiver positions are accurately calculated (error < 0.20 m) with the aid of GPS antennas attached to the end of each of three 24-channel streamers. Two telescopic booms hold the streamers at a distance of 7.5 m from each other. With a receiver spacing of 2.5 m, the bin dimension is 1.25 m in inline and 3.75 m in crossline direction. To test the system, we conducted a 3D survey of about 1 km(2) in Lake Geneva, Switzerland, over a complex fault zone. A 5-m shot spacing resulted in a nominal fold of 6. A double-chamber bubble-cancelling 15/15 in(3) air gun (40-650 Hz) operated at 80 bars and 1 m depth gave a signal penetration of 300 m below water bottom and a best vertical resolution of 1.1 m. Processing followed a conventional scheme, but had to be adapted to the high sampling rates, and our unconventional navigation data needed conversion to industry standards. The high-quality data enabled us to construct maps of seismic horizons and fault surfaces in three dimensions. The system proves to be well adapted to investigate complex structures by providing non-aliased images of reflectors with dips up to 30 degrees.

Reproducibility of 3D free-breathing magnetic resonance coronary vessel wall imaging.

AIMS: Although the coronary artery vessel wall can be imaged non-invasively using magnetic resonance imaging (MRI), the in vivo reproducibility of wall thickness measures has not been previously investigated. Using a refined magnetization preparation scheme, we sought to assess the reproducibility of three-dimensional (3D) free-breathing black-blood coronary MRI in vivo. METHODS AND RESULTS: MRI vessel wall scans parallel to the right coronary artery (RCA) were obtained in 18 healthy individuals (age range 25-43, six women), with no known history of coronary artery disease, using a 3D dual-inversion navigator-gated black-blood spiral imaging sequence. Vessel wall scans were repeated 1 month later in eight subjects. The visible vessel wall segment and the wall thickness were quantitatively assessed using a semi-automatic tool and the intra-observer, inter-observer, and inter-scan reproducibilities were determined. The average imaged length of the RCA vessel wall was 44.5+/-7 mm. The average wall thickness was 1.6+/-0.2 mm. There was a highly significant intra-observer (r=0.97), inter-observer (r=0.94), and inter-scan (r=0.90) correlation for wall thickness (all P<0.001). There was also a significant agreement for intra-observer, inter-observer, and inter-scan measurements on Bland-Altman analysis. The intra-class correlation coefficients for intra-observer (r=0.97), inter-observer (r=0.92), and inter-scan (r=0.86) analyses were also excellent. CONCLUSION: The use of black-blood free-breathing 3D MRI in conjunction with semi-automated analysis software allows for reproducible measurements of right coronary arterial vessel-wall thickness. This technique may be well-suited for non-invasive longitudinal studies of coronary atherosclerosis.

3D virtual rendering before laparoscopic partial splenectomy in children.

INTRODUCTION: Partial splenectomy in children is a good surgical option for hematological diseases and focal splenic tumors because it allows the preservation of the spleen's immunological function. Furthermore, it can be performed by laparoscopy in children as it is a safe procedure, offering the benefits of a minimally invasive approach. MATERIALS AND METHODS: The software VR-render LE version 0.81 is a system that enables the visualization of bidimentional 3D images with magnification of anatomical details. We have applied this system to five cases of non-parasitic splenic cysts before laparoscopic partial splenectomy. RESULTS: The images obtained with VR rendering software permitted the preoperative reconstruction of the vascularization of the splenic hilum, allowing the surgeon safe vessel control during laparoscopic procedures. All five partial splenectomies were carried out with no complications or major blood loss. CONCLUSIONS: Laparoscopic partial splenectomy should be a first choice procedure because it is feasible, reproducible, and safe for children; furthermore, it preserves enough splenic tissue thereby preventing post-splenectomy infections. Volume rendering provides high anatomical resolution and can be useful in guiding the surgical procedure.

3D reconstruction and comparison of shapes of DNA minicircles observed by cryo-electron microscopy.

We use cryo-electron microscopy to compare 3D shapes of 158 bp long DNA minicircles that differ only in the sequence within an 18 bp block containing either a TATA box or a catabolite activator protein binding site. We present a sorting algorithm that correlates the reconstructed shapes and groups them into distinct categories. We conclude that the presence of the TATA box sequence, which is believed to be easily bent, does not significantly affect the observed shapes.

Comparison of 3D segmented gradient-echo and steady-state free precession coronary MRI sequences in patients with coronary artery disease.

OBJECTIVE: Our objective was to compare two state-of-the-art coronary MRI (CMRI) sequences with regard to image quality and diagnostic accuracy for the detection of coronary artery disease (CAD). SUBJECTS AND METHODS: Twenty patients with known CAD were examined with a navigator-gated and corrected free-breathing 3D segmented gradient-echo (turbo field-echo) CMRI sequence and a steady-state free precession sequence (balanced turbo field-echo). CMRI was performed in a transverse plane for the left coronary artery and a double-oblique plane for the right coronary artery system. Subjective image quality (1- to 4-point scale, with 1 indicating excellent quality) and objective image quality parameters were independently determined for both sequences. Sensitivity, specificity, and accuracy for the detection of significant (> or = 50% diameter) coronary artery stenoses were determined as defined in invasive catheter X-ray coronary angiography. RESULTS: Subjective image quality was superior for the balanced turbo field-echo approach (1.8 +/- 0.9 vs 2.3 +/- 1.0 for turbo field-echo; p < 0.001). Vessel sharpness, signal-to-noise ratio, and contrast-to-noise ratio were all superior for the balanced turbo field-echo approach (p < 0.01 for signal-to-noise ratio and contrast-to-noise ratio). Of the 103 segments, 18% of turbo field-echo segments and 9% of balanced turbo field-echo segments had to be excluded from disease evaluation because of insufficient image quality. Sensitivity, specificity, and accuracy for the detection of significant coronary artery stenoses in the evaluated segments were 92%, 67%, 85%, respectively, for turbo field-echo and 82%, 82%, 81%, respectively, for balanced turbo field-echo. CONCLUSION: Balanced turbo field-echo offers improved image quality with significantly fewer nondiagnostic segments when compared with turbo field-echo. For the detection of CAD, both sequences showed comparable accuracy for the visualized segments.