TY - JOUR
T1 - A novel constrained reconstruction model towards high-resolution sub-tomogram averaging.
AU - Han, Renmin
AU - Li, Lun
AU - Yang, Peng
AU - Zhang, Fa
AU - Gao, Xin
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): FCC/1/1976-18-01, FCC/1/1976-23-01, FCC/1/1976-25-01, FCC/1/1976-26-01
Acknowledgements: The research reported in this paper was supported by funding from the National Key Research and Development Program of China (No. 2017YFE0103900 and 2017YFA0504702), the King Abdullah University of Science and Technology (KAUST) Office of Sponsored Research (OSR) under Awards No. FCC/1/1976-18-01, FCC/1/1976-23-01, FCC/1/1976-25-01, FCC/1/1976-26-01 and FCS/1/4102-02-01, and FCS/1/4102-02-01, and the NSFC projects Grant (No. U1611263, U1611261, 61932018, and 61672493).
PY - 2019/10/17
Y1 - 2019/10/17
N2 - MOTIVATION:Electron tomography (ET) offers a unique capacity to image biological structures in situ. However, the resolution of ET reconstructed tomograms is not comparable to that of the single-particle cryo-EM. If many copies of the object of interest are present in the tomograms, their structures can be reconstructed in the tomogram, picked, aligned and averaged to increase the signal-to-noise ratio and improve the resolution, which is known as the subtomogram averaging (STA). To date, the resolution improvement of the subtomogram averaging is still limited because each reconstructed subtomogram is of low reconstruction quality due to the missing wedge issue.% issue in the tilt series images. RESULTS:In this paper, we propose a novel computational model, the constrained reconstruction model (CRM), to better recover the information from the multiple subtomograms and compensate for the missing wedge issue in each of them. CRM is supposed to produce a refined reconstruction in the final turn of subtomogram averaging after alignment, instead of directly taking the average. We first formulate the averaging method and our CRM as linear systems, and prove that the solution space of CRM is no larger, and in practice much smaller, than that of the averaging method. We then propose a sparse Kaczmarz algorithm to solve the formulated CRM, and further extend the solution to the simultaneous algebraic reconstruction technique (SART). Experimental results demonstrate that CRM can significantly alleviate the missing wedge issue and improve the final reconstruction quality. In addition, our model is robust to the number of images in each tilt series, the tilt range, and the noise level. AVAILABILITY:The codes of CRM-SIRT and CRM-SART are available at https://github.com/icthrm/CRM.
AB - MOTIVATION:Electron tomography (ET) offers a unique capacity to image biological structures in situ. However, the resolution of ET reconstructed tomograms is not comparable to that of the single-particle cryo-EM. If many copies of the object of interest are present in the tomograms, their structures can be reconstructed in the tomogram, picked, aligned and averaged to increase the signal-to-noise ratio and improve the resolution, which is known as the subtomogram averaging (STA). To date, the resolution improvement of the subtomogram averaging is still limited because each reconstructed subtomogram is of low reconstruction quality due to the missing wedge issue.% issue in the tilt series images. RESULTS:In this paper, we propose a novel computational model, the constrained reconstruction model (CRM), to better recover the information from the multiple subtomograms and compensate for the missing wedge issue in each of them. CRM is supposed to produce a refined reconstruction in the final turn of subtomogram averaging after alignment, instead of directly taking the average. We first formulate the averaging method and our CRM as linear systems, and prove that the solution space of CRM is no larger, and in practice much smaller, than that of the averaging method. We then propose a sparse Kaczmarz algorithm to solve the formulated CRM, and further extend the solution to the simultaneous algebraic reconstruction technique (SART). Experimental results demonstrate that CRM can significantly alleviate the missing wedge issue and improve the final reconstruction quality. In addition, our model is robust to the number of images in each tilt series, the tilt range, and the noise level. AVAILABILITY:The codes of CRM-SIRT and CRM-SART are available at https://github.com/icthrm/CRM.
UR - http://hdl.handle.net/10754/659965
UR - https://academic.oup.com/bioinformatics/advance-article/doi/10.1093/bioinformatics/btz787/5588412
U2 - 10.1093/bioinformatics/btz787
DO - 10.1093/bioinformatics/btz787
M3 - Article
C2 - 31617571
SN - 1367-4803
JO - Bioinformatics (Oxford, England)
JF - Bioinformatics (Oxford, England)
ER -