GPU机群系统上加速EMAN(英文)

近年来,高质量的三维重构对计算机系统的计算能力提出越来越高的需求。EMAN是对冷冻电子显微镜拍摄出来的图像进行单颗粒三维重构的流行软件包。由于GPU显示出优于CPU的性能水平,在GPU机群系统上对EMAN进行加速将有可能取得较好的效果。首先,有必要调查EMAN在计算方面的特点,包括计算密度,访存情况和并行潜力。结果显示其在GPU体系结构上有获得较高性能的潜力。然后,根据此特点完成一款并行EMAN软件,使用CUDA和MPI对EMAN实施不同粒度的并行处理。在一个拥有16颗Intel6核处理器的机群系统上,用8颗FermiGPU加速后的EMAN(CUDA-EMAN)相比于仅使用CPU的并行EMAN程序获得了2.5倍的加速比,也就是1颗FermiGPU相比一个Intel6核处理器可对EMAN加速5倍。     

封面故事

<正>电子显微镜三维重构技术近年来有了飞速的发展,已逐渐成为研究生物大分子复合体三维结构的重要手段。电子晶体学、单颗粒三维重构技术和电子断层技术是电子显微     

图形处理器在实时光学相干断层成像中的应用

光学相干断层成像(optical coherence tomography,OCT)技术在成像过程中具有极大的数据量和计算量,传统的基于中央处理器(central processing unit,CPU)的计算平台难以满足OCT实时成像的需求。图形处理器(graphics processing unit,GPU)在通用计算方面具有强大的并行处理能力和数值计算能力,可以突破OCT实时成像的瓶颈。本文对GPU做了简要介绍并阐述了GPU在OCT实时成像及功能成像中的应用及研究进展。     

生物三维电子显微学进入全新高速发展时期

生物三维电子显微学主要由三个部分组成——电子晶体学、单颗粒技术和电子断层成像术,其结构解析对象的尺度范围介于x射线晶体学与光学显微镜之间,适合从蛋白质分子结构到细胞和组织结构的解析。以冷冻电镜技术与三维重构技术为基础的低温电子显微学代表了生物电子显微学的前沿。低温单颗粒技术对于高度对称的病毒颗粒的解析最近已达到3．8A分辨率,正在成为解析分子量很大的蛋白质复合体高分辨结构的有效技术手段。低温电子断层成像技术目前对于真核细胞样品的结构解析已达到约40A的分辨率,在今后5年有望达到20A。这样,把x射线晶体学、NMR以及电镜三维重构获得的蛋白质分子及复合体的高分辨率的结构,锚定到较低分辨率的电子断层成像图像中,从而在细胞水平上获得高精确的蛋白质空间定位和原子分辨率的蛋白质相互作用的结构信息。这将成为把分子水平的结构研究与细胞水平的生命活动衔接起来的可行途径。     

冷冻电子断层成像技术及其在生物研究领域的应用

冷冻电子断层成像可以在纳米级尺度上研究那些结构不具有均一性的分子、病毒、细胞器以及它们之间组成的复合体的三维结构。在过去的十年中,电子显微镜硬件、冷冻制样设备和技术,以及自动化断层数据收集方法的进步使得本研究领域得到快速发展。本文对冷冻电子断层成像的方法,包括基本原理、样品制备、断层数据采集和图像处理、三维重构以及重建信息的理解和展示、近年来在生物样品领域的一些典型应用以及前景作一简单介绍。     

始于努力,续于真诚,臻于共赢

2010年第7期,<学报>为在生物物理研究所举办的电镜三维重构学习班出版了一期专刊,其中除了两篇与电镜三维重构相关的研究论文之外,在朝仡夕拾栏目发表了北京大学医学部尹长城老师撰写的一篇介绍运用电子显微镜研究生物人分子起源及进展的文章,另外四篇综述分别讨论了电镜三维重构、单颗粒电镜、冷冻电镜断层成像技术及三维电镜自动化...     

头部组织三维核磁共振电阻抗成像算法的仿真研究

文章给出了一种基于核磁共振技术的三维阻抗成像（电导率分布）重构算法,并将该方法应用于人体头部组织电导率分布重构上。该代数重构方法是利用高分辨率的核磁共振成像系统对成像物体进行三维构建和不同组织的边界区分,根据核磁共振系统中测量得到的磁感应强度Bx和By分量并结合有限元数值计算得到的电流密度分布J组成非线性矩阵,通过迭代求解此非线性矩阵,来解决三维电导率分布的重构问题。在三层球头模型（包括头皮、颅骨和大脑）上分别进行的仿真实验结果表明,该算法具有较强的抗噪声能力和较好的收敛性,重构的头部电导率分布图像具有较高的精确性。     

电子断层成像技术及其在细胞生物学中的应用

电子断层成像技术（electrontomography）是近年来发展起来一项三维成像技术,可以在纳米分辨率（2-10nm）水平上获得生物大分子及其复合物或聚集体、细胞器、细胞以及组织的三维结构,而且可以用于研究生物大分子在细胞中的定位、排列、分布以及相互作用,已逐渐成为细胞生物学领域中的一项重要技术手段。该文针对这项技术及其在细胞生物学中的应用作一简要介绍。     

使用电子断层技术对IgG1抗体逐个分子的三维重构与结构分析

抗体(antibody)又称免疫球蛋白(immunoglobulin,Ig),是人体免疫反应的重要参与者.了解抗体的结构和结构动态特征,是理解人体免疫作用机理、修复或提高免疫能力、定向设计抗体以治疗各种疾病的基础.本文以人体IgG1抗体为对象,综述了使用透射电子显微学方法研究IgG1抗体结构方向的最新进展.详细介绍了使用逐个分子的电子断层三维重构技术(individual-particle electron tomography,IPET)对抗体进行结构研究的方法,包括样品制备、图像处理和数据分析等.并描述了利用该技术,在研究抗体结合肽分子后的结构形变和通过收集不同构象来研究抗体动态结构特征方面所取得的阶段性成果.最后,对尚待解决的关键问题与该技术未来的发展方向进行了讨论与展望.     

光学透明技术在植物多尺度成像中的应用

光学透明技术是一种通过各种化学试剂, 将原本不透明的生物样本实现透明化, 并在光学显微镜下深度成像的技术。结合多种光学显微成像新技术, 光学透明技术可对整个组织进行成像和三维重建, 深度剖析生物体内部空间特征与形成机制。近年来, 多种植物光学透明技术和多尺度成像技术被陆续研发, 并取得了丰硕的研究成果。该文综述了生物体光学透明技术的基本原理和一些新技术, 重点介绍基于光学透明技术开发的新型成像方法及其在植物成像与细胞生物学中的应用, 为后续植物整体、组织或器官的透明、成像与三维重构及功能研究提供理论依据和技术支持。     

Compensation of Missing Wedge Effects with Sequential Statistical Reconstruction in Electron Tomography

Electron tomography (ET) of biological samples is used to study the organization and the structure of the whole cell and subcellular complexes in great detail. However, projections cannot be acquired over full tilt angle range with biological samples in electron microscopy. ET image reconstruction can be considered an ill-posed problem because of this missing information. This results in artifacts, seen as the loss of three-dimensional (3D) resolution in the reconstructed images. The goal of this study was to achieve isotropic resolution with a statistical reconstruction method, sequential maximum a posteriori expectation maximization (sMAP-EM), using no prior morphological knowledge about the specimen. The missing wedge effects on sMAP-EM were examined with a synthetic cell phantom to assess the effects of noise. An experimental dataset of a multivesicular body was evaluated with a number of gold particles. An ellipsoid fitting based method was developed to realize the quantitative measures elongation and contrast in an automated, objective, and reliable way. The method statistically evaluates the sub-volumes containing gold particles randomly located in various parts of the whole volume, thus giving information about the robustness of the volume reconstruction. The quantitative results were also compared with reconstructions made with widely-used weighted backprojection and simultaneous iterative reconstruction technique methods. The results showed that the proposed sMAP-EM method significantly suppresses the effects of the missing information producing isotropic resolution. Furthermore, this method improves the contrast ratio, enhancing the applicability of further automatic and semi-automatic analysis. These improvements in ET reconstruction by sMAP-EM enable analysis of subcellular structures with higher three-dimensional resolution and contrast than conventional methods.     

Three-dimensional reconstruction using an adaptive simultaneous algebraic reconstruction technique in electron tomography

Three-dimensional (3D) reconstruction of electron tomography (ET) has emerged as an important technique in analyzing structures of complex biological samples. However most of existing reconstruction methods are not suitable for extremely noisy and incomplete data conditions. We present an adaptive simultaneous algebraic reconstruction technique (ASART) in which a modified multilevel access scheme and an adaptive relaxation parameter adjustment method are developed to improve the quality of the reconstructed 3D structure. The reconstruction process is facilitated by using a column-sum substitution approach. This modified multilevel access scheme is adopted to arrange the order of projections so as to minimize the correlations between consecutive views within a limited angle range. In the adaptive relaxation parameter adjustment method, not only the weight matrix (as in the existing methods) but the gray levels of the pixels are employed to adjust the relaxation parameters so that the quality of the reconstruction is improved while the convergence process of the reconstruction is accelerated. In the column-sum substitution approach, the computation to obtain the reciprocal of the sum for the columns in each view is avoided so that the needed computations for each iteration can be reduced. Experimental results show that the proposed technique ASART is better based on objective quality measures than other methods, especially when data is noisy and limited in tilt angles. At the same time, the reconstruction by ASART outperforms that of simultaneous algebraic reconstruction technique (SART) in speed.     

UCSF tomography: an integrated software suite for real-time electron microscopic tomographic data collection, alignment, and reconstruction

A real-time alignment and reconstruction scheme for electron microscopic tomography (EMT) has been developed and integrated within our UCSF tomography data collection software. This newly integrated software suite provides full automation from data collection to real-time reconstruction by which the three-dimensional (3D) reconstructed volume is immediately made available at the end of each data collection. Real-time reconstruction is achieved by calculating a weighted back-projection on a small Linux cluster (five dual-processor compute nodes) concurrently with the UCSF tomography data collection running on the microscope's computer, and using the fiducial-marker free alignment data generated during the data collection process. The real-time reconstructed 3D volume provides users with immediate feedback to fully asses all aspects of the experiment ranging from sample choice, ice thickness, experimental parameters to the quality of specimen preparation. This information can be used to guide subsequent data collections. Access to the reconstruction is especially useful in low-dose cryo EMT where such information is very difficult to obtain due to extraordinary low signal to noise ratio in each 2D image. In our environment, we generally collect 2048 x 2048 pixel images which are subsequently computationally binned four-fold for the on-line reconstruction. Based upon experiments performed with thick and cryo specimens at various CCD magnifications (50000x-80000x), alignment accuracy is sufficient to support this reduced resolution but should be refined before calculating a full resolution reconstruction. The reduced resolution has proven to be quite adequate to assess sample quality, or to screen for the best data set for full-resolution reconstruction, significantly improving both productivity and efficiency of system resources. The total time from start of data collection to a final reconstructed volume (512 x 512 x 256 pixels) is about 50 min for a +/-70 degrees 2k x 2k pixel tilt series acquired at every 1 degrees.     

Altered role of smooth muscle endothelin receptors in coronary endothelin-1 and alpha1-adrenoceptor-mediated vasoconstriction in Type 2 diabetes

Regulation of vascular tone and blood flow involves interactions between numerous local and systemic vascular control signals, many of which are altered by Type 2 diabetes (T2D). Vascular responses to endothelin-1 (ET-1) are mediated by endothelin type A (ET(A)) and type B (ET(B)) receptors that have been implicated in cross talk with alpha(1)-adrenoceptors (alpha(1)-AR). ET(A) and ET(B) receptor expression and plasma ET-1 levels are elevated in T2D; however, whether this influences coronary alpha(1)-AR function has not been examined. Therefore, we examined the effect of ET(A) and ET(B) receptor inhibition on coronary vasoconstriction to ET-1 and alpha(1)-AR activation in a mouse model of T2D. Coronary vascular responses were examined in isolated mouse hearts from control and diet-induced T2D C57BL/6J mice. Responses to ET-1 and the selective alpha(1)-AR agonist phenylephrine (PE) were examined alone and in the presence of the nitric oxide synthase inhibitor N(omega)-nitro-l-arginine methyl ester (l-NAME) alone or in combination with selective ET(A) or ET(B) receptor inhibitors BQ-123 and BQ-788, respectively. Vasoconstriction to ET-1 was enhanced, whereas ET(B), but not ET(A), receptor blockade reduced basal coronary tone in T2D hearts. In the presence of l-NAME, ET(A) receptor inhibition attenuated ET-1 vasoconstriction in both groups, whereas ET(B) inhibition abolished this response only in control hearts. In addition, ET(A) inhibition enhanced alpha(1)-AR-mediated vasoconstriction in T2D, but not control, hearts following l-NAME treatment. Therefore, in this model, enhanced coronary ET-1 responsiveness is mediated primarily through smooth muscle ET(B) receptors, whereas the interaction with alpha(1)-ARs is mediated solely through the ET(A) receptor subtype.     

Evaluation of a commercial Model Based Iterative reconstruction algorithm in computed tomography

IntroductionIterative reconstruction algorithms have been introduced in clinical practice to obtain dose reduction without compromising the diagnostic performance.PurposeTo investigate the commercial Model Based IMR algorithm by means of patient dose and image quality, with standard Fourier and alternative metrics.Materials and methodsA Catphan phantom, a commercial density phantom and a cylindrical water filled phantom were scanned both varying CTDI_vol and reconstruction thickness. Images were then reconstructed with Filtered Back Projection and both statistical (iDose) and Model Based (IMR) Iterative reconstruction algorithms.Spatial resolution was evaluated with Modulation Transfer Function and Target Transfer Function. Noise reduction was investigated with Standard Deviation. Furthermore, its behaviour was analysed with 3D and 2D Noise Power Spectrum. Blur and Low Contrast Detectability were investigated.Patient dose indexes were collected and analysed.ResultsAll results, related to image quality, have been compared to FBP standard reconstructions.Model Based IMR significantly improves Modulation Transfer Function with an increase between 12% and 64%. Target Transfer Function curves confirm this trend for high density objects, while Blur presents a sharpness reduction for low density details.Model Based IMR underlines a noise reduction between 44% and 66% and a variation in noise power spectrum behaviour. Low Contrast Detectability curves underline an averaged improvement of 35–45%; these results are compatible with an achievable reduction of 50% of CTDI_vol.A dose reduction between 25% and 35% is confirmed by median values of CTDI_vol.ConclusionIMR produces an improvement in image quality and dose reduction.     

TOM software toolbox: acquisition and analysis for electron tomography

Automated data acquisition procedures have changed the perspectives of electron tomography (ET) in a profound manner. Elaborate data acquisition schemes with autotuning functions minimize exposure of the specimen to the electron beam and sophisticated image analysis routines retrieve a maximum of information from noisy data sets. "TOM software toolbox" integrates established algorithms and new concepts tailored to the special needs of low dose ET. It provides a user-friendly unified platform for all processing steps: acquisition, alignment, reconstruction, and analysis. Designed as a collection of computational procedures it is a complete software solution within a highly flexible framework. TOM represents a new way of working with the electron microscope and can serve as the basis for future high-throughput applications.     

Iterative reconstruction of cryo-electron tomograms using nonuniform fast Fourier transforms

Algorithms for three-dimensional (3D) reconstruction of objects based on their projections are essential in various biological and medical imaging modalities. In cryo-electron tomography (CET) a major challenge for reconstruction is the limited range of projection angles, which manifests itself as a “missing wedge” of data in Fourier space making the reconstruction problem ill-posed. Here, we apply an iterative reconstruction method that makes use of nonuniform fast Fourier transform (NUFFT) to the reconstruction of cryo-electron tomograms. According to several measures the reconstructions are superior to those obtained using conventional methods, most notably weighted backprojection. Most importantly, we show that it is possible to fill in partially the unsampled region in Fourier space with meaningful information without making assumptions about the data or applying prior knowledge. As a consequence, particles of known structure can be localized with higher confidence in cryotomograms and subtomogram averaging yields higher resolution densities.     

Development of a 3D optical scanner for evaluating patient-specific dose distributions

PurposeThis paper describes the hardware and software characteristics of a 3D optical scanner (P3DS) developed in-house. The P3DS consists of an LED light source, diffuse screen, step motor, CCD camera, and scanner management software with 3D reconstructed software.Materials and methodWe performed optical simulation, 2D and 3D reconstruction image testing, and pre-clinical testing for the P3DS. We developed the optical scanner with three key characteristics in mind. First, we developed a continuous scanning method to expand possible clinical applications. Second, we manufactured a collimator to improve image quality by reducing scattering from the light source. Third, we developed an optical scanner with changeable camera positioning to enable acquisition of optimal images according to the size of the gel dosimeter.ResultsWe confirmed ray-tracing in P3DS with optic simulation and found that 2D projection and 3D reconstructed images were qualitatively similar to the phantom images. For pre-clinical tests, the dose distribution and profile showed good agreement among RTP, optical CT, and external beam radiotherapy film data for the axial and coronal views. The P3DS has shown that it can scan and reconstruct for evaluation of the gel dosimeter within 1 min. We confirmed that the P3DS system is a useful tool for the measurement of 3D dose distributions for 3D radiation therapy QA. Further experiments are needed to investigate quantitative analysis for 3D dose distribution.     

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.