# CITEsq loading in RNA and ADT data

**URL:** <https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423>\
**Category:** scanpy\
**Created:** [May 12, 2023, 3:44pm UTC](https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423 "2023-05-12T15:44:41Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![pouria](https://avatars.discourse-cdn.com/v4/letter/p/e56c9b/32.png) [@pouria](https://discourse.scverse.org/u/pouria)\
**Post date:** [May 12, 2023, 3:44pm UTC](https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423/1 "2023-05-12T15:44:42Z")

</div>

Hello everyone,

I’m in to process of learning single cell analysis. Thus far I’ve within my first week made good progress in analysing single cell data by using Scanpy.

However I’ve hade som trouble analysing a publicly available CITE-seq dataset. My problem is very basic, the data is uploaded in this manner:

ADT\_barcodes.tsv.gz|  
ADT\_counts.txt.gz  
ADT\_tags.tsv.gz|  
barcodes.tsv.gz|  
features.tsv.gz|  
hashtags.csv.gz|  
matrix.mtx.gz|

I’ve only found tutorials showing how to load h5.files containing the above, but I’ve had a hard time figuring out how to create an AnnData containing all of the above.

I would be very happy for any of your suggestions.

---

<div class="post-metadata">

**Author:** ![bennjmo](https://avatars.discourse-cdn.com/v4/letter/b/edb3f5/32.png) [@bennjmo](https://discourse.scverse.org/u/bennjmo)\
**Post date:** [May 19, 2023, 7:44pm UTC](https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423/2 "2023-05-19T19:44:14Z")

</div>

I just realized I made a post with a very similar problem as yours. Did you figure out a solution?

---

<div class="post-metadata">

**Author:** ![pouria](https://avatars.discourse-cdn.com/v4/letter/p/e56c9b/32.png) [@pouria](https://discourse.scverse.org/u/pouria)\
**Post date:** [July 5, 2023, 4:23pm UTC](https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423/3 "2023-07-05T16:23:20Z")

</div>

Sorry, for late response. yet to find a solution, haven’t had time to work on it since my last message. Please update me if you’ve found a way.

---

<div class="post-metadata">

**Author:** ![bennjmo](https://avatars.discourse-cdn.com/v4/letter/b/edb3f5/32.png) [@bennjmo](https://discourse.scverse.org/u/bennjmo)\
**Post date:** [July 5, 2023, 4:35pm UTC](https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423/4 "2023-07-05T16:35:45Z")

</div>

> [@Help concatenating var for cite seq](https://discourse.scverse.org/t/help-concatenating-var-for-cite-seq/1447):
>
> I have a CITE seq experiment where the GEX and antibody capture were mapped separately, so that I have to load in two different matrices. adata\_adt = sc.read\_10x\_mtx( "data/CITESeq\_381/", gex\_only = False ) adata\_gex = sc.read\_10x\_mtx( "data/GEX\_381/filtered\_feature\_bc\_matrix/", gex\_only = True ) I subset to only cells with GEX and ADT common\_cells = np.intersect1d(adata\_adt.obs\_names, adata\_gex.obs\_names) adata\_adt = adata\_adt[common\_cells] adata\_gex = adata\_gex[common\_cells…

---

<div class="post-metadata">

**Author:** ![ddiez](https://yyz1.discourse-cdn.com/flex035/user_avatar/discourse.scverse.org/ddiez/32/590_2.png) [@ddiez](https://discourse.scverse.org/u/ddiez)\
**Post date:** [July 6, 2023, 2:17pm UTC](https://discourse.scverse.org/t/citesq-loading-in-rna-and-adt-data/1423/5 "2023-07-06T14:17:40Z")

</div>

Basically you have a multiomics experiment so why not use MuData and muon for this? You can try:

```auto
import muon as mu

mdata = mu.MuData({
  "rna": adata_gex,
  "prot": adata_adt
})

```

(Of course, use any names you like in place of “rna” and “prot” in the dict).

It is possible that you do not even need to subset the objects with the intersection of barcodes but depending on your specific problem that might be what you want to do anyway.

You can check the muon website on how to use them:

[Muon](https://muon.readthedocs.io/en/latest/index.html)

And the muon/mdata categories in this site for more specific help on those tools.

Hope this helps!