RNA velocity in spermatogenesis#

RNA velocity analysis with the EM and steady-state models using data preprocessed with kallisto_bustools_prepref_isoseparate_include.

Requires

  • adata_generation.ipynb

  • velocyto_var_names.csv from velocyto_vi.ipynb

Output

  • DATA_DIR/spermatogenesis/velocities/kallisto_isoseparate_include_steady_state.pickle

  • DATA_DIR/spermatogenesis/velocities/kallisto_isoseparate_include_em.pickle

Library imports#

import sys

import pandas as pd

import scanpy as sc
import scvelo as scv

sys.path.insert(0, "../../../")
from paths import DATA_DIR
sc.logging.print_version_and_date()
Running Scanpy 1.9.1, on 2022-07-21 16:46.

General settings#

# set verbosity levels
sc.settings.verbosity = 2
scv.settings.verbosity = 3
scv.settings.set_figure_params('scvelo', dpi_save=400, dpi=80, transparent=True, fontsize=20, color_map='viridis')
scv.settings.plot_prefix = ""

Constants#

N_JOBS = 8
VELOCYTO_VAR_NAMES = pd.read_csv('velocyto_var_names.csv', index_col=0, header=None).index.tolist()

Data loading#

adata = sc.read(
    DATA_DIR / 'spermatogenesis' / "kallisto_bustools_prepref_isoseparate_include.h5ad"
)
adata
AnnData object with n_obs × n_vars = 1829 × 54144
    obs: 'cell_index', 'clusters_coarse', 'clusters'
    layers: 'spliced', 'unspliced'
scv.pl.proportions(adata)
../../_images/b11671d6e5a42593871d62ad8a936b3493f8c425f0b90543fd7bc1a9f2915df9.png

Data pre-processing#

scv.pp.filter_and_normalize(adata, min_shared_counts=20, n_top_genes=2000, retain_genes=VELOCYTO_VAR_NAMES)
scv.pp.moments(adata, n_pcs=30, n_neighbors=30)
Filtered out 42338 genes that are detected 20 counts (shared).
Normalized count data: X, spliced, unspliced.
Extracted 2366 highly variable genes.
Logarithmized X.
computing PCA
    on highly variable genes
    with n_comps=30
    finished (0:00:03)
computing neighbors
    finished (0:00:27) --> added 
    'distances' and 'connectivities', weighted adjacency matrices (adata.obsp)
computing moments based on connectivities
    finished (0:00:03) --> added 
    'Ms' and 'Mu', moments of un/spliced abundances (adata.layers)

Model fitting#

Steady state#

scv.tl.velocity(adata, mode="steady_state")
adata = adata[:, VELOCYTO_VAR_NAMES].copy()
computing velocities
    finished (0:00:00) --> added 
    'velocity', velocity vectors for each individual cell (adata.layers)
pd.DataFrame(
    adata.layers['velocity'],
    index=adata.obs_names,
    columns=adata.var_names
).to_pickle(
    DATA_DIR / 'spermatogenesis' / 'velocities' / 'kallisto_isoseparate_include_steady_state.pickle'
)

EM#

scv.tl.recover_dynamics(adata, var_names='all', n_jobs=N_JOBS)
scv.tl.velocity(adata, mode='dynamical')
recovering dynamics (using 8/64 cores)
    finished (0:03:16) --> added 
    'fit_pars', fitted parameters for splicing dynamics (adata.var)
computing velocities
    finished (0:00:01) --> added 
    'velocity', velocity vectors for each individual cell (adata.layers)
pd.DataFrame(
    adata[:, ~adata.var['fit_alpha'].isnull()].layers['velocity'],
    index=adata.obs_names,
    columns=adata.var_names[~adata.var['fit_alpha'].isnull()]
).to_pickle(
    DATA_DIR / 'spermatogenesis' / 'velocities' / 'kallisto_isoseparate_include_em.pickle'
)
WARNING: 1700095A21Rik not recoverable due to insufficient samples.
WARNING: Gm15904 not recoverable due to insufficient samples.
WARNING: AI507597 not recoverable due to insufficient samples.
WARNING: Ube4bos2 not recoverable due to insufficient samples.
WARNING: Speer4b not recoverable due to insufficient samples.
WARNING: 1700010H22Rik not recoverable due to insufficient samples.
WARNING: Tmem132cos not recoverable due to insufficient samples.
WARNING: Speer3 not recoverable due to insufficient samples.
WARNING: Prkag2os1 not recoverable due to insufficient samples.
WARNING: Gm16226 not recoverable due to insufficient samples.
WARNING: Pard3bos1 not recoverable due to insufficient samples.
WARNING: Tnp1 not recoverable due to insufficient samples.
WARNING: Lypla1 not recoverable due to insufficient samples.
WARNING: Gm7568 not recoverable due to insufficient samples.
WARNING: 4930556I23Rik not recoverable due to insufficient samples.
WARNING: 4930455C13Rik not recoverable due to insufficient samples.
WARNING: 4930545H06Rik not recoverable due to insufficient samples.
WARNING: 4930404N11Rik not recoverable due to insufficient samples.
WARNING: Pln not recoverable due to insufficient samples.
WARNING: Sycp3 not recoverable due to insufficient samples.
WARNING: Gm45604 not recoverable due to insufficient samples.
WARNING: Gm29585 not recoverable due to insufficient samples.
WARNING: Ube2b not recoverable due to insufficient samples.
WARNING: AC154538.1 not recoverable due to insufficient samples.
WARNING: Rab40c not recoverable due to insufficient samples.
WARNING: Gm16555 not recoverable due to insufficient samples.
WARNING: AC165153.2 not recoverable due to insufficient samples.
WARNING: AC126942.1 not recoverable due to insufficient samples.
WARNING: Gm36368 not recoverable due to insufficient samples.
WARNING: AC164010.1 not recoverable due to insufficient samples.
WARNING: AC142100.4 not recoverable due to insufficient samples.
WARNING: 4930471G03Rik not recoverable due to insufficient samples.
WARNING: AC100726.1 not recoverable due to insufficient samples.
WARNING: Rbm4b not recoverable due to insufficient samples.
WARNING: 1700080N15Rik not recoverable due to insufficient samples.
WARNING: 1700042O10Rik not recoverable due to insufficient samples.
WARNING: Lrrc75aos1 not recoverable due to insufficient samples.
WARNING: 4930544D05Rik not recoverable due to insufficient samples.
WARNING: Dbil5 not recoverable due to insufficient samples.
WARNING: Gm47376 not recoverable due to insufficient samples.
WARNING: 4930549C15Rik not recoverable due to insufficient samples.
WARNING: Gm36756 not recoverable due to insufficient samples.
WARNING: Spata31d1d not recoverable due to insufficient samples.
WARNING: Gm48837 not recoverable due to insufficient samples.
WARNING: Gm15935 not recoverable due to insufficient samples.
WARNING: Lrrc18 not recoverable due to insufficient samples.
WARNING: 1700006F04Rik not recoverable due to insufficient samples.
WARNING: AC117700.1 not recoverable due to insufficient samples.
WARNING: 1700015H07Rik not recoverable due to insufficient samples.
WARNING: Gm6569 not recoverable due to insufficient samples.
WARNING: Gm15743 not recoverable due to insufficient samples.
WARNING: Gm15657 not recoverable due to insufficient samples.
WARNING: Efcab1 not recoverable due to insufficient samples.
WARNING: Trp53tg5 not recoverable due to insufficient samples.
WARNING: Alkbh3os1 not recoverable due to insufficient samples.
WARNING: Gm29010 not recoverable due to insufficient samples.
WARNING: 1700003F12Rik not recoverable due to insufficient samples.
WARNING: Gm42609 not recoverable due to insufficient samples.
WARNING: 4930509J09Rik not recoverable due to insufficient samples.
WARNING: 4930503B20Rik not recoverable due to insufficient samples.
WARNING: Fhl5 not recoverable due to insufficient samples.
WARNING: Sox5os3 not recoverable due to insufficient samples.
WARNING: Lipe not recoverable due to insufficient samples.
WARNING: Gm15412 not recoverable due to insufficient samples.
WARNING: Nsmce4a not recoverable due to insufficient samples.
WARNING: Gm29638 not recoverable due to insufficient samples.
WARNING: Gm26891 not recoverable due to insufficient samples.
WARNING: Etfb not recoverable due to insufficient samples.
WARNING: Klk1b7-ps not recoverable due to insufficient samples.
WARNING: Klk1b8 not recoverable due to insufficient samples.
WARNING: Gm15517 not recoverable due to insufficient samples.
WARNING: Prss23os not recoverable due to insufficient samples.
WARNING: Gm5901 not recoverable due to insufficient samples.
WARNING: Gm4275 not recoverable due to insufficient samples.
WARNING: 1700029B22Rik not recoverable due to insufficient samples.
WARNING: Fam24a not recoverable due to insufficient samples.
WARNING: Acadsb not recoverable due to insufficient samples.
WARNING: Gm16351 not recoverable due to insufficient samples.
WARNING: Lncbate1 not recoverable due to insufficient samples.
WARNING: Ankrd11 not recoverable due to insufficient samples.
WARNING: Ccsap not recoverable due to insufficient samples.
WARNING: 1700047A11Rik not recoverable due to insufficient samples.
WARNING: 1700030M09Rik not recoverable due to insufficient samples.
WARNING: 4921509A06Rik not recoverable due to insufficient samples.
WARNING: Gm4665 not recoverable due to insufficient samples.
WARNING: Gm31497 not recoverable due to insufficient samples.
WARNING: Iqcf5 not recoverable due to insufficient samples.
WARNING: Gm27002 not recoverable due to insufficient samples.