- 1Mathematical Intelligence Application LAB, Institute for Mathematical Sciences, Renmin University of China, Beijing, China
- 2Beijing Academy of Intelligence, Beijing, China
Cryo-electron tomography (cryo-ET) plays a critical role in imaging microorganisms in situ in terms of further analyzing the working mechanisms of viruses and drug exploitation, among others. A data processing workflow for cryo-ET has been developed to reconstruct three-dimensional density maps and further build atomic models from a tilt series of two-dimensional projections. Low signal-to-noise ratio (SNR) and missing wedge are two major factors that make the reconstruction procedure challenging. Because only few near-atomic resolution structures have been reconstructed in cryo-ET, there is still much room to design new approaches to improve universal reconstruction resolutions. This review summarizes classical mathematical models and deep learning methods among general reconstruction steps. Moreover, we also discuss current limitations and prospects. This review can provide software and methods for each step of the entire procedure from tilt series by cryo-ET to 3D atomic structures. In addition, it can also help more experts in various fields comprehend a recent research trend in cryo-ET. Furthermore, we hope that more researchers can collaborate in developing computational methods and mathematical models for high-resolution three-dimensional structures from cryo-ET datasets.
1 Introduction
Cryo-electron tomography (cryo-ET) is an important imaging technique that provides high-resolution three-dimensional (3D) structures for biological specimens in situ. Observing the structure of biological specimens in situ can reflect the location of spatial interaction, and further help to analyze the mechanism and function of macromolecules. Single-particle analysis (SPA) in cryogenic electron microscopy (cryo-EM) is a mature approach to obtaining high-resolution 3D reconstructed structures, even imaging nearly native biomolecules in cells or organelles in situ. Indeed, the in situ single-particle analysis (isSPA) (Cheng et al., 2021) method is a potential competitor of cryo-ET. However, a complex system of biological macromolecules is still challenging to construct. Cryo-ET can play a key bridging role in forming high-resolution cellular landscapes in situ among multiple scales. Cryo-ET has been increasingly favored by structural biologists in the past few years. There are several advantages of cryo-ET. Firstly, some artificial alterations or modifications are averted in sample preparation, like particle isolation and purification. The native state can be cryo-preserved by vitrification (Plitzko and Baumeister, 2019). Secondly, cryo-ET allows the sample thickness of biological specimens below 500 nm (Khanna and Villa, 2022; Zuber and Lučić, 2022), such as viruses and bacterial cells. However, a thinner sample means higher resolution recorded in some sense. Thirdly, cryo-ET can analyze heterogeneous (Xu et al., 2021) and pleomorphic (Borgnia et al., 2008; Obr and Schur, 2019) structures like cells, viruses, and tissues (Borgnia et al., 2008). Different from cryo-EM, cryo-ET does not require many copies of the same specimen. Here, we focus on cryo-ET in microorganisms.
The high-resolution structures of microorganisms in cryo-ET can provide reliable 3D structures for drug development, target discovery, vaccine research, and so on. They can promote the development of science and technology. There are some actual examples of microorganisms using cryo-ET, like Bdellovibrio bacteriovorus (Borgnia et al., 2008), Mycoplasma genitalium (Seybert et al., 2018), Candida glabrata (Jiménez-Ortigosa et al., 2021), Mycoplasma pneumoniae (Xue et al., 2022), and Escherichia coli (Chen et al., 2022). Recently, cryo-ET has resolved a giant virus (Fang et al., 2019) at the near-atomic level, namely, a 3.5-Å paramecium bursaria chlorella virus 1 (PBCV-1). The reorganization and host–pathogen interaction (Nans et al., 2014; López-Jiménez and Mostowy, 2021) are also resolved in cryo-ET. The structures further help to analyze the action mechanism of infected cells, which helps researchers to discover anti-infective strategies. Some viruses are also resolved in cryo-ET like severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Klein et al., 2020), human immunodeficiency virus (HIV) (Mendonҫa et al., 2021), and a bullet-shaped vesicular stomatitis virus (Jenni et al., 2022). As coronavirus disease 2019 (COVID-19) has caused great damage and loss to humans, the analysis of SARS-CoV-2 has become an important application of cryo-ET. The 3D structures of SARS-CoV-2 (Turoňová et al., 2020; Calder et al., 2022) by cryo-ET have helped researchers to explore the functional mechanism and promote the development of vaccines and drugs. Critical SARS-CoV-2 structures are resolved in Zhang et al. (2021), such as the virus egress pathway and native virus spike structures. The spike protein plays an important role in virus–host interaction, membrane fusion, and so on. The spike protein is a main target for vaccine development and antigenic analysis (Li, 2022). The structural analysis for alphacoronavirus spike protein in situ revealed the motion of D0 domain (Huang et al., 2022). However, there are still a large number of structures that have not been resolved. It is significant and urgent to develop methods for each step of reconstruction processing from tilt series in cryo-ET to restore high-resolution structures.
There are two main limiting factors for 3D structure reconstruction steps. On the one hand, a signal-to-noise ratio (SNR) (Bepler et al., 2020) is the ratio of signal power to noise power, which is often expressed in decibels (dB). Higher SNR means that the signal is stronger and the noise is weaker; that is, the quality of the image is better. On the other hand, the tilt angles are generally from −70° to 70° in cryo-ET. Thus, the missing angles lead to a lack of information in two continuous, symmetrical wedge-shaped areas in Fourier space (Liu et al., 2022b), which is referred to the missing wedge of data. For convenience, we use “missing wedge” to mean “missing wedge of data” in the rest of the paper. Low SNR and missing wedge are two factors that make the reconstruction procedure difficult. In addition, inter-structure heterogeneity can also negatively affect the reconstruction procedure. Structural heterogeneity refers to the fact that macromolecules have many different discrete or continuous conformations. Heterogeneity can affect the steps of 3D particle picking, classification of 3D particles, and so on.
Data websites and software for cryo-ET have been developed in the past few years. The Electron Microscopy Data Bank (EMDB) (Lawson et al., 2016) and Electron Microscopy Public Image Archive (EMPIAR) (Iudin et al., 2016) are two popular databases. The EMPIAR provides tilt series projections, and EMDB provides 3D density maps. Subtomogram averaging (STA) is a technique of averaging 3D biological samples within tomograms (Bharat and Scheres, 2016b). The STA is a popular image restoration technique in cryo-ET to obtain the 3D structure of biological samples in their native environments (Scaramuzza and Castaño-Díez, 2021). The STA technique is increasingly important in structural biology. The concepts of each reconstruction steps will be described in specific sections. To demonstrate the current ability of cryo-ET to resolve structures, we show the resolution distribution and software rank by the STA from the EMDB before 28 December 2022, as Figure 1 shows. The resolution of most biological samples analyzed by cryo-ET is low, and only a tiny part can reach more than 10 Å resolution, which is far from enough for the study of the internal mechanism of organisms. The top 11 software/methods are IMOD (Mastronarde and Held, 2017), Relion (Zivanov et al., 2022), Particle Estimation for Electron Tomography (PEET) (Nicastro et al., 2006), Dynamo (Scaramuzza and Castaño-Díez, 2021), Ctffind (Rohou and Grigorieff, 2015), Novactf (Turonova et al., 2017), PyTom (Hrabe et al., 2012), Warp (Tegunov and Cramer, 2019), Av3 (Förster et al., 2005), EMAN2 (Murray et al., 2014), and GCTF (Zhang, 2016). The ScipionTomo (Morena et al., 2022) provides a platform containing most software for users applying different software to solve the same task. Recently, with the rapid popularity of artificial intelligence algorithms, many new computational methods have emerged to improve the reconstruction resolution in cryo-ET. This review summarizes some new computational methods and classical mathematical models for 3D structure reconstruction from cryo-ET data. We hope that this review can attract more researchers to collaborate in developing reconstruction methods in cryo-ET and provide guidelines for other researchers.
Figure 1 Left: the numbers of 3D density maps as a function of resolutions by the subtomogram averaging (STA) in cryo-ET. Here, we only capture numbers of the resolution from 1 to 51 Å. Right: the numbers of times the software is used by the STA. Specific data are before 28 December 2022 from the EMDB website: https://www.ebi.ac.uk/emdb/.
Further improvement in cryo-ET data processing is still significant and urgent for more precise analysis and higher throughput. Here are some examples of the current problems. The limitation of the number of fiducial markers may lead to over-fitting in motion correction and tilt series alignment (Fernandez and Li, 2021). The reason for over-fitting is that the data from marked regions would be less than the imaging model’s parameters in proportion. Segmentation and particle picking also rely, to some extent, on manual labeling based on prior knowledge (Pyle and Zanetti, 2021). Image denoising can help particle picking and be the initial steps of an STA analysis. In later steps, since the coordinates are known, particles are extracted and used from noisy data. Furthermore, denoised data can be used to produce an initial model. Full-resolution data can be used for structural analysis to reach the highest resolution as possible (Taylor et al., 2006). While the experimental data in cryo-ET are progressively enriched, computational methods ought to leverage such abundant information for higher-resolution restoration and more accurate analysis.
This review summarizes current methods of artificial intelligence and mathematics for reconstructing 3D density maps and atomic models from two-dimensional (2D) tilt-series projections in cryo-ET. Here, we draw a general flowchart to show about 10 steps of 3D reconstruction in cryo-ET in Figure 2. It is worth noting that the process is subjected to appropriate adjustments and changes owing to differences in biological samples. The three steps in the dashed box in Figure 2, namely, motion correction, contrast transfer function (CTF) correction, and tilt-series alignment, are sometimes chosen to group into a program (Fernandez and Li, 2021). A denoiser is used to remove noise. While removing noise, it also removes the high-resolution signal. Thus, a denoiser is just used for some initial steps. The 2D denoiser is just an option to remove noise in 2D tilt-series projections, and 3D denoisers can remove noise in 3D subtomograms. The 3D denoiser can help improve the accuracy of 3D particle picking and classification. For the same classification, the subtomogram alignment and averaging (STAA) restores the 3D volumes by averaging to remove noise. Post-processing, also called refinement, corrects the 3D density map. Finally, the last step is building an atomic model from a 3D density map. In addition, denoising can be used for particle picking, initial averaging, and alignments, and classification. For higher-resolution information, the full-resolution volumes should be used for STA, alignments, and 3D density map refinement. Moreover, for real high resolution in STA analysis, sub-tilt optimization will take place. The methods for each step are shown in detail in the following sections, arranged as follows.
Figure 2 The general flowchart of the 3D reconstruction from tilt-series projections by cryo-ET. From top to bottom, the first frame shows the imaging process by cryo-ET. Then, the tilt series are the inputs of the second frame to reconstruct a 3D structure. The three parts in the dashed box are sometimes handled together. The last frame displays some applications to the microorganism.
1.1 Organization of this paper
The 3D reconstruction in cryo-ET roughly contains 10 steps. Sections 2 to Section 11 show motion correction, correction of CTF, tilt series alignment, 2D and 3D denoisers, 3D reconstruction, 3D particle picking, classification, subtomograms alignment and averaging, post-processing for 3D density map, and atomic structure building in turn, respectively. We discuss the difficulties of each small task, the methods, and the prospects for the future. Section 12 discusses cryo-ET reconstruction research’s current status, difficulties, and future research prospects.
2 Motion correction
Motion happens because the samples are exposed to the electron beam, which will cause unsatisfactory blurs and high-resolution informational losses. This motion is called beam-induced motion. The electronic beam is produced through the acceleration of electrons at high voltage (approximately 300 kV). In cryo-electron tomography, the beam-induced motion is usually attributed to multiple factors, including electron beam interaction with the specimen, stage vibration, sample deformation, intense ionization in materials, and measurement error (Naydenova et al., 2020; Thorne, 2020). The local sample motion induced by the electron beam can be observed as a “doming” motion (Ganguly et al., 2022), where the sample suffers an upward deformation perpendicular to the sample plane. The intermolecular correlation shifts and the measurement error of the projections should be considered. Therefore, these motions must be corrected to reduce their impacts on low SNR and artifacts.
In cryo-electron tomography, a frame is a single image recorded by the direct-electron detectors. A movie is a sequence of frames during the electron beam exposure of the biological specimens (Li et al., 2013). Simulating the beam-induced movement and correcting the movie/frame data are an inverse problems in a mathematical sense. As the complex electron–sample interactions presented above, some mathematical models are developed based on the anisotropies. The two-dimensional or polynomial spline interpolation and stochastic processes are capable of fitting. In TomoAlign (Fernandez et al., 2018; Fernandez et al., 2019), the algebraic equation includes the 3D deformation of the sample as the motion correction through the different tilts as homogeneous polynomial vectors in the 3D-axis of a microscope, , where rj= (xj,yj,zj) are the jth fiducial marker’s coordinates in the axis of the sample stage, and i is the index of the image of the tilt series. In RELION 4.0 (Zivanov et al., 2022), the stochastic process estimation of the electron beam-induced motions is derived within Gaussian processes, which is an analogous method in RELION 3.0. Specifically, the motion of the interest particle is modeled on the path through the tilt series in 3-dimensional space. The optimization of parameters comprises the maximum likelihood of the posterior estimation determining the motion pattern and the irradiate damage accumulation. With the assumption that the sparse fiducial movements locally estimate the whole motion of the sample, it follows the limitation of the dependence of fiducial markers. The number of markers may contribute to an over-fitted approximation of the parameters introduced by polynomials.
For the particles of growing interest in marker-free alignments, such as cryo-lamella, Aretomo (Zheng et al., 2022) provides an integrated marker-free solution in a fully acquired tilt series for correcting the beam-induced motion. Aretomo regards the drift paths of feature in patches as the beam-induced local deformation throughout the whole tilt series. The rotation of the sample stage in cryo-electron microscopy can essentially be described by 3D rigid-body motion. By matching of this motion with the features, residuals of the original and matched samples represent the local deformation. Informed about the measured inter-subframe shifts, MotionCor2 (Zheng et al., 2017) describes these non-coherent and idiosyncratic shifts in the projection plane as smooth two-dimensional polynomials, gives least-squared estimates of correlated shifts between adjacent subframes, and maps backward to the damage-weighted image. The patches are artificially divided acquired images and the features are respectively extracted from them. MotionCor2 v1.5.0 is used to match the centers of these patches in the tilt series and detect features for local motion correction. Warp (Tegunov and Cramer, 2019) takes a similar process strategy as MotionCor2 but with no assumption on the 2D polynomials, where two grids decompose the global and local motion sources and perform motion correction separately. The dose-fraction frames are aligned and corrected in Fourier space. Instead of manually or semi-manually labeling marker locations, SparseAlign (Ganguly et al., 2022) is a mathematical method for marker localization and deformation estimation simultaneously. The SparseAlign solves a few deformation parameters based on an image-based loss between the forward projection of markers and the observed projections. Software such as TomoAlign (Fernandez and Li, 2021), Aretomo (Zheng et al., 2022), MotionCor2 (Zheng et al., 2017), Ximpp (Strelak et al., 2021), Warp (Tegunov and Cramer, 2019), RELION (Zivanov et al., 2022), and Unblur (Grant and Grigorieff, 2015) all provide solutions for motion correction.
Motion correction is indeed a vital step in ensuring the successful analysis of cryo-ET data. Most methods of motion correction in cryo-ET are based on those in cryo-EM. Consequently, many of these methods may not be specifically optimized for cryo-ET. In contrast, the different sample preparation, electron dose, and relative issues are supposed to be in consideration of image processing. It is a valuable recognition that algorithms correct beam-induced motion in the electron microscope or specimen stage rather than in the projection plane. More accurate tilt series alignments, tomography reconstruction, and subtomogram analysis will be carried out with considerable motion correction. However, the calculating limitations should be taken into account for the experimental approaches of raw image acquisition, such as the number of fiducial markers or intrinsic features and more accurate but omitted factors of beam-induced deformation beyond doming motion. Sometimes, the performance of motion correction is also limited by the low SNR and the high-tilt angles.
3 Determination, correction, and refinement of CTF
During the imaging process in cryo-ET, the low electron dose, the thickness of the specimen, and the tilt angles result in intricately noisy and blurred projections. The aberration of the lens and the defocusing give rise to the CTF, which restricts the recording of the high-resolution information (Marabini et al., 2015). The CTF causes periodic signal inversion over a range of frequencies and completely loses information on signals at zero-crossings (Turonova et al., 2017). Furthermore, the tilted images in cryo-ET have defocus gradients perpendicular to the tilt axis, resulting in different defocus values for different angles of projections. Therefore, CTF determination and correction are very significant steps before further reconstruction processing.
The CTF is an oscillatory function, which is more severe at high frequencies and defoci. These oscillations cause changes in contrast and spectrum amplitudes for projections in Fourier space, which weakens high-resolution information. The CTF is affected by the defocus, spherical aberration, the wavelength of electrons, etc. (Voortman, 2014). From a physical point of view, the CTF (Zhang, 2016; Kunz and Frangakis, 2017) is defined as
where γ(t) = πλt2(Δz − 0.5λ2t2Ct), t is the modulus of the spatial frequency t, λ is the electron wavelength, Ct is the spherical aberration constant, Δz is the defocus, and A is the amplitude contrast coefficient. The Fourier transform of the projected tomography is multiplied by the CTF to obtain the Fourier transform of the observed image, i.e.,
where f, I are the observed image and the underlying projection, and is the Fourier transform. The effect of the CTF on the images’ magnitude spectrum is shown in Figure 3. The Fourier inverse transformation of the CTF is called the point spread function (PSF), denoted as A, i.e., −1(CTF(t)) = A(x) (Croxford et al., 2021). The PSF in real space blurs projection images as Figure 3 shows.
Figure 3 An illustration of the effects of the CTF and PSF on the 2D projected images. First row: blurred images; second row: the magnitude spectrum. From left to right, the first ones are original, and the others are increasingly blurry images and their corresponding magnitude spectrum.
Most existing literature mainly focuses on the methods of CTF determination and correction in Fourier space, which contains 2D-CTF and 3D-CTF determination and correction. Furthermore, these CTF correction methods usually combine with other steps like tilt series alignment and 3D refinement (Fernandez and Li, 2021). The 2D CTF correction approaches include e2ctf in EMAN2 (Tang et al., 2007), CTER (Penczek et al., 2014), CTFFIND4 (Rohou and Grigorieff, 2015), Gctf (Zhang, 2016), and goCTF (Su, 2019). The CTER method (Penczek et al., 2014) can estimate the parameters of CTF rapidly and accurately. These advantages imply that the CTER can spend less time analyzing high-throughput data and obtain more accurate results than other methods. The CTER can also provide an error assessment which is a cross-resolution correlation 1D function, which evaluates the similarity of the 3D structures with CTF correction or not. In 2015, the effects of 2D CTF parameters and correction methods were tested (Marabini et al., 2015) among ctfind, xmipp, appion, e2ctf, spider, and eman, but not for goCTF and Gctf. It shows that most CTF correction methods have minor differences even though each approach has its unique advantage. In Marabini et al. (2015), relatively speaking, the CTFFIND4 (Rohou and Grigorieff, 2015) has more stable and better results. CTFFIND4 (Rohou and Grigorieff, 2015) defined a specific formula for the spatial frequency to estimate CTF. Gctf (Zhang, 2016) is a GPU-accelerated global and local CTF estimation program based on astigmatism rotational averaging and self-consistency verification. In Gctf, Δz = z(θ) = zucos2(θ −θast)+zvsin2(θ −θast), where zu, zv are maximum and minimum defocus, and θast is the angle between zu axis and x-axis. The goal in Gctf is
where F(t) is the amplitude spectrum, Bg is the calculated background, CTFsim(t) is the simulated CTF, B is the B-factor, and CC is the cross-correlation. Another method, GoCTF (Su, 2019), is a geometrically optimized approach for estimating the global focus gradient and local focus refinement per particle. The typical cryo-EM SPA approaches for the CTF correction do not take the defocus gradient into consideration.
There are some 2D correction CTF methods designed for cryo-electron tomography. The method in Fernández et al. (2006) is a classical method for CTF correction. Only a part of this method is implemented in the software CTFFIND3. To be specific, the defocus D(x) at any images of tilt series can be determined by the simple formula, i.e., D(x) = D0 + d(x)tan(θ), where D0 is the defocus at the untitled plane, d(x) is the distance to the tilt axis and θ is the tilt angle. Then, the CTF is estimated by the defocus, and its inverse Fourier transform is used to correct the particular tilted projection. An extended image acquisition method (Eibauer et al., 2012) is proposed to determine CTF by recording two high-dose images along the tilt axis for each tilt angle. Galaz-Montoya et al. (2016) proposed a per-particle CTF correction method for a single particle in cryo-ET, which automatically computes CTF for each particle and tilt-series images. This method proposes a linear regression strategy to fit the defocus gradient and further perform CTF correction. The Bsoft (Heymann, 2018) software is a tool including 2D-CTF correction for reconstructing macromolecular structures. The other software and methods for 2D CTF correction are shown in Table 1.
3D CTF correction methods for cryo-ET also exist, like NovaCTF (Turonova et al., 2017), TomoAlign (Fernandez and Li, 2021), 3D CTF (Kunz and Frangakis, 2017), and Entropy-regularized deconvolution (ERDC) (Croxford et al., 2021). The NovaCTF (Turonova et al., 2017) is a popular procedure for CTF correction in cryo-ET, which realizes CTF correction in each image with different defocus many times and improves reconstruction resolution. TomoAlign (Fernandez and Li, 2021) is also famous. The TomoAlign procedure realizes CTF correction and tilt-series alignment simultaneously, where the results of 2D CTF correction and tilt-series alignment in IMOD can be the initial inputs. To be specific, the defocus of the ith image in tilt series is Di(p) = pxsin(θi) + pzcos(θi) + Di(o), where p = (px, py, pz) is the particle coordinates with the center o of the tomogram, Di(p) and Di(o) are the averaged defocus and the defocus of the ith image, and θi is the tilt angle. The estimated CTF is determined by the defocus. Three-dimensional CTF (Kunz and Frangakis, 2017) is performed on each slice of the reconstruction instead of the micrographs, which is based on the Jensen-Kornberg scheme. The ERDC approach by Croxford et al. (2021) is a deconvolution process based on entropy regularization in real space, where the regularization plays a role of robustness for noise.
Both 2D and 3D CTF correction approaches, the accurate estimation of parameters in CTF is crucial and challenging. With the increase in cryo-ET data and the development of artificial intelligence algorithms, it may be possible to estimate the parameters of physical formulas for CTF by network methods in the future. At present, there are still many gaps before the network approach is applied to CTF. For instance, there is no benchmark dataset and no ground truth to verify the accuracy of CTF correction methods.
4 Tilt series alignment
The probability of misalignment of the projection is induced by the distortions during tilt series acquisition, including shifts, rotations, translations, and magnifications. The misalignment of the projections would hamper the SNR enhancement and 3D restoration in the following workflow. Tilt series alignment is an essential component of procedures modifying the impacts of several kinds of distortions by modeling the projections onto 2D images. In addition, full tracking in every tilt angle allows more complete data collection at high magnification, which can be applied to various tilt strategies, including the dose symmetric scheme. Specifically, the dose-symmetric tilt series scheme (Hagen et al., 2017) starts at zero degrees’ tilt and moves up to the highest tilt in both tilt directions simultaneously. This method concentrates high-resolution information in the lower tilts where the sample is thinner, thus providing maximum information transfer. With the assistance of full tracking, the dose-symmetric tilt series scheme avoids resolution limitations, improves information transfer, and minimizes the effect of sample distortions.
After motion correction, CTF correction and coarse tilt series alignment, a more precise projection model is needed for fiducial-based alignment and feature-based alignment. For fiducial alignment and calibration, “gold beads” in sample preparation are tiny gold nanoparticles used as fiducial markers. Geometrically, once the gold beads are placed as common fiducial markers in sample preparation for precise alignment, it follows the coordinate systems of the specimen, microscope, and each projection image, S = {Os,s1,s2,s3}, M = {Om, m1, m2, m3}, Bi= {Oi, ui, vi}. In Figure 4, the geometric projection formula of a position in the sample onto the ith image (Hanssen, 2018) can be represented as:
where N is the number of images of tilt series, Nm is the number of fiducial markers, pi,j= (ui,j, vi,j) is the coordinate of the projection of the jth fiducial marker in the ith view, rj is the coordinate in M of the jth fiducial marker, di= (Δui,j,Δvi,j)T is the shift relative to Om, P denotes the matrix for projection onto (u, v) plane, Di denotes the matrix representing specimen changes with respect to the scaling factor mi, the additive scaling factor si, the skew rotation δi and multiplicative factor ti, and R(γi),R(βi),R(αi) are rotation matrices around the z-, y-, and the x-axis of M. By minimizing the total error relative to measured projection and estimated projection pi,j(Hanssen, 2018),
The over-fitting problem is common in fiducial-based alignments because the number of parameters is relatively large for the numbers of fiducial markers as well as the limited number of tilt series, albeit there are techniques for reducing parameters and making the model more robust. In TomoAlign (Fernandez and Li, 2021), the projection matrix is reduced and the motion-aware polynomials are induced as a term on the right side of Eq. (4), pi,j= Mirj+ Si,j+ di, where Mi= miPRi, Ri including the rotation matrix around three axes in M, denoting the beam-induced motion undergone by the sample in the ith image, which corresponds to a number of parameters by giving a different type of the doming deformation polynomials or splines. Specifically, the doming motion causes an upward rise at the center of the carbon hole in the electron microscope grid while the periphery remains relatively unchanged.
Figure 4 The imaging of cryo-ET tilt series. (A) A sample diagram of obtaining tilt series by projecting specimens in cryo-ET, where the tilt angles range from −70◦ to 70◦. (B) A three-coordinate system of obtaining tilt series in cryo-ET. The relative deviation between theoretical and measured coordinates pi,jand . The different colors, namely, red, green, and yellow, represent different particles of interest. In the projection plane, there are theoretical and measured coordinates of features/fiducials.
Feature-based alignment is a classical marker-free alignment approach suitable for lamellae, which is recognized for practicality due to the fewer requirements for sample preparation and the broader scenarios, such as the prohibition of the fiducials. The difficulties in feature localization, extraction, and matching exceed similar procedures in fiducial-based alignments like ClusterAlign (Seifer and Elbaum, 2022). In contrast, the substantial number of virtual markers contributes a statistically robust and accurate alignment with relatively expensive computation. Satisfactorily, several automatic platforms and programs show their performance on feature-based alignments and joint alignments like AuTom (Mastronarde and Held, 2017), AuTom-dualx (Han et al., 2019b), and a joint method (Han et al., 2019a). The AuTom package contains marker-free (Han et al., 2014)and marker-based (Han et al., 2015) alignment. The marker-free method (Han et al., 2014) is based on Scale-Invariant Feature Transfrom (SIFT), which realizes more accurate feature extraction, more stable feature description, and more robust feature matching and tracking. The marker-based method is a fully automatic alignment of data with fiducial markers (Han et al., 2015). AuTom-dualx (Han et al., 2019b) is a tool for marker-based alignment of cryo-ET dual-axis tilt series. At first, the two tilt series are dealt independently with detecting and tracking of fiducial marks, and bundle adjustment for single-axis data. Then, the two tilt-series results are combined by transforming their axes to obtain global projection parameters and tracks by solving a non-linear least square problem. Finally, the last stage is global parameter optimization of tilt series based on a simultaneous reconstruction. The standard method in Han et al. (2019a) is an iterative approach for marker-free alignment, which combines track-based and intensity-based alignment. Some software also provide automated tilt-served alignment like IMOD (Mastronarde and Held, 2017), EMAN2 (Murray et al., 2014), and Dynamo (Scaramuzza and Castaño-Díez, 2021).
The traditional approaches in marker-free alignments are dedicated to calculating the shift for the highest cross-correlation in the Fourier space, determining the common line of each 2D projection based on the Fourier slice theorem and the projection matching by constructing a 3D model as an intermediate reference. Sometimes, the traditional approaches serve as coarse alignment and preprocessing schemes. Aretomo (Zheng et al., 2022) utilizes the traditional marker-free methods with stages of tilt-angle offset determining planarity and recursion, in-plane rotational alignments, translational alignments, and local alignments with patches, which reduces the impacts of the distortions above progressively at sufficient accuracy.
5 2D and 3D denoisers
The electron dose is strictly limited due to radiation damage to the biological samples (Han et al., 2022), which is an essential factor for noise. Electron beam intensity and exposure time are factors to control the final exposure dose on the sample. During the imaging process, tilt angles usually vary from −70◦ to 70◦, which causes a missing wedge (Moebel and Kervrann, 2018). In other words, the parts of projections will be missing since the specimen cannot be rotated to plus/minus 90◦. As the tilt angle increases, so does the ice thickness that the electrons have to penetrate. Low electron dose and thick ice are two significant factors to lower the SNR. Figure 5 shows the effect of noise on an image.
Figure 5 An illustration of the effect of noise on an image. From left to right, the left one is an underlying image; the others are noisy images as the noise becomes larger and larger.
Noise in cryo-ET is a critical factor in limiting the resolution of three-dimensional reconstructed structures (Frangakis, 2021). In the reconstruction processing, Aligning and averaging of the same particles is a well-known approach to remove noise and enhance the signal. However, the STAA step is appropriate for a number of homogeneous macromolecules, not for intrinsic molecular flexibility and dynamics with a continuum or discrete changes of conformation, such as proteins in folding, antibodies, and lipoproteins (Han et al., 2022). The intrinsic molecular flexibility and dynamics needs some denoisers for resolution enhancement. Denoisers usually remove noise and enhance the image contrast in the low frequencies. High noise, complicated noise types, and no ground truths are three factors that make the denoising task in cryo-ET complex (Yang et al., 2021). These three factors make most existing denoisers unsuitable for tomographic images. Next, we introduce 2D and 3D image denoisers in detail.
2D denoisers are applied to tomographic projections, where we focus on relevant mathematical methods and deep learning methods. Natural images are taken by the camera according to the sensitive unit, like images taken by mobile phones. Here, 2D projection images by cryo-ET are unnatural images. The existing mathematical methods for denoising natural images include total variation (TV)-based variants (Pang et al., 2020), block matching and 3D collaborative filtering (BM3D) (Maggioni et al., 2013), k-singular value decomposition (k-SVD) (Aharon et al., 2005), and traditional Wiener filter (TWF). However, the types and levels of noise in the tomography of cryo-ET are different from those of natural images, which cannot obtain a satisfying result on tomographic projections (Gu et al., 2022). The imaging principle of cryo-EM and cryo-ET is the same. However, compared with cryo-EM, cryo-ET needs to project specimens under different view angles and assign electronic dose for limiting radiation damage. Thus, 2D projection images in the tilt series by cryo-ET have a lower SNR than 2D projections by cryo-EM. To date, there are few 2D denoisers specifically designed for cryo-ET, while many 2D denoisers are designed for cryo-EM. Some researchers applied the methods of denoising images for cryo-EM to 2D projection images in cryo-ET as a preprocessing step. Mathematical models have been developed for denoising images in cryo-EM in the past few years. The noise of 2D projections is usually assumed to follow Gaussian distribution, while the actual noise is complicated. Although these assumptions for mathematical modeling are simplified relative to the actual situation, the experimental results are acceptable. A mathematical approach named Covariance Wiener Filtering (CWF) (Bhamre et al., 2016) was proposed for denoising images and correcting the CTF of 2D projection images simultaneously. In the method, the imaging model in Fourier space is Yi= AiXi+ ξi for i = 1,2,…,n, where Yi,Ai,Xi,ξi are the 2D observed image, spread point function, 2D clean image, and noise. The CWF approach estimates the mean µ and covariance Σ of the underlying clean images, chooses Fourier Bessel basis to represent the images, and then obtains restored images by the Wiener filter based on and instead of estimating the spectral SNR in TWF. The CWF method takes projection directions, ice thickness, and structural variability into account in Σ. The estimation can be calculated by
where λ > 0 is the regularization parameter. The estimation can be calculated by
where σ2 is the variance of Gaussian noise. After obtaining and , the restored images can be solved by the Wiener filter as follows:
where is the Wiener filter (Frappart et al., 2016). In 2022, an improved version of the CWF method named Ab-initio Contrast Estimation (ACE) is proposed (Shi and Singer, 2022), which improves contrast estimation. The approach is derived under the imaging model Yi= ciAiXi+ ξi, where ci is the unknown contrast of the clean image. In the ACE method, ci makes the estimation of more accurate and the resolution of the repaired images higher. Although the CWF and ACE methods can realize an improved result, this issue still has much room for improvement due to their strict assumptions.
Next, we turn to deep learning methods for 2D denoising. Many deep learning-based methods have a satisfying performance on denoising in cryo-ET (Frangakis, 2021). The Generative Adversarial Network (GAN) is a network for generating images. GAN consists of a generator and a discriminator, which are neural networks. The generator tries to generate fake data and cheat the discriminator. The discriminator distinguishes the real and fake data. The generator and discriminator are neural networks that run in competition during training. In 2018, the β−GAN approach with l1 and l2 loss (Gu et al., 2022) was proposed for denoising 2D cryo-EM images. To overcome the lack of ground truths, Noise-Transfer2Clean (NT2C) (Li et al., 2022a) applies the simulation software InSilicoTEM (Vulović et al., 2013) to produce simulated 2D projection images in cryo-EM. There are three parts in NT2C: noise extraction, noise modeling, and denoiser training. There is the coarse convolutional neural network (CNN) denoiser and the refined CNN denoiser, which are connected by fusing the noise in real noisy images to the simulated clean images through the contrast-guided noise and signal re-weighting step. NT2C can realize many denoising results, but its operation is tedious, and its results seriously depend on simulated clean images. Furthermore, some deep learning-based methods do not rely on ground truths, such as Noise2Noise (N2N) (Lehtinen et al., 2018) and Noise2void (N2V) (Krull et al., 2019). The Topaz method (Bepler et al., 2020) is a version of N2N, which is for 2D denoising. This method contains four network structures, including affine, fully connected neural network (FCNN), small U-net convolutional network, and U-net convolutional network. The loss function is a measure of evaluating the difference between the predicted results by the algorithm and the real dataset. The loss function is usually reduced during training the network. Its loss can be l0, l1, and l2 as follows:
where i = 0,1,2, and xa, xb are noisy neighboring odd and even frames, f is the network, X represents the distribution of images, and represents mathematical expectation. At the overall results, the U-net can realize a better result. However, the model’s training data cannot cover all cases of noisy images in cryo-EM. From another point of view, the losses of the Topaz model are simple, and maybe one can replace the losses with an adaptive loss based on real data.
3D denoisers can be applied to tomograms. The popular methods include Block-Matching and 4D filtering (BM4D) (Maggioni et al., 2013), a Monte Carlo Framework (Moebel and Kervrann, 2018), 3D Noise2Noise (3D-Topaz) (Bepler et al., 2020), 3D Noise2void (Krull et al., 2019), and the Sparsity Constrained Network (SC-Net) (Yang et al., 2021). BM4D was designed for denoising a single 3D natural image. One can apply BM4D directly to a single noisy volume image of cryo-ET. The Monte Carlo Framework is a general denoising framework to overcome noise in cryo-ET. One can integrate any 3D Gaussian denoiser like BM4D into the framework. The Topaz method proposes two network frameworks, including U-net-3d-10a and U-net-3d-20a, where the difference is the average pixel size of training data 10 Å and 20 Å, respectively. The training data of the 3D Topaz method are based on the even/odd frames of 32 aligned cryo-ET tilt series and the training time is over 1 month. The 3D Topaz can usually realize a good denoising performance and is regarded as a benchmark 3D volume image denoiser in cryo-ET. The 3D Noise2void method is trained in the CNN on noisy images, which does not need noisy image pairs or ground truths. SC-Net is a self-supervised learning network with a sparsity constraint. The SC-Net is based on the 2D filter to produce the 2D reconstructed images and then apply the 3D reconstruction method to obtain smoothed 3D volumes to guide the network training. The inputs of SC-Net are a linear combination of four parts, including raw noisy volumes and smoothed volumes. The loss function of SC-Net contains volumetric reconstruction, sparsity-guided smoothing, expectancy constraint, and regularization. There are some comparisons among SC-Net, BM4D, N2N, and N2V on simulated and real-world data (Yang et al., 2021). SC-Net can realize competitive performances with Topaz.
Although the existing methods can achieve acceptable denoising results for some cases in cryo-ET, there are still room to improve resolutions. For 2D denoisers, most denoising methods are specially designed for cryo-EM, not cryo-ET. Because of the different image acquisition strategy between cryo-EM and cryo-ET, the noise in projections is different. On one hand, due to the constraints of the total electron dose, noise is higher in the per-tilt projection image recorded for cryo-ET. In comparison to cryo-EM acquisitions, the full dose is used to record only one image in zero tilt. On the other hand, as the tilt angle increases, so does the thickness of the ice and sample the electrons need to penetrate, prohibiting the recording of the high-resolution signal to these projection images. Hence, the 2D projections by cryo-ET is very noisy. We should be careful to apply 2D denoisers in cryo-EM, which are directly applied to denoise 2D projection images of tilt series in cryo-ET. The factors such as lower electron dose and changes in ice thickness should be considered. Designing more accurate 2D tomography denoisers according to different imaging conditions in cryo-ET is an urgent task in the future. For 3D denoisers, there are several popular methods for denoising 3D volume images in cryo-ET, including BM4D, Topaz, N2V, and SC-Net. BM4D is for a single 3D volume image. Topaz, N2V, and SC-Net are networks based on training data. The difficulties of denoising 3D volume images are no ground truths, high noise levels, and complex noise types. To address this, some imaging simulation algorithms for cryo-ET can be designed to fill up the shortage of no ground truths, or some novel 3D denoising network structures that do not need ground truths can be constructed. Alternatively, some adaptive 3D denoising mathematical models can be proposed in the future. Restoring 3D volume images in cryo-ET is a tough but important task for improving resolution.
6 3D reconstruction and reduction of missing wedge
The 3D reconstruction in cryo-ET aims to restore 3D volume images from 2D aligned images (Frank, 2005). Owing to the limitation of electron dose and constraint of tilt angles, low SNR and artifacts caused by the missing wedge are two difficult points to handle in cryo-ET (Ding et al., 2019); thus, many existing 3D reconstruction methods are not suitable for cryo-ET. Figure 6 shows a simple illustration of 3D reconstruction step from the tilt series in cryo-ET, where the complete volume suffers from the missing data. In electron tomography (ET), the most common strategy to acquire data is the single-axis tomography (Han et al., 2017), although dual-axis tomography (Tong et al., 2006) is also used. Single-axis imaging rotates around the same axis during the imaging process, while dual-axis imaging rotates around two axes that are perpendicular to each other, which can obtain two different tilt-series images. The latter is designed to reduce the missing information to a “missing pyramid” experimentally. For the single-axis imaging process, the missing wedge (MW) information makes the reconstruction problem ill-posed. In some sense, MW can be restored by 3D reconstruction methods (Zhai et al., 2021). Then, the restored 3D volumes can be transported to 3D denoisers or 3D subtomogram image alignment and averaging to obtain high-resolution results. For the dual-axis imaging process (Han et al., 2019b), combining information from two tilt-series images is a difficult task.
Figure 6 This is an illustration of the 3D reconstruction step from tilt series in cryo-ET, where the final volume suffers from the missing wedge.
Over the past few decades, many approaches have risen to deal with 3D reconstruction for the single-axis imaging process in cryo-ET. We firstly summarize the existing literature on 3D reconstructions proposed for the single-axis imaging process. The mathematical models can be divided into two categories: real space and Fourier space. The Fourier projection slice theorem (also known as the central section theorem) (Frank, 2005) is the key point of Fourier-space-based methods. The Fourier projection slice theorem tells us that the Fourier transform of 2D projection images equals the slice operator acting on the Fourier transform of 3D volume, whose mathematical expression is
where V is a 3D volume, Fi is the i−D Fourier transformation for i = 2,3, P is the 2D projection operator that integrates along the z-axis, and S is the slice operator that sets the component of the z-axis equal to 0. The Fourier transformation does not lose any image information. The Fourier-based approaches contain filtered back-projection (FBP) (Zeng, 2012), non-uniform fast Fourier transforms (NUFFT) (Chen and Förster, 2014), and the iterative compressed-sensing optimized non-uniform fast Fourier transform reconstruction (ICON) (Deng et al., 2016). The FBP method uses the same filter function for all projections. The ICON is a popular 3D reconstruction software based on NUFFT with compressed sensing and sparsity assumption in cryo-ET. In ICON, the imaging process is expressed as f = Ax + N in Fourier space, where f and x are the discrete Fourier transform of projections and the density of specimen, respectively, A is a non-uniform Fourier sampling matrix, and N is noise. Then, the mathematical model of ICON (Deng et al., 2016) is
where P is a sparse transform matrix, W is the weight of nonuniform sampling, AH is the conjugate transpose of A, and ϵ is a constant.
In addition, there are many mathematical models in real space for 3D reconstruction in cryo-ET. Because of the missing wedge, the mathematical problems for 3D reconstruction are ill-posed. The most common methods include weighted back-projection (WBP) (Carazo et al., 2006), a regularization-based method (Albarqouni et al., 2015), the Monte Carlo approach (Moebel and Kervrann, 2020), sMAP-EM (Paavolainen et al., 2014), an iteration-based method (Zhai et al., 2021), and an image-inpainting-based method (Trampert et al., 2018). The WBP method applies different filters for different projections according to their geometric characters. For the iteration-based approaches, the imaging formation process can be simplified as f = Ax + N. Specifically, there are the algebraic reconstruction technique (ART) (Gordon et al., 1970), an adaptive simultaneous algebraic reconstruction technique (ASART) (Wan et al., 2011), the simultaneous iterative reconstruction technique (SIRT) by Trampert and Leveque (1990), the progressive stochastic reconstruction technique (PSRT) by Turoňová et al. (2015), direct iterative reconstruction of computed tomography trajectories (DIRECTT) algorithm (Kupsch et al., 2015), model-based iterative reconstruction (MBIR) (Deng et al., 2016), low-tilt tomographic reconstruction (LoTToR) (Zhai et al., 2021), and the constrained reconstruction model (CRM) (Han et al., 2021). The ART method aims to solve a large and sparse linear equation system, where one can obtain an approximate solution by Kaczmarz’s iterative method. There are some variants of the SIRT method that contain modified simultaneous iterative reconstruction (MSIRT) (Wan et al., 2009) and weighted simultaneous iterative reconstruction (WSIRT) (Wolf et al., 2014). The WSIRT method is a combination of the WBP and the SIRT methods. The CRM introduces a constrained term based on SART and SIRT. The minimization formulation of CRM (Han et al., 2021) is
where bi and gi are sample images and their relevant backgrounds, Di and Bi are the image formation operators of bi and gi, respectively, and fi is their measured projection. The sparsity Kaczmarz’s algorithm solves the CRM problem. Another alternative approach is to regard missing wedge as an image inpainting problem (Trampert et al., 2018), which repairs the missing part with some general information about the images.
There are several neural network methods (Ding et al., 2019; Liu et al., 2022b) for 3D reconstruction and the reduction of missing wedge in the single-axis imaging process of cryo-ET. IsoNet (Liu et al., 2022b) is a neural network method for reconstructing the missing wedge information and increasing SNR. The IsoNet contains five steps: deconvolve CTF, generate mask, extract, refine, and predict steps. The refine step plays a crucial role in IsoNet. The refine step needs a pair of tomograms to train. However, the extracted subtomograms have missing wedge. The IsoNet adopts an iterative way to update the pairs with the predicted subtomograms, which takes advantage of deep learning to fill up the missing wedge by predicting information. The joint deep learning model (JDLMID) (Ding et al., 2019) is for the inpainting of the missing wedge and de-artifact. JDLMID tries to fill up the missing wedge by image inpainting. The JDLMID method contains two parts: GAN frameworks with residual in a residual dense block for inpainting and a U-net for eliminating artifacts. In the future, more neural networks for 3D natural image reconstruction can be revised to be suitable for 3D reconstruction in cryo-ET.
Besides the single-axis imaging process, the dual-axis imaging process (Mastronarde, 1997; Tong et al., 2006; Han et al., 2019b) is an acquisition strategy to collect more data and reduce the missing information to a “missing pyramid” in the experiment (Guesdon et al., 2013). However, the projections of the dual-axis tomograms have a lower SNR than the single-axis type, and matching two orthogonally oriented tilt series is a new and difficult task. For the dual-axis imaging process in cryo-ET, the Etomo (Mastronarde and Held, 2017) in the software IMOD and PEET can deal with the dual-axis data. Moreover, Han et al. proposed a new toolkit called AuTom-dualx (Han et al., 2019b) for the 3D reconstruction of dual-axis tomograms. The combination of both two tilt series in AuTom-dualx is based on fiducial markers. They designed a new algorithm to match the two tilt-series coordinates and then embedded it into the SART and SIRT methods for 3D reconstruction. The dual-axis tomogram (Winter and Chlanda, 2021) has analyzed the interactions of the Ebola virus and action-VP40.
In brief, there are two imaging processes: singe-axis and dual-axis. As shown above, many mathematical models and neural networks have been developed for 3D reconstruction in cryo-ET for the singe-axis imaging process. They aim to reconstruct 3D volumes from 2D projections and try to fill up the missing wedge by prediction. Some restored methods depend on prior knowledge and are not data-driven, like MTV (Albarqouni et al., 2015). Others rely on the known information of observed projections to restore the missing wedge (Liu et al., 2022b). Although the existing methods have realized some acceptable 3D density maps of specimens in some cases, there are some areas for improvement. The methods based on prior constraints may generate elongation, streaking, and ghost tail artifacts due to the lack of using the information of the acquired projections. The existing methods depending on the acquired projections may not obtain ideal reconstructed results from incomplete projection orientations. In the future, we hope that more well-performed restored-processing methods can be developed for reconstructing 3D density maps in cryo-ET image processing.
7 3D particle picking
Since a 3D reconstructed tomogram contains multiple copies of specimens, the 3D particle picking task focuses on detecting and picking the particles of interest to biologists. As Figure 7 shows, there are many noisy particles of different orientations. The higher the noise, the harder to detect the particles of interest. There are four factors that make 3D particle detection difficult: the low SNR in 3D reconstructed tomograms, missing wedge, the small sizes of particles, and the crowded environment (Castaño-Díez et al., 2017; Mo et al., 2021). Figure 8 shows simulations of 2D projection images by rotating the biological sample, where the particles are located in an overcrowded environment. The overcrowded environment and the difficulties of picking particles in situ are shown in Figure 9, that is, for the case of particles inside the cells. The 3D particle picking task contains two tasks: detection and segmentation. The low SNR makes it impossible to detect particles by eyes and requires computational methods. The existing methods of particle picking usually focus on coordinate identification and extraction of 3D targeted particles in volumes. There are segmentation-based methods like instance segmentation (Liu et al., 2022a), CASSPER (George et al., 2021), and detection-based methods like one-stage detection (Liu et al., 2022a). These methods for picking 3D particles in cryo-ET are done in a manual or semi-automated pattern. Manual marking particles are empirical and time-consuming. In recent years, many methods based on artificial intelligence have been developed.
Figure 7 An illustration of the 3D particle packing task. For the first row, from left to right, the first one are the underlying particles. The second and the third rows are noisy tomograms as the noise becomes larger, so does the left one in the second row. The last one detects and marks particles in 3D particle picking.
Figure 8 An illustration of simulation of 2D projection images in different tilt angles by rotating the biological sample. The top picture shows the particles located in an overcrowded environment. The other pictures are 2D projection images as the tilt angle increases from left to right and from top to bottom. The cubes stand for particles.
Figure 9 This is the biological environment in situ for picking macromolecular complexes (represented by green balls) inside a cell.
Recently, many new neural network methods have been developed, like those based on 3D-CNN, manifold, U-net, active learning, and so forth. It can be divided into template-based and template-free methods. For the former, the Dynamo provides a basic geometrical library, and users can select a specific one to do 3D particle picking. Here, we only focus on the latter. We collect some methods for detecting and segmenting particles. The MemBrain (Lamm et al., 2022) and 3D-UCaps (Hajarolasvadi et al., 2022) are designed for detecting membrane–protein complexes. The MemBrain is a deep learning frame combined with using CNN for scoring and means shift clustering. The 3D-UCaps combines 3D Capsule Network (CapsNet) and CNN to decrease the computing cost and improve accuracy. A learning method by integrating dense CNN and autoencoder into U-net (Zhou et al., 2023) is for the cell membrane segmentation. In addition, Pyle et al. (2022) propose several strategies for picking and cleaning membrane–protein complexes. The Deepfinder (Moebel et al., 2021) and TomoTwin (Rice et al., 2023) are designed for detecting macromolecules. The Deepfinder can identify multiple macromolecular complexes simultaneously by combining U-net, a 3D CNN-based model, and a clustering algorithm. TomoTwin maps particles into a high-dimensional space that is called the embedding space. It measures particles by embedding their representations. To save computing time, the YOLOv3 (Wu et al., 2022) regards particle picking as an object detection task using one-stage detection. To combine 2D detection and segmentation methods into 3D cases, a novel 2D-to-3D framework (Zeng et al., 2022) for detection is used by detecting objects firstly in 2D projections and then localizing objects in 3D structures. Active learning (Mo et al., 2021) is an active-learning-based approach designed with region proposal and Resnet3d under a new class label. There are some learning-based approaches, like the cryo-ET one-shot network (COS-Net) (Zhou et al., 2021). In addition, there is a new method to simulate crowded cells in cryo-ET (Pei et al., 2016) to check the effects of 3D particle picking methods in a crowded particle environment. Some software contains 3D particle picking like Dynamo (Castaño-Díez et al., 2017), EMAN2 (Murray et al., 2014), Bsoft (Heymann, 2018), TomoTwin (Rice et al., 2023), Deepfinder (Moebel et al., 2021), and VP-Detector (Hao et al., 2022).
In fact, 3D particle picking is a time-consuming and computing-consuming task. Because the sizes and shapes of particles vary, the existing methods have their cases for the application. How to select an appropriate method is a question that is worth exploring. In the future, embedding biological samples in manifolds can be further developed. However, finding a suitable embedding space and representation needs to be explored and experimented. Some novel neural networks can also be extended to 3D particle detection or segmentation in cryo-ET.
8 Classification
Biomolecules constantly change their structures and conformations to function in the cellular environment. Cryo-ET enables the discovery of these diverse conformations in situ (Zhang, 2019). Yet, classifying these conformations accurately is challenging due to low SNR, uncertain classification numbers, and imbalanced class proportions (Taylor et al., 2006). Additionally, the missing wedge introduces variations in the structural information, depending on the orientations of the biomolecules (Heumann et al., 2011).
Classification of the subtomogram is used to extract consistent structural features from heterogeneous conformations (Wan and Briggs, 2016). Specifically, the fundamental purpose of classification is to find multiple conformations and compositions of known complexes, and discover unknown complexes. Similar to tilt-series alignment and subtomogram alignment, classification also relies on cross-correlation optimization in the tomograms (Frangakis et al., 2002; Heumann et al., 2011; Hrabe et al., 2012). In addition, pattern recognition schemes, based on a priori knowledge, are commonly used in research to classify the particles of interest. These schemes based on a priori knowledge require high-resolution experimental structures to generate clusters that are subjected to imaging experimental parameters. At the same time, low SNR and missing wedge also limit the feature representation of clusters.
The current approaches to perform classification are simultaneous within subtomogram alignment or are based on a priori knowledge. On the one hand, when the classification is decoupled from the subtomogram alignment, the constrained principal component analysis (PCA) method based on K-means or hierarchical clustering is often used, which brings higher computational consumption as its voxel-by-voxel comparison. For the interested complex and its specific functional regions, the PCA method has limitations when dealing with data distributions that exhibit higher levels of heterogeneity. Generally, the binning method is used to reduce the dimensions of data. On the other hand, the comparison method with multiple references is often used in simultaneous alignment and classification. The multiple references approach (MRA) relies on a priori knowledge of existing high-resolution structures. Thus, it can achieve higher cross-correlation coefficients by minimizing the probability of maximum likelihood estimation. The noise outside regions of structural change can reduce classification accuracy. To address this limitation, MRA was developed for inter-reference structural differences (Chen et al., 2014). Although the MRA is theoretically computationally simpler, it requires iterative operations to achieve classification convergence, which requires more demanding multi-reference sampling. Software packages for classification are available in Dynamo (Scaramuzza and Castaño-Díez, 2021), RELION (Bharat and Scheres, 2016b), emClarity (Ni et al., 2022), and EMAN2 (Murray et al., 2014; Chen et al., 2019).
In recent years, bioinformatic deep learning methods have significantly improved the classification of heterogeneous subtomogram of interest particles. The tomoDRGN (Powell and Davis, 2023) can reconstruct a heterogeneous ensemble of structures by learning a new neural network including a variational autoencoder (VAE) and decoder. This network can learn a continuous low-dimensional representation of structural heterogeneity in cryo-ET data (Powell and Davis, 2023). A learning-based approach for sub-break classification is proposed in Xu et al. (2017), addressing the challenges of complex data distributions and scalability. This approach combines deep learning for supervised feature extraction with unsupervised clustering and template-free classification, offering a complementary solution to traditional methods. Thus, it improves the performance of classification based on deep learning of three designed CNN models for large-scale systematic macromolecular structures. The Respond-weighted Class Activation Mapping (Respond-CAM) algorithm presents a method for extracting and interpreting significant structural features by CNNs (Zhao et al., 2018). An adversarial domain adaptation method is developed to enhance the discovery of new structures and facilitate data recovery in the context of imbalanced classification data (Luo et al., 2019). A combined data processing strategy addresses the challenges of imbalanced classification data and the composition of real molecules, which demonstrates significant improvements in classification and segmentation performance (Che et al., 2018). SHREC is an online 3D Shape Retrieval Challenge that has seen the development of novel deep learning-based methods for classifying cryo-electron tomograms. SHERC 2020 (Zhao et al., 2018) and 2021 (Gubins et al., 2021) evaluated classification results from six and seven groups, respectively. In SHERC 2021, methods like URFinder, DeepFinder, U-CLSTM, MC DS DS Net, YOPO, CFN, and TM-T/TM-F were proposed, with MC DS Net showing superior performance overall.
Classification is an essential step in the whole procedure, which provides input for subtomograms alignment and averaging. Although there are satisfying results for classification, better methods are still needed for higher accuracy and precision. In addition, ContinuousFlex (Harastani et al., 2022) is specifically designed to analyze continuous conformational variability. However, the continuous conformational variability of macromolecules still poses a challenge in this field.
9 Subtomogram alignment and averaging
After the classification of 3D particles (subtomograms), the next step is usually aligning and averaging 3D particles of the same class to denoise (Castaño-Díez and Zanetti, 2019; Zhang, 2019). It is worth noting that STAA is suitable for static conformational variability of macromolecules not for heterogeneous structures. The heterogeneous structures mean that there exists discrete or continuous conformational variability of macromolecules. The heterogeneous macromolecules can be denoised by the denoising algorithms in Section 5. Random orientation and low SNR are two factors that make the STAA task difficult in cryo-ET. There are two keys in subtomograms aligning which are 3D translations and rotations (Leigh et al., 2019). Because of the large shapes of 3D tomograms, computing cost and time should be considered. In this section, we summarize the existing approaches to the STAA problem.
The existing methods for STAA can be divided into three classes: maximum likelihood estimation (MLE) (Stölken et al., 2011) in real space, maximum cross-correlation function (MCCF) (Leigh et al., 2019), and neural networks (Zeng and Xu, 2020). Because the MLE method needs to go through all pixels of 3D subtomograms, the computational cost and time of the MLE method are relatively high. For speeding up the calculation of the MCCF-based methods in real space, there are some fast aligning approaches based on common gradient (Xu and Alber, 2012) and stochastic average gradient (SAG) (Lü et al., 2019) by parallel optimization. The multi-reference alignment (Zhao et al., 2022) is a simple mathematical model for particles alignment and averaging. Furthermore, some fast rotational correction-based aligning methods have been developed by the translation-variant rotational matching in Fourier space under some constraints (Xu et al., 2012). A dissimilarity function of 3D volumes (Förster et al., 2008) has been defined to obtain the 3D translations and rotations iteratively. A constrained cross-correction (CCC) was proposed to reduce the missing information to a “missing pyramid”. The spherical cross-correlation function of Chen et al. (2013) can be regarded as an extension of the MCCF method. There is a Fourier space-constrained fast volumetric matching based on Fourier space equivalence of constrained correlation measure (Xu et al., 2012). There is an improved STA method based on biophysical analysis and supramolecular context (Metskas et al., 2022). Some software can realize the task of STAA, such as Relion (Bharat and Scheres, 2016a), Dynamo (Scaramuzza and Castaño-Díez, 2021), Emclarity (Ni et al., 2022), STOPGAP (Wan et al., 2020), AItom (Zeng and Xu, 2020), and TomoMiner and TomoMinerCloud (Frazier et al., 2017).
In recent years, some methods based on neural networks have arisen for the task of STAA, like Gum-Net (Zeng and Xu, 2020) and Jim-Net (Zeng et al., 2021). The Gum-Net is a geometric unsupervised matching network only for 3D subtomogram alignment. The Gum-Net is based on 3D discrete cosine transform (DCT) for feature extraction, spatial transformer network for alignment, and Pearson’s correlation as the loss to train (Zeng and Xu, 2020). The Jim-Net is an end-to-end unsupervised CNN for the two tasks of cluster and alignment simultaneously (Zeng et al., 2021). There are two parts of the Jim-Net, namely, an aligning process based on constrained cross-correlation loss and a clustering process based on the Gaussian mixture model (GMM). The neural network methods can be better realized when the training sets are large.
In the STAA task, 3D subtomogram aligning is challenging and highly time-consuming. The existing methods have applied parallel optimization in real space and variant translation in Fourier space to speed up the computing process. The neural networks can realize a good result for the STAA that is based on the train data. Some methods for 2D particle alignment and averaging in cryo-EM can be generalized to the 3D STAA in cryo-ET, such as invariant-rotation features: first-, second-, and third-order autocorrelation (Singer, 2019; Marshall et al., 2020). Some manifold learning methods by Kileel et al. (2021) for dimensionality reduction may be considered to solve the STAA task.
10 Post-processing for 3D density map
Owing to the imaging noise, radiation damage, the motion of the samples, and biases in the above reconstruction steps, it is necessary to correct the details in the reconstructed 3D density maps (Rosenthal and Henderson, 2003). The post-processing procedure is used to modify the loss of contrast and fill in the details in the restored 3D density maps to visualize high-resolution features as the final reconstructed step. IsoNet (Liu et al., 2022b) is a deep learning method containing the refinement part for 3D reconstructed density maps, which has already been introduced in Section 5. The post-processing step is used to restore details of 3D density maps, in both cryo-EM and cryo-ET. Here, we summarize post-processing methods into three kinds, namely, map sharpening (Rosenthal and Henderson, 2003), density modification (Terwilliger et al., 2019), and refinement (Urzhumtsev et al., 2022).
Firstly, map sharpening can be roughly divided into global and local approaches. The global sharpening methods for the full volumes are almost based on B-factor correction (Rosenthal and Henderson, 2003). There is also an adjusted surface area measure (Terwilliger et al., 2018) for global sharpening by combining details and connectivity of the maps. The resolutions of the whole reconstructed volumes are usually not the same, which means some parts of 3D volumes are restored well while others are flawed. The local sharpening approaches for raw subtomograms are designed for local restoration, including LocalDeblur (Ramírez-Aportela et al., 2019), LocalScale (Jakobi et al., 2017), LocSpiral (Kaur et al., 2021), an adjustment of density maps in Phenix (Terwilliger et al., 2018), and DeepEMhancer (Sanchez-Garcia et al., 2021). The LocalDeblur is a Wiener filter for image deblurring based on local resolutions of the original maps. The LocSpiral provides a method to compute the local B-factor by the spiral phase transform. The LocalScale and Phenix need atomic maps as references, which can help to obtain good results. However, the requirement of atomic maps is a limitation to use these two approaches. The training data of the DeepEMhancer are pairs of experimental maps and post-processed maps by the LocScale. The well-trained DeepEMhancer is applied to reconstructed maps and realizes better results than the global sharpening methods. Secondly, the density modification (Terwilliger, 2000; Terwilliger, 2001) utilizes the maximum likelihood estimation based on prior knowledge and other known macromolecular structures. The improved procedure by Terwilliger et al. (2019); Terwilliger et al. (2020) needs two independent half-maps and some prior information of the reconstructed maps as inputs which available in the Phenix software. This method is suitable for those maps with preliminary knowledge or similar known macromolecular structures. Thirdly, there are some refinement methods for post-processing, like Servalcat (Yamashita et al., 2021), Relion (Zivanov et al., 2022), and M (Tegunov et al., 2021). The Servalcat computes a weighted Fo− Fc difference map, which is a part of REFMAC5 in the Collaborative Computational Project for electron cryo-microscopy (CCP-EM), where Fo, Fc is the Fourier transform of the observed map and calculated map from atoms. The weighted Fo−Fc difference map can help to enhance weak features like hydrogens. The refinement in Relion (Zivanov et al., 2022) contains CTF refinement and tilt-series refinement, which correspond to optical aberration and geometrical refinement separately. The M (Tegunov et al., 2021) software realizes refinement by a CNN denoiser for 2D images before 3D construction to improve local resolution. Furthermore, there is a framework by Burt et al. (2021) for three-dimensional reconstruction to multi-particle refinement in M, where it provides a link between the Dynamo and the Warp-Relion-M.
In short, the existing post-processing methods further optimize 3D volumes by calculating the whole or local B-factor, deep learning-based density modification, de-blurring, and denoising. In future, some representation methods and measures can be designed to describe the difference between density maps and atomic maps. In addition, as the 3D atomic structures increase, more new neural network methods can be applied to post-processing.
11 Atomic structure building
Biologists can analyze biological functions and conduct medical research through three-dimensional atomic structures (Xu et al., 2019). Atomic structure reconstruction (ASR) from density maps is significant. With the increase of atomic structure data and the improvement of computing equipment, artificial intelligence methods (Giri et al., 2022) have been applied to build atomic models automatically instead of by hand. Some deep learning methods have been summarized in Giri et al. (2022) for building atomic structures from density maps. To avoid repetition, we briefly introduce novel methods from two aspects of building and refining atomic models.
Some deep learning methods build atomic models using the training data of 3D density maps and the corresponding atomic structures. Amino Acid Network (A2-net) (Xu et al., 2019) built an Amino Acid dataset of density maps. The A2-Net contains two stages for detecting 3D amino acids and searching the leading chains of amino acids. It can determine a molecular structure rapidly and accurately. Recently, as the predicting structure methods with sequences develop, there appears some novel ASR methods that apply structure prediction as one of their parts. The inputs are usually 3D density maps and sequences. The Auto-DRRAFTER method by Ma et al. (2022) is designed to build an RNA model automatically based on the structure prediction method Rosetta. Based on the protein structure prediction I-TASSER method, the CR-I-TASSER (Zhang et al., 2022) identifies the templates from the Protein Data Bank (PDB) to promote atomic assembly. The CR-I-TASSER combines the I-TASSER folding simulations and the Cα atom model by deep residual convolutional neural networks and molecular dynamics (MD) simulation to build the atomic structures. There are some similar methods, like domain-enhanced modeling using cryo-electron microscopy (DEMO-EM) (Zhou et al., 2020). The constructed atomic models can be refined further by fitting them to density maps such as an accurate real-space refinement (Urzhumtsev et al., 2022) and refinement in Phenix (Afonine et al., 2018). For the atomic structures, one can also perform interface refinement by HADDOCK2.4 (Neijenhuis et al., 2022) and refine further by correlation-driven molecular dynamics (CDMD) (Igaev et al., 2019). In turn, the CERES method (Liebschner et al., 2021) can refine 3D density maps by the constructed atomic models.
Above all, the existing methods can build or refine atomic models automatically. Although the existing methods can realize good results, there are some limitations like insufficient training data and the low resolution of density maps. Furthermore, as the structure prediction field develops rapidly, novel methods for the ASR can combine them into building atomic model procedures in the future.
12 Discussion
Structural biologists increasingly favor cryo-ET because of its ability to image large biological specimens in situ. The atomic structure of microorganisms is important to make drug discoveries, host–pathogen interactions, and so on. Publicly available databases like EMDB (Lawson et al., 2016) and EMPIAR (Iudin et al., 2016) separately provide 3D density maps and original 2D tilt-series projections. It involves a small amount of resolved density maps and corresponding tilt series applied in the microorganism in nature. However, the amount of data occupies a vast storage, which is a challenge for computer equipment. To the best of our knowledge, there is no unified benchmark of cryo-ET for deep learning method training, testing, and comparing. In future, computer hardware development for massive data storage and uniform standard training and test datasets will be needed.
With the development of computing technology and artificial intelligence, the resolution of resolved structures is increasing. Although many reconstruction software shown in Table 1 have been developed, such as IMOD, RELION, DYNAMO, and EMAN2, each has its characteristics. However, there are several problems: (1) For structural biologists, after obtaining the cryo-ET tilt-series projections, selecting suitable software quickly and accurately remains a difficult problem; (2) Files of different software are in different formats, which cannot be used among different software. The ScipionTomo (Morena et al., 2022) provides an integrated platform where all software is available, including IMOD, Relion, Dynamo, motioncorr2, Gctf, Tomo3d, Aretomo, and NovaCTF. Users can compare different methods for solving the same task instead of wasting time installing software. However, ScipionTomo does not solve the problem of different input/output file types between different software. Furthermore, another new platform can be developed for converting files from different software to uniform formats.
Mathematically, the reconstruction process from a two-dimensional tilt series to a three-dimensional electron density map involves various inverse problems. Specifically, these inverse problems are included in the standard processing discussed above: beam-induced motion fitting, parameter determination of CTF, tilt-series acquisition, noise modeling, and multi-view geometry of three-dimensional reconstruction. In addition, the inherent problems of cryo-ET are the limiting electron dose, the algorithmic enhancement, and the non-unique and unstable numerical solution due to the missing wedge. Thus, a bundle of inverse problems could lead to multiple paths in the pursuit of highly accurate resolution of cells in situ. In general, there should be a wide range of mathematical applications in cryo-ET, which may continue to produce outstanding data processing strategies and methods.
With the rapid development of artificial intelligence (AI), many new deep learning-based methods have been designed for cryo-ET like U-net, Resnet, manifold learning, active learning, automatic encoder, and CNN. However, there are several problems. One is that the disunity of training and testing datasets makes the results of comparison among approaches of the same task less reliable. It also confuses users in choosing a suitable one when they want to analyze a new specimen. In addition, input–output gaps also exist with existing software, making some methods less usable. There is much room for improvement in accuracy and speed of calculation. Furthermore, it is certain that more deep learning methods will emerge for 3D reconstruction in cryo-ET. Many new artificial intelligence approaches may be further revised to 3D reconstruction cryo-ET such as methods based on a transformer, manifold learning, and representation theory. Thus, it is necessary to establish a uniform dataset, evaluation criteria, and document format requirements for each task. Cryo-ET will play an increasingly important role in biological science, and may be applied further to medicine, pharmacy, and other fields to promote human progress and the development of science.
Above all, popular deep learning methods and classical mathematical methods have improved the resolution of cryo-ET reconstruction. However, as far as we know, few methods are designed for cryo-ET reconstruction by combining mathematical theory and AI algorithms to complement each other. These kinds of approaches have been used for natural image processing (Wang et al., 2020; Zhao et al., 2021; Li et al., 2022b). The MoAMN method (Zhao et al., 2021) combined the EM algorithm and a denoiser for Gaussian noise to remove a mixture of additive and multiplicative noise. A learnable regularizer provides a good image prior to denoising images in Li et al. (2022b) by integrating the variational method into the architecture of denoiser. Similarly, an image segmentation model with adaptive similarity was developed in Wang et al. (2020). Mathematical knowledge can increase the prior knowledge of AI methods and enrich their interpretability. Meanwhile, the AI approaches improve the mathematical methods to be more flexible and fit to actual data. Therefore, combining mathematical theory and AI algorithms can be an alternative development trend for cryo-ET in the future.
AI professionals and mathematicians have devoted themselves to developing methods for cryo-ET reconstruction procedures. There is still much room to develop new ways to improve resolution besides deep learning methods and mathematical methods (as mentioned above). The combination of the two, and even the integration of biology, physics, and other fields of knowledge. We hope that this review will attract more interdisciplinary researchers to join in developing cryo-ET approaches. In particular, experts from different fields can collaborate to take full advantage of their strength.
Author contributions
This work is organized by XG. This manuscript is the collaboration of all the authors who are responsible for reading literature, reporting, and discussing the results. The Abstract and Sections 1, 3, 5, 6, and 8–12 are written by CZ. Sections 2, 4, and 7 are written by DL. Details, such as checking language and making images, are completed by CZ. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the National Natural Science Foundation of China (No. 31670725), the Beijing Advanced Innovation Center for Structural Biology of Tsinghua University, and the Public Computing Cloud of Renmin University of China, and the Beijing Academy of Intelligence.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Afonine, P., Poon, B., Read, R., Sobolev, O., Terwilliger, T., Urzhumtsev, A., et al. (2018). Real-space refinement in phenix for cryo-em and crystallography. Acta Crystallographica Section D Struct. Biol. 74, 531–544. doi: 10.1107/S2059798318006551
Aharon, M., Elad, M., Bruckstein, A. M., Katz, Y. (2005). “K-svd: An algorithm for designing of overcomplete dictionaries for sparse representation,” in IEEE transactions on signal processing (Piscataway, NJ, United States: IEEE), 4311–4322.
Albarqouni, S., Lasser, T., Alkhaldi, W., Al-Amoudi, A., Navab, N. (2015). Gradient projection for regularized cryo-electron tomographic reconstruction Vol. chap. 1 (Cham: Springer International Publishing), 43–51. doi: 10.1007/978-3-319-18431-95
Bepler, T., Kelley, K., Noble, A. J., Berger, B. (2020). Topaz-denoise: general deep denoising models for cryoem and cryoet. Nat. Commun. 11, 5208. doi: 10.1038/s41467-020-18952-1
Bhamre, T., Zhang, T., Singer, A. (2016). Denoising and covariance estimation of single particle cryo-em images. J. Struct. Biol. 195, 72–81. doi: 10.1016/j.jsb.2016.04.013
Bharat, T., Russo, C., Löwe, J., Passmore, L., Scheres, S. (2015). Advances in single-particle electron cryomicroscopy structure determination applied to sub-tomogram averaging. Structure 23, 1743–1753. doi: 10.1016/j.str.2015.06.026
Bharat, T. A., Scheres, S. H. (2016a). Sub-tomogram averaging in relion. Microscopy Microanal. 22, 1124–1125. doi: 10.1017/S1431927616006462
Bharat, T. A. M., Scheres, S. H. W. (2016b). Resolving macromolecular structures from electron cryo-tomography data using subtomogram averaging in relion. Nat. Protoc. 11, 2054–2065. doi: 10.1038/nprot.2016.124
Borgnia, M., Subramaniam, S., Milne, J. (2008). Three-dimensional imaging of the highly bent architecture of bdellovibrio bacteriovorus by using cryo-electron tomography. J. bacteriology 190, 2588–2596. doi: 10.1128/JB.01538-07
Burt, A., Gaifas, L., Dendooven, T., Gutsche, I. (2021). A flexible framework for multi-particle refinement in cryo-electron tomography. PloS Biol. 19, e3001319. doi: 10.1371/journal.pbio.3001319
Calder, L. J., Calcraft, T., Hussain, S., Harvey, R., Rosenthal, P. B. (2022). Electron cryotomography of sars-cov-2 virions reveals cylinder-shaped particles with a double layer rnp assembly. Commun. Biol. 5, 1210. doi: 10.1038/s42003-022-04183-1
Carazo, J.-M., Herman, G. T., Sorzano, C. O. S., Marabini, R. (2006). Algorithms for threedimensional reconstruction from the imperfect projection data provided by electron microscopy Vol. chap. 7 (New York, NY: Springer New York), 217–243. doi: 10.1007/978-0-387-69008-78
Castaño-Díez, D., Kudryashev, M., Arheit, M., Stahlberg, H. (2012). Dynamo: A flexible, userfriendly development tool for subtomogram averaging of cryo-em data in high-performance computing environments. J. Struct. Biol. 178, 139–151. doi: 10.1016/j.jsb.2011.12.017. Special Issue: Electron Tomography.
Castaño-Díez, D., Kudryashev, M., Stahlberg, H. (2017). Dynamo catalogue: Geometrical tools and data management for particle picking in subtomogram averaging of cryo-electron tomograms. J. Struct. Biol. 197, 135–144. doi: 10.1016/j.sbi.2019.05.011
Castaño-Díez, D., Zanetti, G. (2019). In situ structure determination by subtomogram averaging. Curr. Opin. Struct. Biol. 58, 68–75. doi: 10.1016/j.sbi.2019.05.011
Che, C., Lin, R., Zeng, X., Elmaaroufi, K., Galeotti, J., Xu, M. (2018). Improved deep learningbased macromolecules structure classification from electron cryo-tomograms. Mach. Vision Appl. 29, 1227–1236. doi: 10.1007/s00138-018-0949-4
Chen, M., Bell, J. M., Shi, X., Sun, S. Y., Wang, Z., Ludtke, S. J. (2019). A complete data processing workflow for cryo-et and subtomogram averaging. Nat. Methods 16, 1161—1168. doi: 10.1038/s41592-019-0591-8
Chen, Y., Förster, F. (2014). Iterative reconstruction of cryo-electron tomograms using nonuniform fast fourier transforms. J. Struct. Biol. 185, 309–316. doi: 10.1016/j.jsb.2013.12.001
Chen, Y., Pfeffer, S., Fernández, J., Sorzano, C., Förster, F. (2014). Autofocused 3d classification of cryoelectron subtomograms. Structure 22, 1528–1537. doi: 10.1016/j.str.2014.08.007
Chen, Y., Pfeffer, S., Hrabe, T., Schuller, J. M., Förster, F. (2013). Fast and accurate reference-free alignment of subtomograms. J. Struct. Biol. 182, 235–245. doi: 10.1016/j.jsb.2013.03.002
Chen, M., Shi, X., Yu, Z., Fan, G., Serysheva, I. I., Baker, M. L., et al. (2022). In situ structure of the acrab-tolc efflux pump at subnanometer resolution. Structure 30, 107–113.e3. doi: 10.1016/j.str.2021.08.008
Cheng, J., Li, B., Si, L., Zhang, X. (2021). Determining structures in a native environment using single-particle cryoelectron microscopy images. Innovation 2, 100166. doi: 10.1016/j.xinn.2021.100166
Croxford, M., Elbaum, M., Arigovindan, M., Kam, Z., Agard, D., Villa, E., et al. (2021). Entropyregularized deconvolution of cellular cryotransmission electron tomograms. Proc. Natl. Acad. Sci. 118, e2108738118. doi: 10.1073/pnas.2108738118
Deng, Y., Chen, Y., Zhang, Y., Wang, S., Zhang, F., Sun, F. (2016). ICON: 3d reconstruction with’missing-information’ restoration in biological electron tomography. J. Struct. Biol. 195, 100–112. doi: 10.1016/j.jsb.2016.04.004
Ding, G., Liu, Y., Zhang, R., Xin, H. L. (2019). A joint deep learning model to recover information and reduce artifacts in missing-wedge sinograms for electron tomography and beyond. Sci. Rep. 9, 12803. doi: 10.1038/s41598-019-49267-x
Ding, H. J., Oikonomou, C. M., Jensen, G. J. (2015). The caltech tomography database and automatic processing pipeline. J. Struct. Biol. 192, 279–286. doi: 10.1016/j.jsb.2015.06.016. Recent Advances in Detector Technologies and Applications for Molecular TEM.
Eibauer, M., Hoffmann, C., Plitzko, J. M., Baumeister, W., Nickell, S., Engelhardt, H. (2012). Unraveling the structure of membrane proteins in situ by transfer function corrected cryo-electron tomography. J. Struct. Biol. 180, 488–496. doi: 10.1016/j.jsb.2012.09.008
Fang, Q., Zhu, D., Agarkova, I. V., Adhikari, J., Klose, T., Liu, Y., et al. (2019). Near-atomic structure of a giant virus. Nat. Commun. 10, 1–11. doi: 10.1038/s41467-019-08319-6
Fernandez, J.-J., Li, S. (2021). Tomoalign: A novel approach to correcting sample motion and 3d ctf in cryoet. J. Struct. Biol. 213, 107778. doi: 10.1016/j.jsb.2021.107778
Fernandez, J.-J., Li, S., Agard, D. (2019). Consideration of sample motion in cryo-tomography based on alignment residual interpolation. J. Struct. Biol. 205, 1–6. doi: 10.1016/j.jsb.2019.01.005
Fernandez, J.-J., Li, S., Bharat, T. A., Agard, D. A. (2018). Cryo-tomography tilt-series alignment with consideration of the beam-induced sample motion. J. Struct. Biol. 202, 200–209. doi: 10.1016/j.jsb.2018.02.001
Fernández, J., Li, S., Crowther, R. (2006). Ctf determination and correction in electron cryotomography. Ultramicroscopy 106, 587–596. doi: 10.1016/j.ultramic.2006.02.004
Förster, F., Medalia, O., Zauberman, N., Baumeister, W., Fass, D. (2005). Retrovirus envelope protein complex structure in situ studied by cryo-electron tomography. Proc. Natl. Acad. Sci. United States America 102, 4729–4734. doi: 10.1073/pnas.0409178102
Förster, F. G., Pruggnaller, S., Seybert, A., Frangakis, A. S. (2008). Classification of cryo-electron sub-tomograms using constrained correlation. J. Struct. Biol. 161, 276–286. doi: 10.1016/j.jsb.2007.07.006
Frangakis, A. S. (2021). It’s noisy out there! a review of denoising techniques in cryo-electron tomography. J. Struct. Biol. 213, 10127–10134. doi: 10.1016/j.jsb.2021.107804
Frangakis, A. S., Boh, J., Förster, F., Nickell, S., Nicastro, D., Typke, D., et al. (2002). Identification of macromolecular complexes in cryoelectron tomograms of phantom cells. Proc. Natl. Acad. Sci. 99, 14153–14158. doi: 10.1073/pnas.172520299
Frank, J. (2005). Electron tomography: methods for three-dimensional visualization of structures in the cell (New York: Springer Science). doi: 10.1007/978-0-387-69008-7
Frappart, F., Ramillien, G., Seoane, L. (2016). “8 - monitoring water mass redistributions on land and polar ice sheets using the grace gravimetry from space mission,” in Land surface remote sensing in continental hydrology. Eds. Baghdadi, N., Zribi, M. (London, UK: Elsevier), 255–279. doi: 10.1016/B978-1-78548-104-8.50008-5
Frazier, Z., Xu, M., Alber, F. (2017). Tomominer and tomominercloud: A software platform for large-scale subtomogram structural analysis. Structure 25, 951–961.e2. doi: 10.1016/j.str.2017.04.016
Galaz-Montoya, J. G., Hecksel, C. W., Baldwin, P. R., Wang, E., Weaver, S. C., Schmid, M. F., et al. (2016). Alignment algorithms and per-particle ctf correction for single particle cryo-electron tomography. J. Struct. Biol. 194, 383–394. doi: 10.1016/j.jsb.2016.03.018
Ganguly, P. S., Lucka, F., Kohr, H., Franken, E., Hupkes, H. J., Batenburg, K. J. (2022). Sparsealign: A grid-free algorithm for automatic marker localization and deformation estimation in cryo-electron tomography. IEEE Trans. Comput. Imaging 8, 651–665. doi: 10.1109/TCI.2022.3194719
George, B., Assaiya, A., Jacob Roy, R., Kembhavi, A., Chauhan, R., Paul, G., et al. (2021). Cassper is a semantic segmentation-based particle picking algorithm for single-particle cryo-electron microscopy. Commun. Biol. 4, 1–12. doi: 10.1038/s42003-021-01721-1
Giri, N., Roy, R., Cheng, J. (2022). Deep learning for reconstructing protein structures from cryo-em density maps: recent advances and future directions. Curr. Opin. Struct. Biol. 79 (102536), 1–9. doi: 10.1016/j.sbi.2023.102536
Gordon, R., Bender, R., Herman, G. T. (1970). Algebraic reconstruction techniques (art) for threedimensional electron microscopy and x-ray photography. J. Theor. Biol. 29, 471–481. doi: 10.1016/0022-5193(70)90109-8
Grant, T., Grigorieff, N. (2015). Measuring the optimal exposure for single particle cryo-em using a 2.6 a˚ reconstruction of rotavirus vp6. elife 4, e06980. doi: 10.7554/eLife.06980
Gu, H., Xian, Y., Unarta, I. C., Yao, Y. (2023). “Generative adversarial networks for robust cryo-EM image denoising,” in Handbook of mathematical models and algorithms in computer vision and imaging: mathematical imaging and vision. Eds. Chen, K., Schönlieb, C.-B., Tai, X.-C., Younes, L. (Cham, Switzerland: Springer International Publishing), 969–1000. doi: 10.1007/978-3-030-98661-2_126
Gubins, I., Chaillet, M. L., Schot, G., Trueba, M. C., Veltkamp, R. C., Förster, F., et al. (2021). “SHREC 2021: classification in cryo-electron tomograms,” in Eurographics workshop on 3D object retrieval. Eds. Biasotti, S., Dyke, R. M., Lai, Y., Rosin, P. L., Veltkamp, R. C. (Eindhoven, Netherlands: The Eurographics Association), 5–17. doi: 10.2312/3dor.20211307
Guesdon, A., Blestel, S., Kervrann, C., Chrétien, D. (2013). Single versus dual-axis cryo-electron tomography of microtubules assembled in vitro: Limits and perspectives. J. Struct. Biol. 181, 169–178. doi: 10.1016/j.jsb.2012.11.004
Hagen, W. J., Wan, W., Briggs, J. A. (2017). Implementation of a cryo-electron tomography tilt-scheme optimized for high resolution subtomogram averaging. J. Struct. Biol. 197, 191–198. doi: 10.1016/j.jsb.2016.06.007
Hajarolasvadi, N., Sunkara, V., Khavnekar, S., Beck, F., Brandt, R., Baum, D. (2022). Volumetric macromolecule identification in cryo-electron tomograms using capsule networks. BMC Bioinf. 23, 1–16. doi: 10.1186/s12859-022-04901-w
Han, R., Bao, Z., Zeng, X., Niu, T., Zhang, F., Xu, M., et al. (2019a). A joint method for markerfree alignment of tilt series in electron tomography. Bioinf. (Oxford England) 35, i249–i259. doi: 10.1093/bioinformatics/btz323
Han, R., Li, L., Yang, P., Zhang, F., Gao, X. (2021). A novel constrained reconstruction model towards high-resolution subtomogram averaging. Bioinformatics 37, 1616–1626. doi: 10.1093/bioinformatics/btz787
Han, R., Wan, X., Li, L., Lawrence, A. F., Yang, P., Li, Y., et al. (2019b). Autom-dualx: a toolkit for fully automatic fiducial marker-based alignment of dual-axis tilt series with simultaneous reconstruction. Bioinformatics 35, 319–328. doi: 10.1093/bioinformatics/bty620
Han, R., Wan, X., Wang, Z., Hao, Y., Zhang, J., Chen, Y., et al. (2017). Autom: A novel automatic platform for electron tomography reconstruction. J. Struct. Biol. 199, 196–208. doi: 10.1016/j.jsb.2017.07.008
Han, R., Wang, L., Liu, Z., Sun, F., Zhang, F. (2015). A novel fully automatic scheme for fiducial marker-based alignment in electron tomography. J. Struct. Biol. 192, 403–417. doi: 10.1016/j.jsb.2015.09.022
Han, X., Zhang, M., Liu, J., Wang, J., Ren, G. (2022). Cryo-electron tomography related radiationdamage parameters for individual-molecule 3d structure determination. Front. Chem. 10. doi: 10.3389/fchem.2022.889203
Han, R., Zhang, F., Wan, X., Fernandez, J.-J., Sun, F., Liu, Z. (2014). A marker-free automatic alignment method based on scale-invariant features. J. Struct. Biol. 186, 167–180. doi: 10.1016/j.jsb.2014.02.011
Hanssen, E. (2018). Cellular imaging, electron tomography and related techniques (Switzerland: Springer Nature). Available at: https://link.springer.com/content/pdf/bfm. doi: 10.1007/978-3-319-68997-5
Hao, Y., Wan, X., Yan, R., Liu, Z., Li, J., Zhang, S., et al. (2022). Vp-detector: A 3d multi-scale dense convolutional neural network for macromolecule localization and classification in cryo-electron tomograms. Comput. Methods Programs Biomedicine 221, 106871. doi: 10.1016/j.cmpb.2022.106871
Harastani, M., Vuillemot, R., Hamitouche, I., Moghadam, N. B., Jonic, S. (2022). Continuousflex: Software package for analyzing continuous conformational variability of macromolecules in cryo electron microscopy and tomography data. J. Struct. Biol. 214, 107906. doi: 10.1016/j.jsb.2022.107906
Henderson, R. (2013). Avoiding the pitfalls of single particle cryo-electron microscopy: Einstein from noise. Proc. Natl. Acad. Sci. U. S. A. 110, 18037–18041. doi: 10.1073/pnas.1314449110
Heumann, J. M., Hoenger, A., Mastronarde, D. N. (2011). Clustering and variance maps for cryoelectron tomography using wedge-masked differences. J. Struct. Biol. 175, 288–299. doi: 10.1016/j.jsb.2011.05.011
Heymann, J. B. (2018). Guidelines for using bsoft for high resolution reconstruction and validation of biomolecular structures from electron micrographs. Protein Sci. 27, 159–171. doi: 10.1002/pro.3293
Hrabe, T., Chen, Y., Pfeffer, S., Kuhn Cuellar, L., Mangold, A.-V., Förster, F. (2012). Pytom: A python-based toolbox for localization of macromolecules in cryo-electron tomograms and subtomogram analysis. J. Struct. Biol. 178, 177–188. doi: 10.1016/j.jsb.2011.12.003
Huang, C.-Y., Draczkowski, P., Wang, Y.-S., Chia-Yu, C., Chien, Y.-C., Yun-han, C., et al. (2022). In situ structure and dynamics of an alphacoronavirus spike protein by cryo-et and cryo-em. Nat. Commun. 13, 4877. doi: 10.1038/s41467-022-32588-3
Igaev, M., Kutzner, C., Bock, L. V., Vaiana, A. C., Grubmüller, H. (2019). Automated cryo-em structure refinement using correlation-driven molecular dynamics. eLife 8, e43542. doi: 10.7554/eLife.43542
Iudin, A., Korir, P. K., Salavert-Torres, J., Kleywegt, G. J., Patwardhan, A. (2016). Empiar: a public archive for raw electron microscopy image data. Nat. Methods 13, 387—388. doi: 10.1038/nmeth.3806
Jakobi, A. J., Wilmanns, M., Sachse, C. (2017). Model-based local density sharpening of cryo-em maps. eLife 6, 1–26. doi: 10.7554/eLife.27131.036
Jenni, S., Horwitz, J. A., Bloyet, L.-M., Whelan, S. P. J., Harrison, S. C. (2022). Visualizing molecular interactions that determine assembly of a bullet-shaped vesicular stomatitis virus particle. Nat. Commun. 13, 4802. doi: 10.1038/s41467-022-32223-1
Jensen, G., Kornberg, R. (2000). Defocus-gradient corrected back-projection. Ultramicroscopy 84, 57–64. doi: 10.1016/S0304-3991(00)00005-X
Jiménez-Ortigosa, C., Jiang, J., Chen, M., Xuyuan, K., Healey, K., Castellano, P., et al. (2021). Preliminary structural elucidation of β-(1,3)-glucan synthase from candida glabrata using cryo-electron tomography. J. Fungi 7, 120. doi: 10.3390/jof7020120
Kaur, S., Gomez-Blanco, J., Khalifa, A. A. Z., Adinarayanan, S., Sanchez-Garcia, R., Wrapp, D., et al. (2021). Local computational methods to improve the interpretability and analysis of cryo-em maps. Nat. Commun. 12, 1240. doi: 10.1038/s41467-021-21509-5
Khanna, K., Villa, E. (2022). Revealing bacterial cell biology using cryo-electron tomography. Curr. Opin. Struct. Biol. 75, 102419. doi: 10.1016/j.sbi.2022.102419
Kileel, J., Moscovich, A., Zelesko, N., Singer, A. (2021). Manifold learning with arbitrary norms. J. Fourier Anal. Appl. 27, 1–56. doi: 10.1007/s00041-021-09879-2
Klein, S., Cortese, M., Winter, S., Wachsmuth-Melm, M., Neufeldt, C., Cerikan, B., et al. (2020). Sars-cov-2 structure and replication characterized by in situ cryo-electron tomography. Nat. Commun. 11, 5885. doi: 10.1038/s41467-020-19619-7
Krull, A., Buchholz, T.-O., Jug, F. (2019). “Noise2void - learning denoising from single noisy images,” in 2019 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) (Long Beach, CA, United States: IEEE). 2124–2132. doi: 10.1109/CVPR.2019.00223
Kunz, M., Frangakis, A. S. (2017). Three-dimensional ctf correction improves the resolution of electron tomograms. J. Struct. Biol. 197, 114–122. doi: 10.1016/j.jsb.2016.06.016. Electron Tomography.
Kupsch, A., Lange, A., Hentschel, M., Lück, S., Schmidt, V., Grothausmann, R., et al. (2015). Missing wedge computed tomography by iterative algorithm DIRECTT. J. microscopy 261, 36–45. doi: 10.1111/jmi.12313
Lamm, L., Righetto, R. D., Wietrzynski, W., Pöge, M., Martinez-Sanchez, A., Peng, T., et al. (2022). Membrain: A deep learning-aided pipeline for detection of membrane proteins in cryoelectron tomograms. Comput. Methods Programs Biomedicine 224, 106990. doi: 10.1016/j.cmpb.2022.106990
Lawson, C. L., Patwardhan, A., Baker, M. L., Hryc, C. F., Sanz-García, E., Hudson, B. P., et al. (2016). Emdatabank: unified data resource for 3dem. Nucleic Acids Res. 44, D396–D403. doi: 10.1093/nar/gkv1126
Lehtinen, J., Munkberg, J., Hasselgren, J., Laine, S., Karras, T., Aittala, M., et al. (2018). Noise2noise: Learning image restoration without clean data (Stockholm, Sweden: ICML), 2965–2974.
Leigh, K. E., Navarro, P. P., Scaramuzza, S., Chen, W., Zhang, Y., Castaño-Díez, D., et al. (2019). “Chapter 11 - subtomogram averaging from cryo-electron tomograms,” in Three-dimensional electron microscopy, vol. 152 . Eds. Müller-Reichert, T., Pigino, G. (San Diego, CA, United States: Academic Press), 217–259. of Methods in Cell Biology. doi: 10.1016/bs.mcb.2019.04.003
Li, S. (2022). Cryo-electron tomography of enveloped viruses. Trends Biochem. Sci. 47, 173–186. doi: 10.1016/j.tibs.2021.08.005. Special Issue: Pushing boundaries of cryo-EM.
Li, X., Mooney, P., Zheng, S., Booth, C. R., Braunfeld, M. B., Gubbens, S., et al. (2013). Electron counting and beam-induced motion correction enable near-atomic-resolution single-particle cryo-em. Nat. Methods 10, 584–590. doi: 10.1038/nmeth.2472
Li, Z., Wang, F., Cui, L., Liu, J. (2022b). Dual mixture model based cnn for image denoising. IEEE Trans. Image Process. 31, 3618–3629. doi: 10.1109/TIP.2022.3173814
Li, H., Zhang, H., Wan, X., Yang, Z., Li, C., Li, J., et al. (2022a). Noise-Transfer2Clean: denoising cryo-EM images based on noise modeling and transfer. Bioinformatics 38, 2022–2029. doi: 10.1093/bioinformatics/btac052
Liebschner, D., Afonine, P., Moriarty, N., Poon, B., Chen, V., Adams, P. (2021). Ceres: A cryo-em re-refinement system for continuous improvement of deposited models. Acta Crystallographica Section D Struct. Biol. 77, 48–61. doi: 10.1107/S2059798320015879
Liu, G., Luo, Y., Yang, G. (2022a). “3d particle picking in cryo-electron tomograms using instance segmentation,” in 2022 IEEE International Conference on Image Processing (ICIP) (Bordeaux, France: IEEE). 2157–2161. doi: 10.1109/ICIP46576.2022.9897829
Liu, Y.-T., Zhang, H., Wang, H., Tao, C., Bi, G.-Q., Zhou, H. (2022b). Isotropic reconstruction of electron tomograms with deep learning. Nat. Commun. 13, 6482. doi: 10.1101/2021.07.17.452128
López-Jiménez, A. T., Mostowy, S. (2021). Emerging technologies and infection models in cellular microbiology. Nat. Commun. 12, 6764. doi: 10.1038/s41467-021-26641-w
Lü, Y., Zeng, X., Zhao, X., Li, S., Li, H., Gao, X., et al. (2019). Fine-grained alignment of cryoelectron subtomograms based on mpi parallel optimization. BMC Bioinf. 20, 443. doi: 10.1186/s12859-019-3003-2
Luo, Z., Zeng, X., Bao, Z., Xu, M. (2019). “Deep learning-based strategy for macromolecules classification with imbalanced data from cellular electron cryotomography,” in 2019 International Joint Conference on Neural Networks (IJCNN) (Budapest, Hungary: IEEE). 1–8.
Ma, H., Pham, P., Luo, B., Rangan, R., Kappel, K., Su, Z., et al. (2022). Auto-DRRAFTER: automated RNA modeling based on cryo-EM density Vol. chap. 2568 (New York: Springer US), 193–211. doi: 10.1007/978-1-0716-2687-013
Maggioni, M., Katkovnik, V., Egiazarian, K., Foi, A. (2013). Nonlocal transform-domain filter for volumetric data denoising and reconstruction. IEEE Trans. Image Process. 22, 119–133. doi: 10.1109/TIP.2012.2210725
Marabini, R., Carragher, B., Chen, S., Chen, J., Cheng, A., Downing, K. H., et al. (2015). Ctf challenge: Result summary. J. Struct. Biol. 190, 348–359. doi: 10.1016/j.jsb.2015.04.003
Marshall, N. F., Lan, T.-Y., Bendory, T., Singer, A. (2020). “Image recovery from rotational and translational invariants,” in ICASSP 2020 - 2020 IEEE international conference on acoustics, speech and signal processing (ICASSP) (Barcelona, Spain: IEEE), 5780–5784. doi: 10.1109/ICASSP40776.2020.9053932
Mastronarde, D. N. (1997). Dual-axis tomography: An approach with alignment methods that preserve resolution. J. Struct. Biol. 120, 343–352. doi: 10.1006/jsbi.1997.3919
Mastronarde, D. N., Held, S. R. (2017). Automated tilt series alignment and tomographic reconstruction in imod. J. Struct. Biol. 197, 102–113. doi: 10.1016/j.jsb.2016.07.011. Electron Tomography.
McLeod, R. A., Kowal, J., Ringler, P., Stahlberg, H. (2017). Robust image alignment for cryogenic transmission electron microscopy. J. Struct. Biol. 197, 279–293. doi: 10.1016/j.jsb.2016.12.006
Mendonҫa, L., Sun, D., Ning, J., Liu, J., Kotecha, A., Olek, M., et al. (2021). Cryoet structures of immature hiv gag reveal six-helix bundle. Commun. Biol. 4, 481. doi: 10.1038/s42003-021-01999-1
Metskas, L. A., Wilfong, R., Jensen, G. J. (2022). Subtomogram averaging for biophysical analysis and supramolecular context. J. Struct. Biology: X 6, 100076. doi: 10.1016/j.yjsbx.2022.100076
Mo, M., Kong, F., Liu, Q. (2021). “Particle picking method for cryo electron tomography image based on active learning,” in Web information systems and applications (Kaifeng, China: Springer-Verlag), 468–479. doi: 10.1007/978-3-030-87571-840
Moebel, E., Kervrann, C. (2018). “A monte carlo framework for denoising and missing wedge reconstruction in cryo-electron tomography,” in Patch-based techniques in medical imaging (Cham: Springer International Publishing), 28–35.
Moebel, E., Kervrann, C. (2020). A monte carlo framework for missing wedge restoration and noise removal in cryo-electron tomography. J. Struct. Biology: X 4, 100013. doi: 10.1016/j.yjsbx.2019.100013
Moebel, E., Martínez-Sánchez, A., Lamm, L., Righetto, R., Wietrzynski, W., Albert, S., et al. (2021). Deep learning improves macromolecule identification in 3d cellular cryo-electron tomograms. Nat. Methods 18, 1–9. doi: 10.1038/s41592-021-01275-4
Morena, J., Conesa, P., Fonseca-Reyna, Y., Isidro-Gómez, F., Herreros, D., Gimenez Fernandez, E., et al. (2022). Scipiontomo: Towards cryo-electron tomography software integration, reproducibility, and validation. J. Struct. Biol. 214, 107872. doi: 10.1016/j.jsb.2022.107872
Murray, S., Galaz-Montoya, J. G., Tang, G., Flanagan, J., Ludtke, S. J. (2014). Eman2.1 - a new generation of software for validated single particle analysis and single particle tomography. Microscopy Microanalysis 20, 832–833. doi: 10.1017/S1431927614005881
Nans, A., Saibil, H. R., Hayward, R. D. (2014). Pathogen-host reorganization during chlamydia invasion revealed by cryo-electron tomography. Cell. Microbiol. 16, 1457—1472. doi: 10.1111/cmi.12310
Naydenova, K., Jia, P., Russo, C. J. (2020). Cryo-em with sub–1 a˚ specimen movement. Science 370, 223–226. doi: 10.1126/science.abb7927
Neijenhuis, T., van Keulen, S. C., Bonvin, A. M. J. (2022). Interface refinement of low- to mediumresolution cryo-em complexes using haddock2.4. Structure 30, 476–484.e3. doi: 10.1016/j.str.2022.02.001
Ni, T., Frosio, T., Mendonҫa, L., Sheng, Y., Clare, D., Himes, B. A., et al. (2022). High-resolution in situ structure determination by cryo-electron tomography and subtomogram averaging using emclarity. Nat. Protoc. 17, 421—444. doi: 10.1038/s41596-021-00648-5
Nicastro, D., Schwartz, C., Pierson, J., Gaudette, R., Porter, M. E., McIntosh, J. R. (2006). The molecular architecture of axonemes revealed by cryoelectron tomography. Science 313, 944–948. doi: 10.1126/science.1128618
Obr, M., Schur, F. K. (2019). “Chapter five - structural analysis of pleomorphic and asymmetric viruses using cryo-electron tomography and subtomogram averaging,” in Complementary strategies to understand virus structure and function, vol. 105 . Ed. Rey, F. A. (San Diego, CA, United States: Academic Press), 117–159. Advances in Virus Research. doi: 10.1016/bs.aivir.2019.07.008
Paavolainen, L., Acar, E., Tuna, U., Peltonen, S., Moriya, T., Soonsawad, P., et al. (2014). Compensation of missing wedge effects with sequential statistical reconstruction in electron tomography. PloS One 9, 1–14. doi: 10.1371/journal.pone.0108978
Pang, Z.-F., Zhang, H.-L., Luo, S., Zeng, T. (2020). Image denoising based on the adaptive weighted tvp regularization. Signal Process. 167, 107325. doi: 10.1016/j.sigpro.2019.107325
Pei, L., Xu, M., Frazier, Z., Alber, F. (2016). Simulating cryo electron tomograms of crowded cell cytoplasm for assessment of automated particle picking. BMC Bioinf. 17, 405. doi: 10.1186/s12859-016-1283-3
Penczek, P. A., Fang, J., Li, X., Cheng, Y., Loerke, J., Spahn, C. M. (2014). CTER—rapid estimation of ctf parameters with error assessment. Ultramicroscopy 140, 9–19. doi: 10.1016/j.ultramic.2014.01.009
Plitzko, J., Baumeister, W. P. (2019). “Cryo-electron yomography,” in Springer handbook of microscopy. Eds. Hawkes, P. W., Spence, J. C. H. (Springer, Cham: Springer Handbooks), 189–228. doi: 10.1007/978-3-030-00069-1_4
Powell, B. M., Davis, J. (2023). Learning structural heterogeneity from cryo-electron sub-tomograms with tomodrgn. bioRxiv, 1-29. doi: 10.1101/2023.05.31.542975
Pyle, E., Hutchings, J., Zanetti, G. (2022). Strategies for picking membrane-associated particles within subtomogram averaging workflows. Faraday Discussions 240, 101–113. doi: 10.1039/D2FD00022A
Pyle, E., Zanetti, G. (2021). Current data processing strategies for cryo-electron tomography and subtomogram averaging. Biochem. J. 478, 1827–1845. doi: 10.1042/BCJ20200715
Ramírez-Aportela, E., Vilas, J. L., Glukhova, A., Melero, R., Conesa, P., Martínez, M., et al. (2019). Automatic local resolution-based sharpening of cryo-em maps. Bioinformatics 36, 765–772. doi: 10.1093/bioinformatics/btz671
Rice, G., Wagner, T., Stabrin, M., Raunser, S. (2023). Tomotwin: Generalized 3d localization of macromolecules in cryo-electron tomograms with structural data mining. Nat. Methods 20, 871—880. doi: 10.1038/s41592-023-01878-z
Rohou, A., Grigorieff, N. (2015). Ctffind4: Fast and accurate defocus estimation from electron micrographs. J. Struct. Biol. 192, 216–221. doi: 10.1016/j.jsb.2015.08.008. Recent Advances in Detector Technologies and Applications for Molecular TEM.
Rosenthal, P. B., Henderson, R. (2003). Optimal determination of particle orientation, absolute hand, and contrast loss in single-particle electron cryomicroscopy. J. Mol. Biol. 333, 721–745. doi: 10.1016/j.jmb.2003.07.013
Rubinstein, J. L., Brubaker, M. A. (2015). Alignment of cryo-em movies of individual particles by optimization of image translations. J. Struct. Biol. 192, 188–195. doi: 10.1016/j.jsb.2015.08.007
Sanchez-Garcia, R., Gomez-Blanco, J., Cuervo, A., Carazo, J., Sorzano, C., Vargas, J. (2021). Deepemhancer: A deep learning solution for cryo-em volume post-processing. Commun. Biol. 4, 874. doi: 10.1101/2020.06.12.148296
Sanchez-Garcia, R., Segura, J., Maluenda, D., Carazo, J. M., Sorzano, C. O. S. (2018). Deep consensus, a deep learning-based approach for particle pruning in cryo-electron microscopy. IUCrJ 5, 854–865. doi: 10.1107/S2052252518014392
Scaramuzza, S., Castaño-Díez, D. (2021). Step-by-step guide to efficient subtomogram averaging of virus-like particles with dynamo. PloS Biol. 19, 1–22. doi: 10.1371/journal.pbio.3001318
Seifer, S., Elbaum, M. (2022). Clusteralign: A fiducial tracking and tilt series alignment tool for thick sample tomography. Biol. Imaging 2, e7. doi: 10.1017/S2633903X22000071
Seybert, A., Gonzalez-Gonzalez, L., Scheffer, M., Lluch, M., Martinez Mariscal, A., Querol, E., et al. (2018). Cryo-electron tomography analyses of terminal organelle mutants suggest the motility mechanism of mycoplasma genitalium. Mol. Microbiol. 108, 57. doi: 10.1111/mmi.13938
Shi, Y., Singer, A. (2022). Ab-initio contrast estimation and denoising of cryo-em images. Comput. Methods Prog. Biomed. 224, 107018. doi: 10.1016/j.cmpb.2022.107018
Singer, A. (2019). “Mathematics for cryo-electron microscopy,” in 2018 international congress of mathematicians, ICM 2018 (Singapore: World Scientific Publishing Co. Pte Ltd), 3995–4014. doi: 10.1142/97898132728800209
Stölken, M., Beck, F., Haller, T., Hegerl, R., Gutsche, I., Carazo, J.-M., et al. (2011). Maximum likelihood based classification of electron tomographic data. J. Struct. Biol. 173, 77–85. doi: 10.1016/j.jsb.2010.08.005
Strelak, D., Jiménez-Moreno, A., Vilas, J. L., Ramírez-Aportela, E., Sanchez-García, R., Maluenda, D., et al. (2021). Advances in xmipp for cryo–electron microscopy: From xmipp to scipion. Molecules 26, 6224. doi: 10.3390/molecules26206224
Su, M. (2019). goctf: Geometrically optimized ctf determination for single-particle cryo-em. J. Struct. Biol. 205, 22–29. doi: 10.1016/j.jsb.2018.11.012
Tang, G., Peng, L., Baldwin, P. R., Mann, D. S., Jiang, W., Rees, I., et al. (2007). Eman2: An extensible image processing suite for electron microscopy. J. Struct. Biol. 157, 38–46. doi: 10.1016/j.jsb.2006.05.009
Taylor, K., Liu, J., Winkler, H. (2006). Localization and classification of repetitive structures in electron tomograms of paracrystalline assemblies Vol. chap. 7 (New York: Springer), 417–439.
Tegunov, D., Cramer, P. (2019). Real-time cryo-electron microscopy data preprocessing with warp. Nat. Methods 16, 1146–1152. doi: 10.1038/s41592-019-0580-y
Tegunov, D., Xue, L., Dienemann, C., Cramer, P., Mahamid, J. (2021). Multi-particle cryo-em refinement with m visualizes ribosome-antibiotic complex at 3.5 a˚ in cells. Nat. Methods 18, 186–193. doi: 10.1038/s41592-020-01054-7
Terwilliger, T. C. (2000). Maximum-likelihood density modification. Acta Crystallographica Section D: Biol. Crystallogr. 56, 965 –972. doi: 10.1107/S0907444900005072
Terwilliger, T. (2001). Maximum-likelihood density modification using pattern recognition of structural motifs. Acta crystallographica. Section D Biol. Crystallogr. 57, 1755—1762. doi: 10.1107/s0907444901013737
Terwilliger, T. C., Ludtke, S. J., Read, R. J., Adams, P. D., Afonine, P. V. (2019). Improvement of cryo-em maps by density modification. Nat. Methods 17, 923-927. doi: 10.1101/845032
Terwilliger, T. C., Sobolev, O. V., Afonine, P. V., Adams, P. D. (2018). Automated map sharpening by maximization of detail and connectivity. Acta Crystallographica. Section D Struct. Biol. 74, 545–559. doi: 10.1107/S2059798318004655
Terwilliger, T. C., Sobolev, O. V., Afonine, P. V., Adams, P. D., Read, R. J. (2020). Density modification of cryo-em maps. Acta Crystallographica. Section D. Struct. Biol. 76, 912–925. doi: 10.1107/S205979832001061X
Thorne, R. E. (2020). Hypothesis for a mechanism of beam-induced motion in cryo-electron microscopy. IUCrJ 7, 416–421. doi: 10.1107/S2052252520002560
Tong, J., Arslan, I., Midgley, P. (2006). A novel dual-axis iterative algorithm for electron tomography. J. Struct. Biol. 153, 55—63. doi: 10.1016/j.jsb.2005.10.005
Trampert, J., Leveque, J.-J. (1990). Simultaneous iterative reconstruction technique: Physical interpretation based on the generalized least squares solution. J. Geophysical Research: Solid Earth 95, 12553–12559. doi: 10.1029/JB095iB08p12553
Trampert, P., Wang, W., Chen, D., Ravelli, R. B., Dahmen, T., Peters, P. J., et al. (2018). Exemplar-based inpainting as a solution to the missing wedge problem in electron tomography. Ultramicroscopy 191, 1–10. doi: 10.1016/j.ultramic.2018.04.001
Turoňová, B., Marsalek, L., Davidovic, T., Slusallek, P. (2015). Progressive stochastic reconstruction technique (psrt) for cryo electron tomography. J. Struct. Biol. 189, 195–206. doi: 10.1016/j.jsb.2015.01.011
Turonova, B., Schur, F., Wan, W., Briggs, J. (2017). Efficient 3d-ctf correction for cryo-electron tomography using novactf improves subtomogram averaging resolution to 3.4 A.˚. J. Struct. Biol. 199, 187–195. doi: 10.1016/j.jsb.2017.07.007
Turoňová, B., Sikora, M., Schürmann, C., Hagen, W. J. H., Welsch, S., Blanc, F. E. C., et al. (2020). In situ structural analysis of sars-cov-2 spike reveals flexibility mediated by three hinges. Sci. (New York N.y.) 370, 203–208. doi: 10.1101/2020.06.26.173476
Urzhumtsev, A. G., Urzhumtseva, L. M., Lunin, V. Y. (2022). Direct calculation of cryo-em and crystallographic model maps for real-space refinement. Acta Crystallographica Section D 78, 1451–1468. doi: 10.1107/S2059798322010907
Voortman, L. M. (2014). “Ctf correction in cryo-electron tomography,” in Ridderprint. Ed. Ridderkerk, B. V. (Netherlands: Technische Universiteit Delft), 1–135.
Vulović, M., Ravelli, R. B., van Vliet, L. J., Koster, A. J., Lazić, I., Lücken, U., et al. (2013). Image formation modeling in cryo-electron microscopy. J. Struct. Biol. 183, 19–32. doi: 10.1016/j.jsb.2013.05.008
Wan, W., Briggs, J. (2016). “Chapter thirteen - cryo-electron tomography and subtomogram averaging,” in The resolution revolution: recent advances in cryoEM, vol. 579 . Ed. Crowther, R. (San Diego, CA, United States: Academic Press), 329–367. Methods in Enzymology. doi: 10.1016/bs.mie.2016.04.014
Wan, W., Khavnekar, S., Wagner, J., Erdmann, P. S., Baumeister, W. (2020). Stopgap: A software package for subtomogram averaging and refinement. Microscopy Microanalysis 26, 2516–2516. doi: 10.1017/S143192762002187X
Wan, X., Zhang, F., Chu, Q., Zhang, K., Sun, F., Yuan, B., et al. (2011). Three-dimensional reconstruction using an adaptive simultaneous algebraic reconstruction technique in electron tomography. J. Struct. Biol. 175, 277–287. doi: 10.1016/j.jsb.2011.06.002
Wan, X., Zhang, F., Liu, Z. (2009). “Modified simultaneous algebraic reconstruction technique and its parallelization in cryo-electron tomography,” in 2009 15th International Conference on Parallel and Distributed Systems (Shenzhen, China: IEEE). 384–390. doi: 10.1109/ICPADS.2009.21
Wang, F., Zhao, C., Liu, J., Huang, H. (2020). A variational image segmentation model based on norMalized cut with adaptive similarity and spatial regularization. SIAM J. Imaging Sci. 13, 651–684. doi: 10.1137/18M1192366
Winter, S. L., Chlanda, P. (2021). Dual-axis volta phase plate cryo-electron tomography of ebola virus-like particles reveals actin-vp40 interactions. J. Struct. Biol. 213, 107742. doi: 10.1016/j.jsb.2021.107742
Wolf, D., Lubk, A., Lichte, H. (2014). Weighted simultaneous iterative reconstruction technique for single-axis tomography. Ultramicroscopy 136, 15–25. doi: 10.1016/j.ultramic.2013.07
Wu, S., Liu, G., Yang, G. (2022). “Fast particle picking for cryo-electron tomography using onestage detection,” in 2022 IEEE 19th international symposium on biomedical imaging (ISBI) (Kolkata, India: IEEE), 1–5. doi: 10.1109/ISBI52829.2022.9761580
Xu, M., Alber, F. (2012). High precision alignment of cryo-electron subtomograms through gradientbased parallel optimization. BMC Syst. Biol. 6 Suppl 1, S18. doi: 10.1186/1752-0509-6-S1-S18
Xu, M., Beck, M., Alber, F. (2012). High-throughput subtomogram alignment and classification by fourier space constrained fast volumetric matching. J. Struct. Biol. 178, 152–164. doi: 10.1016/j.jsb.2012.02.014. Special Issue: Electron Tomography.
Xu, M., Chai, X., Muthakana, H., Liang, X., Yang, G., Zeev-Ben-Mordehai, T., et al. (2017). Deep learning-based subdivision approach for large scale macromolecules structure recovery from electron cryo tomograms. Bioinformatics 33, i13–i22. doi: 10.1093/bioinformatics/btx230
Xu, C. S., Pang, S., Shtengel, G., Müller, A., Ritter, A. T., Hoffman, H. K., et al. (2021). An open-access volume electron microscopy atlas of whole cells and tissues. Nature 599, 147–151. doi: 10.1101/2020.11.13.382457
Xu, K., Wang, Z., Shi, J., Li, H., Zhang, Q. C. (2019). A2-net: Molecular structure estimation from cryo-em density volumes. Proc. AAAI Conf. Artif. Intell. 33, 1230–1237. doi: 10.1609/aaai.v33i01.33011230
Xue, L., Lenz, S., Zimmermann-Kogadeeva, M., Tegunov, D., Cramer, P., Bork, P., et al. (2022). Visualizing translation dynamics at atomic detail inside a bacterial cell. Nature 610, 205–211. doi: 10.1038/s41586-022-05255-2
Yamashita, K., Palmer, C. M., Burnley, T., Murshudov, G. N. (2021). Cryo-em single-particle structure refinement and map calculation using servalcat. Acta Crystallographica. Section D Struct. Biol. 77, 1282–1291. doi: 10.1107/S2059798321009475
Yang, Z., Zhang, F., Han, R. (2021). “Self-supervised cryo-electron tomography volumetric image restoration from single noisy volume with sparsity constraint,” in 2021 IEEE/CVF International Conference on Computer Vision (ICCV) (Montreal, QC, Canada: IEEE). 4036–4045. doi: 10.1109/ICCV48922.2021.00402
Zeng, G. (2012). A filtered backprojection algorithm with characteristics of the iterative landweber algorithm. Med. Phys. 39, 603–607. doi: 10.1118/1.3673956
Zeng, X., Howe, G., Xu, M. (2021). “End-to-end robust joint unsupervised image alignment and clustering,” in 2021 IEEE/CVF International Conference on Computer Vision (ICCV) (Montreal, QC, Canada: IEEE). 3834–3846. doi: 10.1109/ICCV48922.2021.00383
Zeng, X., Lin, Z., Uddin, M., Zhou, B., Cheng, C., Zhang, J., et al. (2022). Structure detection in three-dimensional cellular cryoelectron tomograms by reconstructing two-dimensional annotated tilt series. J. Comput. Biol. 29, 932–941. doi: 10.1089/cmb.2021.0606
Zeng, X., Xu, M. (2020). “Gum-net: Unsupervised geometric matching for fast and accurate 3d subtomogram image alignment and averaging,” in 2020 IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) (Seattle, WA, United States: IEEE). 4072–4083. doi: 10.1109/CVPR42600.2020.00413
Zhai, X., Lei, D., Zhang, M., Liu, J., Wu, H., Yu, Y., et al. (2021). Lottor: An algorithm for missing-wedge correction of the low-tilt tomographic 3d reconstruction of a single-molecule structure (lottor) v1.0. Sci. Rep. 10, 1–17. doi: 10.11578/dc.20211122.1
Zhang, K. (2016). Gctf: Real-time ctf determination and correction. J. Struct. Biol. 193, 1–12. doi: 10.1016/j.jsb.2015.11.003
Zhang, P. (2019). Advances in cryo-electron tomography and subtomogram averaging and classification. Curr. Opin. Struct. Biol. 58, 249–258. doi: 10.1016/j.sbi.2019.05.021
Zhang, P., Mendonҫa, L., Howe, A., Gilchrist, J., Sun, D., Knight, M., et al. (2021). Correlative multi-scale cryo-imaging unveils sars-cov-2 assembly and egress. Nat. Commun. 12, 4629. doi: 10.21203/rs.3.rs-134794/v1
Zhang, X., Zhang, B., Freddolino, P., Zhang, Y. (2022). Cr-i-tasser: assemble protein structures from cryo-em density maps using deep convolutional neural networks. Nat. Methods 19, 1–10. doi: 10.1038/s41592-021-01389-9
Zhao, C., Liu, J., Gong, X. (2022). “An adaptive variational model for multireference alignment with mixed noise,” in 2022 IEEE international conference on bioinformatics and biomedicine (BIBM) (Las Vegas, NV, USA: IEEE), 692–699. doi: 10.1109/BIBM55620.2022.9995092
Zhao, C., Liu, J., Zhang, J. (2021). A dual model for restoring images corrupted by mixture of additive and multiplicative noise. IEEE Access PP, 168869–168888. doi: 10.1109/ACCESS.2021.3137995
Zhao, G., Zhou, B., Wang, K., Jiang, R., Xu, M. (2018). “Respond-cam: Analyzing deep models for 3d imaging data by visualizations,” in Medical image computing and computer assisted intervention – MICCAI 2018 (Cham: Springer International Publishing), 485–492.
Zheng, S. Q., Keszthelyi, B., Branlund, E., Lyle, J. M., Braunfeld, M. B., Sedat, J. W., et al. (2007). UCSF tomography: an integrated software suite for real-time electron microscopic tomographic data collection, alignment, and reconstruction. J. Struct. Biol. 157, 138–147. doi: 10.1016/j.jsb.2006.06.005
Zheng, S. Q., Palovcak, E., Armache, J.-P., Verba, K. A., Cheng, Y., Agard, D. A. (2017). Motioncor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy. Nat. Methods 14, 331–332. doi: 10.1038/nmeth.4193
Zheng, S., Wolff, G., Greenan, G., Chen, Z., Faas, F. G., Bárcena, M., et al. (2022). Aretomo: An integrated software package for automated marker-free, motion-corrected cryo-electron tomographic alignment and reconstruction. J. Struct. Biology: X 6, 100068. doi: 10.1016/j.yjsbx.2022.100068
Zhou, X., Li, Y., Zhang, C., Zheng, W., Guijun, Z., Zhang, Y. (2020). Progressive and accurate assembly of multi-domain protein structures from cryo-em density maps. Nat. Comput. Sci. 2, 265–275. doi: 10.1101/2020.10.15.340455
Zhou, L., Yang, C., Gao, W., Perciano, T., Davies, K. M., Sauter, N. K. (2023). A machine learning pipeline for membrane segmentation of cryo-electron tomograms. J. Comput. Sci. 66, 101904. doi: 10.1016/j.jocs.2022.101904
Zhou, B., Yu, H., Zeng, X., Yang, X., Zhang, J., Xu, M. (2021). One-shot learning with attentionguided segmentation in cryo-electron tomography. Front. Mol. Biosci. 7. doi: 10.3389/fmolb.2020.613347
Zivanov, J., Nakane, T., Scheres, S. H. W. (2019). A bayesian approach to beam-induced motion correction in cryo-em single-particle analysis. IUCrJ 6, 5–17. doi: 10.1107/S205225251801463X
Zivanov, J., Otón, J., Ke, Z., von Kügelgen, A., Pyle, E., Qu, K., et al. (2022). A bayesian approach to single-particle electron cryo-tomography in relion-4.0. eLife 11, e83724. doi: 10.7554/eLife.83724
Keywords: cryo-electron tomography (cryo-ET), subtomogram averaging (STA), microorganism in situ, 3D reconstruction, deep learning, mathematical models
Citation: Zhao C, Lu D, Zhao Q, Ren C, Zhang H, Zhai J, Gou J, Zhu S, Zhang Y and Gong X (2023) Computational methods for in situ structural studies with cryogenic electron tomography. Front. Cell. Infect. Microbiol. 13:1135013. doi: 10.3389/fcimb.2023.1135013
Received: 31 December 2022; Accepted: 29 August 2023;
Published: 04 October 2023.
Edited by:
Feng Chen, Peking University, ChinaReviewed by:
Anastasia D. Gazi, Institut Pasteur, FranceGeng Meng, China Agricultural University, China
Copyright © 2023 Zhao, Lu, Zhao, Ren, Zhang, Zhai, Gou, Zhu, Zhang and Gong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xinqi Gong, xinqigong@ruc.edu.cn
†These authors share first authorship