biopipen 0.31.7__py3-none-any.whl → 0.32.0__py3-none-any.whl

biopipen/scripts/tcr/GIANA/query.py

@@ -5,220 +5,222 @@ import subprocess as sp
 import pandas as pd
 from GIANA4 import *
 
-def CreateReference(rFile, outdir='./', Vgene=True, ST=3):
+
+def CreateReference(rFile, outdir="./", Vgene=True, ST=3):
     ## convert input reference file into a python workplace
-    h=open(rFile)
-    alines=h.readlines()
-    ww=alines[0].strip().split('\t')
-    if not ww[0].startswith('C'):
+    h = open(rFile)
+    alines = h.readlines()
+    ww = alines[0].strip().split("\t")
+    if not ww[0].startswith("C"):
         ## header line
-        hline=alines[0]
-        alines=alines[1:]        
-    elif 'CDR3' in ww[0]:
-        hline=alines[0]
-        alines=alines[1:]
+        hline = alines[0]
+        alines = alines[1:]
+    elif "CDR3" in ww[0]:
+        hline = alines[0]
+        alines = alines[1:]
     else:
-        hline='CDR3\t'+'\t'.join(['Info'+str(x) for x in range(len(ww)-1)])        
-    seqs=[]
-    vgs=[]
-    infoList=[]
-    count=0
+        hline = "CDR3\t" + "\t".join(["Info" + str(x) for x in range(len(ww) - 1)])
+    seqs = []
+    vgs = []
+    infoList = []
+    count = 0
     for ll in alines:
-        ww=ll.strip().split('\t')
-        cdr3=ww[0]
-        if '*' in cdr3:
+        ww = ll.strip().split("\t")
+        cdr3 = ww[0]
+        if "*" in cdr3:
             continue
-        if '_' in cdr3:
+        if "_" in cdr3:
             continue
         seqs.append(ww[0])
         if Vgene:
             vgs.append(ww[1])
-            infoList.append('\t'.join(ww[1:]))
+            infoList.append("\t".join(ww[1:]))
         else:
-            infoList.append('\t'.join(ww[1:]))
-        count+=1
-    LD,VD, ID,SD= BuildLengthDict(seqs, vGene=vgs,INFO=infoList,sIDs=[x for x in range(len(seqs))])
+            infoList.append("\t".join(ww[1:]))
+        count += 1
+    LD, VD, ID, SD = BuildLengthDict(
+        seqs, vGene=vgs, INFO=infoList, sIDs=[x for x in range(len(seqs))]
+    )
     LDu_r, VDu_r, IDu_r, SDu_r = CollapseUnique(LD, VD, ID, SD)
-    flagLD_r={}
-    dMD_r={}
+    flagLD_r = {}
+    dMD_r = {}
     for kk in LDu_r:
-        vss=SDu_r[kk]
-        vInfo=IDu_r[kk]
-        flagL=[len(x)-1 for x in vInfo]
-        flagLD_r[kk]=flagL
-        dM=np.array([EncodingCDR3(x[ST:-2], M6, n0) for x in vss])
-        dM=dM.astype("float32")        
-        dMD_r[kk]=dM
-##    ff0=re.sub('.txt','',rFile)
-##    outfile=outdir+ff0+'_giana_ref.shelve'
-##    giana_shelf = shelve.open(outfile, 'n')
-##    giana_shelf['flagLD']=flagLD_r
-##    giana_shelf['dMD']=dMD_r
-##    giana_shelf['LDu']=LDu_r
-##    giana_shelf['VDu']=VDu_r
-##    giana_shelf['IDu']=IDu_r
-##    giana_shelf['SDu']=SDu_r
-##    giana_shelf.close()
+        vss = SDu_r[kk]
+        vInfo = IDu_r[kk]
+        flagL = [len(x) - 1 for x in vInfo]
+        flagLD_r[kk] = flagL
+        dM = np.array([EncodingCDR3(x[ST:-2], M6, n0) for x in vss])
+        dM = dM.astype("float32")
+        dMD_r[kk] = dM
+    ##    ff0=re.sub('.txt','',rFile)
+    ##    outfile=outdir+ff0+'_giana_ref.shelve'
+    ##    giana_shelf = shelve.open(outfile, 'n')
+    ##    giana_shelf['flagLD']=flagLD_r
+    ##    giana_shelf['dMD']=dMD_r
+    ##    giana_shelf['LDu']=LDu_r
+    ##    giana_shelf['VDu']=VDu_r
+    ##    giana_shelf['IDu']=IDu_r
+    ##    giana_shelf['SDu']=SDu_r
+    ##    giana_shelf.close()
     return [LDu_r, VDu_r, IDu_r, SDu_r, dMD_r]
 
-def MakeQuery(qFile, rData=[],dbFile=None, Vgene=True, thr=7, ST=3, thr_s=3.3):
+
+def MakeQuery(qFile, rData=[], dbFile=None, Vgene=True, thr=7, ST=3, thr_s=3.3):
     if dbFile is not None:
         with shelve.open(dbFile) as db:
             for key in db:
-                globals()[key]=db[key]
+                globals()[key] = db[key]
     else:
-        if len(rData)==0:
-            raise("Need to provide either a reference file or a shelve")
-        LDu_r=rData[0]
-        VDu_r=rData[1]
-        IDu_r=rData[2]
-        SDu_r=rData[3]
-        dMD_r=rData[4]
-    h=open(qFile)
-    alines=h.readlines()
-    ww=alines[0].strip().split('\t')
-    if not ww[0].startswith('C'):
+        if len(rData) == 0:
+            raise ("Need to provide either a reference file or a shelve")
+        LDu_r = rData[0]
+        VDu_r = rData[1]
+        IDu_r = rData[2]
+        SDu_r = rData[3]
+        dMD_r = rData[4]
+    h = open(qFile)
+    alines = h.readlines()
+    ww = alines[0].strip().split("\t")
+    if not ww[0].startswith("C"):
         ## header line
-        hline=alines[0]
-        alines=alines[1:]        
-    elif 'CDR3' in ww[0]:
-        hline=alines[0]
-        alines=alines[1:]
+        hline = alines[0]
+        alines = alines[1:]
+    elif "CDR3" in ww[0]:
+        hline = alines[0]
+        alines = alines[1:]
     else:
-        hline='CDR3\t'+'\t'.join(['Info'+str(x) for x in range(len(ww)-1)])        
-    seqs=[]
-    vgs=[]
-    infoList=[]
-    count=0
+        hline = "CDR3\t" + "\t".join(["Info" + str(x) for x in range(len(ww) - 1)])
+    seqs = []
+    vgs = []
+    infoList = []
+    count = 0
     for ll in alines:
-        ww=ll.strip().split('\t')
-        cdr3=ww[0]
-        if '*' in cdr3:
+        ww = ll.strip().split("\t")
+        cdr3 = ww[0]
+        if "*" in cdr3:
             continue
-        if '_' in cdr3:
+        if "_" in cdr3:
             continue
         seqs.append(ww[0])
         if Vgene:
             vgs.append(ww[1])
-            infoList.append('\t'.join(ww[1:]))
+            infoList.append("\t".join(ww[1:]))
         else:
-            infoList.append('\t'.join(ww[1:]))
-        count+=1
-    LD,VD, ID,SD= BuildLengthDict(seqs, vGene=vgs,INFO=infoList,sIDs=[x for x in range(len(seqs))])
+            infoList.append("\t".join(ww[1:]))
+        count += 1
+    LD, VD, ID, SD = BuildLengthDict(
+        seqs, vGene=vgs, INFO=infoList, sIDs=[x for x in range(len(seqs))]
+    )
     LDu, VDu, IDu, SDu = CollapseUnique(LD, VD, ID, SD)
-    tmpFile='tmp_query.txt'
-    g=open(tmpFile,'w')
+    tmpFile = "tmp_query.txt"
+    g = open(tmpFile, "w")
     for kk in LDu:
-        vss=SDu[kk]
-        vInfo=IDu[kk]
-        vss_r=SDu_r[kk]
-        vInfo_r=IDu_r[kk]
-        flagL=[len(x)-1 for x in vInfo]
-        dM_r=dMD_r[kk]
-        dM=np.array([EncodingCDR3(x[ST:-2], M6, n0) for x in vss])
-        dM=dM.astype("float32")
-        nq=dM.shape[0]
-        nr=dM_r.shape[0]
-        vssc=vss+vss_r
-        vInfoc=vInfo+vInfo_r
-        dMc=np.concatenate((dM, dM_r))
-        index = faiss.IndexFlatL2(Ndim*6)
+        vss = SDu[kk]
+        vInfo = IDu[kk]
+        vss_r = SDu_r[kk]
+        vInfo_r = IDu_r[kk]
+        flagL = [len(x) - 1 for x in vInfo]
+        dM_r = dMD_r[kk]
+        dM = np.array([EncodingCDR3(x[ST:-2], M6, n0) for x in vss])
+        dM = dM.astype("float32")
+        nq = dM.shape[0]
+        nr = dM_r.shape[0]
+        vssc = vss + vss_r
+        vInfoc = vInfo + vInfo_r
+        dMc = np.concatenate((dM, dM_r))
+        index = faiss.IndexFlatL2(Ndim * 6)
         index.add(dMc)
         D, I = index.search(dM, 2)
-        vv=np.where((D[0:nq,1]<=thr))[0]
-        flagL=np.array(flagL)
-        vv0=np.where((D[0:nq,1]>thr) & (flagL>0))[0]
-        curList=[]
+        vv = np.where((D[0:nq, 1] <= thr))[0]
+        flagL = np.array(flagL)
+        vv0 = np.where((D[0:nq, 1] > thr) & (flagL > 0))[0]
+        curList = []
         for v in vv0:
             for ii in range(len(vInfoc[v])):
-                line=vssc[v]+'\t'+vInfoc[v][ii]+'\t'+'query\n'
-                _=g.write(line)
+                line = vssc[v] + "\t" + vInfoc[v][ii] + "\t" + "query\n"
+                _ = g.write(line)
         for v in vv:
-            tmpI=I[v,]
+            tmpI = I[v,]
             if v not in tmpI:
-                tmpI[0]=v
-            idx1=tmpI[0]
-            idx2=tmpI[1]
-            c1=vssc[idx1]
-            c2=vssc[idx2]
-            info1=vInfoc[idx1]
-            info2=vInfoc[idx2]
+                tmpI[0] = v
+            idx1 = tmpI[0]
+            idx2 = tmpI[1]
+            c1 = vssc[idx1]
+            c2 = vssc[idx2]
+            info1 = vInfoc[idx1]
+            info2 = vInfoc[idx2]
             for tmpInfo in info1:
-                tup1=(c1, tmpInfo)
+                tup1 = (c1, tmpInfo)
                 if tup1 not in curList:
-                    if idx1<nq:
-                        line1=c1+'\t'+tmpInfo+'\t'+'query\n'
+                    if idx1 < nq:
+                        line1 = c1 + "\t" + tmpInfo + "\t" + "query\n"
                     else:
-                        line1=c1+'\t'+tmpInfo+'\t'+'ref\n'
-                    _=g.write(line1)
+                        line1 = c1 + "\t" + tmpInfo + "\t" + "ref\n"
+                    _ = g.write(line1)
                     curList.append(tup1)
             for tmpInfo in info2:
-                tup2=(c2, tmpInfo)
+                tup2 = (c2, tmpInfo)
                 if tup2 not in curList:
-                    if idx2<nq:
-                        line2=c2+'\t'+tmpInfo+'\t'+'query\n'
+                    if idx2 < nq:
+                        line2 = c2 + "\t" + tmpInfo + "\t" + "query\n"
                     else:
-                        line2=c2+'\t'+tmpInfo+'\t'+'ref\n'
-                    _=g.write(line2)
+                        line2 = c2 + "\t" + tmpInfo + "\t" + "ref\n"
+                    _ = g.write(line2)
                     curList.append(tup2)
     g.close()
-    cmd='python3 GIANA4.1.py -f tmp_query.txt -S '+str(thr_s) ## updated to GIANA4.1
-    p=sp.run(cmd, shell=True)
-
-def MergeExist(refClusterFile, outFile='queryFinal.txt',queryClusterFile='tmp_query--RotationEncodingBL62.txt', direction='q'):
+    cmd = "python3 GIANA4.1.py -f tmp_query.txt -S " + str(
+        thr_s
+    )  ## updated to GIANA4.1
+    p = sp.run(cmd, shell=True)
+
+
+def MergeExist(
+    refClusterFile,
+    outFile="queryFinal.txt",
+    queryClusterFile="tmp_query--RotationEncodingBL62.txt",
+    direction="q",
+):
     ## This function compare the query file with ref cluster file and merge the two based on shared TCRs
     ## If direction is 'q', the overlapping clusters will be added to the query file
     ## If direction is 'r', the overlapping and non-overlapping clusters will be added to the reference file
-    refT=pd.read_table(refClusterFile, skiprows=2, delimiter='\t', header=None)
-    queryT=pd.read_table(queryClusterFile, skiprows=2, delimiter='\t', header=None)
-    nq=queryT.shape[1]
-    nr=refT.shape[1]
-    if nr != nq-1:
+    refT = pd.read_table(refClusterFile, skiprows=2, delimiter="\t", header=None)
+    queryT = pd.read_table(queryClusterFile, skiprows=2, delimiter="\t", header=None)
+    nq = queryT.shape[1]
+    nr = refT.shape[1]
+    if nr != nq - 1:
         print("ERROR: Make sure reference and the query samples have the same columns!")
         print("No query file is generated.")
         return
-    gn=np.unique(queryT[1])
-    queryTs=pd.DataFrame([], columns=queryT.columns)
+    gn = np.unique(queryT[1])
+    queryTs = pd.DataFrame([], columns=queryT.columns)
     for nn in gn:
-            tmp_ddq=queryT.loc[np.where(queryT[1]==nn)[0],:]
-            cls_lab=np.unique(tmp_ddq[nq-1])
-            if len(cls_lab)==1:
-                if cls_lab[0]=='ref':
-                    continue
-            queryTs=queryTs.append(tmp_ddq)
-    queryTs.index=range(queryTs.shape[0])
-    keyr=refT[0]+'_'+refT[2]
-    keyq=queryTs[0]+'_'+queryTs[2]
-    vvr=np.where(queryTs[nq-1]=='ref')[0]
-    vvr_in=np.where(keyr.isin(keyq[vvr]))[0]
-    gn_r=list(refT.loc[vvr_in,1].drop_duplicates())
-    ddo=pd.DataFrame([], columns=refT.columns)
+        tmp_ddq = queryT.loc[np.where(queryT[1] == nn)[0], :]
+        cls_lab = np.unique(tmp_ddq[nq - 1])
+        if len(cls_lab) == 1:
+            if cls_lab[0] == "ref":
+                continue
+        queryTs = queryTs.append(tmp_ddq)
+    queryTs.index = range(queryTs.shape[0])
+    keyr = refT[0] + "_" + refT[2]
+    keyq = queryTs[0] + "_" + queryTs[2]
+    vvr = np.where(queryTs[nq - 1] == "ref")[0]
+    vvr_in = np.where(keyr.isin(keyq[vvr]))[0]
+    gn_r = list(refT.loc[vvr_in, 1].drop_duplicates())
+    ddo = pd.DataFrame([], columns=refT.columns)
     for nn in gn_r:
-        tmp_dd=refT.loc[np.where(refT[1]==nn)[0],:]
-        tmpkey=tmp_dd[0]+'_'+tmp_dd[2]
-        vv=np.where(keyq.isin(tmpkey))[0][0]
-        gq=queryTs[1][vv]
-        tmp_dd[1]=gq
-        ddo=ddo.append(tmp_dd)
-    if direction=='q':
-        ddo[nq-1]='ref'
+        tmp_dd = refT.loc[np.where(refT[1] == nn)[0], :]
+        tmpkey = tmp_dd[0] + "_" + tmp_dd[2]
+        vv = np.where(keyq.isin(tmpkey))[0][0]
+        gq = queryTs[1][vv]
+        tmp_dd[1] = gq
+        ddo = ddo.append(tmp_dd)
+    if direction == "q":
+        ddo[nq - 1] = "ref"
         ## remove groups that contain only ref group
-        queryTs=queryTs.append(ddo)
-        queryTs=queryTs.drop_duplicates()
-        queryTs.to_csv(outFile, sep='\t',header=False,index=False)
-#    queryTs.index=range(queryTs.shape[0])
-    if direction=='r':
+        queryTs = queryTs.append(ddo)
+        queryTs = queryTs.drop_duplicates()
+        queryTs.to_csv(outFile, sep="\t", header=False, index=False)
+    #    queryTs.index=range(queryTs.shape[0])
+    if direction == "r":
         ## to be developed
         pass
-
-
-
-
-
-
-
-
-
-
-
-

biopipen/scripts/tcr/TCRClustering.R

@@ -60,10 +60,12 @@ prepare_clustcr = function(clustcr_dir) {
     }
     clustcr_source = '
 import sys
-import pandas as pd
-import clustcr
 import atexit
 
+import pandas as pd
+from scipy import sparse as scipy_sparse
+
+
 @atexit.register
 def clustcr_exit():
     import pandas as pd
@@ -78,13 +80,32 @@ def clustcr_exit():
     sys.stderr.write(f"- sklearn: {sklearn.__version__}\\n")
     sys.stderr.write(f"- matplotlib: {matplotlib.__version__}\\n")
 
+
+# Monkey-patch scipy.sparse.isspmatrix to adopt latest scipy v1.14
+# If not, an error is raised:
+#   numpy.linalg.LinAlgError: 0-dimensional array given.
+#   Array must be at least two-dimensional
+scipy_sparse.isspmatrix = lambda x: isinstance(
+    x,
+    (
+        scipy_sparse.spmatrix,
+        scipy_sparse.csr_array,
+        scipy_sparse.csr_matrix,
+        scipy_sparse.csc_array,
+        scipy_sparse.csc_matrix,
+    ),
+)
+
+
+import clustcr  # noqa: #402
+
 clustcr_dir, clustcr_infile = sys.argv[1:3]
 cdr3df = pd.read_csv(clustcr_infile, index_col=None)
 cdr3 = cdr3df.iloc[:, 0]
 
-clustering = clustcr.Clustering(%s)
+clustering = clustcr.Clustering()
 output = clustering.fit(cdr3)
-output.clusters_df.to_csv(clustcr_dir + "/clusters.txt", sep="\\t", index=False)
+output.clusters_df.to_csv(clustcr_dir + "/clusters.txt", sep="\t", index=False)
 '
     clustcr_file = file.path(clustcr_dir, "_clustcr.py")
     cat(sprintf(clustcr_source, clustering_args), file=clustcr_file)

biopipen/utils/single_cell.R

@@ -1,7 +1,7 @@
 suppressPackageStartupMessages(library(rlang))
 suppressPackageStartupMessages(library(dplyr))
 suppressPackageStartupMessages(library(tidyr))
-suppressPackageStartupMessages(library(immunarch))
+try(suppressPackageStartupMessages(library(immunarch)))
 
 #' Expand a Immunarch object into cell-level
 #'
@@ -114,3 +114,94 @@ immdata_from_expanded <- function(
     )
     out
 }
+
+#' Convert Seurat object to Anndata
+#'
+#' @param sobjfile Seurat object file
+#' @param outfile Output file
+#' @param assay Assay to be used
+#'
+#' @export
+seurat_to_anndata <- function(sobjfile, outfile, assay = NULL, log_info, tmpdir = NULL, log_indent = "") {
+    library(Seurat)
+    library(SeuratDisk)
+    library(hdf5r)
+    if (endsWith(sobjfile, ".rds") || endsWith(sobjfile, ".RDS")) {
+        library(digest)
+
+        dig <- digest::digest(sobjfile, algo = "md5")
+        dig <- substr(dig, 1, 8)
+        assay_name <- ifelse(is.null(assay), "", paste0("_", assay))
+        tmpdir <- tmpdir %||% dirname(outfile)
+        dir.create(tmpdir, showWarnings = FALSE)
+        h5seurat_file <- file.path(
+            tmpdir,
+            paste0(
+                tools::file_path_sans_ext(basename(outfile)),
+                assay_name, ".", dig, ".h5seurat"
+            )
+        )
+        if (file.exists(h5seurat_file) &&
+            (file.mtime(h5seurat_file) < file.mtime(sobjfile))) {
+            file.remove(h5seurat_file)
+        }
+        if (!file.exists(h5seurat_file)) {
+            log_info("{log_indent}Reading RDS file ...")
+            sobj <- readRDS(sobjfile)
+            assay <- assay %||% DefaultAssay(sobj)
+            # In order to convert to h5ad
+            # https://github.com/satijalab/seurat/issues/8220#issuecomment-1871874649
+            sobj$RNAv3 <- as(object = sobj[[assay]], Class = "Assay")
+            DefaultAssay(sobj) <- "RNAv3"
+            sobj$RNA <- NULL
+            sobj <- RenameAssays(sobj, RNAv3 = "RNA")
+
+            log_info("{log_indent}Saving to H5Seurat file ...")
+            SaveH5Seurat(sobj, h5seurat_file)
+            rm(sobj)
+            gc()
+            sobjfile <- h5seurat_file
+        } else {
+            log_info("{log_indent}Using existing H5Seurat file ...")
+        }
+    }
+
+    if (!endsWith(sobjfile, ".h5seurat")) {
+        stop(paste0("Unknown input file format: ",
+            tools::file_ext(sobjfile),
+            ". Supported formats: .rds, .RDS, .h5seurat"))
+    }
+
+    log_info("{log_indent}Converting to Anndata ...")
+    Convert(sobjfile, dest = outfile, assay = assay %||% "RNA", overwrite = TRUE)
+
+    log_info("{log_indent}Fixing categorical data ...")
+    # See: https://github.com/mojaveazure/seurat-disk/issues/183
+    H5.create_reference <- function(self, ...) {
+        space <- self$get_space()
+        do.call("[", c(list(space), list(...)))
+        ref_type <- hdf5r::h5const$H5R_OBJECT
+        ref_obj <- hdf5r::H5R_OBJECT$new(1, self)
+        res <- .Call("R_H5Rcreate", ref_obj$ref, self$id, ".", ref_type,
+                    space$id, FALSE, PACKAGE = "hdf5r")
+        if (res$return_val < 0) {
+            stop("Error creating object reference")
+        }
+        ref_obj$ref <- res$ref
+        return(ref_obj)
+    }
+
+    h5ad <- H5File$new(outfile, "r+")
+    cats <- names(h5ad[["obs/__categories"]])
+    for (cat in cats) {
+        catname <- paste0("obs/__categories/", cat)
+        obsname <- paste0("obs/", cat)
+        ref <- H5.create_reference(h5ad[[catname]])
+        h5ad[[obsname]]$create_attr(
+            attr_name = "categories",
+            robj = ref,
+            space = H5S$new(type = "scalar")
+        )
+    }
+    h5ad$close()
+}

{biopipen-0.31.7.dist-info → biopipen-0.32.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: biopipen
-Version: 0.31.7
+Version: 0.32.0
 Summary: Bioinformatics processes/pipelines that can be run from `pipen run`
 License: MIT
 Author: pwwang

{biopipen-0.31.7.dist-info → biopipen-0.32.0.dist-info}/RECORD

@@ -1,4 +1,4 @@
-biopipen/__init__.py,sha256=APQVRwZptBFPacKCHqg_tW4g4--qdUtMJoA6GprSuSI,23
+biopipen/__init__.py,sha256=s1FApnNsg-mmWP3shY6ep6EAjZrNFXHlyBXWcP2KV6Q,23
 biopipen/core/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
 biopipen/core/config.py,sha256=edK5xnDhM8j27srDzsxubi934NMrglLoKrdcC8qsEPk,1069
 biopipen/core/config.toml,sha256=7IXvviRicZ2D1h6x3BVgbLJ96nsh-ikvZ0sVlQepqFE,1944
@@ -22,7 +22,7 @@ biopipen/ns/plot.py,sha256=XzLq0A8qCIQRbxhPEdWhEfbRZ8g3e4KriVz0RP8enNY,18078
 biopipen/ns/protein.py,sha256=33pzM-gvBTw0jH60mvfqnriM6uw2zj3katZ82nC9owI,3309
 biopipen/ns/regulatory.py,sha256=gJjGVpJrdv-rg2t5UjK4AGuvtLNymaNYNvoD8PhlbvE,15929
 biopipen/ns/rnaseq.py,sha256=bKAa6friFWof4yDTWZQahm1MS-lrdetO1GqDKdfxXYc,7708
-biopipen/ns/scrna.py,sha256=fXP_h7gchcuk_Jwos0IgY_P8ON6Q995OgKHgdrxfvAY,112868
+biopipen/ns/scrna.py,sha256=viIfqTH_0f7kSpHpKfrUf7HUr_w-V_oRYRj1TwuXYNs,120878
 biopipen/ns/scrna_metabolic_landscape.py,sha256=6AhaynGG3lNRi96N2tReVT46BJMuEwooSSd2irBoN80,28347
 biopipen/ns/snp.py,sha256=-Jx5Hsv_7KV7TqLU0nHCaPkMEN0CFdi4tNVlyq0rUZ4,27259
 biopipen/ns/stats.py,sha256=DlPyK5Vsg6ZEkV9SDS3aAw21eXzvOHgqeZDkXPhg7go,20509
@@ -46,6 +46,7 @@ biopipen/reports/delim/SampleInfo.svelte,sha256=Vpjrdd3AXNDNhx2pzDaFA61xPMYAFKvv
 biopipen/reports/gsea/FGSEA.svelte,sha256=3gNbZovmRWK6QjtxfgZgmft6LUurVZfQyHBRqyGMosk,405
 biopipen/reports/gsea/GSEA.svelte,sha256=lYHf8h8RLx3i-jNCEGu_LM-dbYm9ZJDzyAEadsZXzmQ,417
 biopipen/reports/protein/ProdigySummary.svelte,sha256=WEYPwmcmgtSqpFpTvNPV40yZR-7ERa5LgZni2RXxUZQ,348
+biopipen/reports/scrna/CellCellCommunicationPlots.svelte,sha256=McLQtgzbJGnb9o3GC17sRee96e-naPNPIF8SQyRp-Kc,346
 biopipen/reports/scrna/CellsDistribution.svelte,sha256=Mg6P0tazpzOxsOAtc-0LGEim5KprPt5KIgoW3TzrxxA,413
 biopipen/reports/scrna/DimPlots.svelte,sha256=ZLbtN0ioevRyEky4jb_DkDGAcpy_jAhaHfFym5ELEPM,479
 biopipen/reports/scrna/MarkersFinder.svelte,sha256=77rD1psj0VJykPDhfwS-B8mubvaasREAE6RYR2atTN4,444
@@ -53,7 +54,7 @@ biopipen/reports/scrna/MetaMarkers.svelte,sha256=iIFRKjvVYrM1AtDWqq8UfeS8q23R8FK
 biopipen/reports/scrna/RadarPlots.svelte,sha256=g_fp9d3vdnzk-egXPhkhhfWXOeG569Rj8rYLRIKmlLc,396
 biopipen/reports/scrna/ScFGSEA.svelte,sha256=Gqt-XjqsB3XgdR3XukvphwyMExZpScwqgEo7AD-gK6g,491
 biopipen/reports/scrna/SeuratClusterStats.svelte,sha256=rTsoOCiVm7lEyBYgu2ZVcKNGi4g6PynYoOSx0DHdTHo,429
-biopipen/reports/scrna/SeuratMap2Ref.svelte,sha256=SYuoKJR1PuBf9c7WvUXXjmokVo-JxFLqZqiu9kTtLuY,542
+biopipen/reports/scrna/SeuratMap2Ref.svelte,sha256=1g-FEmvR2YMn0OF66qESNBWmGdHQY3u8QuG4M-T7Y4E,659
 biopipen/reports/scrna/SeuratPreparing.svelte,sha256=CCQ3qodjfoQeMZhSZp7irvSqBt10m7urzcWEOHsZa8M,365
 biopipen/reports/scrna/TopExpressingGenes.svelte,sha256=h2Or2cKCzuY2cypq1PKrqiHFVMZmFUgFOXoQumvfGSk,471
 biopipen/reports/scrna_metabolic_landscape/MetabolicFeatures.svelte,sha256=4M8V6jDv1_ZviTsRdoYNQixoTnKiD06oPHzVabUTq6k,978
@@ -138,7 +139,10 @@ biopipen/scripts/rnaseq/Simulation-ESCO.R,sha256=68cEHDdJclX8P8Q7ey9yBOfK09M_kxl
 biopipen/scripts/rnaseq/Simulation-RUVcorr.R,sha256=6C6Ke5RLF0fC2V9WQPoFEdqoDabCnhslZBIyB6zhIxc,1155
 biopipen/scripts/rnaseq/Simulation.R,sha256=PK9tZS88AcBPStcFalZlMU0KE0gSqFSQvhUoQ-8eg90,871
 biopipen/scripts/rnaseq/UnitConversion.R,sha256=pr1uZuzjapZKwIupE-xtE2HW-OwdVCGBVaHQ_3ipIDw,11377
-biopipen/scripts/scrna/AnnData2Seurat.R,sha256=8VEXHnknoEPWRQ0cj80so6O5FAd1_YlwQ0CP0TwmZVo,2555
+biopipen/scripts/scrna/AnnData2Seurat.R,sha256=PohvXtUziu0_A_HL5HaOYG1wSVdLgMaV8yhzEKDQT84,2814
+biopipen/scripts/scrna/CCPlotR-patch.R,sha256=KpB8fwacBaWaUNjIidcLUkMShLjS4Gq9UY8LUgIITB0,8369
+biopipen/scripts/scrna/CellCellCommunication.py,sha256=1txufcxFsniB9UET_L2RA2qaQZfjidjUbUsese_hoGE,2899
+biopipen/scripts/scrna/CellCellCommunicationPlots.R,sha256=KE-yy7TkMjPz4opfQgosr6NfMQQYeg4mgzfDod-RJrQ,6104
 biopipen/scripts/scrna/CellTypeAnnotation-celltypist.R,sha256=2lkIwtSw2A21RCUBdo3BVUdAwsLTIq65NDA-ZR1YWE0,9644
 biopipen/scripts/scrna/CellTypeAnnotation-common.R,sha256=s2exGdJ3WM8mDDeBXGUxuKO1DAGVgRiirHaJD0lPgDw,336
 biopipen/scripts/scrna/CellTypeAnnotation-direct.R,sha256=hmuk3G4IbF66Sv4RmrcrKn8b4J7ZEBD5PzJHa5EKYB8,2125
@@ -159,7 +163,7 @@ biopipen/scripts/scrna/RadarPlots.R,sha256=4zs0hAm7yq1Ls62f_29koPLqAKCeKbYiztNM-
 biopipen/scripts/scrna/SCImpute.R,sha256=dSJOHhmJ3x_72LBRXT72dbCti5oiB85CJ-OjWtqONbk,2958
 biopipen/scripts/scrna/ScFGSEA.R,sha256=MFoJ3i3LFBfsPCxwLPnTh141ZJyrzwnrTuCZIFwvYjU,6318
 biopipen/scripts/scrna/ScSimulation.R,sha256=b2LtL68ucxLoI57tSEDD0hOSbVHUA_x88Y96eK07N-s,1712
-biopipen/scripts/scrna/Seurat2AnnData.R,sha256=G7bcHGffdNlz6Uuy98tQdlahXiPkTDokflp1yTUgcSQ,1578
+biopipen/scripts/scrna/Seurat2AnnData.R,sha256=zXkdHEqv6H8RCw1bEPStvYHVpu5oRnioWxz164MgHxQ,276
 biopipen/scripts/scrna/SeuratClusterStats-clustree.R,sha256=FkbniQMHiZGrFAIuS8nUfPVJKFofSL-ZnpxTqIhTa54,3058
 biopipen/scripts/scrna/SeuratClusterStats-dimplots.R,sha256=NEdlJHNXnJZfF7YkefYVWTPO8Z_KAppRAs9rNvB8TXs,2360
 biopipen/scripts/scrna/SeuratClusterStats-features.R,sha256=DeGo7AkBRq0V3Y3JDaifId6rrr5dwawTzcSAJ3W1lxE,15614
@@ -171,7 +175,7 @@ biopipen/scripts/scrna/SeuratClustering-common.R,sha256=JX4Cn2FC6GOcBqaVyGDD3MM5
 biopipen/scripts/scrna/SeuratClustering.R,sha256=0OKRBQ5rFuupK7c03_sSt2HMwMdMnCYFqTvkRXFKchs,1706
 biopipen/scripts/scrna/SeuratFilter.R,sha256=BrYK0MLdaTtQvInMaQsmOt7oH_hlks0M1zykkJtg2lM,509
 biopipen/scripts/scrna/SeuratLoading.R,sha256=ekWKnHIqtQb3kHVQiVymAHXXqiUxs6KKefjZKjaykmk,900
-biopipen/scripts/scrna/SeuratMap2Ref.R,sha256=B0RZ2k1IUKdhe34SjU-8CxvYGj7ub-z7JJaSE0snCok,12040
+biopipen/scripts/scrna/SeuratMap2Ref.R,sha256=MP6ttYv8bK5kzWNoP3DXEad06snotVOqCqdvBHsOasc,12776
 biopipen/scripts/scrna/SeuratMetadataMutater.R,sha256=PMwG0Xvl_EEVKkicfrIi4arEqpY948PkYLkb59kTAXI,1135
 biopipen/scripts/scrna/SeuratPreparing-common.R,sha256=WuD7lGS17eAUQWSiIdAoV0EIeqS3Tnkkx-7PbP6Q3tc,16279
 biopipen/scripts/scrna/SeuratPreparing-doublet_detection.R,sha256=TNN2lfFjpnnO0rguMsG38JYCP1nFUhcPLJ1LqGj-Sc8,6674
@@ -212,10 +216,10 @@ biopipen/scripts/tcr/Attach2Seurat.R,sha256=0KZaBkuPvqOBXq4ZG3pzIIua5HL-161K5dVX
 biopipen/scripts/tcr/CDR3AAPhyschem.R,sha256=ZxHyn4vJF-1PDbnbozovOipET_QZKqeFrKUibRpLgzw,16825
 biopipen/scripts/tcr/CloneResidency.R,sha256=KAcFB39vTYsk8IEj44s8oSFVhKjpAdJ8hkpKxtdWzRA,21540
 biopipen/scripts/tcr/CloneSizeQQPlot.R,sha256=hds1C80Q_W40Ikp-BrFfvh_aBf_V61lz-1YAvkDESyk,4569
-biopipen/scripts/tcr/GIANA/GIANA.py,sha256=0qLhgCWxT8K-4JvORA03CzBPTT5pd4Di5B_DgrHXbFA,47198
-biopipen/scripts/tcr/GIANA/GIANA4.py,sha256=Z7Q3cUr1Pvmy4CFADN0P7i9g1-HbzWROMqk5HvL_F1Q,45762
+biopipen/scripts/tcr/GIANA/GIANA.py,sha256=jo0d58K57CF4W6mc2Q-hQn9rLl6oLHTsr5JceP8xqN0,54874
+biopipen/scripts/tcr/GIANA/GIANA4.py,sha256=X_iXdgKZ9V5T87jXBdml0H6bwiDjsAqoPKhsdF72lj8,53904
 biopipen/scripts/tcr/GIANA/Imgt_Human_TRBV.fasta,sha256=XUwDPXJxVH5O9Q0gCL6EILKXEwVyiAZXm4VS2vkPcnQ,15371
-biopipen/scripts/tcr/GIANA/query.py,sha256=5NWSEDNrJomMt48tzLGpRwJTZB0zQxvMVTilyG8osX8,7298
+biopipen/scripts/tcr/GIANA/query.py,sha256=YJeCk24Hab7XfQxZ5VNesdjr_LVxPH0KH9rwGm1ePf4,7680
 biopipen/scripts/tcr/Immunarch-basic.R,sha256=2UAQtGEDKgf5LmPeGlNkkAfQOCtdYpCiBMAn-JHOV4A,3036
 biopipen/scripts/tcr/Immunarch-clonality.R,sha256=NN6WikD-9BKodgKRAbvM2oMLrhn9FUWjqIb69DpHnAk,3744
 biopipen/scripts/tcr/Immunarch-diversity.R,sha256=6T8-IXIMyy4YXgIA4TXJThy91FqgLVsljqxd3PxLDys,27136
@@ -232,7 +236,7 @@ biopipen/scripts/tcr/ImmunarchLoading.R,sha256=Vw2oIza3mDJzg9kuo-w5jvwdivk4AtDA6
 biopipen/scripts/tcr/ImmunarchSplitIdents.R,sha256=FGCeGV0uSmFU91lKkldUAeV4A2m3hHw5X4GNi8ffGzI,1873
 biopipen/scripts/tcr/SampleDiversity.R,sha256=oipN4-2nQZe8bYjI0lZ0SvZ7T8GZ_FWkpkobi1cwmWE,2664
 biopipen/scripts/tcr/TCRClusterStats.R,sha256=QhXgfKSh27VHO901hDizyGYTXMYFJxW22StchQUq3uE,12906
-biopipen/scripts/tcr/TCRClustering.R,sha256=MRCaFbgMhhW9uimVBrE1Q5VqqxQpnKv5W2ILWmHdkU8,9593
+biopipen/scripts/tcr/TCRClustering.R,sha256=Lt_e32c52gJ9fz1abboGXLIR--Lx2dqW0oFlSOCrCVs,10079
 biopipen/scripts/tcr/TCRDock.py,sha256=jjzxMWp-hs0LDtA1mVbiWDvUieSO7X-F9yeKGy1LSTM,3026
 biopipen/scripts/tcr/TESSA.R,sha256=XFC2P_e_Gm83jG5EjzVIW6KcyG5IieAaK7sXDA3_oZ0,6864
 biopipen/scripts/tcr/TESSA_source/Atchley_factors.csv,sha256=SumqDOqP67P54uM7Cuc5_O_rySTWcGo7eX3psMSPX9s,763
@@ -284,9 +288,9 @@ biopipen/utils/mutate_helpers.R,sha256=Bqy6Oi4rrPEPJw0Jq32bVAwwBfZv7JJL9jFcK5x-c
 biopipen/utils/plot.R,sha256=fmWnCv6EpOU8NvHFvShbdPRqB659vHcDlJAqWIXM8XQ,4415
 biopipen/utils/reference.py,sha256=oi5evicLwHxF0KAIPNZohBeHJLJQNWFJH0cr2y5pgcg,5873
 biopipen/utils/rnaseq.R,sha256=Ro2B2dG-Z2oVaT5tkwp9RHBz4dp_RF-JcizlM5GYXFs,1298
-biopipen/utils/single_cell.R,sha256=pJjYP8bIZpNAtTQ32rOXhZxaM1Y-6D-xUcK3pql9tbk,4316
+biopipen/utils/single_cell.R,sha256=3jmTV9-kHTXyBA6GRlRf_Ig0jx_5jRr6E1THoa1869k,7564
 biopipen/utils/vcf.py,sha256=ajXs0M_QghEctlvUlSRjWQIABVF02wPdYd-0LP4mIsU,9377
-biopipen-0.31.7.dist-info/METADATA,sha256=dmsDb7Q7iTWvkZjChqUgVqhb2CiQP4j8OA0jBzRYgGY,882
-biopipen-0.31.7.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
-biopipen-0.31.7.dist-info/entry_points.txt,sha256=BYqHGBQJxyFDNLYqgH64ycI5PYwnlqwYcCFsMvJgzAU,653
-biopipen-0.31.7.dist-info/RECORD,,
+biopipen-0.32.0.dist-info/METADATA,sha256=1kWhdtPEVLqTrBh-9eK1I-Fb7KVCycmWJTA4_HXQLlw,882
+biopipen-0.32.0.dist-info/WHEEL,sha256=sP946D7jFCHeNz5Iq4fL4Lu-PrWrFsgfLXbbkciIZwg,88
+biopipen-0.32.0.dist-info/entry_points.txt,sha256=BYqHGBQJxyFDNLYqgH64ycI5PYwnlqwYcCFsMvJgzAU,653
+biopipen-0.32.0.dist-info/RECORD,,