Skip to content

Code to reproduce analyses of mouse cortical gradients

Notifications You must be signed in to change notification settings

benfulcher/mouseGradients

Repository files navigation

mouseGradients

This repository contains code to reproduce all figures from our paper:

Fulcher, B. D., Murray, J. D., Zerbi, V., & Wang, X.-J. (2019). Multimodal gradients across mouse cortex. PNAS, 201814144.

Downloading processed data

You will need to fill the Data directory with some processed data that you can download from this figshare repository.

This paper would not have been possible without the generous sharing of open data (thank you!). Please consult each original data source for detailed licensing information for each dataset.

Allen Institute gene-expression data was retrieved using custom code that queries the Allen SDK.

Figures

Figure 1A generated manually mouse-brain imaging software (by Valerio Zerbi).

Figure 1B: T1w:T2w across all brain areas

FunctionalClasses()

Results

Figure 1C-F: Scatter plots versus T1w:T2w

T1T2Corrs()

Results

Figure 2A,B: Zooming in on some manually-selected genes

SpecificGeneFigure()

Note that B was removed from the final paper; which only plots the Grin3a scatter Results

Figure 2C: Low-dimensional projection

MassAction('brainExpressed','scaledSigmoid','benCombo','bayesian')

Note that additional outputs are also produced, visualizing analytics about the principal components, including Fig.S4. Results

Figure 3: Correlations across individual cortical layers

LaminarPlot()

Results

Figure 4: Combination of diverse measurements

DiverseMatrixT1T2()
PlotBigDamnMatrix()

Results

Figure 5: Interspecies comparison

G = LoadMeG('cortexAll');
HumanMouseComparisonAnal(G,'brainRelated','benCombo',true)

Results

Supplementary Figures

Figs S1, S2, S3 generated manually using mouse-brain imaging software (by Valerio Zerbi).

Figure S4: Correlations of T1w:T2w with Principal Component projections of brain-expressed genes

MassAction('brainExpressed','scaledSigmoid','benCombo','bayesian')

Results

Figure S5: Correlations between T1w:T2w values measured in individual cortical layers

T1T2ByLayer

Results

Figure S6: Correlations to Ero et al. (2019) cell-density estimates

MyelinCorrs('CellAtlas')

Results

Figure S7: Layer-specific gene-expression gradients

TheBestLayers('genes',true,[]);

Results

Figure S8: Comparison of Harris hierarchy and T1w:T2w

T1T2_or_hierarchy

Results

Figure S9: Comparison of mouse-human correspondence as a function of expression processing and gene-set filtering

whatGeneSets = {'all','brainExpressed','brainRelated','CahoyNeuron','CahoyOgligodendrocyte','CahoyAstrocyte','myelinSetOf999'};
whatSectionFilters = {'sagittal','coronal','combZ','benCombo','replicated'};
HumanMouseSectionCompare(whatGeneSets,whatSectionFilters);

Results

About

Code to reproduce analyses of mouse cortical gradients

Topics

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published