R

This notebook uses the rOMERO-gateway and EBImage to process an Image associated to the paper 'Timing of gene expression in a cell-fate decision system'.

The Image "Pos22" is taken from the dataset idr0040-aymoz-singlecell/experimentA/YDA306_AGA1y_PRM1r_Mating. It is a timelapse Image with 42 timepoints separated by 5 minutes. This Image is used to fit a model for the growth of the yeast cells. The notebook does not replicate any of the analysis of the above mentioned paper.

Its purpose is mainly to demonstrate the use of Jupyter, rOMERO-gateway and EBimage.

What it does:

• For each time point of one movie:
  • Read the image for this time point  from the IDR
  • Threshold the images and count the cells using EBimage functions
• Fit an exponential model to the count of cells against time to get a coefficient of grow (exponential factor)

This workflow can be ran with data from 3D-SIM showing the centrosomes in order to compare the distribution of diameters of rings (or toroids) of different proteins from the centrioles or the peri centriolar material. It aims to reproduce the results of the Nature Cell Biology Paper Subdiffraction imaging of centrosomes reveals higher-order organizational features of pericentriolar material  from the same data set but with a different analysis method.

It is slightly different from the methods described in the paper itself, where the method was to work on a maximum intensity projection of a 3D-SIM stack, and then to fit circle to the centrioles to estimate the diameters of the toroids.

In this workflow, the images are read from the IDR , then process by thresholding (Maximum entropy auto thresholding with Image J), and processed by Analyze Particles  with different measurement sets, including the bouding box. Then the analysis of diameters and the statistical test are performed using R. All the code and data sets are available, and in the case of this paper have shown a layered organisation of the proteins.

Bioconductor provides tools for the analysis and comprehension of high-throughput genomic data. Bioconductor uses the R statistical programming language, and is open source and open development. It has two releases each year, 1560 software packages, and an active user community. Bioconductor is also available as an AMI (Amazon Machine Image) and a series of Docker images.

Neuroconductor is an open-source platform for rapid testing and dissemination of reproducible computational imaging software, specialized in brain medical imaging (MRI, fMRI, DTI, etc...) but that could be used on a wider range of images. The goals of the project are to:

• provide a centralized repository of R software dedicated to image analysis;
• disseminate quickly software updates;
• educate a large, diverse community of scientists using detailed tutorials and short courses;
• ensure quality via automatic and manual quality controls; and
• promote reproducibility of image data analysis.

Based on the programming language R, Neuroconductor starts with 68 inter-operable packages that cover multiple areas of imaging including visualization, data processing and storage, and statistical inference. Neuroconductor accepts new R package submissions, which are subject to a formal review and continuous automated testing.

This R package implements the NBLAST neuron similarity algorithm described in a preprint available at http://dx.doi.org/10.1101/006346. In addition to basic pairwise comparison, the package implements search of databases of neurons. There is also suport for all x all comparison for a group of neurons. This can produce a distance matrix suitable for hierarchical clustering, which is also implemented in the package.