Standardizing workflows in imaging transcriptomics with the abagen toolbox

  1. Ross D Markello  Is a corresponding author
  2. Aurina Arnatkevičiūtė
  3. Jean-Baptiste Poline
  4. Ben D Fulcher
  5. Alex Fornito
  6. Bratislav Misic  Is a corresponding author
  1. McGill University, Canada
  2. Monash University, Australia
  3. University of Sydney, Australia

Abstract

Gene expression fundamentally shapes the structural and functional architecture of the human brain. Open-access transcriptomic datasets like the Allen Human Brain Atlas provide an unprecedented ability to examine these mechanisms in vivo; however, a lack of standardization across research groups has given rise to myriad processing pipelines for using these data. Here, we develop the abagen toolbox, an open-access software package for working with transcriptomic data, and use it to examine how methodological variability influences the outcomes of research using the Allen Human Brain Atlas. Applying three prototypical analyses to the outputs of 750,000 unique processing pipelines, we find that choice of pipeline has a large impact on research findings, with parameters commonly varied in the literature influencing correlations between derived gene expression and other imaging phenotypes by as much as p ≥ 1:0. Our results further reveal an ordering of parameter importance, with processing steps that influence gene normalization yielding the greatest impact on downstream statistical inferences and conclusions. The presented work and the development of the abagen toolbox lay the foundation for more standardized and systematic research in imaging transcriptomics, and will help to advance future understanding of the influence of gene expression in the human brain.

Data availability

All datasets used in this study are publicly available. Detailed information about the datasets and how to access them are described in the manuscript.

The following previously published data sets were used

Article and author information

Author details

  1. Ross D Markello

    McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University, Montreal, Canada
    For correspondence
    ross.markello@mail.mcgill.ca
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-1057-1336
  2. Aurina Arnatkevičiūtė

    Monash University, Clayton, Australia
    Competing interests
    No competing interests declared.
  3. Jean-Baptiste Poline

    Montreal Neurological Institute & Hospital,, McGill University, Montreal, QC, Canada, Canada
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9794-749X
  4. Ben D Fulcher

    School of Physics, University of Sydney, Sydney, Australia
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3003-4055
  5. Alex Fornito

    Monash University, Clayton, Australia
    Competing interests
    Alex Fornito, Reviewing editor, eLife.
  6. Bratislav Misic

    McConnell Brain Imaging Centre, Montreal Neurological Institute, McGill University, Montreal, Canada
    For correspondence
    bratislav.misic@mcgill.ca
    Competing interests
    No competing interests declared.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0307-2862

Funding

Natural Sciences and Engineering Research Council of Canada (017-04265)

  • Bratislav Misic

National Health and Medical Research Council (3274306)

  • Alex Fornito

National Institutes of Health (NIH-NIBIB P41 EB019936)

  • Jean-Baptiste Poline

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Copyright

© 2021, Markello et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,401
    views
  • 634
    downloads
  • 233
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ross D Markello
  2. Aurina Arnatkevičiūtė
  3. Jean-Baptiste Poline
  4. Ben D Fulcher
  5. Alex Fornito
  6. Bratislav Misic
(2021)
Standardizing workflows in imaging transcriptomics with the abagen toolbox
eLife 10:e72129.
https://doi.org/10.7554/eLife.72129

Share this article

https://doi.org/10.7554/eLife.72129

Further reading

    1. Neuroscience
    Magdalena Solyga, Georg B Keller
    Research Article

    Our movements result in predictable sensory feedback that is often multimodal. Based on deviations between predictions and actual sensory input, primary sensory areas of cortex have been shown to compute sensorimotor prediction errors. How prediction errors in one sensory modality influence the computation of prediction errors in another modality is still unclear. To investigate multimodal prediction errors in mouse auditory cortex, we used a virtual environment to experimentally couple running to both self-generated auditory and visual feedback. Using two-photon microscopy, we first characterized responses of layer 2/3 (L2/3) neurons to sounds, visual stimuli, and running onsets and found responses to all three stimuli. Probing responses evoked by audiomotor (AM) mismatches, we found that they closely resemble visuomotor (VM) mismatch responses in visual cortex (V1). Finally, testing for cross modal influence on AM mismatch responses by coupling both sound amplitude and visual flow speed to the speed of running, we found that AM mismatch responses were amplified when paired with concurrent VM mismatches. Our results demonstrate that multimodal and non-hierarchical interactions shape prediction error responses in cortical L2/3.

    1. Neuroscience
    Gyeong Hee Pyeon, Hyewon Cho ... Yong Sang Jo
    Research Article Updated

    Recent studies suggest that calcitonin gene-related peptide (CGRP) neurons in the parabrachial nucleus (PBN) represent aversive information and signal a general alarm to the forebrain. If CGRP neurons serve as a true general alarm, their activation would modulate both passive nad active defensive behaviors depending on the magnitude and context of the threat. However, most prior research has focused on the role of CGRP neurons in passive freezing responses, with limited exploration of their involvement in active defensive behaviors. To address this, we examined the role of CGRP neurons in active defensive behavior using a predator-like robot programmed to chase mice. Our electrophysiological results revealed that CGRP neurons encode the intensity of aversive stimuli through variations in firing durations and amplitudes. Optogenetic activation of CGRP neurons during robot chasing elevated flight responses in both conditioning and retention tests, presumably by amplifying the perception of the threat as more imminent and dangerous. In contrast, animals with inactivated CGRP neurons exhibited reduced flight responses, even when the robot was programmed to appear highly threatening during conditioning. These findings expand the understanding of CGRP neurons in the PBN as a critical alarm system, capable of dynamically regulating active defensive behaviors by amplifying threat perception, and ensuring adaptive responses to varying levels of danger.