Preprocessing - RELION Documentation
Preprocessing - RELION Documentation
Getting organised
We recommend to create a single directory per project, i.e. per structure you want to
determine. We’ll call this the project directory. It is important to always launch
the RELION graphical user-interface (GUI) from the project directory. In-
side the project directory you should make a separate directory to store all your raw mi-
crographs or micrograph movies in MRC, TIFF or EER format. We like to call this di-
rectory Movies/ if all movies are in one directory, or for example
Movies/15jan16/ and Movies/23jan16/ if they are in different directories (e.g.
because they were collected on different dates). If for some reason you do not want to
place your movies inside the RELION project directory, then inside the project directory
you can also make a symbolic link to the directory where your movies are stored.
Warning:
Symbolic links must be made by an absolute path (e.g.
/storage/data/15jan16 ). Use of relative paths (e.g.
../../storage/data/15jan16 ) cancause problems in later steps. Our
precalculated example contains a symbolic
v: latest link from Movies to ../Tutori-
al4.0/Movies/ but please do not follow this practice. Because we don’t know
where you decompress the archive, we cannot include a link by an absolute
path.
Single-image micrographs should have a .mrc extension, while movies can have a
.mrc , .mrcs , .tif , .tiff or .eer extension. For EER movies, see this for the
details. RELION 5.0 can also read MRC movies compressed by bzip2, xz, zstd or gzip.
When you unpacked the tutorial test data, the ( Movies/ ) directory was created. It
should contain 24 movies in compressed TIFF format, a gain-reference file
( gain.mrc ) and a NOTES file with information about the experiment.
We will start by launching the RELION GUI. As said before, this GUI always needs to be
launched from the project directory. To prevent errors with this, the GUI will ask for
confirmation the first time you launch it in a new directory. Make sure you are inside
the project directory, and launch the GUI by typing:
relion &
and answer Yes when prompted to set up a new RELION project here.
The first thing to do is to import the set of recorded micrograph movies into the pipe-
line. Select Import from the job-type browser on the left, and fill in the following pa-
rameters on the Movies/mics tab:
The MTF file can be obtained from the Gatan Web site. If you are working offline and
cannot obtain the file, you can ignore it. The MTF correction does not change the final
resolution (but changes the B factor). You can also apply it in PostProcessing.
You may provide a meaningful alias (for example: movies) for this job in the white field
named Current job: Give_alias_here . Clicking the Run! button will launch
the job. A directory called Import/job001/ will be created, together with a symbolic
link to this directory that is called Import/movies . Inside the newly created directory
a STAR file with all the movies is created. Have a look at it using:
less Import/job001/movies.star
If you had extracted your particles in a different software package, then instead of going
through the Preprocessing steps below, you would use the same Import job-type to
import particles STAR file, 3D references, 3D masks, etc. Note that this is NOT the rec-
ommended way to run RELION, and that the user is responsible for generating correct
STAR files.
Bayesian polishing (for per-particle motion-correction; see this section) can only
read local motion tracks from our own implementation;
The MOTIONCOR2 program performs outlier-pixel detection on-the-fly, and this infor-
mation is not conveyed to Bayesian polishing , which may result in unexpectedly
bad particles after polishing;
Our own implementation can write out the sum of power spectra over several movie
frames, which can be passed directly into CTFFIND 4.1 for faster CTF-estimation.
For these three reasons, we now favour running our own implementation.
files.)
First frame for 1
corrected sum::
Last frame for -1
corrected sum::
(This will result in using all movie frames.)
Dose per frame 1.277
(e/A2):
Pre-exposure 0
(e/A2):
EER fractionation: 32
Bfactor:: 150
Number of MPI 1
procs::
(Assuming you’re running this tutorial on a local computer)
Number of threads:: 12
(As these movies are 24 frames, each thread will do two
movie frames)
Submit to queue?: No
(Again, assuming you’re running this tutorial on a local
computer)
Executing this program takes approximately 5 minutes when using 12 threads on a rea-
sonably modern machine. Note that our own implementation of the MOTIONCOR2 algo-
rithm does not use a GPU. This program is multi-threaded. As each thread will work
independently on a movie frame, it is optimal to use a number of threads such that the
number of movie frames divided by the number threads is an integer number. As these
movies have 24 frames, using 12 threads will result in 2 frames being processed by each
thread. You can look at the estimated beam-induced shifts, and their statistics over the
entire data set, by selecting the out: logfile.pdf from the Display: button below
the run buttons, or you can look at the summed micrographs by selecting out: correct-
ed_micrographs.star. Depending on the size of your screen, you should probably
downscale the micrographs ( Scale: 0.3 ) and use Sigma contrast: 3 and few
columns (something like Number of columns: 3 ) for convenient visualisation.
Note that you cannot select any micrographs from this display. If you want to exclude
micrographs at this point (which we will not do, because they are all fine), you could use
the Subset selection job-type.
CTF estimation
Next, we will estimate the CTF parameters for each corrected micrograph. You can use
the CTF estimation job-type as a wrapper to Alexis Rohou and Niko Grigorieff’s
CTFFIND 4.1 to execute efficiently on the CPU. We now prefer CTFFIND 4.1, as it is the only
open-source option, and because it allows reading in the movie-averaged power spectra
calculation by RELION’s own implementation of the MOTIONCOR2 algorithm. Support for
GCTF was dropped in RELION 5.0. Fill in the settings as follows:
On the I/O :
Input Motioncorr/job002/corrected_micrographs.star
micrograp
hs STAR
file:: (You can again use the Browse button to select the corrected_micro-
graphs.star file of the Motion correction job.)
Use No
micrograp
h without (These may have better Thon rings than the dose-weighted ones, but we
dose-
decided in the previous step not to write these out)
weighting?
:
Estimate No
phase
shifts?: (This is only useful for phase-plate data)
Amount of 100
astigmatis
m (A):: (Assuming your scope was reasonably well aligned, this value will be
suitable for many data sets.)
On the Running tab, use six MPI processes to process the 24 micrographs in parallel.
This took less than 10 seconds on our machine. Once the job finishes there are addition-
al files for each micrograph inside the output CtfFind/job003/Movies directory:
the .ctf file contains an image in MRC format with the computed power spectrum
and the fitted CTF model; the .log file contains the output from CTFFIND; the .com file
contains the script that was used to launch CTFFIND.
You can visualise all the Thon-ring images using the Display button, selecting out:
micrographs_ctf.star . The zeros between the Thon rings in the experimental im-
ages should coincide with the ones in the model. Note that you can sort the display in
order of defocus, maximum resolution, figure-of-merit, etc. The logfile.pdf file
contains plots of useful parameters, such as defocus, astigmatism, estimated resolution,
etc for all micrographs, and histograms of these values over the entire data set.
Analysing these plots may be useful to spot problems in your data acquisition.
If you see CTF models that are not a satisfactory fit to the experimental Thon rings, you
can delete the .log files for those micrographs, select the CtfFind/job003 entry
from the Finished jobs list, alter the parameters in the parameter-panel, and then re-
run the job by clicking the Continue! button. Only those micrographs for which a
.log file does not exist will be re-processed. You can do this until all CTF models are
satisfactory. If this is not possible, or if you decide to discard micrographs because they
have unsatisfactory Thon rins, you can use the Subset selection job-type to do this.
Picking particles manually is a personal experience! If you don’t like to pick particles in
RELION, we also support coordinate file formats for Jude Short’s ximdisp [Smi99] (with
any extension); for xmipp-2.4 [SNRS+08] (with any extension); and for Steven Ludtke’s
e2boxer.py [TPB+07] (with a .box extension). If you use any of these, make sure to
save the coordinate files as a text file in the same directory as from where you imported
the micrographs (or movies), and with the same micrograph rootname, but a different
(suffix+) extension as the micrograph, e.g. Movies/006.box or
Movies/006_pick.star for micrograph Movies/006.mrc . You should then use
the Import job-type and set Node type: to 2D/3D particle coordinates .
Make sure that the Input Files: field contains a linux wildcard, followed by the co-
ordinate-file suffix, e.g. for the examples above you have to give Movies/*.box or
Movies/*_pick.star , respectively.
On the I/O tab of the Manual picking job-type, use the Browse button to select the
micrographs_ctf.star file that was created in CtfFind/job003 , ignore the
Colors tab, and fill in the Display tab as follows:
(Use this to manually set which value will be white. For this
to work, Sigma contrast should be set to 0)
Black value:: 0
(Use this to manually set which value will be black. For this
to work, Sigma contrast should be set to 0)
Lowpass filter (A):: -1
(Playing with this may help you to see particles better in
very noisy micrographs)
Highpass filter (A):: -1
Note:
As of RELION-5.0, Topaz comes pre-installed in the relion-5 conda environment,
which should be picked up automatically by the GUI.
Run the job by clicking the Run! button and click on a few particles if you want to.
However, as we will use the LoG-based autopicking in the next section, you do not
need to pick any if you don’t want to. If you were going to use manually picked
particles for an initial 2D classification job, then you would need approximately 500-
1,000 particles in order to calculate reasonable class averages. Left-mouse click for
picking, middle-mouse click for deleting a picked particle, right-mouse click for a pop-
up menu in which you will need to save the coordinates!. Note that you can al-
ways come back to pick more from where you left it (provided you saved the STAR files
with the coordinates throught the pop-up menu), by selecting ManualPick/job004
from the Finished jobs and clicking the Continue! button.
Ad by EthicalAds ·