Back

Strategies for Integrating Single-Cell RNA Sequencing Results With Multiple Species

Hart, R. P.

2019-06-14 bioinformatics
10.1101/671115 bioRxiv
Show abstract

Single-cell RNA sequencing (scRNAseq) is a robust technology for parsing gene expression in individual cells from a tissue or other complex source. One application involves experiments where cells from multiple species are recovered from a single sample, such as when human cells are transplanted into an animal model. We transplanted microglial precursor cells into newborn mouse brain and then recovered unenriched cortical tissue six months later. Dissociated cells were assessed by scRNAseq. The default method for analyzing these results begins by aligning sequencing reads with a mixture of both mouse and human reference genomes. While this clearly identifies the human cells as a distinct cluster, the clustering is artificially driven by expression from non-comparable gene identifiers from different species. We devised a method for translating expression counts from human to mouse and evaluated four algorithms for parsing mixed-species scRNAseq data. Our optimal approach split raw sequencing reads according to the best alignment score in each genome, and then re-aligned reads only with the appropriate genome. After gene symbol translation, pooled results indicate that cell types are more appropriately clustered and that differential expression analysis identifies species-specific patterns. This method should be applicable to any mixed-species scRNAseq experiment.\n\nSummary of optimal strategyO_LIMixed-species scRNAseq data are aligned with mixture of mouse and human reference genomes\nC_LIO_LIThe BAM file is scanned to find the best alignment score for each sequencing read identifier; these are used to split the paired FASTQ files into two sets of files\nC_LIO_LIEach set of species-specific, paired FASTQ files is re-aligned with only the appropriate reference genome\nC_LIO_LIRaw counts imported into Seurat\nC_LIO_LIThe human counts table is translated to mouse gene symbols using a custom HomoloGene translation table\nC_LIO_LIResults are merged and analyzed\nC_LI

Matching journals

The top 6 journals account for 50% of the predicted probability mass.

1
BMC Bioinformatics
383 papers in training set
Top 0.7%
12.0%
2
Bioinformatics
1061 papers in training set
Top 2%
12.0%
3
Nature Protocols
30 papers in training set
Top 0.1%
8.2%
4
Genome Research
409 papers in training set
Top 0.3%
7.0%
5
Genome Biology
555 papers in training set
Top 1.0%
6.6%
6
Cell Reports Methods
141 papers in training set
Top 0.4%
6.1%
50% of probability mass above
7
STAR Protocols
15 papers in training set
Top 0.1%
4.2%
8
Bioinformatics Advances
184 papers in training set
Top 2%
3.5%
9
Nature Methods
336 papers in training set
Top 3%
2.6%
10
PLOS ONE
4510 papers in training set
Top 49%
2.0%
11
eneuro
389 papers in training set
Top 4%
2.0%
12
Scientific Reports
3102 papers in training set
Top 52%
2.0%
13
F1000Research
79 papers in training set
Top 1%
2.0%
14
Nature Communications
4913 papers in training set
Top 47%
2.0%
15
Cell Systems
167 papers in training set
Top 6%
2.0%
16
BMC Genomics
328 papers in training set
Top 2%
1.8%
17
Nature Biotechnology
147 papers in training set
Top 5%
1.6%
18
iScience
1063 papers in training set
Top 16%
1.6%
19
Briefings in Bioinformatics
326 papers in training set
Top 4%
1.6%
20
Frontiers in Genetics
197 papers in training set
Top 6%
1.3%
21
Nucleic Acids Research
1128 papers in training set
Top 13%
1.3%
22
RNA
169 papers in training set
Top 0.3%
1.2%
23
Life Science Alliance
263 papers in training set
Top 0.9%
0.9%
24
PeerJ
261 papers in training set
Top 12%
0.9%
25
eLife
5422 papers in training set
Top 59%
0.7%
26
GigaScience
172 papers in training set
Top 3%
0.7%
27
Journal of Open Source Software
22 papers in training set
Top 0.3%
0.6%
28
NAR Genomics and Bioinformatics
214 papers in training set
Top 4%
0.6%
29
PLOS Computational Biology
1633 papers in training set
Top 28%
0.6%