This repository contains the source code to reproduce the figures of the manuscript entitled The RNA binding proteome of axonal mRNAs in sympathetic neurons by R Luisier, C Andreassi, L Fournier and A Riccio. The repository contains:
Bioconductor version 3.16 (BiocManager 1.30.19), R 4.2.2 (2022-10-31)
The following R packages should be installed:
GenomicRanges_1.50.2
Rsamtools_2.14.0
rtracklayer_1.58.0
IRanges_2.32.0
geneplotter_1.76.0
multtest_2.54.0
mclust_6.0.0
knitr_1.42
edgeR_3.40.2
topGO_2.50.0
SparseM_1.81
graph_1.76.0
plotly_4.10.1
fitdistrplus_1.1-8
GO.db_3.16.0
The python requirements are listed in python_requirements.txt.
To run the jupyter notebook, create a python environment with conda by running the command conda create --name <env_name> --file python_requirements.txt (tested on osx-64). You can then activate the environment by running conda activate <env_name>.
- data: folder containing the data for examples matrix of gene expression; etc. Raw sequencing data will be deposited publicly.
- Scripts:
R,PythonandBashcustome code
The attached code can be run with data that are publicly available on Zenodo (https://zenodo.org/record/8047412).
3’ end sequencing of RNA isolated from axons and cell bodies of sympathetic neurons exposed to either Nerve Growth factor (NGF) or Neurotrophin 3 (NT3). To identify mRNAs localised in sympathetic cell bodies and axons, we combined compartmentalized cultures and 3’UTR sequencing (Andreassi et al. 2021). Compartmentalized chambers allow the physical separation of cell bodies from distal axons and are especially suited for sympathetic neurons because these grow in culture as a highly homogeneous population without glial cells. Neurons were seeded in the central compartment with NGF (100 ng/ml) and after 5 days, the NGF was lowered to 10ng/ml in the cell bodies, and NGF or NT3 were supplied at high concentration to the peripheral compartments to stimulate axon growth. After 7 additional days in vitro, RNA was isolated from either cell bodies or axons, subject to 2 rounds of linear amplification and sequenced by using 3’end-Seq . This technique allowed the sequencing of transcripts 3’ends independently of the length of the transcript (Andreassi et al. 2021). There are two technical replicates for each neurotrophin condition.
The analysis of the read count and samples is presented in the preprocessing analysis that relates to Supplementary Figure 1.
The analysis of the differential gene expression analysis between NT3 and NGF conditions together with transcription factor analysis is presented in DGE analysis. This relates to first half of Figure 1.
The analysis of differential APA between NGF and NT3 and related results are presented in APA analysis, that relates to the second half of Figure 1.
The RBPome study that underlies APA in developing sympathetic neurons was performed by integrating the 3'end sequencing data obtained from sympathetic neurons with publicly available CLIP-sequencing data obtained from human cell lines and lifted over the rat genome. The analysis and results are presented in APA RBPome and relates to Figure 2.
The analysis of the compartment-specific mRNA pools is presented in compartment analysis, that relates to the first part of the Figure 3.
The modeling of the axonal localisation alongside analysis of differential localisation between NGF and NT3 condition are presented in the LS analysis analysis, that relates to the second part of the Figure 3 and Figure 4.
The analysis of the synergistic regulatory potential of RBPs in axonal localisation in presented synergistic analysis, that relates to the second part of Figure 5.
The analysis of the axonal remodeling and related RBPome in presented axonal remodelling analysis analysis, that relates to the second part of Figure 6.
This work has been published in Genome Research (doi: 10.1101/gr.277804.123).