<figure><img src="image_1.png" /></figure>

HAL
open science

# FAIR_Bioinfo: a turnkey training course and protocol for reproducible computational biology

Thomas Denecker, Claire Toffano-Nioche

► To cite this version:

Thomas Denecker, Claire Toffano-Nioche. FAIR_Bioinfo: a turnkey training course and protocol for reproducible computational biology. JOSE - Journal of Open Source Education, 2021, 10.21105/jose.00068. hal-02880655v2

HAL Id: hal-02880655

https://hal.science/hal-02880655v2

Submitted on 10 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# FAIR_Bioinfo: a turnkey training course and protocol for reproducible computational biology

Thomas Denecker¹ and Claire Toffano-Nioche¹

1 CEA, CNRS, Univ. Paris-Sud, Institute for Integrative Biology of the Cell (I2BC), Gif-sur-Yvette, France.

DOI: 10.21105/jose.00068

## Software

- Review ⚙️

- Repository ⚙️

- Archive ⚙️

Submitted: 18 July 2019

Published: 26 August 2021

## License

Authors of papers retain copyright and release the work under a Creative Commons Attribution 4.0 International License (CC BY 4.0).

## Summary

Reproducibility plays an essential part in the success of a bioinformatics project. Indeed, Reproducibility makes it possible to guarantee the validity of scientific results and to simplify the dissemination of projects. To help disseminate Reproducibility principles among bioinformatics students, engineers and scientists, we created the FAIR_Bioinfo course, which presents a set of features we consider necessary to make a complete bioinformatics analysis reproducible. To illustrate the theoretical concepts of reproducibility, we use as an example a classic bioinformatics analysis (differential gene expression analysis from RNA-seq data). In short, we retrieve the data from public databases (ENA/SRA), we perform a reproducible analysis using a workflow management system (snakemake) in a virtual environment (Docker). The entire versioned (git) code is open source (Github https://github.com/thomasdenecker/FAIR_Bioinfo and dockerhub https://hub.docker.com/r/tdenecker/fair_bioinfo). The visualization of the results is dynamic (Shiny app) and the PDF or HTML report (R markdown) provides the results of the analysis and lists all user-selected parameters.

## How to use

The whole example proposed in FAIR_Bioinfo can be simply reproduced with the following 4 command lines executed in a terminal (note that docker must be installed beforehand - https://www.docker.com/products/docker-desktop):

```
$ git clone https://github.com/thomasdenecker/FAIR_Bioinfo
```

```
$ cd FAIR_Bioinfo
```

```
$ sudo docker run --rm -d -p 8888:8888 --name fair_bioinfo -v ${PWD}:/home/rstudio
```

```
$ sudo docker exec -it fair_bioinfo bash ./FAIR_script.sh
```

Once the analysis is complete, the web application can be launched with the following command:

```
$ sudo docker exec -it fair_bioinfo bash ./FAIR_app.sh
```

and used by copying http://localhost:8888/rstudio/p/4444/ into your browser.

Note: The use of the sudo command is not necessary for Mac OS and Windows. The demonstrations of the course have been made on a Linux environment (Ubuntu 18.04).
<figure><img src="image_3.png" /><figcaption>Figure 1: 7-steps</figcaption></figure>

To illustrate FAIR principles for bioinformatics workflows, we used a classic RNA-seq data analysis workflow. Starting from raw data (FASTQ files), this analysis aims to identify genes that are differentially expressed between different conditions. The underlying biology is not detailed, as our focus is on the reproducibility aspects rather than the specific bioinformatics protocol or biological question. Eventually, the RNA-seq analysis protocol used as an example here could be replaced by any other type of bioinformatics analysis.

Instructor notes

The training course is divided into 8 sessions, each one bringing an additional level of reproducibility to the global workflow. Each session lasted an hour and a half; they were carried out at the rate of one session per month. Table 1 shows the detailed program and learning objectives of each session.

<table><caption>Table 1: Overview of sessions</caption><thead><tr><td>Title</td><td>Session</td><td>Description</td><td>Tools</td></tr></thead><tbody><tr><td>This is not magic</td><td>1</td><td>Open a terminal and retrieve data for the analysis pipeline (loop & variable concepts, bash command)</td><td>shell, wget, md5sum</td></tr><tr><td>The code memory</td><td>2</td><td>Initiation to code versioning</td><td>Git & Github</td></tr><tr><td>Play with analysis tools</td><td>3</td><td>Implementation of a first part of the analysis pipeline from the fastq file to the count table</td><td>conda, third-party tools used for bioinformatics analysis (FastQC, bowtie2, samtools, HTseq-counts)</td></tr></tbody></table>
<table><thead><tr><td>Title</td><td>Session</td><td>Description</td><td>Tools</td></tr></thead><tbody><tr><td>A trip to the sea</td><td>4</td><td>Control the computing environment using a container and compute a cloud platform</td><td>Docker, ssh</td></tr><tr><td>I've got the power!</td><td>5</td><td>Parallel computing and use of a computer cluster</td><td>Snakemake, slurm, Singularity, ascp</td></tr><tr><td>LoveR</td><td>6</td><td>Upgrade to dynamic rendering using a web application and presentation of a few powerful R packages</td><td>Shiny, R packages</td></tr><tr><td>Sharing results and protocols using notebooks</td><td>7</td><td>Creating an analysis report and sharing it</td><td>Jupyter &amp; R-markdown</td></tr><tr><td>Disseminate your project</td><td>8</td><td>Disseminate a reproducible project. Perspective for improvement</td><td>Githubpages, License, Release, Zenodo</td></tr></tbody></table>

The slides of each session are available on the FAIR_Bioinfo github in French. They have been translated into English and inserted as figures in the Gitbook. English speaking instructors can very easily create slides from the Gitbook and these figures.

## Conclusion and Perspectives

We proposed a training course to help scientists improve their computational practice by ensuring reproducibility of bioinformatic analysis. The training course is divided into 8 sessions each introducing a core competence for reproducible analysis. We also address possible future developments (last session), including: - Installing all tools in a Conda environment itself residing in a Docker image. The goal is to isolate the installation layer (Grüning et al., 2018) - Setting up a Virtual machine to run the Docker container. The goal is to stand the test of time by also setting the Operating System. - Implement continuous integration to ensure that code changes do not cause changes in results.

Through the FAIR_Bioinfo training, we offer users a solution to make Bioinformatics analyses entirely reproducible. All the tools used can be replaced by other equivalent solutions (docker by singularity, snakemake by nextflow,...). Thus, with FAIR data and protocols according to the principles of FAIR_Bioinfo, we obtain FAIR-processed data. This same data can be used again and again to enter a virtuous circle of reproducible analyses.

## References

* Baker, M. (2016). 1,500 scientists lift the lid on reproducibility: Nature News & Comment. https://www.nature.com/news/1-500-scientists-lift-the-lid-on-reproducibility-1.19970

* Denecker, T., & Toffano-Nioche, C. (2019). 7-steps.png. https://doi.org/10.6084/m9.fishshare.8950505.v1
Grüning, B., Chilton, J., Köster, J., Dale, R., Soranzo, N., Beek, M. van den, Goecks, J., Backofen, R., Nekrutenko, A., & Taylor, J. (2018). Practical Computational Reproducibility in the Life Sciences. *Cell Systems, 6*(6), 631–635. https://doi.org/10.1016/j.cels.2018.03.014

Warren, A. (2015). *Repeatability and Benefaction in Computer Systems Research — A Study and a Modest Proposal*.

Wilkinson, M. D., Dumontier, M., Aalbersberg, Ij. J., Appleton, G., Axton, M., Baak, A., Blomberg, N., Boiten, J.-W., Silva Santos, L. B. da, Bourne, P. E., Bouwman, J., Brookes, A. J., Clark, T., Crosas, M., Dillo, I., Dumon, O., Edmunds, S., Evelo, C. T., Finkers, R., ... Mons, B. (2016). The FAIR Guiding Principles for scientific data management and stewardship. *Scientific Data, 3*, 160018. https://doi.org/10.1038/sdata.2016.18