Advancing brain barriers RNA sequencing: guidelines from experimental design to publication

Abstract

Background: RNA sequencing (RNA-Seq) in its varied forms has become an indispensable tool for analyzing differential gene expression and thus characterization of specific tissues. Aiming to understand the brain barriers genetic signature, RNA seq has also been introduced in brain barriers research. This has led to availability of both, bulk and single-cell RNA-Seq datasets over the last few years. If appropriately performed, the RNA-Seq studies provide powerful datasets that allow for significant deepening of knowledge on the molecular mechanisms that establish the brain barriers. However, RNA-Seq studies comprise complex workflows that require to consider many options and variables before, during and after the proper sequencing process. Main body: In the current manuscript, we build on the interdisciplinary experience of the European PhD Training Network BtRAIN (https ://www.btrai n-2020.eu/) where bioinformaticians and brain barriers researchers collaborated to analyze and establish RNA-Seq datasets on vertebrate brain barriers. The obstacles BtRAIN has identified in this process have been integrated into the present manuscript. It provides guidelines along the entire workflow of brain barriers RNA-Seq studies starting from the overall experimental design to interpretation of results. Focusing on the vertebrate endothelial blood-brain barrier (BBB) and epithelial blood-cerebrospinal-fluid barrier (BCSFB) of the choroid plexus, we provide a step-by-step description of the workflow, highlighting the decisions to be made at each step of the workflow and explaining the strengths and weaknesses of individual choices made. Finally, we propose recommendations for accurate data interpretation and on the information to be included into a publication to ensure appropriate accessibility of the data and reproducibility of the observations by the scientific community. Conclusion: Next generation transcriptomic profiling of the brain barriers provides a novel resource for understanding the development, function and pathology of these barrier cells, which is essential for understanding CNS homeostasis and disease. Continuous advancement and sophistication of RNA-Seq will require interdisciplinary approaches between brain barrier researchers and bioinformaticians as successfully performed in BtRAIN. The present guidelines are built on the BtRAIN interdisciplinary experience and aim to facilitate collaboration of brain barriers researchers with bioinformaticians to advance RNA-Seq study design in the brain barriers community.

Both the Complexity of Tight Junctions and Endothelial Transcytosis Are Increased During BBB Postnatal Development in Rats

“…The BBB is located at the level of microvascular endothelial cells in the CNS, which includes capillaries, pre-capillary arterioles and post-capillary venules ( Northrop and Yamamoto, 2015 ; Francisco et al, 2020 ). Previous gene profiling studies of brain vessels had several limitations, including the use of whole vessel fractions that contain many different cell types, the lack of description of the cerebral endothelial isolation procedures, and separation of the entire brain instead of specific areas ( Enerson and Drewes, 2006 ; Richard et al, 2010 ; Munji et al, 2019 ).…”

Brain microvascular endothelial cell dysfunction in an isogenic juvenile iPSC model of Huntington’s disease

“…Two biological replicates were analyzed of: HD180 iPSCs, HD180 iBMECs, HD-corrected iPSCs, and HD-corrected iBMECs. Given recent guidelines for RNA sequencing which suggest a minimum of three replicates [ 18 ], a limitation is that our analysis may be underpowered. iPSCs were harvested prior to differentiation in UM/F- media, while iBMECs were harvested as confluent monolayers 2 days following sub-culture on collagen IV and fibronectin-coated tissue-culture plates.…”

Brief Developmental Exposure to Fluoxetine Causes Life-Long Alteration of the Brain Transcriptome in Zebrafish

“…This function also reported an adjusted p-value (padj) based on the Benjamini-Hochberg procedure, which controls the false discovery rate (FDR). The threshold to identify DEGs was an adjusted p value ≤ 0.05 and a fold change (FC) ≥1.2, which has been reported as the optimal FC for nervous system transcriptomic analysis ( 45 ). A hierarchal clustering ( Figure 2 ) and a principle component analysis (PCA) ( Supplementary Figure 3 ) was performed on the differentially expressed set of genes to further evaluate and compare the transcription patterns between treatment conditions.…”