Name: analyzing multifactorial rna-seq experiments with dicoexpress
Uploaded: 2022-07-29T12:30:00
Description: Watch this Scientific Journal Video about Analyzing Multifactorial RNA-Seq Experiments with DiCoExpress at JoVE.com

This work was mainly supported by the ANR PSYCHE (ANR-16-CE20-0009). The authors thank F. Desprez for the construction of the container of DiCoExpress. KB work is supported by the Investment for the Future ANR-10-BTBR-01-01 Amaizing program. The GQE and IPS2 laboratories benefit from the support of Saclay Plant Sciences-SPS (ANR-17-EUR-0007).

1. DiCoExpress
<ol>
	<li>Open a R studio session and set directory to Template_scripts.</li>
	<li>Open the DiCoExpress_Tutorial.R script in R studio.</li>
	<li>Load DiCoExpress functions in the R session with the following commands: 
	&#62; source(&#34;../Sources/Load_Functions.R&#34;) 
	&#62; Load_Functions() 
	&#62; Data_Directory = &#34;../Data&#34; 
	&#62; Results_Directory = &#34;../Results/&#34;</li>
	<li>Load data files in the R session with the following commands: 
	&#62; Project_Name =...

<ol>
	<li>Wang, Z., Gerstein, M., Snyder, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RNA-Seq:+a+revolutionary+tool+for+transcriptomics.">RNA-Seq: a revolutionary tool for transcriptomics.</a> Nature reviews. Genetics. 10 (1), 57-63 (2009).</li><li>Yang, I. S., Kim, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Analysis+of+Whole+Transcriptome+Sequencing+Data:+Workflow+and+Software.">Analysis of Whole Transcriptome Sequencing Data: Workflow and Software.</a> Genomics &amp; Informatics. 13 (4), 119-125 (2015).</li><li>R Core Team. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=R:+A+language+and+environment+for+statistical+computing.">R: A language and environment for statistical computing.</a> R Foundation for Statistical Computing. , (2020).</li><li>Huber, W., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Orchestrating+high-throughput+genomic+analysis+with+Bioconductor.">Orchestrating high-throughput genomic analysis with Bioconductor.</a> Nature Methods. 12 (2), 115-121 (2015).</li><li>Smith, D. R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+battle+for+user-friendly+bioinformatics.">The battle for user-friendly bioinformatics.</a> Frontiers in Genetics. 4, 187 (2013).</li><li>Pavelin, K., Cham, J. A., de Matos, P., Brooksbank, C., Cameron, G., Steinbeck, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Bioinformatics+Meets+User-Centred+Design:+A+Perspective.">Bioinformatics Meets User-Centred Design: A Perspective.</a> PLoS Computational Biology. 8 (7), 1002554 (2012).</li><li>. Shiny: web application framework Available from: <a target="_blank" href="https://rdrr.io/cran/shiny/">https://rdrr.io/cran/shiny/</a> (2021)</li><li>Lambert, I., Roux, C. P. -. L., Colella, S., Martin-Magniette, M. -. L. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=DiCoExpress:+a+tool+to+process+multifactorial+RNAseq+experiments+from+quality+controls+to+co-expression+analysis+through+differential+analysis+based+on+contrasts+inside+GLM+models.">DiCoExpress: a tool to process multifactorial RNAseq experiments from quality controls to co-expression analysis through differential analysis based on contrasts inside GLM models.</a> Plant methods. 16 (1), 68 (2020).</li><li>Dillies, M. -. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+comprehensive+evaluation+of+normalization+methods+for+Illumina+high-throughput+RNA+sequencing+data+analysis.">A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis.</a> Briefings in bioinformatics. 14 (6), 671-683 (2012).</li><li>Rigaill, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Synthetic+data+sets+for+the+identification+of+key+ingredients+for+RNA-seq+differential+analysis.">Synthetic data sets for the identification of key ingredients for RNA-seq differential analysis.</a> Briefings in Bioinformatics. 19 (1), (2016).</li><li>Rau, A., Maugis-Rabusseau, C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Transformation+and+model+choice+for+RNA-seq+co-expression+analysis.">Transformation and model choice for RNA-seq co-expression analysis.</a> Briefings in Bioinformatics. 19 (3), (2017).</li><li>Robinson, M. D., McCarthy, D. J., Smyth, G. K. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=edgeR:+a+Bioconductor+package+for+differential+expression+analysis+of+digital+gene+expression+data.">edgeR: a Bioconductor package for differential expression analysis of digital gene expression data.</a> Bioinformatics. 26 (1), 139-140 (2009).</li><li>Wilkinson, M. D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+FAIR+Guiding+Principles+for+scientific+data+management+and+stewardship.">The FAIR Guiding Principles for scientific data management and stewardship.</a> Scientific Data. 3 (1), 160018 (2016).</li><li>Stark, R., Grzelak, M., Hadfield, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=RNA+sequencing:+the+teenage+years.">RNA sequencing: the teenage years.</a> Nature Reviews Genetics. 20 (11), 631-656 (2019).</li></ol>

Because RNA-Seq has become a ubiquitous method in biological studies, there is a constant need to develop versatile and user-friendly analytical tools. A critical step within most of the analytical workflows is often to identify with confidence the genes differentially expressed between biological conditions and/or treatments15. The production of reliable results requires proper statistical modeling, which has been the motivation for the development of DiCoExpress.
DiCo...

Next-generation RNA sequencing (RNA-Seq) technology is now the gold standard of transcriptome analysis1. Since the early days of the technology, the combined efforts of bioinformaticians and biostatisticians have resulted in the development of numerous methods tackling all the essential steps of transcriptomic analyses, from mapping to transcript quantification2. Most of the tools available today to the biologist are developed within the R software environment for statistical computing and graphs3, and many packages for biological data analysis are available in the Bioconductor repository4. These packages offer total control and customization of the analysis, but they come at the cost of extensive use of a command-line interface. Because many biologists are more comfortable with a &#34;point and click&#34; approach5, the democratization of RNA-Seq analyses requires the development of more user-friendly interfaces or protocols6. For example, it is possible to build web interfaces of R packages using Shiny7, and command-line data analysis is made more intuitive with the R-studio8 interface. The development of dedicated, step-by-step tutorials can also help the novel user. In particular, a video tutorial supplements a classic text one, leading to a deeper understanding of all the procedure steps.
We recently developed DiCoExpress9, a tool for analyzing multifactorial RNA-Seq experiments in R using methods considered to be the best ones based on neutral comparison studies10,11,12. Starting from a count table, DiCoExpress proposes a data quality control step followed by a differential gene expression analysis (edgeR package13) using a generalized linear model (GLM) and the generation of co-expression clusters using Gaussian mixture models (coseq package12). DiCoExpress handles complete and unbalanced design up to 2 biological factors (i.e., genotype and treatment) and one technical factor (i.e., replicate). The originality of DiCoExpress lies in its directory architecture storing and organizing data, scripts, and results and in the automation of the writing of the contrasts allowing the user to investigate numerous questions within the same statistical model. An effort was also made to provide graphical outputs illustrating the statistical results.
The DiCoExpress workspace is available at https://forgemia.inra.fr/GNet/dicoexpress.&#160;It contains four directories, two pdf, and two text files. The Data/ directory contains the input datasets; for this protocol, we will use the &#34;tutorial&#34; dataset. The Sources/ directory contains seven R functions necessary to perform the analysis, and must not be modified by the user. The analysis is run using scripts stored in the Template_scripts/ directory. The one used in this protocol is called DiCoExpress_Tutorial_JoVE.R and can be easily adapted to any transcriptomic project. All the results are written in the Results/ directory and stored in a subdirectory named according to the project. The README.md file contains useful installation information, and any specific details concerning the method and its use can be found in the DiCoExpress_Reference_Manual.pdf file.
This video tutorial guides the user through the different features of DiCoExpress with the aim to overcome the reluctance felt by biologists using command-line-based tools. We present here the analysis of an artificial RNA-Seq dataset describing gene expression in three biological replicates of four genotypes, with or without treatment. We will now go through the different steps of the DiCoExpress workflow illustrated in Figure 1. The script described in the Protocol section and input files are available on the site: https://forgemia.inra.fr/GNet/dicoexpress
Prepare data files 
The four csv files stored in the Data/ directory should be named according to the project name. In our example, all the names, therefore, begin with &#34;Tutorial&#34;, and we will set Project_Name = &#34;Tutorial&#34; in Step 4 of the protocol. The separator used in the csv files must be indicated in the Sep variable in Step 4. In our &#34;tutorial&#34; dataset, the separator is a tabulation. For advanced users the full dataset can be reduced to a subset by providing a list of instructions and a new Project_Name through the Filter variable. This option avoids redundant copies of the input files and verifies FAIR principles14.
Among the four csv files, only the COUNTS and TARGET files are mandatory. They contain the raw counts for every gene (here Tutorial_COUNTS.csv) and the experimental design description (here Tutorial_TARGET.csv). The TARGET.csv file describes every sample (one sample per row) with a modality for each biological or technical factor (in the columns). We strongly recommend that the names chosen for the modalities start with a letter, not a number. The name of the last column (&#34;Replicate&#34;) cannot be changed. Finally, the sample names (first column) must match the names in the headings of the COUNTS.csv file (Genotype1_control_rep1 in our example). The Enrichment.csv file in which every line contains one Gene_ID and one annotation term is only required if the user plans to run the enrichment analysis. If one gene has several annotations, they will have to be written on different lines. The Annotation.csv file is optional and is used to add a short description of every gene in the output files. The best way to get an annotation file is to retrieve the information from dedicated databases (e.g., Thalemine: https://bar.utoronto.ca/thalemine/begin.do for Arabidopsis).
Installation of DiCoExpress 
DiCoExpress requires specific R packages. Use the command line source(&#34;../Sources/Install_Packages.R&#34;) in the R console to check the required package installation status. For users on Linux, another solution is to install the container dedicated to DiCoExpress and available at https://forgemia.inra.fr/GNet/dicoexpress/container_registry. By definition, this container contains DiCoExpress with all of the parts needed, such as libraries and other dependencies.

analyzing multifactorial rna-seq experiments with dicoexpress

The proper use of statistical modeling in NGS data analysis requires an advanced level of expertise. There has recently been a growing consensus on using generalized linear models for differential analysis of RNA-Seq data and the advantage of mixture models to perform co-expression analysis. To offer a managed setting to use these modeling approaches, we developed DiCoExpress that provides a standardized R pipeline to perform an RNA-Seq analysis. Without any particular knowledge in statistics or R programming, beginners can perform a complete RNA-Seq analysis from quality controls to co-expression through differential analysis based on contrasts inside a generalized linear model. An enrichment analysis is proposed both on the lists of differentially expressed genes, and the co-expressed gene clusters. This video tutorial is conceived as a step-by-step protocol to help users take full advantage of DiCoExpress and its potential in empowering the biological interpretation of an RNA-Seq experiment.

All the DiCoExpress outputs are saved in the Tutorial/ directory, itself placed within the Results/ directory. We provide here some guidance for assessing the overall quality of the analysis.
Quality Control 
The quality control output, located in the Quality_Control/ directory, is essential to verify that the RNA-Seq analysis results are reliable. The Data_Quality_Control.pdf file contains several plots obtained with raw and normalized data that can be used to identify a...

Watch this Scientific Journal Video about Analyzing Multifactorial RNA-Seq Experiments with DiCoExpress at JoVE.com

Analyzing Multifactorial RNA-Seq Experiments with DiCoExpress

DiCoExpress is a script-based tool implemented in R to perform an RNA-Seq analysis from quality control to co-expression. DiCoExpress handles complete and unbalanced design up to 2 biological factors. This video tutorial guides the user through the different features of DiCoExpress.

The proper use of statistical modeling in NGS data analysis requires an advanced level of expertise. There has recently been a growing consensus on using ...

Research

JoVE Journal

Engineering

Characterization of Surface Modifications by White Light Interferometry: Applications in Ion Sputtering, Laser Ablation, and Tribology Experiments

In materials science and engineering it is often necessary to obtain quantitative measurements of surface topography with micrometer lateral resolution. ...

Flame Experiments at the Advanced Light Source: New Insights into Soot Formation Processes

The following experimental protocols and the accompanying video are concerned with the flame experiments that are performed at the Chemical Dynamics ...

Analyzing the Movement of the Nauplius 'Artemia salina' by Optical Tracking of Plasmonic Nanoparticles

Analyzing the Movement of the Nauplius 'Artemia salina' by Optical Tracking of Plasmonic Nanoparticles

We demonstrate how optical tweezers may provide a sensitive tool to analyze the fluidic vibrations generated by the movement of small aquatic organisms. A ...

Experiments on Ultrasonic Lubrication Using a Piezoelectrically-assisted Tribometer and Optical Profilometer

Friction and wear are detrimental to engineered systems. Ultrasonic lubrication is achieved when the interface between two sliding surfaces is vibrated at ...

Visualizing Hyporheic Flow Through Bedforms Using Dye Experiments and Simulation

Advective exchange between the pore space of sediments and the overlying water column, called hyporheic exchange in fluvial environments, drives solute ...

Wind Tunnel Experiments to Study Chaparral Crown Fires

The present protocol presents a laboratory technique designed to study chaparral crown fire ignition and spread. Experiments were conducted in a low ...

Conducting Elevated Temperature Normal and Combined Pressure-Shear Plate Impact Experiments Via a Breech-end Sabot Heater System

A novel approach for conducting normal and/or combined pressure-shear plate impact experiments at test temperatures up to 1000 &#176;C is presented. The ...

Enhanced Electron Injection and Exciton Confinement for Pure Blue Quantum-Dot Light-Emitting Diodes by Introducing Partially Oxidized Aluminum Cathode

Stable and efficient red (R), green (G), and blue (B) light sources based on solution-processed quantum dots (QDs) play important roles in next-generation ...