Microglia single cell

Questions to address from Lot 10X: How many MG (microglia) nuclei did we find in the total nuclei (% of total)?

10X: How do these MG nuclei compare to the fresh microglia in terms of QC (number of transcripts per cell; number of genes per cell) and gene expression (we could look at the Patir core genes first?)

10X us versus 10X others: How do our QC metrics related to reported QC metrics of others?

Setup

#load libaries
library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

library(Seurat)
library(patchwork)
reticulate::py_install(packages ='umap-learn')

library(plotly)

## Loading required package: ggplot2

## 
## Attaching package: 'plotly'

## The following object is masked from 'package:ggplot2':
## 
##     last_plot

## The following object is masked from 'package:stats':
## 
##     filter

## The following object is masked from 'package:graphics':
## 
##     layout

library(future)

print(paste("Seurat ", packageVersion("Seurat")))

## [1] "Seurat  3.2.3"

#Helper function from Brian
createDT <- function(DF, caption="", scrollY=500){
  data <- DT::datatable(DF, caption=caption,
    extensions =  'Buttons',
    options = list( dom = 'Bfrtip', 
                    buttons = c('copy', 'csv', 'excel', 'pdf', 'print'), 
                    scrollY = scrollY, scrollX=T, scrollCollapse = T, paging = F,  
                      columnDefs = list(list(className = 'dt-center', targets = "_all"))
    )
  ) 
   return(data)
}
# Need to wrap inside tagList when rendering DTs within a for loop
createDT_html <- function(DF, caption="", scrollY=400){
  print( htmltools::tagList( createDT(DF, caption, scrollY)) ) 
}

#setwd("/sc/arion/projects/ad-omics/emily/Microglia/fresh")
#setwd("/sc/arion/projects/ad-omics/emily/Microglia")
setwd("~/Raj_Lab/Microglia")
#MG22 <- readRDS("~/hpc/users/kozike01/test/MG22_combined_filtered.rds")
MG22 <- readRDS("MG22_combined_filtered.rds")
NUC <- readRDS("NUC_filtered.rds")

Pre-processing

Normalizing the Data

Next we use a global-scaling normalization method “LogNormalize” that normalizes the feature expression measurements for each cell by the total expression, multiplies this by a scale factor (10,000 by default), and log-transforms the result.

MG22 <- NormalizeData(MG22, normalization.method = "LogNormalize", scale.factor = 10000)
NUC <- NormalizeData(NUC, normalization.method = "LogNormalize", scale.factor = 10000)

levels(Idents(MG22))

## [1] "MG-22-HIPP" "MG-22-MFG"

levels(Idents(NUC))

## [1] "17-016-MFG"

Feature Selection

Highly variable features are those transcripts that have high cell-cell variation in the dataset. Identifying the highly variable transcripts makes it easier in downstream analysis to extract biologically significant information.

MG22 <- FindVariableFeatures(MG22, selection.method = "vst", nfeatures = 2000)
NUC <- FindVariableFeatures(NUC, selection.method = "vst", nfeatures = 2000)


var.genes.MG22 <- Seurat::VariableFeatures(MG22)
appendedGenes.MG22 <- c(var.genes.MG22, "CD68","CCL6", "LYZ2")


var.genes.NUC <- Seurat::VariableFeatures(NUC)
appendedGenes.NUC <- c(var.genes.NUC , "CD68","CCL6", "LYZ2")

# Identify the 10 most highly variable genes
top10_MG22 <- head(VariableFeatures(MG22), 10)
top10_NUC <- head(VariableFeatures(NUC), 10)


highly_var_feats <- data.frame(Top10Features = c(1:10),
                               MG22 = top10_MG22,
                               NUC = top10_NUC)
highly_var_feats

##    Top10Features    MG22    NUC
## 1              1  S100A9 CXCL14
## 2              2  S100A8  TAGLN
## 3              3   CXCL8   PLP1
## 4              4     LYZ PTPRZ1
## 5              5     LTF   CNR1
## 6              6    CCL4 SLC1A2
## 7              7  CCL4L2   BCAN
## 8              8    SPP1  GPR17
## 9              9    LCN2   MYL9
## 10            10 S100A12  CLDN5

# plot variable features with and without labels
plot1_MG22 <- VariableFeaturePlot(MG22)
plot2_MG22<- LabelPoints(plot = plot1_MG22, points = top10_MG22, repel = TRUE)

## Warning: Using `as.character()` on a quosure is deprecated as of rlang 0.3.0.
## Please use `as_label()` or `as_name()` instead.
## This warning is displayed once per session.

## When using repel, set xnudge and ynudge to 0 for optimal results

plot1_MG22 + plot2_MG22

## Warning: Transformation introduced infinite values in continuous x-axis

## Warning: Transformation introduced infinite values in continuous x-axis

plot1_NUC<- VariableFeaturePlot(NUC)
plot2_NUC <- LabelPoints(plot = plot1_NUC, points = top10_NUC, repel = TRUE)

## When using repel, set xnudge and ynudge to 0 for optimal results

plot1_NUC + plot2_NUC

## Warning: Transformation introduced infinite values in continuous x-axis

## Warning: Transformation introduced infinite values in continuous x-axis

# Dimension Reduction ## Scaling the data

Standard pre-processing for dimensional reduction. The ScaleData function:

• Shifts the expression of each gene, so that the mean expression across cells is 0
• Scales the expression of each gene, so that the variance across cells is 1
• This step gives equal weight in downstream analyses, so that highly-expressed genes do not dominate The results of this are stored in pbmc[[“RNA”]]@scale.data

#use if making heatmaps
all.genes_MG22 <- rownames(MG22) 
MG22 <- ScaleData(MG22, 
                      features = all.genes_MG22,
                      vars.to.regress = c("percent.mt","nCount_RNA"))

## Regressing out percent.mt, nCount_RNA

## Centering and scaling data matrix

all.genes_NUC <- rownames(NUC)
NUC <- ScaleData(NUC,features = all.genes_NUC,
                     vars.to.regress = c("percent.mt","nCount_RNA"))

## Regressing out percent.mt, nCount_RNA
## Centering and scaling data matrix

#MG22 <- saveRDS("MG22_scaled.rds")
#NUC <- saveRDS("NUC_scaled.rds")

PCA

MG22 <- RunPCA(MG22, features = VariableFeatures(object = MG22))

## PC_ 1 
## Positive:  FOSB, RHOB, JUN, EGR1, DNAJB1, HSPA1B, HSP90AA1, PPP1R15A, HSPH1, KLF2 
##     FOS, ATF3, IER5, HSPB1, HSPA1A, GADD45B, NR4A2, BTG1, DNAJB4, HSPD1 
##     DUSP1, JUND, IER2, IER3, JUNB, IFRD1, DDIT4, NR4A1, KLF4, TUBB4B 
## Negative:  APOE, RPL21, RPS12, LIPA, RPS26, HLA-DRB1, APOC1, RPL10, NPC2, GLDN 
##     OLFML3, NAP1L1, S100A11, IFI44L, CTSS, HLA-DRA, CACNB4, FTL, RTTN, RPL41 
##     APOC2, SPP1, DAB2, FOLR2, GPNMB, LPL, MRC1, CTSD, DBI, PTMS 
## PC_ 2 
## Positive:  FCN1, S100A9, S100A8, CD44, S100A6, S100A12, TREM1, S100P, ANXA1, LYZ 
##     CD52, CD48, CDA, VIM, FGR, VCAN, AQP9, G0S2, IFITM1, APOBEC3A 
##     LITAF, CRIP1, IFITM2, FLNA, CXCL8, S100A10, BCL2A1, RIPOR2, TAGLN2, CYP4F3 
## Negative:  HSPB1, APOE, PLK2, RGS16, SPP1, CSKMT, TRNAU1AP, DDIT4, DNAJB4, AL031777.3 
##     SGK1, HSD17B7, LINC02256, HIST1H2BG, AC025164.1, IER5L, AL451085.1, IER3, IDI1, HSPA2 
##     MAFB, AL118516.1, KCNB1, KCNQ1OT1, HLA-DRA, ID2, HLA-DRB1, HIST3H2BB, TM2D2, AL021453.1 
## PC_ 3 
## Positive:  CCL5, NKG7, KLRD1, CST7, GZMB, FYN, IL32, PRF1, GZMA, CD247 
##     CTSW, SAMD3, LINC01871, SPON2, IL2RB, IFNG, IFITM1, LINC00861, GNLY, SLAMF7 
##     PYHIN1, MYBL1, RORA, CRIP1, CD8A, DUSP2, TRGC2, SLC38A1, DOK2, SYNE2 
## Negative:  S100A9, S100A8, S100P, S100A12, FCN1, CYP4F3, CDA, CXCL8, FCGR3B, TREM1 
##     G0S2, AC245128.3, AQP9, VNN2, LYZ, NAMPT, ACSL1, PGLYRP1, HCAR3, CMTM2 
##     BCL2A1, BASP1, PADI4, PROK2, SOD2, C5AR1, CXCR1, SAT1, FPR1, CD55 
## PC_ 4 
## Positive:  S100A8, S100A12, S100A9, CSKMT, CYP4F3, S100P, PGLYRP1, FCN1, LTF, PLK2 
##     CDA, AC016745.2, LCN2, CD24, CRISP3, MMP8, PADI4, AC025164.1, MIR222HG, ARG1 
##     AC024257.3, RETN, HP, AC243960.1, ARHGAP25, CAMP, PDK4, HSD17B7, CACNB4, AC016596.1 
## Negative:  HLA-DRA, FTH1, HLA-DRB1, SRGN, RPL41, FTL, F13A1, HLA-DPB1, RGS1, RPS12 
##     CD83, RPL21, RPL10, SELENOP, GPR183, MARCO, CEBPD, RPS26, CD163, DUSP1 
##     TGFBI, CXCR4, NR4A2, RNASE1, NR4A1, FCGR2B, INSIG1, FOS, SIGLEC1, JUNB 
## PC_ 5 
## Positive:  LYVE1, MARCO, RNASE1, MRC1, SELENOP, F13A1, CD163, FOLR2, LILRB5, COLEC12 
##     SIGLEC1, CD14, TGFBI, LGMN, CCDC152, ADGRG6, DAB2, FCGR2B, CH25H, CCL2 
##     LINC00996, CCDC170, CD28, PLTP, UBAC2, NRP1, VMP1, LDLR, AC103591.3, RND3 
## Negative:  APOC1, LGALS1, NUPR1, NAP1L1, S100A11, FTL, SPP1, GPNMB, CD9, APOE 
##     S100A6, GLDN, ALCAM, VIM, TUBB4B, RPS26, SH3BGRL3, RPS12, S100A4, RPL10 
##     RPL41, LSP1, RPL21, HSP90AA1, FABP5, DBI, CXCR4, FTH1, LPL, S100A1

NUC <- RunPCA(NUC, features = VariableFeatures(object = NUC))

## PC_ 1 
## Positive:  MEG3, MIAT, SYT1, MEF2C, CELF2, GRIN2B, ATP2B1, RIMS2, CAMK2A, FBXW7 
##     PRKCB, SNAP25, KCNIP4-IT1, NTRK3, UCHL1, CHN1, R3HDM1, CNKSR2, SLITRK5, CAMK1D 
##     NRGN, MKL2, OLFM1, PCLO, SLC8A1, BCL11A, ST8SIA3, STXBP5L, KCTD16, SCN1A 
## Negative:  PLP1, TMEM144, CNP, MBP, PTGDS, TF, MOBP, AL359091.1, ST18, QKI 
##     ERMN, MAG, MOG, MYRF, GPRC5B, CLDND1, CLDN11, SPP1, HHIP, SELENOP 
##     ENPP2, LAMP2, SCD, SLCO1A2, PLEKHH1, FAM107B, GLUL, GPR37, CERCAM, AMER2 
## PC_ 2 
## Positive:  MBP, TF, PLP1, MOBP, TMEM144, SLC44A1, CLDND1, ST18, CNP, MAG 
##     ENPP2, PLEKHH1, AL359091.1, MOG, HHIP, MYRF, MEG3, EDIL3, SCD, SLCO1A2 
##     ERMN, CERCAM, CLDN11, FAM107B, GPR37, CNTN2, SPP1, CALM1, PIP4K2A, QDPR 
## Negative:  ATP1A2, GJA1, AQP4, SLC1A3, SLC1A2, APOE, PTPRZ1, HIF3A, CST3, SLC14A1 
##     FGFR3, ETNPPL, ADGRV1, SLCO1C1, AGT, GFAP, C1orf61, F3, ATP1B2, RANBP3L 
##     TPD52L1, SLC25A18, MRVI1, WLS, CLU, RFX4, SPARCL1, DTNA, CPE, COL5A3 
## PC_ 3 
## Positive:  SLC6A1, GAD1, DLX6-AS1, ERBB4, GAD2, ZNF385D, MAF, NRIP3, KCNQ1OT1, ANK1 
##     CXCL14, IGF1, GRIK1, VWC2, ARL4C, NXPH1, DNER, NRXN3, DLX1, GRIP2 
##     PROX1, DGKD, CNTNAP2, ZNF536, ADARB2, BTBD11, KIT, TAC1, ARX, MAFB 
## Negative:  CHN1, NRGN, SLC17A7, PSD3, ENC1, AC009041.2, CAMK2A, KCNIP4-IT1, NPTX1, FBXW7 
##     MKL2, LMO4, PTK2B, NECAB1, LMO3, SV2B, TBR1, ARPP21, ARPP19, LGI4 
##     TESPA1, CAMK1D, PRKCB, EPB41L2, HOPX, CNKSR2, ONECUT2, R3HDM1, ANO3, SERPINE2 
## PC_ 4 
## Positive:  PDGFRA, VCAN, SMOC1, GPNMB, TNR, FERMT1, BCAN, CSPG4, MEGF11, SPRY4 
##     PLEKHH2, PCDH15, SULF2, BEST3, MYT1, PHLDA1, COL9A1, LHFPL3, PTPRZ1, SOX4 
##     SNX22, SEMA5A, OLIG1, BCHE, GPR17, MARCKS, APOD, EPN2, HLA-B, FGFR1 
## Negative:  AQP4, GJA1, ADGRV1, SLC14A1, AGT, DTNA, GLUL, ETNPPL, SLCO1C1, NEAT1 
##     APOE, SLC25A18, TPD52L1, GFAP, GAD2, ID4, WLS, CNTNAP3B, PAX6, RYR3 
##     COL5A3, SLC4A4, FGFR3, FAM189A2, LINC00499, ACSBG1, ERBB4, AHCYL1, MGST1, DAAM2 
## PC_ 5 
## Positive:  A2M, COL1A2, VIM, DCN, CFH, FN1, CD74, MGP, ITIH5, CAVIN1 
##     IFITM3, MYL9, HLA-DRB1, TAGLN, HLA-E, SLCO2B1, HLA-DRA, IGFBP2, TFPI, S100A6 
##     EPAS1, C3, MS4A6A, CSF1R, GGT5, PRELP, ID3, ADH1B, MYOF, NCKAP1L 
## Negative:  VCAN, PTPRZ1, BCAN, SMOC1, TNR, FERMT1, C1orf61, PDGFRA, PCDH15, OLIG1 
##     NKAIN4, PLEKHH2, EPN2, SPRY4, MYT1, PHLDA1, KANK1, SLC1A2, CSPG4, LHFPL3 
##     MEGF11, GPR37L1, BEST3, BCHE, HIF3A, APCDD1, LUZP2, SERPINE2, COL9A1, SNX22

DimPlot(MG22, reduction = "pca")

DimPlot(NUC, reduction = "pca")

DimHeatmap(MG22, dims = 1:15, cells = 500, balanced = TRUE)

DimHeatmap(NUC, dims = 1:15, cells = 500, balanced = TRUE)

MG22 <- JackStraw(MG22, num.replicate = 100, dims=20)
MG22 <- ScoreJackStraw(MG22, dims = 1:20)

NUC <- JackStraw(NUC, num.replicate = 100, dims = 20)
NUC <- ScoreJackStraw(NUC, dims = 1:20)

JackStrawPlot(MG22, dims = 1:20)

## Warning: Removed 28772 rows containing missing values (geom_point).

JackStrawPlot(NUC, dims = 1:20)

## Warning: Removed 28156 rows containing missing values (geom_point).

ElbowPlot(MG22)

ElbowPlot(NUC)

Based on Elbow plots, we lose most variation after PC 11 or 12.

MG22HIPP use PCs 1 through 13 MG22 MFG use PCs 1 through 13 17_016 MFG use PCs 1 through 18

Clustering

MG_22

Cluster

## Computing nearest neighbor graph

## Computing SNN

## Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
## 
## Number of nodes: 5990
## Number of edges: 198131
## 
## Running Louvain algorithm...
## Maximum modularity in 10 random starts: 0.9220
## Number of communities: 5
## Elapsed time: 0 seconds

## Warning: The default method for RunUMAP has changed from calling Python UMAP via reticulate to the R-native UWOT using the cosine metric
## To use Python UMAP via reticulate, set umap.method to 'umap-learn' and metric to 'correlation'
## This message will be shown once per session

## 23:31:58 UMAP embedding parameters a = 0.9922 b = 1.112

## 23:31:58 Read 5990 rows and found 13 numeric columns

## 23:31:58 Using Annoy for neighbor search, n_neighbors = 30

## 23:31:58 Building Annoy index with metric = cosine, n_trees = 50

## 0%   10   20   30   40   50   60   70   80   90   100%

## [----|----|----|----|----|----|----|----|----|----|

## **************************************************|
## 23:31:58 Writing NN index file to temp file /var/folders/b8/csskzk797yqdzbrxw6stvzhh0000gn/T//RtmprdwnBz/filefc9539a2b0dc
## 23:31:58 Searching Annoy index using 1 thread, search_k = 3000
## 23:31:59 Annoy recall = 100%
## 23:32:00 Commencing smooth kNN distance calibration using 1 thread
## 23:32:00 Initializing from normalized Laplacian + noise
## 23:32:00 Commencing optimization for 500 epochs, with 252582 positive edges
## 23:32:07 Optimization finished

MG22 <- RunTSNE(object=MG22,  reduction.use = "pca", 
               dims.use = 1:13, 
               do.fast = T)
DimPlot(MG22 , reduction = "tsne")

DimPlot(MG22 , reduction = "tsne", group.by = "orig.ident")

Save RDS file

saveRDS(MG22, file="MG22_final.RDS")

Differential Expression with FindAllMarkers

https://satijalab.org/seurat/archive/v3.2/de_vignette.html This vignette describes the different test that can be used and the ways to use differential expression in seurat

##MG22H <-readRDS("output/MG22_seurat_final.rds")
MG22.markers <- FindAllMarkers(MG22, only.pos = TRUE, min.pct = 0.25, logfc.threshold = 0.25)

## Calculating cluster 0

## Calculating cluster 1

## Calculating cluster 2

## Calculating cluster 3

## Calculating cluster 4

top10markers <- MG22.markers %>% group_by(cluster) %>% top_n(n = 10, wt = avg_logFC)
top10markers_table <-data.frame("Cluster0" = top10markers[1:10,7], 
                                "Cluster1" = top10markers[11:20,7],
                                "Cluster2" = top10markers[21:30,7], 
                                "Cluster3" = top10markers[31:40,7],
                                "Cluster4" = top10markers[41:50,7], 
                                "Cluster5" = top10markers[51:60,7],
                                "Cluster6" = top10markers[61:70,7], 
                                "Cluster7" = top10markers[71:80,7],
                                "Cluster8" = top10markers[81:90,7])

## Warning: The `i` argument of ``[.tbl_df`()` must lie in [0, rows] if positive, as of tibble 3.0.0.
## Use `NA_integer_` as row index to obtain a row full of `NA` values.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.

top10markers_table

##         gene gene.1  gene.2  gene.3 gene.4 gene.5 gene.6 gene.7 gene.8
## 1     CX3CR1  HSPB1 S100A12  RNASE1   CCL5   <NA>   <NA>   <NA>   <NA>
## 2     P2RY12 DNAJB1   S100P   MARCO   NKG7   <NA>   <NA>   <NA>   <NA>
## 3        A2M    JUN    FCN1   F13A1  KLRD1   <NA>   <NA>   <NA>   <NA>
## 4      SOCS6   RGS1  S100A8   CD163   IL32   <NA>   <NA>   <NA>   <NA>
## 5     IFNGR1   EGR1  S100A9 SELENOP   IFNG   <NA>   <NA>   <NA>   <NA>
## 6    SLCO2B1  NR4A2    G0S2 CCDC152 IFITM1   <NA>   <NA>   <NA>   <NA>
## 7     P2RY13   ATF3   NAMPT   LYVE1   GZMB   <NA>   <NA>   <NA>   <NA>
## 8      MKNK1   FOSB   CXCL8    DAB2  DUSP2   <NA>   <NA>   <NA>   <NA>
## 9  LINC01268 TUBB4B  IFITM2    MRC1  ARL4C   <NA>   <NA>   <NA>   <NA>
## 10    CACNB4 DNAJB4     LYZ    CD14   GNLY   <NA>   <NA>   <NA>   <NA>

top10.MG22 <- MG22.markers %>% group_by(cluster) %>% top_n(n = 10, wt = avg_logFC)
DoHeatmap(MG22, features = top10.MG22$gene) + NoLegend()

DotPlot(MG22, features = c("CX3CR1", "P2RY12", "SOCS6", "TUBB48", "RGS1", "HSPB1", "CXCR4", "JUN", "FOSB", "CCL3", "PPP1R15A", "S100P", "FCGR3B","S100A12", "S100A9", "IFITM2", "RNASE1", "MARCO", "F13A1", "CD163", "CCL5", "NKG7", "KLRD1", "IL32", "IFITM1", "DUSP2", "VCAN", "EREG", "FCN1", "THBS1", "CD52", "CD44")) + RotatedAxis()

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: TUBB48

top5.MG22 <- MG22.markers %>% group_by(cluster) %>% top_n(n = 5, wt = avg_logFC)
DotPlot(MG22, features = top5.MG22$gene) + RotatedAxis()

Gene enrichment with gprofiler

Although a useful tool, in this case gprofiler was not very useful because there is not as much literature defining the different subtypes of microglia. Most hits come up as realting to immune function.

##   
## #### Cluster 0   
##  [1] "Biomarker: CX3CR1"    "Biomarker: P2RY12"    "Biomarker: A2M"      
##  [4] "Biomarker: SOCS6"     "Biomarker: IFNGR1"    "Biomarker: SLCO2B1"  
##  [7] "Biomarker: P2RY13"    "Biomarker: MKNK1"     "Biomarker: LINC01268"
## [10] "Biomarker: CACNB4"

## Warning: `group_by_()` is deprecated as of dplyr 0.7.0.
## Please use `group_by()` instead.
## See vignette('programming') for more help
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.

## [1] "No significant GO terms at Bonferonni-corrected p-value < 0.05."
##   
##   
## #### Cluster 1   
##  [1] "Biomarker: HSPB1"  "Biomarker: DNAJB1" "Biomarker: JUN"   
##  [4] "Biomarker: RGS1"   "Biomarker: EGR1"   "Biomarker: NR4A2" 
##  [7] "Biomarker: ATF3"   "Biomarker: FOSB"   "Biomarker: TUBB4B"
## [10] "Biomarker: DNAJB4"
## [1] "No significant GO terms at Bonferonni-corrected p-value < 0.05."
##   
##   
## #### Cluster 2   
##  [1] "Biomarker: S100A12" "Biomarker: S100P"   "Biomarker: FCN1"   
##  [4] "Biomarker: S100A8"  "Biomarker: S100A9"  "Biomarker: G0S2"   
##  [7] "Biomarker: NAMPT"   "Biomarker: CXCL8"   "Biomarker: IFITM2" 
## [10] "Biomarker: LYZ"    
## <div id="htmlwidget-e9d49c972d73e761db5f" style="width:960px;height:500px;" class="datatables html-widget"></div>
## <script type="application/json" data-for="htmlwidget-e9d49c972d73e761db5f">{"x":{"filter":"none","extensions":["Buttons"],"caption":"<caption><\/caption>","data":[["1","4","5","6","7","657","658","659","730","731"],["query_1","query_1","query_1","query_1","query_1","query_1","query_1","query_1","query_1","query_1"],[true,true,true,true,true,true,true,true,true,true],[6.8615226007384e-06,1.95120261377797e-07,1.7357725883558e-06,1.91232084666745e-06,9.56710736103201e-06,4.54505699958536e-07,4.89604366829837e-07,5.07979664777956e-07,1.3512585850603e-09,7.40930004293025e-06],[2,673,500,507,144,320,324,326,11,83],[2,10,10,10,10,10,10,10,9,9],[2,8,7,7,5,6,6,6,4,4],[1,0.8,0.7,0.7,0.5,0.6,0.6,0.6,0.444444444444444,0.444444444444444],[1,0.011887072808321,0.014,0.0138067061143984,0.0347222222222222,0.01875,0.0185185185185185,0.0184049079754601,0.363636363636364,0.0481927710843374],["CORUM:6826","GO:0002274","GO:0042119","GO:0036230","GO:0019730","GO:0034774","GO:0060205","GO:0031983","GO:0050786","GO:0048306"],["CORUM","GO:BP","GO:BP","GO:BP","GO:BP","GO:CC","GO:CC","GO:CC","GO:MF","GO:MF"],["Calprotectin heterotetramer","myeloid leukocyte activation","neutrophil activation","granulocyte activation","antimicrobial humoral response","secretory granule lumen","cytoplasmic vesicle lumen","vesicle lumen","RAGE receptor binding","calcium-dependent protein binding"],[3624,18092,18092,18092,18092,18963,18963,18963,18694,18694],[2399,942,11592,10993,6972,1811,2613,1463,9105,8432],["CORUM:0000000","GO:0045321","GO:0036230","GO:0002274",["GO:0006959","GO:0051707"],["GO:0030141","GO:0060205"],["GO:0031410","GO:0031983"],["GO:0031982","GO:0043233"],"GO:0005102","GO:0005515"],[0.0325647862631044,0.000926040760499024,0.00823797670433662,0.00907587473828372,0.0454054915354579,0.00215708405200321,0.00232366232497441,0.00241087148903618,6.41307324469619e-06,0.035164538003747]],"container":"<table class=\"display\">\n  <thead>\n    <tr>\n      <th> <\/th>\n      <th>query<\/th>\n      <th>significant<\/th>\n      <th>p_value<\/th>\n      <th>term_size<\/th>\n      <th>query_size<\/th>\n      <th>intersection_size<\/th>\n      <th>precision<\/th>\n      <th>recall<\/th>\n      <th>term_id<\/th>\n      <th>source<\/th>\n      <th>term_name<\/th>\n      <th>effective_domain_size<\/th>\n      <th>source_order<\/th>\n      <th>parents<\/th>\n      <th>q_value<\/th>\n    <\/tr>\n  <\/thead>\n<\/table>","options":{"dom":"Bfrtip","buttons":["copy","csv","excel","pdf","print"],"scrollY":400,"scrollX":true,"scrollCollapse":true,"paging":false,"columnDefs":[{"className":"dt-center","targets":"_all"},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false}},"evals":[],"jsHooks":[]}</script>
##   
##   
## #### Cluster 3   
##  [1] "Biomarker: RNASE1"  "Biomarker: MARCO"   "Biomarker: F13A1"  
##  [4] "Biomarker: CD163"   "Biomarker: SELENOP" "Biomarker: CCDC152"
##  [7] "Biomarker: LYVE1"   "Biomarker: DAB2"    "Biomarker: MRC1"   
## [10] "Biomarker: CD14"   
## <div id="htmlwidget-be9297ecea7408e1286a" style="width:960px;height:500px;" class="datatables html-widget"></div>
## <script type="application/json" data-for="htmlwidget-be9297ecea7408e1286a">{"x":{"filter":"none","extensions":["Buttons"],"caption":"<caption><\/caption>","data":[["602"],["query_1"],[true],[9.96853215988115e-06],[77],[10],[4],[0.4],[0.051948051948052],["GO:0038024"],["GO:MF"],["cargo receptor activity"],[18694],[6534],["GO:0003674"],[0.0467623843620025]],"container":"<table class=\"display\">\n  <thead>\n    <tr>\n      <th> <\/th>\n      <th>query<\/th>\n      <th>significant<\/th>\n      <th>p_value<\/th>\n      <th>term_size<\/th>\n      <th>query_size<\/th>\n      <th>intersection_size<\/th>\n      <th>precision<\/th>\n      <th>recall<\/th>\n      <th>term_id<\/th>\n      <th>source<\/th>\n      <th>term_name<\/th>\n      <th>effective_domain_size<\/th>\n      <th>source_order<\/th>\n      <th>parents<\/th>\n      <th>q_value<\/th>\n    <\/tr>\n  <\/thead>\n<\/table>","options":{"dom":"Bfrtip","buttons":["copy","csv","excel","pdf","print"],"scrollY":400,"scrollX":true,"scrollCollapse":true,"paging":false,"columnDefs":[{"className":"dt-center","targets":"_all"},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false}},"evals":[],"jsHooks":[]}</script>
##   
##   
## #### Cluster 4   
##  [1] "Biomarker: CCL5"   "Biomarker: NKG7"   "Biomarker: KLRD1" 
##  [4] "Biomarker: IL32"   "Biomarker: IFNG"   "Biomarker: IFITM1"
##  [7] "Biomarker: GZMB"   "Biomarker: DUSP2"  "Biomarker: ARL4C" 
## [10] "Biomarker: GNLY"  
## [1] "No significant GO terms at Bonferonni-corrected p-value < 0.05."
## 

FindMarkers(MG22, ident.1 = "3", ident.2 = "2")

##                    p_val  avg_logFC pct.1 pct.2    p_val_adj
## SELENOP     2.436055e-63  3.1325870 0.992 0.078 4.626312e-59
## S100A9      4.475581e-54 -4.8267407 0.167 0.985 8.499576e-50
## S100A8      2.295506e-52 -4.7855533 0.129 0.971 4.359396e-48
## F13A1       3.808922e-46  1.8518107 0.765 0.015 7.233525e-42
## MARCO       3.658808e-45  2.5896631 0.788 0.054 6.948443e-41
## RNASE1      5.131373e-43  2.8497536 0.727 0.024 9.744991e-39
## DAB2        5.819818e-43  1.9507596 0.848 0.137 1.105242e-38
## LYVE1       5.599124e-42  2.0936778 0.856 0.156 1.063330e-37
## CD163       2.401129e-41  1.6538694 0.758 0.059 4.559984e-37
## MS4A4A      8.363732e-40  1.6002759 0.826 0.132 1.588356e-35
## C1QA        2.057284e-37  1.4886350 0.970 0.517 3.906987e-33
## CD74        2.134694e-33  0.9540210 1.000 0.888 4.053997e-29
## MRC1        2.567385e-33  1.5961971 0.727 0.102 4.875721e-29
## FOLR2       4.367602e-32  1.3244523 0.833 0.215 8.294513e-28
## MXD1        9.623623e-32 -2.0471456 0.159 0.815 1.827622e-27
## C1QB        5.615771e-31  1.2289599 0.962 0.527 1.066491e-26
## C1QC        8.782318e-31  1.2769799 0.879 0.307 1.667850e-26
## S100A6      2.148389e-30 -2.4338983 0.030 0.678 4.080005e-26
## LSP1        9.409605e-30 -1.9161987 0.045 0.678 1.786978e-25
## CXCL8       1.137527e-29 -3.8600250 0.129 0.737 2.160277e-25
## IFITM2      1.971358e-29 -2.2704895 0.568 0.893 3.743806e-25
## CCDC152     2.036838e-29  1.5147658 0.515 0.005 3.868158e-25
## HLA-DRB1    9.213406e-29  1.0781218 0.962 0.698 1.749718e-24
## NAMPT       9.653600e-29 -1.7282862 0.636 0.893 1.833315e-24
## CST3        7.396157e-28  1.1053000 0.970 0.673 1.404604e-23
## MNDA        1.150429e-27 -1.5924674 0.432 0.888 2.184780e-23
## NPC2        1.408782e-27  0.9976115 0.955 0.532 2.675418e-23
## COLEC12     1.515149e-27  1.1735727 0.477 0.000 2.877419e-23
## CD14        1.633384e-27  1.3577508 0.902 0.459 3.101960e-23
## G0S2        1.649246e-26 -2.7013703 0.015 0.605 3.132084e-22
## HLA-DPA1    3.739405e-26  1.0973112 0.939 0.668 7.101504e-22
## FCGR2B      1.488563e-25  1.2480639 0.523 0.039 2.826930e-21
## S100A11     1.552120e-25 -1.2623873 0.652 0.883 2.947632e-21
## SLC4A7      2.009742e-25  1.2458126 0.879 0.415 3.816700e-21
## GYPC        9.136898e-25  1.1051556 0.735 0.185 1.735188e-20
## HLA-DPB1    2.715794e-24  1.0939342 0.909 0.488 5.157564e-20
## RPL18A      3.296044e-24  0.9280475 0.992 0.620 6.259517e-20
## TGFBI       3.619876e-24  1.1523126 0.644 0.146 6.874507e-20
## S100P       4.456678e-24 -2.5262095 0.008 0.556 8.463677e-20
## CD55        7.912614e-24 -1.6787542 0.159 0.683 1.502685e-19
## HLA-DRA     1.015304e-23  0.9155957 0.985 0.810 1.928164e-19
## RPS18       1.296535e-23  0.9304441 0.985 0.727 2.462249e-19
## RPSA        1.734217e-23  0.9639945 0.879 0.405 3.293452e-19
## LST1        2.266878e-23 -1.2349179 0.545 0.859 4.305027e-19
## PSAP        2.370755e-23  0.8174256 0.970 0.785 4.502301e-19
## BCL2A1      2.564631e-23 -1.6622115 0.061 0.605 4.870492e-19
## BASP1       2.770865e-23 -1.6854895 0.311 0.766 5.262150e-19
## CSF3R       3.201941e-23 -1.4645735 0.470 0.849 6.080806e-19
## VMO1        5.057434e-23  1.0774784 0.538 0.059 9.604573e-19
## MSR1        6.041045e-23  1.0630258 0.674 0.171 1.147255e-18
## S100A12     6.500603e-23 -3.1698585 0.000 0.527 1.234530e-18
## SLC2A3      8.025342e-23 -1.4068158 0.379 0.854 1.524093e-18
## HLA-DMA     1.044833e-22  0.9712040 0.871 0.434 1.984241e-18
## RPL5        2.057855e-22  0.9151165 0.962 0.654 3.908073e-18
## EGFL7       2.675608e-22  0.9333116 0.417 0.010 5.081247e-18
## TCF4        2.941623e-22  0.9921608 0.841 0.317 5.586436e-18
## EPS15       3.501651e-22  0.9654580 0.705 0.185 6.649986e-18
## RPL13A      8.248084e-22  0.8232986 0.992 0.790 1.566394e-17
## H3F3A       1.672516e-21 -1.0109777 0.818 0.922 3.176275e-17
## HLA-DMB     2.893068e-21  0.9398401 0.886 0.449 5.494225e-17
## RPL13       3.301116e-21  0.8080337 0.992 0.746 6.269149e-17
## LILRB5      1.139845e-20  1.1824247 0.379 0.005 2.164680e-16
## VIM         1.205615e-20 -2.0147278 0.030 0.517 2.289584e-16
## AXL         1.879633e-20  0.9456403 0.538 0.083 3.569611e-16
## CTSB        2.278984e-20  0.8734538 0.962 0.639 4.328019e-16
## H3F3B       3.682056e-20 -1.5903470 0.803 0.922 6.992593e-16
## FCGR3B      6.052596e-20 -1.9429838 0.015 0.488 1.149448e-15
## FCN1        6.717325e-20 -2.1180136 0.015 0.488 1.275687e-15
## STAB1       9.364421e-20  1.0886445 0.742 0.278 1.778397e-15
## APOE        1.021470e-19  1.1115662 0.909 0.512 1.939874e-15
## NCF1        1.070322e-19 -1.4510543 0.167 0.654 2.032649e-15
## RPL38       1.846908e-19  0.7212148 0.977 0.790 3.507463e-15
## SIGLEC1     2.319465e-19  0.8820835 0.356 0.005 4.404896e-15
## RPS12       2.618721e-19  0.8191358 0.985 0.741 4.973213e-15
## RPS21       3.093893e-19  0.8542134 0.932 0.541 5.875611e-15
## ADGRG6      4.554291e-19  0.8720265 0.371 0.010 8.649054e-15
## NFKBIA      5.294466e-19 -1.6078869 0.288 0.698 1.005472e-14
## MAF         6.972309e-19  0.8616694 0.985 0.615 1.324111e-14
## ARHGAP18    7.044277e-19  0.7998075 0.439 0.044 1.337779e-14
## MS4A6A      7.087263e-19  0.9044689 0.894 0.532 1.345942e-14
## SH3BGRL3    1.308938e-18 -1.2683950 0.402 0.741 2.485805e-14
## LIPA        1.319436e-18  0.9386520 0.765 0.322 2.505741e-14
## RPS6        2.469430e-18  0.7921287 0.955 0.732 4.689694e-14
## RPS5        2.707801e-18  0.8341986 0.856 0.454 5.142385e-14
## FPR1        3.202072e-18 -1.3803572 0.197 0.644 6.081055e-14
## SOD2        3.654761e-18 -1.3026664 0.523 0.829 6.940757e-14
## TREM1       4.104852e-18 -1.4395741 0.015 0.454 7.795524e-14
## MTSS1       4.428110e-18  0.9778846 0.614 0.185 8.409424e-14
## TAGLN2      4.814104e-18 -1.3520475 0.053 0.498 9.142465e-14
## S100A4      5.474265e-18 -1.7860940 0.144 0.571 1.039618e-13
## SELL        7.859894e-18 -1.5033949 0.038 0.483 1.492672e-13
## LCP1        9.048648e-18 -1.1456919 0.523 0.776 1.718429e-13
## SLC40A1     1.206941e-17  0.9452947 0.606 0.156 2.292101e-13
## SLC25A37    1.639258e-17 -1.2747253 0.447 0.761 3.113115e-13
## RPS29       1.709164e-17  0.7489132 0.977 0.751 3.245873e-13
## EEF1A1      2.054213e-17  0.5874464 1.000 0.907 3.901155e-13
## IFI44L      2.085531e-17  1.0667171 0.424 0.049 3.960632e-13
## ITSN1       2.272700e-17  0.9026599 0.386 0.034 4.316085e-13
## NRP1        2.728880e-17  1.0045712 0.455 0.068 5.182417e-13
## LGMN        3.474344e-17  1.2876928 0.538 0.132 6.598126e-13
## LITAF       3.725863e-17 -1.2569665 0.462 0.751 7.075786e-13
## RPS27A      3.801549e-17  0.7710297 0.917 0.678 7.219522e-13
## PLXND1      1.484457e-16  0.8440941 0.561 0.141 2.819133e-12
## RPL35       1.567346e-16  0.7742947 0.909 0.595 2.976546e-12
## RPL12       1.710705e-16  0.6745190 0.977 0.756 3.248800e-12
## BST2        1.793855e-16  0.8798402 0.659 0.259 3.406710e-12
## COTL1       1.977334e-16 -1.1374659 0.364 0.702 3.755155e-12
## VSIG4       2.069519e-16  0.8210353 0.576 0.141 3.930224e-12
## KCNAB1      2.075056e-16  0.6411157 0.311 0.005 3.940739e-12
## RPS8        2.643335e-16  0.6766031 0.977 0.756 5.019957e-12
## JUND        2.822686e-16 -1.2524915 0.606 0.844 5.360563e-12
## RPL37A      4.370078e-16  0.6311398 1.000 0.800 8.299214e-12
## RPL10       4.444343e-16  0.6492625 0.970 0.756 8.440251e-12
## ATP1B1      6.902988e-16  0.8581051 0.583 0.161 1.310946e-11
## RPS19       7.466256e-16  0.6273847 0.970 0.688 1.417917e-11
## RPS2        1.170781e-15  0.7071682 0.924 0.673 2.223429e-11
## RPL19       1.688265e-15  0.6816856 0.962 0.663 3.206183e-11
## RPL21       1.722411e-15  0.7081275 0.970 0.737 3.271031e-11
## RPS23       2.008190e-15  0.7216233 0.955 0.688 3.813754e-11
## MAMDC2      2.168901e-15  0.6746721 0.341 0.024 4.118961e-11
## RPS3A       2.385686e-15  0.7032552 0.970 0.688 4.530656e-11
## CD81        2.412429e-15  0.7326594 0.811 0.385 4.581445e-11
## RPL17       4.305488e-15  0.7068942 0.909 0.571 8.176551e-11
## CSRNP1      4.362671e-15 -1.2057089 0.083 0.483 8.285148e-11
## RPL41       4.862205e-15  0.5958631 0.985 0.893 9.233814e-11
## ARHGAP5     5.671206e-15  0.7695692 0.576 0.166 1.077019e-10
## SNX6        6.027777e-15  0.7575377 0.614 0.205 1.144735e-10
## RPS13       6.680412e-15  0.7474624 0.947 0.620 1.268677e-10
## RPS16       1.091454e-14  0.7146468 0.917 0.620 2.072781e-10
## RPL22       1.176525e-14  0.7334359 0.902 0.566 2.234338e-10
## RPL23       1.318711e-14  0.6232804 0.985 0.756 2.504364e-10
## RPL26       1.335328e-14  0.6347440 0.977 0.712 2.535921e-10
## SESN1       1.337608e-14  0.8510012 0.576 0.171 2.540252e-10
## TTN         1.776478e-14  0.9471525 0.326 0.029 3.373710e-10
## CYP4F3      1.891953e-14 -1.5275656 0.000 0.356 3.593008e-10
## FTH1        1.942095e-14 -0.9913457 1.000 0.956 3.688233e-10
## SEPT11      2.825676e-14  0.7808740 0.326 0.029 5.366242e-10
## RPL36       3.083052e-14  0.6806943 0.924 0.712 5.855025e-10
## RPS24       3.169748e-14  0.6468122 0.985 0.795 6.019669e-10
## RPS28       3.753308e-14  0.5343958 0.992 0.820 7.127906e-10
## WWP1        4.582528e-14  0.5956375 0.364 0.044 8.702680e-10
## FCGR2A      4.773303e-14  0.7199030 0.856 0.615 9.064980e-10
## RPL35A      5.059154e-14  0.6332310 0.939 0.634 9.607839e-10
## OLR1        5.726626e-14 -1.2343198 0.220 0.590 1.087544e-09
## HCAR3       5.753178e-14 -1.6327713 0.008 0.356 1.092586e-09
## FCHO2       6.080061e-14  0.6559626 0.341 0.034 1.154664e-09
## RPL27A      6.727921e-14  0.7112954 0.917 0.610 1.277700e-09
## RPS3        7.975720e-14  0.6769179 0.894 0.585 1.514669e-09
## AQP9        8.449696e-14 -1.1348332 0.000 0.341 1.604682e-09
## RPS27       8.687978e-14  0.6329543 0.985 0.873 1.649934e-09
## RPLP2       9.518604e-14  0.5558461 0.992 0.805 1.807678e-09
## OAZ1        1.254523e-13 -0.7298005 0.697 0.868 2.382464e-09
## RPL10A      1.305148e-13  0.7476230 0.826 0.493 2.478607e-09
## CYBB        1.368300e-13  0.5701493 0.947 0.678 2.598538e-09
## CD68        1.374753e-13  0.6309642 0.788 0.405 2.610794e-09
## MS4A7       1.387846e-13  0.8042071 0.795 0.439 2.635659e-09
## RPL32       1.418813e-13  0.5701025 0.977 0.805 2.694469e-09
## RPL31       1.713257e-13  0.6410560 0.902 0.595 3.253647e-09
## RPS11       2.071397e-13  0.5514221 0.977 0.805 3.933791e-09
## CDA         2.259755e-13 -1.2249322 0.000 0.332 4.291502e-09
## RPL18       2.514972e-13  0.6301228 0.909 0.644 4.776183e-09
## CHCHD10     2.796851e-13  0.8238319 0.538 0.190 5.311500e-09
## RPLP1       2.813567e-13  0.5247345 0.992 0.883 5.343246e-09
## CEBPB       2.849047e-13 -1.0564555 0.682 0.844 5.410626e-09
## ABHD5       2.921770e-13 -1.1888868 0.114 0.498 5.548733e-09
## ZFP36       2.986598e-13 -1.0906726 0.583 0.790 5.671848e-09
## SNX9        3.072185e-13  0.6830658 0.341 0.049 5.834386e-09
## CXCL16      3.550616e-13  0.8367129 0.674 0.317 6.742974e-09
## RPS10       3.606780e-13  0.7431671 0.826 0.498 6.849636e-09
## IVNS1ABP    3.875460e-13 -1.1812352 0.447 0.722 7.359886e-09
## RPS4X       3.958098e-13  0.6404406 0.871 0.551 7.516824e-09
## PEAK1       4.129983e-13  0.6734714 0.470 0.112 7.843250e-09
## CTSZ        4.441556e-13  0.6822274 0.833 0.502 8.434959e-09
## FGR         4.870542e-13 -1.1208971 0.008 0.337 9.249646e-09
## RPL3        5.335411e-13  0.6911925 0.894 0.610 1.013248e-08
## HSPB1       6.561698e-13  0.7845689 0.939 0.649 1.246132e-08
## MYL6        6.565870e-13 -0.9304522 0.667 0.839 1.246924e-08
## RPL36A      6.572944e-13  0.7044836 0.871 0.566 1.248268e-08
## ADAP2       6.739458e-13  0.6047114 0.818 0.400 1.279891e-08
## LRP1        7.072465e-13  0.8073589 0.682 0.327 1.343132e-08
## RPS17       7.544909e-13  0.7202701 0.841 0.556 1.432854e-08
## RPL7A       8.095594e-13  0.6732754 0.864 0.571 1.537434e-08
## RPS26       1.057479e-12  0.5659883 0.970 0.722 2.008258e-08
## TPT1        1.142635e-12  0.4730920 1.000 0.941 2.169977e-08
## SLC11A1     1.186073e-12 -1.0968958 0.220 0.571 2.252471e-08
## RPS14       1.194823e-12  0.5942693 0.985 0.761 2.269088e-08
## RIPOR2      1.392406e-12 -1.0474323 0.030 0.366 2.644318e-08
## RPL14       1.400273e-12  0.6069143 0.955 0.610 2.659258e-08
## EEF1B2      1.403918e-12  0.7318737 0.773 0.420 2.666181e-08
## APOBEC3A    1.564052e-12 -1.2416859 0.000 0.312 2.970290e-08
## A2M         1.742845e-12  0.5484848 0.886 0.459 3.309837e-08
## RPL15       1.746104e-12  0.6012272 0.939 0.693 3.316026e-08
## RPL11       1.797814e-12  0.5815789 0.977 0.761 3.414229e-08
## TRPS1       1.952458e-12  0.8591124 0.455 0.132 3.707913e-08
## RPL27       2.165573e-12  0.6497141 0.856 0.522 4.112640e-08
## RPL34       2.414686e-12  0.5156397 0.985 0.849 4.585730e-08
## SPP1        3.008786e-12 -0.8632846 0.386 0.693 5.713985e-08
## FUCA1       3.150882e-12  0.6616931 0.311 0.039 5.983839e-08
## TUBA4A      3.272230e-12 -1.1235787 0.023 0.341 6.214291e-08
## C20orf194   3.355590e-12  0.7117270 0.379 0.073 6.372600e-08
## SCN1B       3.795007e-12  0.6359639 0.424 0.098 7.207097e-08
## VNN2        4.904012e-12 -1.3273659 0.197 0.527 9.313210e-08
## ACSL1       5.015825e-12 -1.1709444 0.220 0.556 9.525553e-08
## PLTP        5.138540e-12  0.7391133 0.311 0.039 9.758601e-08
## GPCPD1      5.393863e-12 -1.1510006 0.235 0.576 1.024348e-07
## CCNL1       5.575180e-12 -0.9678484 0.477 0.722 1.058783e-07
## AP2A2       5.593322e-12  0.6664412 0.508 0.161 1.062228e-07
## BTG1        5.826645e-12 -0.6877842 0.765 0.907 1.106538e-07
## HLA-DQA1    5.887517e-12  0.8064935 0.758 0.424 1.118098e-07
## KCNQ1OT1    6.220370e-12  1.0593771 0.568 0.229 1.181311e-07
## CD44        6.266881e-12 -1.4033048 0.023 0.332 1.190143e-07
## MGAT4A      6.634537e-12  0.6698013 0.652 0.283 1.259965e-07
## B2M         6.696497e-12 -0.3826355 1.000 0.995 1.271732e-07
## TIMP2       6.885625e-12  0.6541680 0.689 0.302 1.307649e-07
## RPS15       7.138547e-12  0.5730153 0.924 0.737 1.355681e-07
## IGF1        7.614571e-12  0.7626315 0.364 0.068 1.446083e-07
## ANP32A      7.813880e-12 -0.8554134 0.303 0.629 1.483934e-07
## GCA         8.062983e-12 -1.1538062 0.242 0.576 1.531241e-07
## RGL1        8.293526e-12  0.6867043 0.341 0.054 1.575024e-07
## YWHAZ       8.505172e-12 -0.8837127 0.492 0.717 1.615217e-07
## KCTD12      9.896186e-12  0.5358598 0.864 0.561 1.879385e-07
## GPR155      1.025971e-11  0.8763200 0.508 0.176 1.948421e-07
## RPL30       1.133565e-11  0.5467786 0.977 0.785 2.152753e-07
## SLC16A3     1.188030e-11 -1.1237585 0.167 0.502 2.256188e-07
## ITGB5       1.196601e-11  0.5872125 0.356 0.063 2.272464e-07
## PLD3        1.198925e-11  0.7640770 0.508 0.176 2.276878e-07
## TUBA1B      1.436904e-11  0.6624108 0.758 0.385 2.728825e-07
## KLF2        1.659160e-11 -0.9154887 0.477 0.756 3.150910e-07
## CD28        1.667370e-11  0.5398637 0.205 0.000 3.166501e-07
## XAF1        1.707621e-11  0.8591698 0.553 0.205 3.242943e-07
## CDKN2D      2.025882e-11 -1.1107096 0.053 0.371 3.847352e-07
## R3HDM4      2.370174e-11 -0.9735414 0.023 0.327 4.501198e-07
## NEAT1       2.423925e-11 -0.5856231 0.970 0.980 4.603276e-07
## RPL23A      2.630736e-11  0.6575659 0.909 0.668 4.996030e-07
## PNRC1       2.743104e-11 -0.6610182 0.788 0.883 5.209429e-07
## SCARB2      2.876144e-11  0.5481721 0.508 0.161 5.462085e-07
## SCAMP2      3.326885e-11  0.6146148 0.500 0.161 6.318088e-07
## RPS25       3.333403e-11  0.6334020 0.924 0.673 6.330466e-07
## SAMD4A      3.390913e-11  0.6149245 0.258 0.024 6.439682e-07
## RPL37       3.787773e-11  0.5003222 0.985 0.824 7.193360e-07
## ANKH        4.182045e-11  0.6162865 0.402 0.093 7.942121e-07
## TMEM154     4.187358e-11 -1.1366322 0.045 0.346 7.952211e-07
## HLA-DQB1    4.203966e-11  0.6796663 0.833 0.541 7.983752e-07
## CSF1R       4.453380e-11  0.6118193 0.841 0.512 8.457413e-07
## CTSL        4.592798e-11  0.7157795 0.553 0.224 8.722182e-07
## SLC38A6     4.880246e-11  0.6192033 0.280 0.034 9.268075e-07
## CCND1       6.140938e-11  0.6635475 0.485 0.151 1.166226e-06
## RPS20       7.055366e-11  0.5348789 0.932 0.673 1.339885e-06
## MYADM       7.517468e-11 -1.1382650 0.098 0.400 1.427642e-06
## PLEK        7.817834e-11 -1.1036459 0.515 0.702 1.484685e-06
## ITGAX       8.041574e-11 -1.0415529 0.295 0.571 1.527175e-06
## PTGS2       8.461878e-11 -1.4726499 0.083 0.390 1.606995e-06
## SLC46A3     9.527655e-11  0.6249683 0.318 0.054 1.809397e-06
## SCIN        9.619070e-11  0.5745429 0.621 0.244 1.826758e-06
## ADGRG3      1.041577e-10 -0.9328696 0.000 0.268 1.978060e-06
## PROK2       1.041586e-10 -1.2956286 0.000 0.268 1.978076e-06
## CANX        1.042030e-10  0.6323806 0.697 0.376 1.978920e-06
## VASP        1.231128e-10 -0.9694410 0.189 0.512 2.338036e-06
## RPLP0       1.235748e-10  0.6813289 0.765 0.454 2.346809e-06
## PEBP1       1.251007e-10  0.5712867 0.432 0.127 2.375787e-06
## HRH1        1.315835e-10  0.5239320 0.250 0.024 2.498902e-06
## RPL6        1.330729e-10  0.6667723 0.932 0.668 2.527188e-06
## HCST        1.413263e-10  0.6798665 0.788 0.459 2.683928e-06
## CCDC170     1.597576e-10  0.5287898 0.227 0.015 3.033957e-06
## SERPING1    1.692317e-10  0.5544930 0.265 0.029 3.213878e-06
## TREM2       1.762578e-10  0.6036262 0.682 0.288 3.347313e-06
## RPL39       2.304702e-10  0.4787441 0.992 0.883 4.376859e-06
## P2RY14      2.613743e-10  0.5282835 0.197 0.005 4.963760e-06
## CDC42EP3    2.715343e-10 -1.0447156 0.053 0.346 5.156707e-06
## TRIB1       2.958601e-10 -0.9128388 0.288 0.590 5.618679e-06
## AL031056.1  2.964068e-10  0.5311720 0.197 0.005 5.629061e-06
## RB1         2.980962e-10  0.5635778 0.621 0.293 5.661145e-06
## IFITM1      3.152711e-10 -1.2108054 0.015 0.283 5.987314e-06
## PRKACB      3.848542e-10  0.5272330 0.311 0.059 7.308767e-06
## EFR3A       4.166558e-10  0.5704691 0.348 0.083 7.912711e-06
## SMIM25      4.504412e-10 -0.9960988 0.030 0.302 8.554329e-06
## EPB41L3     4.518077e-10  0.5592335 0.568 0.234 8.580280e-06
## FSCN1       5.050003e-10  0.7292698 0.553 0.239 9.590461e-06
## FCGBP       5.172481e-10  1.2150503 0.333 0.078 9.823059e-06
## ETV5        6.402695e-10  0.5930568 0.515 0.185 1.215936e-05
## SLC25A6     6.950414e-10  0.5552488 0.674 0.351 1.319953e-05
## ZFHX3       7.382802e-10  0.5662717 0.803 0.502 1.402068e-05
## HMGB2       7.406968e-10 -1.2380377 0.189 0.483 1.406657e-05
## TRIM14      7.677605e-10  0.5409893 0.326 0.068 1.458054e-05
## RAPH1       7.711748e-10  0.5662522 0.432 0.132 1.464538e-05
## PPIF        7.854341e-10 -1.2779389 0.121 0.410 1.491618e-05
## IL18        7.924604e-10  0.5648788 0.583 0.254 1.504962e-05
## AC245128.3  8.065214e-10 -1.0385799 0.045 0.322 1.531665e-05
## DBI         9.386898e-10  0.6761165 0.523 0.210 1.782666e-05
## ADGRE5      9.688903e-10 -0.9213279 0.061 0.346 1.840019e-05
## ICAM3       9.828275e-10 -0.9726724 0.000 0.244 1.866488e-05
## FNBP1       9.849049e-10  0.5604742 0.606 0.268 1.870433e-05
## MFSD1       1.029772e-09  0.5385314 0.614 0.263 1.955640e-05
## FCGRT       1.079087e-09  0.5647200 0.750 0.429 2.049293e-05
## FCHSD2      1.251659e-09  0.5040052 0.523 0.195 2.377025e-05
## SPON1       1.326753e-09  0.4907824 0.235 0.024 2.519637e-05
## IDH1        1.337942e-09  0.5916942 0.409 0.122 2.540886e-05
## RGS10       1.350794e-09  0.4844662 0.780 0.395 2.565294e-05
## UACA        1.380694e-09  0.5751090 0.235 0.024 2.622077e-05
## SPRED1      1.458921e-09  0.5936223 0.515 0.205 2.770636e-05
## AC020916.1  1.530961e-09 -1.1182767 0.076 0.351 2.907448e-05
## KDM6B       1.591476e-09 -1.0058110 0.212 0.488 3.022372e-05
## LMAN1       1.620474e-09  0.5608442 0.386 0.112 3.077441e-05
## HIST2H2AA4  1.879053e-09 -1.0786957 0.697 0.844 3.568509e-05
## HLA-B       2.038091e-09 -0.4597805 0.947 0.976 3.870538e-05
## RND3        2.274131e-09  0.5992471 0.205 0.015 4.318801e-05
## DUSP1       2.425802e-09 -0.6816089 0.818 0.922 4.606841e-05
## RBPJ        2.483344e-09  0.6014096 0.705 0.385 4.716118e-05
## PIK3R1      2.680923e-09  0.5506230 0.674 0.307 5.091341e-05
## FAM107B     3.016561e-09 -0.8393783 0.053 0.322 5.728752e-05
## STARD13     3.368336e-09  0.4438395 0.189 0.010 6.396807e-05
## CXCR1       3.671588e-09 -1.0125158 0.000 0.229 6.972713e-05
## PTPRC       3.751889e-09 -0.5848500 0.826 0.902 7.125213e-05
## FAM198B     3.835371e-09  0.6791754 0.394 0.117 7.283753e-05
## APPL1       3.923853e-09  0.6009186 0.470 0.176 7.451789e-05
## NRP2        4.093372e-09  0.7755508 0.447 0.156 7.773723e-05
## RHOBTB3     4.237218e-09  0.4741109 0.212 0.020 8.046901e-05
## MIDN        4.293266e-09 -0.9522354 0.265 0.527 8.153342e-05
## SFT2D2      4.534371e-09  0.6260763 0.500 0.205 8.611225e-05
## SSR1        4.815914e-09  0.5258987 0.697 0.429 9.145902e-05
## CXCR2       5.022158e-09 -0.9844009 0.030 0.278 9.537581e-05
## RPS15A      5.226413e-09  0.4433627 0.939 0.761 9.925482e-05
## SULF2       5.670504e-09 -0.7826520 0.000 0.224 1.076885e-04
## VCAN        5.670504e-09 -1.5896853 0.000 0.224 1.076885e-04
## TFEC        5.841067e-09  0.5004366 0.515 0.195 1.109277e-04
## TXNIP       6.447929e-09 -0.7272407 0.735 0.854 1.224526e-04
## SMCHD1      6.597918e-09 -0.9562115 0.280 0.522 1.253011e-04
## CD84        6.910226e-09  0.5861379 0.644 0.337 1.312321e-04
## MAP2K3      7.151537e-09 -0.9425828 0.106 0.376 1.358148e-04
## SLC9A9      7.597779e-09  0.4191685 0.568 0.224 1.442894e-04
## TAGAP       8.244033e-09 -0.8452165 0.174 0.468 1.565624e-04
## FOS         8.656065e-09 -0.7334933 0.833 0.956 1.643873e-04
## RPL8        8.920744e-09  0.5704023 0.856 0.605 1.694138e-04
## MEF2A       8.921630e-09  0.4718181 0.886 0.561 1.694307e-04
## CALHM6      9.101644e-09  0.4659612 0.189 0.015 1.728493e-04
## RAC2        9.393282e-09 -0.9204796 0.083 0.346 1.783878e-04
## MBOAT7      1.016746e-08 -0.8990827 0.152 0.429 1.930903e-04
## HERPUD1     1.137890e-08  0.6816401 0.773 0.522 2.160967e-04
## CREM        1.207916e-08  0.7492131 0.379 0.127 2.293954e-04
## RPL24       1.220705e-08  0.5079050 0.879 0.644 2.318241e-04
## CFP         1.343323e-08 -0.6803045 0.000 0.215 2.551104e-04
## NCF2        1.420051e-08 -0.8594982 0.220 0.493 2.696819e-04
## GLUL        1.495055e-08 -0.7645362 0.674 0.829 2.839259e-04
## IER2        1.783242e-08 -0.8298276 0.515 0.727 3.386554e-04
## PHACTR1     1.853843e-08 -0.9616241 0.144 0.424 3.520633e-04
## RAB11FIP1   2.116547e-08 -0.9510085 0.212 0.478 4.019534e-04
## CYTIP       2.316383e-08 -0.7613734 0.068 0.322 4.399044e-04
## SAMHD1      2.327744e-08  0.5339848 0.636 0.346 4.420618e-04
## RERE        2.350313e-08  0.4909219 0.583 0.273 4.463479e-04
## FRMD4B      2.381738e-08  0.3961182 0.485 0.190 4.523158e-04
## PCDH12      2.384840e-08  0.4688207 0.220 0.029 4.529050e-04
## CCDC144A    2.406541e-08  0.6680354 0.439 0.176 4.570262e-04
## AFF3        2.420148e-08  0.3916040 0.159 0.005 4.596103e-04
## RUFY3       2.663609e-08  0.5070375 0.424 0.151 5.058460e-04
## ITGB1       2.798122e-08  0.5385509 0.583 0.302 5.313914e-04
## RPL9        2.837841e-08  0.4640200 0.932 0.737 5.389343e-04
## WIPI2       2.909241e-08  0.5942691 0.364 0.122 5.524939e-04
## CD4         2.926883e-08  0.4966508 0.545 0.249 5.558444e-04
## PADI4       3.153841e-08 -1.1616741 0.000 0.205 5.989459e-04
## CD63        3.364610e-08  0.4657378 0.848 0.595 6.389731e-04
## SRGN        3.414199e-08 -0.5066817 0.864 0.932 6.483906e-04
## MAFF        3.507420e-08 -0.8109903 0.053 0.293 6.660941e-04
## CH25H       3.535768e-08  0.7055167 0.576 0.278 6.714777e-04
## AHNAK       3.540439e-08  0.5520123 0.470 0.190 6.723647e-04
## AES         3.637143e-08  0.5021027 0.364 0.117 6.907298e-04
## RPN2        3.760388e-08  0.5520187 0.409 0.156 7.141353e-04
## ACTN1       3.889330e-08 -0.7381541 0.038 0.273 7.386226e-04
## TNFRSF11A   4.179189e-08  0.4085734 0.182 0.015 7.936698e-04
## HK2         4.324952e-08 -0.7156031 0.045 0.283 8.213515e-04
## FOSB        4.387242e-08 -0.7984641 0.356 0.639 8.331812e-04
## MCL1        4.617556e-08 -0.7095947 0.561 0.741 8.769200e-04
## HCAR2       4.643677e-08 -0.9052021 0.023 0.249 8.818808e-04
## FLNA        4.655284e-08 -0.7417540 0.015 0.229 8.840849e-04
## PGLYRP1     4.816530e-08 -1.5584168 0.000 0.200 9.147072e-04
## SLCO2B1     4.834478e-08  0.4934554 0.712 0.395 9.181158e-04
## IL2RG       4.940250e-08 -0.9293399 0.053 0.288 9.382029e-04
## PCYOX1      5.117535e-08  0.4103981 0.212 0.029 9.718710e-04
## SERPINB6    5.121267e-08  0.5348446 0.636 0.317 9.725799e-04
## ARHGAP21    5.278738e-08  0.5439713 0.576 0.302 1.002485e-03
## LYZ         5.878878e-08 -2.4704268 0.439 0.585 1.116458e-03
## IGFBP4      5.891151e-08  0.4189033 0.136 0.000 1.118788e-03
## CD209       6.052739e-08  0.4217588 0.152 0.005 1.149476e-03
## S100A10     6.990256e-08 -0.8557643 0.015 0.224 1.327520e-03
## CMKLR1      7.244878e-08  0.5780804 0.333 0.098 1.375875e-03
## GPR34       7.256591e-08  0.4736207 0.689 0.356 1.378099e-03
## MME         7.339791e-08 -0.7955365 0.000 0.195 1.393900e-03
## KLF6        7.750643e-08 -0.4041612 0.659 0.824 1.471925e-03
## RIOK3       8.719141e-08 -0.9234026 0.220 0.468 1.655852e-03
## DSE         8.822547e-08  0.5424514 0.530 0.244 1.675490e-03
## CIITA       8.973477e-08  0.4515923 0.515 0.234 1.704153e-03
## FGD6        9.158485e-08  0.3790354 0.197 0.024 1.739288e-03
## SERPINA1    9.205013e-08 -0.7379701 0.265 0.522 1.748124e-03
## DPYSL2      9.729865e-08  0.4843697 0.795 0.537 1.847799e-03
## QKI         1.015135e-07  0.4436554 0.879 0.610 1.927842e-03
## RNASE6      1.073037e-07  0.5486263 0.576 0.283 2.037804e-03
## TNFAIP6     1.116116e-07 -0.8159430 0.000 0.190 2.119616e-03
## APOL6       1.154915e-07  0.4746959 0.356 0.117 2.193299e-03
## KIF16B      1.163959e-07  0.3663713 0.273 0.059 2.210475e-03
## HELZ        1.220316e-07  0.5260635 0.508 0.229 2.317502e-03
## HEXA        1.235216e-07  0.4865756 0.492 0.210 2.345799e-03
## LY6E        1.290801e-07  0.4268595 0.242 0.049 2.451361e-03
## RPL4        1.296024e-07  0.4367461 0.841 0.522 2.461280e-03
## TSC22D3     1.302392e-07 -0.7306695 0.227 0.522 2.473372e-03
## ZDHHC18     1.577972e-07 -0.9147660 0.030 0.239 2.996727e-03
## RACK1       1.586067e-07  0.3981865 0.795 0.478 3.012100e-03
## CD52        1.598877e-07 -1.1517281 0.015 0.215 3.036428e-03
## CPM         1.671454e-07  0.5912882 0.295 0.083 3.174257e-03
## BMPR2       1.673935e-07  0.4291597 0.273 0.063 3.178969e-03
## MAN1A1      1.686722e-07  0.5231477 0.417 0.156 3.203254e-03
## JUNB        1.726714e-07 -0.4281334 0.712 0.844 3.279203e-03
## SNHG8       1.866326e-07  0.4889567 0.568 0.302 3.544341e-03
## PDGFC       1.893666e-07  0.3868961 0.167 0.015 3.596261e-03
## NCOA7       1.941299e-07  0.4985273 0.508 0.234 3.686721e-03
## LEPROT      1.966799e-07  0.4095128 0.462 0.185 3.735148e-03
## APOC1       2.125839e-07  0.6132457 0.576 0.278 4.037181e-03
## BMP2K       2.193271e-07  0.5286626 0.652 0.341 4.165240e-03
## MZT2B       2.261314e-07  0.3937482 0.295 0.078 4.294461e-03
## CXCL12      2.301487e-07  0.4396187 0.432 0.171 4.370754e-03
## HGF         2.339129e-07  0.5237429 0.242 0.054 4.442241e-03
## SPATS2L     2.489602e-07  0.4282586 0.197 0.029 4.728003e-03
## NDUFA4      2.496343e-07  0.3627533 0.545 0.273 4.740804e-03
## STX3        2.551197e-07 -0.8527571 0.083 0.312 4.844979e-03
## CNRIP1      2.812830e-07  0.4477482 0.227 0.044 5.341846e-03
## TMEM106A    2.834573e-07  0.5535335 0.500 0.234 5.383138e-03
## CPPED1      3.227044e-07 -0.7481961 0.076 0.307 6.128479e-03
## EMP2        3.461607e-07  0.2965498 0.121 0.000 6.573937e-03
## HLA-DRB5    3.515184e-07  0.5743423 0.462 0.210 6.675687e-03
## PRKX        3.655274e-07  0.3956140 0.227 0.044 6.941731e-03
## PARVB       3.821837e-07  0.4004845 0.159 0.015 7.258051e-03
## MMP9        3.875409e-07 -1.6425037 0.000 0.176 7.359789e-03
## WDR36       4.032994e-07  0.4597853 0.242 0.054 7.659059e-03
## RIMKLB      4.165994e-07  0.4702744 0.220 0.044 7.911640e-03
## MYO7A       4.178430e-07  0.3580592 0.288 0.073 7.935256e-03
## IFI6        4.361772e-07  0.5554554 0.242 0.059 8.283442e-03
## FPR2        4.406686e-07 -0.5591325 0.015 0.205 8.368737e-03
## AP1S2       4.493211e-07  0.4103013 0.667 0.341 8.533057e-03
## LAIR1       4.666844e-07  0.4658151 0.621 0.317 8.862803e-03
## INSIG1      4.869688e-07  0.8290722 0.485 0.249 9.248025e-03
## CSTA        5.077080e-07 -0.6139367 0.038 0.244 9.641882e-03
## SDCBP       5.178801e-07 -0.5026857 0.735 0.810 9.835061e-03
## PAPSS2      5.244095e-07  0.4618805 0.182 0.024 9.959061e-03
## FXYD5       5.264346e-07  0.3911447 0.568 0.302 9.997520e-03
## SLFN11      5.622633e-07  0.3978971 0.333 0.107 1.067794e-02
## VNN3        6.096431e-07 -0.6832490 0.008 0.185 1.157773e-02
## SSPN        6.215350e-07  0.4258163 0.242 0.054 1.180357e-02
## CNN2        6.227891e-07 -0.6952194 0.008 0.185 1.182739e-02
## KCNJ5       6.409745e-07  0.3164886 0.144 0.010 1.217275e-02
## STK17B      6.589304e-07 -0.7094214 0.189 0.439 1.251375e-02
## EFHD2       6.759399e-07 -0.7081405 0.212 0.449 1.283677e-02
## CLTC        6.810195e-07  0.4671451 0.561 0.317 1.293324e-02
## IL1B        6.909791e-07 -1.6164613 0.114 0.337 1.312238e-02
## MFHAS1      7.379327e-07  0.4144026 0.144 0.010 1.401408e-02
## CREG1       7.459077e-07  0.4303847 0.576 0.293 1.416553e-02
## SCD         7.487106e-07  0.4507874 0.227 0.049 1.421876e-02
## PCBP1       7.682381e-07 -0.4988768 0.545 0.717 1.458961e-02
## RPS7        7.885860e-07  0.4597252 0.841 0.590 1.497604e-02
## GPRIN3      7.945480e-07  0.5250661 0.417 0.176 1.508926e-02
## EIF1        8.227188e-07 -0.4367273 0.856 0.868 1.562425e-02
## DDIT3       8.436933e-07 -0.8463440 0.129 0.361 1.602258e-02
## DSC2        8.508511e-07  0.4341046 0.242 0.059 1.615851e-02
## FAM213A     8.931438e-07  0.3084665 0.235 0.054 1.696169e-02
## ATP1B3      9.343455e-07  0.5047168 0.508 0.254 1.774416e-02
## CMTM2       9.493871e-07 -0.9728358 0.098 0.307 1.802981e-02
## CD300E      9.656165e-07 -0.5490360 0.008 0.180 1.833802e-02
## RAB27A      1.000067e-06 -0.6862017 0.038 0.234 1.899227e-02
## GRK3        1.001973e-06  0.4098851 0.371 0.141 1.902848e-02
## TCF12       1.014594e-06  0.5099114 0.371 0.146 1.926816e-02
## PRDX1       1.018445e-06  0.5146241 0.538 0.302 1.934129e-02
## TTC3        1.027688e-06  0.3935508 0.561 0.263 1.951682e-02
## CYTL1       1.071598e-06  0.4215282 0.295 0.088 2.035071e-02
## ABHD12      1.083133e-06  0.4369955 0.235 0.054 2.056977e-02
## LDHB        1.104349e-06  0.5442833 0.280 0.088 2.097269e-02
## CTSS        1.134683e-06 -0.5526506 0.886 0.893 2.154876e-02
## TPCN1       1.145101e-06  0.4255354 0.386 0.146 2.174662e-02
## CCND2       1.190814e-06  0.4106735 0.379 0.141 2.261474e-02
## RPL29       1.226670e-06  0.4504012 0.909 0.639 2.329569e-02
## TTL         1.228710e-06  0.5005367 0.288 0.088 2.333443e-02
## MSN         1.239259e-06 -0.7505226 0.364 0.537 2.353477e-02
## MRPS34      1.265503e-06  0.3877103 0.235 0.059 2.403316e-02
## CRIM1       1.275267e-06  0.3451496 0.152 0.015 2.421860e-02
## PTMS        1.291810e-06  0.3079309 0.280 0.078 2.453276e-02
## GRID2       1.296055e-06 -0.8252066 0.068 0.273 2.461338e-02
## BST1        1.321334e-06 -0.5457606 0.000 0.161 2.509346e-02
## HNMT        1.435721e-06  0.4720061 0.402 0.180 2.726578e-02
## ZFAND5      1.438952e-06 -0.7337860 0.341 0.527 2.732713e-02
## RNF24       1.516694e-06 -0.6562405 0.030 0.215 2.880354e-02
## CPD         1.558781e-06 -0.6552576 0.053 0.259 2.960281e-02
## LGALS3      1.591995e-06 -0.6215280 0.023 0.205 3.023358e-02
## TMEM70      1.594714e-06  0.3557463 0.280 0.083 3.028521e-02
## FAM91A1     1.613978e-06  0.4463193 0.432 0.185 3.065106e-02
## DNMT1       1.623574e-06  0.4984002 0.303 0.102 3.083329e-02
## CD59        1.637444e-06  0.3608326 0.311 0.098 3.109670e-02
## BCL2        1.655956e-06  0.5121894 0.379 0.146 3.144827e-02
## DRAM1       1.675224e-06  0.3774616 0.288 0.088 3.181418e-02
## EPHX1       1.749758e-06  0.3603716 0.182 0.029 3.322965e-02
## NFIL3       1.767491e-06 -0.7011362 0.152 0.385 3.356642e-02
## SNX13       1.802641e-06  0.4146396 0.417 0.171 3.423395e-02
## SKP2        1.926790e-06  0.3939449 0.250 0.063 3.659166e-02
## AAK1        1.963742e-06  0.5120176 0.576 0.337 3.729343e-02
## MCEMP1      1.981006e-06 -0.6346240 0.000 0.156 3.762129e-02
## CCDC141     2.007074e-06  0.3492003 0.106 0.000 3.811634e-02
## TKT         2.030436e-06 -0.7262461 0.197 0.424 3.856001e-02
## TNFAIP3     2.057979e-06 -1.1131286 0.212 0.415 3.908308e-02
## BCL6        2.069492e-06 -0.8485903 0.250 0.478 3.930173e-02
## CD48        2.102711e-06 -0.6446364 0.038 0.224 3.993259e-02
## EHD1        2.161691e-06 -0.7114254 0.053 0.244 4.105267e-02
## TCEAL4      2.173732e-06  0.4893556 0.402 0.166 4.128135e-02
## RPL28       2.300794e-06  0.3573793 0.977 0.849 4.369437e-02
## SET         2.337666e-06  0.4354577 0.636 0.395 4.439462e-02
## FAM198B-AS1 2.368150e-06  0.4298903 0.159 0.020 4.497354e-02
## UBE2D1      2.469832e-06 -0.7234864 0.106 0.317 4.690458e-02
## DAAM1       2.541859e-06  0.3702175 0.280 0.088 4.827244e-02
## UQCRH       2.583932e-06  0.3599469 0.462 0.220 4.907145e-02
## OGFRL1      2.700758e-06  0.4268781 0.811 0.600 5.129010e-02
## MPEG1       2.714912e-06  0.4695016 0.553 0.298 5.155889e-02
## ABCC5       2.804480e-06  0.4056877 0.273 0.083 5.325989e-02
## EPS8        2.819469e-06  0.3197245 0.538 0.273 5.354453e-02
## LGALS3BP    2.917726e-06  0.3192314 0.189 0.034 5.541054e-02
## CKAP4       2.918968e-06 -0.8389188 0.098 0.307 5.543412e-02
## CRIP1       2.964417e-06 -0.6136325 0.000 0.151 5.629724e-02
## SEL1L       2.992412e-06  0.3112804 0.447 0.200 5.682890e-02
## RDX         3.008952e-06  0.4528664 0.379 0.156 5.714301e-02
## CD164       3.015598e-06  0.4010330 0.742 0.483 5.726923e-02
## TRAF3IP3    3.015821e-06 -0.6746630 0.098 0.312 5.727345e-02
## PLEKHO1     3.062898e-06  0.4377728 0.591 0.337 5.816750e-02
## LTF         3.098028e-06 -3.2901299 0.008 0.166 5.883465e-02
## SYNRG       3.130783e-06  0.2823582 0.280 0.078 5.945670e-02
## FGL2        3.132704e-06  0.4738367 0.848 0.585 5.949318e-02
## NUCKS1      3.143354e-06  0.4703324 0.598 0.346 5.969543e-02
## CMTM3       3.170263e-06  0.5006333 0.379 0.171 6.020647e-02
## CYSTM1      3.256953e-06 -0.8394215 0.045 0.229 6.185279e-02
## FCAR        3.359844e-06 -0.6280777 0.023 0.195 6.380679e-02
## TPP1        3.423646e-06  0.3971504 0.652 0.371 6.501846e-02
## SORBS3      3.439918e-06  0.3162220 0.144 0.015 6.532747e-02
## PPIA        3.477637e-06  0.3727290 0.773 0.468 6.604380e-02
## LTC4S       3.497966e-06  0.4491495 0.523 0.273 6.642986e-02
## GAS6        3.580809e-06  0.4430356 0.333 0.122 6.800314e-02
## ADAM9       3.603271e-06  0.4071728 0.439 0.205 6.842972e-02
## PRRC1       3.739947e-06  0.3373063 0.197 0.039 7.102533e-02
## NCL         3.821514e-06  0.4457017 0.621 0.380 7.257437e-02
## DHRS3       3.879074e-06  0.4172396 0.273 0.083 7.366749e-02
## RHOG        3.890041e-06 -0.7493325 0.250 0.463 7.387578e-02
## RPS4Y1      4.010072e-06  0.5103500 0.470 0.254 7.615527e-02
## MYO1F       4.073503e-06 -0.6100959 0.409 0.605 7.735990e-02
## CDC42       4.254829e-06 -0.5179008 0.614 0.737 8.080345e-02
## TBC1D12     4.260706e-06  0.3992487 0.288 0.093 8.091506e-02
## SERF2       4.415667e-06 -0.4682272 0.735 0.776 8.385793e-02
## MEFV        4.427819e-06 -0.5702096 0.000 0.146 8.408870e-02
## FNDC3A      4.488201e-06  0.4513717 0.386 0.166 8.523542e-02
## CD99        4.651265e-06  0.4343367 0.644 0.400 8.833217e-02
## RNF130      4.682105e-06  0.3937919 0.629 0.371 8.891785e-02
## BIN2        4.700516e-06 -0.6771689 0.167 0.380 8.926750e-02
## PIM3        4.855434e-06 -0.5678796 0.030 0.205 9.220955e-02
## ANXA2       4.959874e-06 -0.6052064 0.023 0.190 9.419297e-02
## DNAJC13     4.962189e-06  0.4104940 0.341 0.132 9.423694e-02
## PGD         5.535413e-06 -0.7415648 0.121 0.327 1.051230e-01
## EEA1        5.626881e-06  0.4567744 0.386 0.180 1.068601e-01
## MSRB1       5.858936e-06 -0.5983620 0.015 0.176 1.112671e-01
## IFITM3      6.037895e-06  0.3788294 0.614 0.346 1.146657e-01
## TSPAN33     6.046204e-06  0.3015407 0.182 0.034 1.148235e-01
## HLA-DOA     6.066301e-06  0.4670690 0.273 0.093 1.152051e-01
## CLTA        6.120219e-06  0.4104373 0.379 0.161 1.162291e-01
## ERCC1       6.185454e-06  0.3450478 0.348 0.132 1.174680e-01
## RTN3        6.416376e-06 -0.7367801 0.280 0.483 1.218534e-01
## EVI2B       6.591799e-06 -0.7061584 0.477 0.620 1.251849e-01
## RNF144B     6.609624e-06  0.4938548 0.371 0.161 1.255234e-01
## PLA2G16     6.708591e-06  0.4250099 0.250 0.078 1.274028e-01
## FAM129A     7.104220e-06 -0.7382880 0.098 0.293 1.349162e-01
## NENF        7.183151e-06  0.3663738 0.174 0.034 1.364152e-01
## VEGFB       7.235963e-06  0.3568352 0.227 0.059 1.374182e-01
## TNFAIP8     7.277640e-06  0.3952046 0.288 0.102 1.382097e-01
## EIF4B       7.377161e-06  0.4221100 0.591 0.341 1.400997e-01
## ARPC1B      7.558891e-06  0.3233119 0.765 0.556 1.435509e-01
## CAP1        7.627513e-06 -0.6842932 0.561 0.683 1.448541e-01
## YWHAE       7.763720e-06  0.3578426 0.750 0.468 1.474408e-01
## GALNT1      7.766292e-06  0.2770712 0.394 0.166 1.474896e-01
## EIF2AK4     8.023097e-06  0.4048577 0.288 0.098 1.523666e-01
## UHRF1BP1    8.305450e-06  0.4127723 0.212 0.054 1.577288e-01
## RHOA        8.351223e-06 -0.5379663 0.644 0.761 1.585981e-01
## GLIPR2      8.495443e-06 -0.5026083 0.023 0.185 1.613370e-01
## SOD1        8.708955e-06  0.4680510 0.409 0.195 1.653918e-01
## KTN1        8.904700e-06  0.3282009 0.682 0.400 1.691092e-01
## CCDC88A     8.952310e-06  0.4083243 0.750 0.502 1.700133e-01
## PRDM1       9.010312e-06  0.4069902 0.356 0.137 1.711148e-01
## CTSC        9.221749e-06  0.4267496 0.659 0.424 1.751302e-01
## PLK3        9.254436e-06 -0.5616372 0.114 0.322 1.757510e-01
## TOMM20      9.281602e-06  0.4361201 0.432 0.220 1.762669e-01
## SPIRE1      9.489403e-06  0.2547447 0.121 0.010 1.802132e-01
## ADAM8       9.825519e-06 -0.5184085 0.000 0.137 1.865964e-01
## RRBP1       9.889279e-06  0.3710204 0.447 0.215 1.878073e-01
## DEDD2       1.034788e-05 -0.7663549 0.098 0.288 1.965167e-01
## SLC18B1     1.037166e-05  0.3745794 0.326 0.122 1.969681e-01
## EREG        1.096234e-05 -1.4474357 0.008 0.151 2.081858e-01
## UBALD2      1.099371e-05 -0.6135665 0.061 0.244 2.087816e-01
## NORAD       1.109485e-05  0.3395539 0.485 0.239 2.107024e-01
## MTDH        1.158844e-05  0.3191593 0.629 0.351 2.200760e-01
## KIDINS220   1.176684e-05  0.3270814 0.477 0.234 2.234641e-01
## P4HA1       1.219963e-05  0.3032136 0.409 0.185 2.316832e-01
## TALDO1      1.231376e-05 -0.6595265 0.212 0.424 2.338506e-01
## PTK2B       1.241991e-05 -0.6988542 0.167 0.361 2.358665e-01
## ZNF281      1.248432e-05 -0.5912909 0.121 0.322 2.370897e-01
## CALM2       1.256228e-05  0.3894622 0.871 0.683 2.385702e-01
## LYN         1.285220e-05 -0.6395778 0.394 0.551 2.440762e-01
## HGSNAT      1.316378e-05  0.3859215 0.303 0.117 2.499934e-01
## HIST1H1C    1.335769e-05 -1.1927972 0.197 0.395 2.536759e-01
## IKBKB       1.335975e-05  0.3378754 0.348 0.137 2.537149e-01
## C6orf48     1.342086e-05  0.3621420 0.348 0.141 2.548756e-01
## BLVRB       1.373218e-05  0.3567735 0.402 0.176 2.607878e-01
## MAP1LC3B    1.374229e-05 -0.7665882 0.371 0.556 2.609799e-01
## TIMP1       1.422706e-05 -1.0179073 0.106 0.288 2.701862e-01
## NOP53       1.436461e-05  0.4544794 0.447 0.234 2.727984e-01
## CRTAP       1.441623e-05  0.3973612 0.439 0.220 2.737786e-01
## MYO1G       1.459857e-05 -0.3838466 0.000 0.132 2.772414e-01
## LILRA5      1.459857e-05 -0.4990805 0.000 0.132 2.772414e-01
## PTX3        1.459857e-05 -0.6479133 0.000 0.132 2.772414e-01
## ARG1        1.459857e-05 -0.9194434 0.000 0.132 2.772414e-01
## LCN2        1.459857e-05 -2.6492150 0.000 0.132 2.772414e-01
## HIST1H2BC   1.504993e-05 -0.7346274 0.061 0.234 2.858132e-01
## IL13RA1     1.513079e-05  0.3548936 0.621 0.361 2.873488e-01
## LGALS1      1.537570e-05 -1.0015887 0.114 0.293 2.919999e-01
## SLC5A3      1.568051e-05  0.6087856 0.402 0.195 2.977886e-01
## SOAT1       1.592458e-05  0.4110260 0.462 0.239 3.024238e-01
## BTF3        1.706140e-05  0.4183080 0.621 0.390 3.240130e-01
## EIF3F       1.740540e-05  0.3627432 0.386 0.176 3.305459e-01
## ZNF358      1.746343e-05  0.2784635 0.129 0.015 3.316479e-01
## ETNK1       1.867482e-05  0.4425047 0.470 0.249 3.546536e-01
## GNB4        1.871544e-05  0.2582502 0.553 0.283 3.554249e-01
## GMFG        1.930386e-05 -0.6670417 0.371 0.541 3.665996e-01
## VCL         1.932066e-05 -0.4581500 0.023 0.176 3.669187e-01
## EBI3        1.939902e-05  0.3664899 0.197 0.049 3.684068e-01
## ETV6        1.976205e-05  0.4423247 0.447 0.220 3.753010e-01
## GOLGA4      2.010839e-05  0.2802198 0.583 0.337 3.818784e-01
## GGCX        2.078543e-05  0.5045288 0.356 0.156 3.947362e-01
## SWAP70      2.138812e-05  0.4094146 0.394 0.180 4.061818e-01
## ARHGDIB     2.153032e-05 -0.5238911 0.470 0.644 4.088823e-01
## QPCT        2.165418e-05 -0.4505200 0.000 0.127 4.112345e-01
## FFAR2       2.165418e-05 -0.5553843 0.000 0.127 4.112345e-01
## CD24        2.165418e-05 -1.2525166 0.000 0.127 4.112345e-01
## SUCLG2      2.166387e-05  0.3410877 0.288 0.107 4.114186e-01
## IRF1        2.339658e-05 -0.8190589 0.205 0.395 4.443245e-01
## FEM1C       2.425129e-05 -0.6282502 0.068 0.244 4.605562e-01
## PRNP        2.481382e-05  0.3833318 0.424 0.205 4.712393e-01
## AC245297.3  2.522097e-05  0.3785256 0.568 0.327 4.789714e-01
## NOP10       2.547360e-05 -0.6396627 0.288 0.468 4.837692e-01
## APMAP       2.548692e-05  0.3474456 0.515 0.268 4.840221e-01
## STX4        2.551928e-05  0.4404581 0.303 0.127 4.846366e-01
## MTR         2.574342e-05  0.3720057 0.235 0.073 4.888933e-01
## HIST2H2AC   2.578447e-05 -0.5510052 0.068 0.244 4.896729e-01
## UBE2R2      2.678206e-05 -0.5971360 0.159 0.356 5.086181e-01
## GLIPR1      2.714898e-05  0.3810928 0.697 0.468 5.155862e-01
## GLCE        2.731254e-05  0.2506215 0.197 0.049 5.186924e-01
## UNC93B1     2.782908e-05  0.4124277 0.402 0.200 5.285020e-01
## PTGER4      2.794203e-05  0.3302713 0.348 0.146 5.306471e-01
## PARP14      2.799136e-05  0.4489899 0.568 0.337 5.315839e-01
## ANXA1       2.803256e-05 -1.0846479 0.174 0.351 5.323664e-01
## TNFRSF25    2.822413e-05  0.2814830 0.265 0.093 5.360045e-01
## LAP3        2.828380e-05  0.4146002 0.356 0.156 5.371377e-01
## THOC1       2.867144e-05  0.3902220 0.250 0.083 5.444994e-01
## RAP2A       2.883180e-05  0.3597609 0.303 0.122 5.475447e-01
## IER5        2.901496e-05 -0.7757696 0.311 0.517 5.510231e-01
## LUCAT1      2.966017e-05 -0.4791749 0.015 0.156 5.632763e-01
## CNEP1R1     3.080464e-05  0.4045307 0.205 0.059 5.850110e-01
## REL         3.093587e-05  0.3958006 0.833 0.644 5.875032e-01
## SNCA        3.114435e-05  0.2734543 0.174 0.039 5.914624e-01
## MXD4        3.138340e-05  0.3241070 0.379 0.171 5.960022e-01
## STAT2       3.226675e-05  0.3111701 0.265 0.093 6.127778e-01
## GRN         3.255043e-05  0.3021245 0.788 0.551 6.181653e-01
## SPI1        3.550284e-05 -0.5675162 0.409 0.576 6.742345e-01
## TAPT1       3.707350e-05  0.3155291 0.295 0.112 7.040629e-01
## SIRPB1      3.715682e-05 -0.6101426 0.030 0.180 7.056451e-01
## STXBP2      3.720874e-05 -0.6454428 0.152 0.341 7.066311e-01
## BTF3L4      3.805922e-05  0.4392643 0.394 0.200 7.227826e-01
## MTUS1       3.892488e-05  0.3183239 0.341 0.141 7.392225e-01
## DOCK7       3.950939e-05  0.2947249 0.220 0.063 7.503228e-01
## CRYL1       3.958355e-05  0.3071971 0.159 0.034 7.517312e-01
## KDSR        3.963491e-05  0.4368144 0.280 0.112 7.527065e-01
## LYST        4.018585e-05 -0.5596114 0.356 0.522 7.631695e-01
## FAM110B     4.025740e-05  0.4152661 0.152 0.029 7.645282e-01
## UQCR10      4.083460e-05  0.2814053 0.591 0.351 7.754899e-01
## C3orf58     4.193463e-05  0.4568640 0.379 0.185 7.963806e-01
## DCP2        4.432950e-05  0.4396232 0.417 0.210 8.418616e-01
## GHDC        4.488840e-05  0.2936141 0.197 0.054 8.524757e-01
## DAZAP2      4.585753e-05 -0.4968211 0.583 0.702 8.708804e-01
## SERINC5     4.596732e-05  0.4006574 0.250 0.093 8.729654e-01
## IRS2        4.600019e-05  0.4434584 0.606 0.424 8.735895e-01
## M6PR        4.619764e-05  0.4429195 0.470 0.278 8.773393e-01
## LGALS2      4.741675e-05 -0.3246539 0.000 0.117 9.004914e-01
## IL1R2       4.741675e-05 -0.4239102 0.000 0.117 9.004914e-01
## ATP2A3      4.741675e-05 -0.4513394 0.000 0.117 9.004914e-01
## TNFRSF10C   4.741675e-05 -0.5008500 0.000 0.117 9.004914e-01
## CRISP3      4.741675e-05 -1.4252664 0.000 0.117 9.004914e-01
## PID1        5.001324e-05  0.3593453 0.258 0.093 9.498014e-01
## HYPK        5.078611e-05  0.3543845 0.189 0.049 9.644789e-01
## ESR1        5.079327e-05  0.2688344 0.106 0.010 9.646149e-01
## PRAM1       5.115059e-05 -0.6278485 0.121 0.293 9.714009e-01
## PTP4A2      5.120478e-05 -0.6802489 0.409 0.532 9.724299e-01
## PTEN        5.178899e-05 -0.6893729 0.295 0.468 9.835247e-01
## TES         5.186778e-05 -0.5009981 0.053 0.215 9.850211e-01
## ZBTB7A      5.304654e-05 -0.5032471 0.152 0.341 1.000000e+00
## JPT1        5.363448e-05 -0.6263182 0.076 0.239 1.000000e+00
## FBXW11      5.368069e-05  0.3267699 0.265 0.098 1.000000e+00
## RETN        5.415617e-05 -1.1484502 0.015 0.146 1.000000e+00
## MGAM        5.443678e-05 -0.6595071 0.008 0.132 1.000000e+00
## CHPT1       5.447182e-05  0.2964031 0.235 0.078 1.000000e+00
## TMCO1       5.614312e-05  0.3506886 0.432 0.224 1.000000e+00
## DPP7        5.719200e-05  0.2992953 0.348 0.151 1.000000e+00
## PLAUR       5.760370e-05 -0.7928442 0.379 0.541 1.000000e+00
## CFAP97      5.771940e-05  0.2910164 0.295 0.117 1.000000e+00
## GPAT3       5.893172e-05 -0.5225008 0.053 0.210 1.000000e+00
## PARP1       5.894906e-05  0.3722744 0.273 0.107 1.000000e+00
## CTDSPL2     5.961572e-05  0.3678109 0.227 0.078 1.000000e+00
## AKIRIN2     6.014869e-05 -0.4656430 0.076 0.249 1.000000e+00
## NDUFB7      6.172797e-05  0.3084155 0.235 0.083 1.000000e+00
## ATG7        6.217874e-05  0.3542692 0.258 0.098 1.000000e+00
## ZFP36L2     6.220482e-05  0.3646749 0.879 0.732 1.000000e+00
## VMP1        6.251810e-05  0.7327747 0.826 0.717 1.000000e+00
## TBC1D4      6.268274e-05  0.2885159 0.386 0.176 1.000000e+00
## FAM20A      6.274244e-05  0.3427425 0.114 0.015 1.000000e+00
## FTX         6.290415e-05  0.3844924 0.447 0.229 1.000000e+00
## CENPC       6.299446e-05  0.3595410 0.273 0.107 1.000000e+00
## DUSP6       6.349080e-05 -0.6891861 0.205 0.380 1.000000e+00
## GTPBP1      6.404147e-05 -0.4945011 0.038 0.185 1.000000e+00
## SELENOM     6.409001e-05  0.3511208 0.136 0.024 1.000000e+00
## CEBPD       6.444105e-05  0.4786419 0.720 0.517 1.000000e+00
## ATF7IP      6.473306e-05  0.3816822 0.500 0.263 1.000000e+00
## FMN1        6.583577e-05  0.3458997 0.273 0.107 1.000000e+00
## IGF2R       6.713256e-05 -0.6107963 0.091 0.259 1.000000e+00
## VAMP8       6.807362e-05  0.3979935 0.439 0.244 1.000000e+00
## GON4L       6.812216e-05  0.3356258 0.250 0.088 1.000000e+00
## RPL7        6.980552e-05  0.3223636 0.826 0.590 1.000000e+00
## CAMP        7.000743e-05 -2.4851260 0.000 0.112 1.000000e+00
## LGALS9      7.016087e-05  0.3511949 0.500 0.273 1.000000e+00
## MYO5A       7.046003e-05  0.2661757 0.318 0.132 1.000000e+00
## GATM        7.124039e-05  0.3032854 0.492 0.244 1.000000e+00
## DMD         7.460601e-05  0.2594031 0.114 0.015 1.000000e+00
## MAST3       7.476642e-05 -0.4621305 0.045 0.195 1.000000e+00
## HACD2       7.523484e-05  0.2900286 0.182 0.049 1.000000e+00
## ATP2B1-AS1  7.529688e-05 -0.6609963 0.197 0.395 1.000000e+00
## TBC1D9      7.601458e-05  0.3749427 0.364 0.166 1.000000e+00
## LMBRD2      7.633710e-05  0.3993502 0.212 0.068 1.000000e+00
## NPM1        7.706726e-05  0.3545762 0.538 0.337 1.000000e+00
## IRAK3       7.775122e-05 -0.5796574 0.152 0.332 1.000000e+00
## SERBP1      8.215622e-05  0.3127421 0.591 0.380 1.000000e+00
## ACTB        8.255745e-05 -0.4251378 0.992 0.976 1.000000e+00
## NDRG1       8.293732e-05 -0.5585839 0.068 0.224 1.000000e+00
## HDAC9       8.481780e-05  0.3215033 0.318 0.146 1.000000e+00
## AP003481.1  8.495199e-05 -0.4609521 0.053 0.210 1.000000e+00
## PDIA6       8.496776e-05  0.3679589 0.364 0.171 1.000000e+00
## VPS13C      8.578765e-05  0.3747987 0.727 0.488 1.000000e+00
## RAB14       8.617782e-05  0.3359481 0.470 0.254 1.000000e+00
## WSB1        8.743889e-05 -0.4749009 0.856 0.859 1.000000e+00
## ZNF331      8.764171e-05  0.6356525 0.295 0.132 1.000000e+00
## ERLEC1      8.794806e-05  0.3480933 0.220 0.073 1.000000e+00
## CLPB        8.861673e-05  0.2576349 0.174 0.044 1.000000e+00
## ENOX2       8.948864e-05  0.3762597 0.121 0.020 1.000000e+00
## SH3PXD2A    8.965634e-05  0.2774430 0.114 0.015 1.000000e+00
## PDGFB       8.973518e-05  0.3747104 0.371 0.180 1.000000e+00
## CALR        9.155877e-05  0.3569876 0.583 0.366 1.000000e+00
## IFI44       9.235486e-05  0.2980127 0.197 0.059 1.000000e+00
## XRN2        9.319946e-05  0.4407787 0.303 0.137 1.000000e+00
## HEATR5A     9.367164e-05  0.3359043 0.152 0.034 1.000000e+00
## FYB1        9.467652e-05 -0.4427255 0.720 0.800 1.000000e+00
## TCIRG1      9.769237e-05 -0.5855150 0.129 0.307 1.000000e+00
## LAPTM4A     9.818372e-05  0.3389551 0.568 0.332 1.000000e+00
## AC015912.3  9.946607e-05 -0.3865775 0.023 0.156 1.000000e+00
## PI4K2A      9.987435e-05  0.2676946 0.250 0.088 1.000000e+00
## EMC1        1.000978e-04  0.2865683 0.189 0.054 1.000000e+00
## ALAS1       1.011957e-04 -0.5268487 0.045 0.190 1.000000e+00
## TMEM59      1.031556e-04  0.3025595 0.477 0.259 1.000000e+00
## WNT5A       1.032147e-04 -0.4109261 0.000 0.107 1.000000e+00
## MMP8        1.032147e-04 -1.0990466 0.000 0.107 1.000000e+00
## SSB         1.050781e-04  0.3491195 0.326 0.146 1.000000e+00
## FAU         1.051019e-04  0.2936723 0.932 0.751 1.000000e+00
## FEZ2        1.053686e-04  0.2637386 0.371 0.176 1.000000e+00
## IFNGR1      1.066621e-04 -0.4539604 0.470 0.620 1.000000e+00
## MGST3       1.067314e-04  0.3494483 0.348 0.166 1.000000e+00
## CCDC112     1.080456e-04  0.3201477 0.212 0.068 1.000000e+00
## LPAR6       1.089204e-04  0.2886740 0.864 0.571 1.000000e+00
## KANTR       1.101403e-04  0.2525999 0.144 0.029 1.000000e+00
## LTA4H       1.130302e-04 -0.8011044 0.152 0.332 1.000000e+00
## LIMK2       1.152156e-04 -0.4873620 0.023 0.151 1.000000e+00
## ALOX5       1.159204e-04 -0.5673357 0.258 0.429 1.000000e+00
## GRK6        1.162518e-04 -0.4468448 0.023 0.151 1.000000e+00
## ATG3        1.177057e-04 -0.5453666 0.144 0.322 1.000000e+00
## STX2        1.199444e-04  0.3653762 0.212 0.073 1.000000e+00
## EIF3G       1.253064e-04  0.4126026 0.379 0.200 1.000000e+00
## YPEL5       1.256630e-04 -0.5814603 0.280 0.449 1.000000e+00
## HOOK3       1.258752e-04  0.2523751 0.462 0.239 1.000000e+00
## GSAP        1.277665e-04  0.3288921 0.364 0.171 1.000000e+00
## DRAM2       1.289395e-04  0.4035242 0.371 0.200 1.000000e+00
## PA2G4       1.329927e-04  0.4599253 0.379 0.210 1.000000e+00
## STAG1       1.333553e-04  0.2540582 0.273 0.107 1.000000e+00
## SSR4        1.358794e-04  0.3670498 0.462 0.263 1.000000e+00
## ARL5B       1.385121e-04 -0.6359814 0.129 0.288 1.000000e+00
## DDX18       1.414258e-04  0.3041113 0.409 0.210 1.000000e+00
## TPM1        1.415059e-04  0.3849994 0.288 0.122 1.000000e+00
## SLC2A5      1.426738e-04  0.3456370 0.402 0.205 1.000000e+00
## ZBTB38      1.427326e-04  0.2501824 0.295 0.122 1.000000e+00
## CREB3L2     1.464714e-04  0.2915342 0.280 0.112 1.000000e+00
## GAL3ST4     1.469274e-04  0.3135602 0.273 0.112 1.000000e+00
## LINC01001   1.469740e-04 -0.6834127 0.076 0.224 1.000000e+00
## ATP2A2      1.478574e-04  0.3562943 0.371 0.200 1.000000e+00
## EFCAB14     1.483931e-04  0.2590762 0.386 0.195 1.000000e+00
## LILRB3      1.490413e-04 -0.5558858 0.227 0.410 1.000000e+00
## COMT        1.507320e-04  0.2971152 0.220 0.078 1.000000e+00
## SYNJ2BP     1.512270e-04  0.3324264 0.197 0.063 1.000000e+00
## PLAC8       1.519702e-04 -0.3225120 0.000 0.102 1.000000e+00
## ANXA3       1.519702e-04 -0.5213474 0.000 0.102 1.000000e+00
## CHI3L1      1.519702e-04 -0.6295371 0.000 0.102 1.000000e+00
## RAP2B       1.519728e-04  0.4133740 0.508 0.307 1.000000e+00
## ST13        1.601940e-04  0.3350638 0.545 0.346 1.000000e+00
## NFKBIZ      1.614950e-04 -0.6072516 0.553 0.649 1.000000e+00
## TMEM236     1.622793e-04  0.3274322 0.174 0.049 1.000000e+00
## HPGDS       1.629155e-04  0.2919568 0.235 0.083 1.000000e+00
## RPS6KA3     1.684489e-04  0.3067470 0.348 0.171 1.000000e+00
## EIF4E3      1.696563e-04 -0.5060128 0.023 0.146 1.000000e+00
## GOLIM4      1.703916e-04  0.3159918 0.212 0.073 1.000000e+00
## FADS1       1.705579e-04  0.4742124 0.220 0.083 1.000000e+00
## EMP3        1.730583e-04 -0.7898814 0.227 0.371 1.000000e+00
## ASH2L       1.792321e-04  0.2713026 0.106 0.015 1.000000e+00
## SEPT2       1.807345e-04  0.3192904 0.447 0.254 1.000000e+00
## NARS        1.825637e-04  0.2664111 0.318 0.141 1.000000e+00
## LINC00528   1.847558e-04 -0.4797668 0.038 0.171 1.000000e+00
## PCM1        1.857772e-04  0.2758965 0.530 0.312 1.000000e+00
## RGCC        1.869667e-04 -0.4813216 0.008 0.117 1.000000e+00
## LY86        1.870176e-04  0.2710923 0.462 0.249 1.000000e+00
## GAA         1.870916e-04  0.3718069 0.326 0.151 1.000000e+00
## GNA13       1.891322e-04  0.3140664 0.545 0.332 1.000000e+00
## GIMAP8      1.906578e-04  0.2744612 0.250 0.098 1.000000e+00
## MBNL2       1.915152e-04  0.3484690 0.250 0.102 1.000000e+00
## LRRK2       1.935263e-04 -0.5829643 0.265 0.420 1.000000e+00
## KCNE1       1.946028e-04 -0.3917402 0.008 0.117 1.000000e+00
## FGD5-AS1    1.947472e-04  0.2620312 0.455 0.234 1.000000e+00
## CAB39       1.980912e-04  0.2950418 0.394 0.215 1.000000e+00
## IQGAP2      1.982871e-04  0.3299057 0.326 0.146 1.000000e+00
## MEGF9       1.986105e-04 -0.6241275 0.098 0.259 1.000000e+00
## SESN3       1.992065e-04  0.3742908 0.326 0.161 1.000000e+00
## SESN2       2.014085e-04 -0.4413866 0.015 0.132 1.000000e+00
## CDK6        2.020150e-04  0.3007034 0.288 0.117 1.000000e+00
## TLN2        2.035850e-04 -0.5061847 0.038 0.171 1.000000e+00
## YPEL3       2.046014e-04 -0.6836018 0.356 0.498 1.000000e+00
## TLR7        2.049940e-04  0.3012949 0.280 0.117 1.000000e+00
## OSCAR       2.051950e-04 -0.3699446 0.015 0.132 1.000000e+00
## TMEM14C     2.053635e-04  0.3270672 0.273 0.112 1.000000e+00
## USP14       2.062798e-04  0.2654852 0.250 0.098 1.000000e+00
## GOLGB1      2.072559e-04  0.3073379 0.508 0.293 1.000000e+00
## ARRDC3      2.097967e-04 -0.3946082 0.326 0.512 1.000000e+00
## SLC9A3R1    2.146804e-04 -0.4495837 0.083 0.234 1.000000e+00
## IFI16       2.147484e-04  0.3012028 0.636 0.410 1.000000e+00
## GAPDH       2.167050e-04 -0.6435029 0.606 0.634 1.000000e+00
## CYTH4       2.217915e-04  0.2577546 0.629 0.400 1.000000e+00
## APH1A       2.223256e-04  0.3365889 0.273 0.117 1.000000e+00
## NF1         2.234741e-04  0.2914041 0.348 0.161 1.000000e+00
## ATRX        2.278246e-04  0.2925664 0.614 0.390 1.000000e+00
## PHF14       2.280682e-04  0.2832804 0.205 0.068 1.000000e+00
## TMX2        2.296570e-04  0.2758590 0.121 0.024 1.000000e+00
## SEM1        2.305242e-04  0.3379112 0.250 0.107 1.000000e+00
## ROCK2       2.308060e-04  0.3709522 0.311 0.151 1.000000e+00
## ACADM       2.362058e-04  0.2605535 0.152 0.039 1.000000e+00
## TM7SF3      2.364368e-04  0.3786928 0.220 0.083 1.000000e+00
## RAB10       2.387229e-04  0.2723702 0.386 0.200 1.000000e+00
## CARD16      2.395382e-04 -0.5387328 0.144 0.312 1.000000e+00
## ZRSR2       2.413891e-04  0.2685509 0.265 0.107 1.000000e+00
## SRGAP1      2.417753e-04  0.3240821 0.273 0.112 1.000000e+00
## ATXN3       2.418675e-04  0.3365013 0.242 0.098 1.000000e+00
## CXCL2       2.439672e-04 -1.0674301 0.061 0.200 1.000000e+00
## CARD19      2.460107e-04 -0.6439015 0.152 0.307 1.000000e+00
## PPP1R15A    2.473018e-04 -0.3565148 0.492 0.673 1.000000e+00
## TRAM1       2.555841e-04  0.3166523 0.439 0.244 1.000000e+00
## ORAI2       2.600193e-04 -0.6042538 0.045 0.176 1.000000e+00
## TOB1        2.605314e-04 -0.5396541 0.152 0.332 1.000000e+00
## MRPL14      2.617837e-04  0.3003324 0.152 0.039 1.000000e+00
## ZFP36L1     2.618140e-04  0.2781587 0.985 0.834 1.000000e+00
## C12orf75    2.619935e-04 -0.5372580 0.068 0.205 1.000000e+00
## LINC01268   2.672109e-04 -0.5563046 0.129 0.283 1.000000e+00
## TMCC3       2.673326e-04 -0.4002706 0.023 0.141 1.000000e+00
## STK17A      2.682014e-04 -0.5113864 0.061 0.200 1.000000e+00
## ADD3        2.710910e-04 -0.3792683 0.030 0.156 1.000000e+00
## XPO1        2.747748e-04  0.3004819 0.462 0.263 1.000000e+00
## VARS        2.782730e-04  0.2847554 0.121 0.024 1.000000e+00
## PTAR1       2.788862e-04  0.3165003 0.326 0.156 1.000000e+00
## PDLIM7      2.793890e-04 -0.5016869 0.015 0.127 1.000000e+00
## BDP1        2.871582e-04  0.3029448 0.379 0.190 1.000000e+00
## VPS41       2.875986e-04  0.3047041 0.212 0.078 1.000000e+00
## UBE2B       2.972301e-04 -0.5699337 0.318 0.468 1.000000e+00
## PIAS2       2.972626e-04  0.3158442 0.258 0.107 1.000000e+00
## HNRNPA1     2.974259e-04  0.2654661 0.621 0.434 1.000000e+00
## SNX24       2.996625e-04  0.2911334 0.121 0.024 1.000000e+00
## YBX1        2.999670e-04  0.2807864 0.697 0.459 1.000000e+00
## CAPRIN1     3.068148e-04  0.3145963 0.288 0.132 1.000000e+00
## ATG2A       3.107698e-04 -0.6022282 0.068 0.205 1.000000e+00
## NDUFC2      3.208774e-04  0.3608584 0.235 0.102 1.000000e+00
## PLIN2       3.216713e-04 -1.2784574 0.114 0.259 1.000000e+00
## NAV3        3.250028e-04 -0.6711538 0.182 0.332 1.000000e+00
## EPSTI1      3.294326e-04  0.3113634 0.235 0.093 1.000000e+00
## MID1IP1     3.325264e-04 -0.5905903 0.083 0.229 1.000000e+00
## SEMA4D      3.361135e-04 -0.5922136 0.167 0.332 1.000000e+00
## ARL5A       3.369132e-04  0.2592735 0.447 0.234 1.000000e+00
## NFATC2      3.381319e-04  0.2603705 0.258 0.102 1.000000e+00
## HLA-C       3.386781e-04 -0.3890991 0.826 0.859 1.000000e+00
## SLC35B4     3.515031e-04  0.2716761 0.129 0.029 1.000000e+00
## ASPH        3.618720e-04  0.2983339 0.258 0.107 1.000000e+00
## MMD         3.699483e-04  0.2958603 0.318 0.151 1.000000e+00
## RNF149      3.716792e-04 -0.5094628 0.500 0.629 1.000000e+00
## WWC2        3.734340e-04  0.2673024 0.121 0.024 1.000000e+00
## SLC12A5     3.767702e-04  0.2819221 0.174 0.054 1.000000e+00
## ULK1        3.782938e-04 -0.6326716 0.091 0.234 1.000000e+00
## MIR22HG     3.813791e-04 -0.3605515 0.023 0.137 1.000000e+00
## FBXL5       3.934347e-04 -0.5871625 0.265 0.429 1.000000e+00
## PRRG4       3.942764e-04 -0.3647622 0.015 0.122 1.000000e+00
## CCDC14      3.992719e-04  0.2772972 0.227 0.088 1.000000e+00
## CSGALNACT2  4.010177e-04 -0.4498305 0.091 0.249 1.000000e+00
## PDIA3       4.022335e-04  0.3934768 0.455 0.283 1.000000e+00
## P2RY13      4.044430e-04 -0.5992640 0.477 0.580 1.000000e+00
## MIR222HG    4.076074e-04 -0.3740012 0.091 0.244 1.000000e+00
## TM9SF2      4.124537e-04  0.2724870 0.439 0.254 1.000000e+00
## LAMP1       4.144160e-04  0.2876537 0.606 0.400 1.000000e+00
## ZC3H12A     4.213243e-04 -0.4422092 0.053 0.185 1.000000e+00
## AKR1A1      4.220299e-04  0.3356557 0.242 0.102 1.000000e+00
## VAMP2       4.246265e-04 -0.5887844 0.212 0.376 1.000000e+00
## CREBL2      4.324422e-04  0.3661268 0.341 0.185 1.000000e+00
## NINJ1       4.399734e-04  0.3030027 0.682 0.459 1.000000e+00
## SUGT1       4.474810e-04  0.2798039 0.242 0.098 1.000000e+00
## TMED10      4.478760e-04  0.3765653 0.439 0.268 1.000000e+00
## ARPIN       4.540062e-04  0.2619530 0.197 0.068 1.000000e+00
## MESD        4.597511e-04  0.2928341 0.341 0.171 1.000000e+00
## TNRC6B      4.655999e-04  0.2851275 0.576 0.337 1.000000e+00
## CALU        4.735588e-04  0.3080580 0.303 0.146 1.000000e+00
## TIMMDC1     4.736047e-04  0.2606709 0.136 0.034 1.000000e+00
## RARS        4.800137e-04  0.3168152 0.174 0.059 1.000000e+00
## CLEC7A      4.851733e-04  0.3051047 0.818 0.634 1.000000e+00
## TMEM176A    4.872666e-04  0.3411326 0.288 0.137 1.000000e+00
## C16orf72    4.922177e-04 -0.4153242 0.303 0.468 1.000000e+00
## PRMT2       4.937040e-04  0.3058824 0.561 0.361 1.000000e+00
## NUP85       4.988683e-04  0.3569175 0.182 0.063 1.000000e+00
## SPOPL       5.095581e-04 -0.5650663 0.144 0.293 1.000000e+00
## POU2F1      5.098254e-04  0.2783362 0.227 0.088 1.000000e+00
## NFIA        5.114060e-04  0.3413578 0.250 0.107 1.000000e+00
## KLHL5       5.146565e-04  0.2648212 0.303 0.137 1.000000e+00
## FMNL1       5.172259e-04 -0.5847789 0.439 0.546 1.000000e+00
## RBM47       5.179641e-04  0.2749138 0.485 0.302 1.000000e+00
## JAK3        5.209327e-04 -0.5016942 0.038 0.156 1.000000e+00
## CDKN1B      5.220644e-04  0.3269688 0.379 0.200 1.000000e+00
## BRCA2       5.258194e-04  0.3272024 0.106 0.020 1.000000e+00
## NR4A2       5.351537e-04  0.5404741 0.515 0.346 1.000000e+00
## PPP1R15B    5.419679e-04 -0.4946359 0.250 0.420 1.000000e+00
## ROCK1       5.533066e-04  0.2683965 0.500 0.302 1.000000e+00
## IPO7        5.617189e-04  0.3912739 0.326 0.176 1.000000e+00
## NECAP2      5.635421e-04  0.3023822 0.265 0.117 1.000000e+00
## ZBTB44      5.754015e-04  0.2601518 0.242 0.098 1.000000e+00
## SNX2        5.812864e-04  0.3704681 0.326 0.176 1.000000e+00
## TOMM7       5.942754e-04  0.3968276 0.636 0.463 1.000000e+00
## SNX10       5.967195e-04 -0.4971168 0.106 0.254 1.000000e+00
## TCEAL8      6.025952e-04  0.3027178 0.212 0.083 1.000000e+00
## SERPINB1    6.048885e-04 -0.5361001 0.205 0.361 1.000000e+00
## RAB5A       6.073278e-04 -0.4337245 0.053 0.176 1.000000e+00
## PUM2        6.187724e-04  0.2887328 0.371 0.185 1.000000e+00
## GNB1        6.197252e-04  0.3055348 0.614 0.439 1.000000e+00
## TCEAL9      6.224103e-04  0.3015750 0.167 0.054 1.000000e+00
## WDFY2       6.324969e-04  0.3207568 0.500 0.298 1.000000e+00
## NACA        6.487531e-04  0.3590252 0.833 0.595 1.000000e+00
## BBC3        6.488731e-04 -0.3494668 0.015 0.117 1.000000e+00
## GIMAP1      6.544613e-04  0.3035901 0.220 0.088 1.000000e+00
## C15orf39    6.709580e-04 -0.4117145 0.076 0.210 1.000000e+00
## UPP1        6.722969e-04 -0.3079617 0.023 0.132 1.000000e+00
## CPSF2       6.776425e-04  0.3107777 0.220 0.088 1.000000e+00
## SOX4        6.850265e-04 -0.3442200 0.174 0.341 1.000000e+00
## UBC         6.882543e-04 -0.4048852 0.894 0.844 1.000000e+00
## NRROS       6.934552e-04  0.2506441 0.174 0.059 1.000000e+00
## MAZ         6.935958e-04  0.2533003 0.189 0.068 1.000000e+00
## EIF3E       6.975264e-04  0.3908874 0.492 0.317 1.000000e+00
## GATAD1      6.980260e-04  0.2669970 0.485 0.278 1.000000e+00
## MTAP        7.027632e-04  0.2982904 0.182 0.063 1.000000e+00
## TXN         7.167192e-04 -0.6166962 0.076 0.210 1.000000e+00
## CXCR4       7.205210e-04 -0.3783761 0.477 0.654 1.000000e+00
## SGPL1       7.214496e-04  0.3148569 0.159 0.049 1.000000e+00
## AC139887.2  7.238144e-04  0.2515666 0.114 0.024 1.000000e+00
## CORO1A      7.275243e-04 -0.5416058 0.379 0.507 1.000000e+00
## TM2D2       7.292043e-04  0.3444935 0.235 0.098 1.000000e+00
## GK          7.380881e-04 -0.5001331 0.083 0.215 1.000000e+00
## ARHGAP9     7.590946e-04 -0.5693091 0.098 0.239 1.000000e+00
## SLCO3A1     7.597269e-04 -0.4300024 0.053 0.176 1.000000e+00
## QPRT        7.683842e-04  0.2776389 0.197 0.073 1.000000e+00
## SRPK1       7.684759e-04 -0.3777552 0.030 0.141 1.000000e+00
## NSA2        7.690491e-04  0.2942104 0.629 0.444 1.000000e+00
## PLEKHA2     7.896107e-04  0.3094188 0.455 0.278 1.000000e+00
## COPS9       7.896750e-04  0.3886844 0.295 0.156 1.000000e+00
## ANKRD17     7.906340e-04  0.3062035 0.258 0.122 1.000000e+00
## RHBDF2      7.932451e-04  0.2532709 0.621 0.405 1.000000e+00
## SHPRH       7.991910e-04  0.3066371 0.356 0.185 1.000000e+00
## SYK         8.082919e-04  0.3158618 0.508 0.317 1.000000e+00
## AC087239.1  8.119223e-04 -0.5856848 0.023 0.127 1.000000e+00
## LCP2        8.210286e-04 -0.4894773 0.348 0.473 1.000000e+00
## KDM4B       8.241442e-04 -0.4424714 0.076 0.205 1.000000e+00
## SAT1        8.288407e-04 -0.3770147 0.970 0.971 1.000000e+00
## FBL         8.311198e-04  0.2869098 0.258 0.122 1.000000e+00
## GALM        8.597345e-04  0.3046329 0.167 0.054 1.000000e+00
## DHRS4L2     8.622048e-04  0.2661651 0.159 0.054 1.000000e+00
## TMIGD3      8.772390e-04 -0.3928726 0.061 0.185 1.000000e+00
## NHLRC2      8.792638e-04  0.2811730 0.258 0.112 1.000000e+00
## TROVE2      8.854272e-04  0.2646956 0.265 0.122 1.000000e+00
## TTC17       8.926929e-04  0.2508152 0.288 0.137 1.000000e+00
## LPP         8.949464e-04  0.3086339 0.470 0.288 1.000000e+00
## IMP4        9.138514e-04  0.2829404 0.121 0.029 1.000000e+00
## RBM7        9.162952e-04 -0.3276124 0.030 0.141 1.000000e+00
## DDAH2       9.220936e-04  0.2934248 0.326 0.171 1.000000e+00
## LINC00891   9.271664e-04  0.3930142 0.159 0.054 1.000000e+00
## PYURF       9.599790e-04  0.2519117 0.227 0.098 1.000000e+00
## SH2D3C      9.629060e-04 -0.2754291 0.015 0.112 1.000000e+00
## DHX34       9.659424e-04 -0.4787685 0.083 0.215 1.000000e+00
## ZNRF2       9.704532e-04  0.2587359 0.333 0.176 1.000000e+00
## FAM126A     9.823174e-04  0.2611614 0.220 0.093 1.000000e+00
## ZMYND8      9.827895e-04  0.2665018 0.265 0.122 1.000000e+00
## TARDBP      9.905890e-04  0.2737718 0.424 0.249 1.000000e+00
## CD300LF     9.966382e-04 -0.3439401 0.030 0.137 1.000000e+00
## SDC3        1.005041e-03  0.2681727 0.144 0.044 1.000000e+00
## RIF1        1.026024e-03  0.3131808 0.318 0.166 1.000000e+00
## CARD8       1.036469e-03  0.2921454 0.439 0.254 1.000000e+00
## FAM219B     1.037475e-03  0.3135669 0.189 0.073 1.000000e+00
## HELZ2       1.048617e-03 -0.4863367 0.023 0.122 1.000000e+00
## FKBP2       1.052749e-03  0.3033378 0.280 0.141 1.000000e+00
## LRRFIP1     1.053473e-03 -0.4415946 0.689 0.702 1.000000e+00
## ARL2        1.053970e-03  0.2600432 0.159 0.054 1.000000e+00
## SRGAP2C     1.081809e-03  0.2551846 0.606 0.376 1.000000e+00
## ABI2        1.083551e-03  0.2889712 0.326 0.166 1.000000e+00
## DYNLL1      1.087006e-03  0.2681187 0.470 0.273 1.000000e+00
## AKAP11      1.089986e-03  0.2947203 0.265 0.122 1.000000e+00
## MKNK2       1.097219e-03 -0.4992864 0.220 0.376 1.000000e+00
## ZBTB20      1.118105e-03  0.3210454 0.386 0.210 1.000000e+00
## SMIM7       1.122609e-03  0.2784678 0.220 0.093 1.000000e+00
## LBR         1.133128e-03 -0.6939796 0.402 0.498 1.000000e+00
## BIRC3       1.141180e-03 -0.4200097 0.030 0.137 1.000000e+00
## FAM8A1      1.190056e-03 -0.4168750 0.023 0.122 1.000000e+00
## PTP4A1      1.191184e-03 -0.5569103 0.356 0.478 1.000000e+00
## USP24       1.191596e-03  0.2540824 0.220 0.093 1.000000e+00
## CHD1        1.204969e-03 -0.4616186 0.303 0.454 1.000000e+00
## PRR13       1.226105e-03 -0.4528993 0.364 0.522 1.000000e+00
## PRDX2       1.241941e-03  0.2895542 0.205 0.083 1.000000e+00
## HIF1A       1.248688e-03 -0.4462099 0.712 0.722 1.000000e+00
## MYDGF       1.260039e-03  0.2765630 0.235 0.107 1.000000e+00
## FAM20C      1.284575e-03  0.2852619 0.106 0.024 1.000000e+00
## MAP1LC3A    1.293922e-03 -0.3478995 0.023 0.122 1.000000e+00
## RIPK2       1.320174e-03 -0.4055164 0.114 0.249 1.000000e+00
## GNAI2       1.320245e-03 -0.4010743 0.621 0.673 1.000000e+00
## NPHP3       1.322276e-03  0.2510238 0.182 0.068 1.000000e+00
## ZNF267      1.332531e-03 -0.3880783 0.091 0.224 1.000000e+00
## PPM1K       1.337688e-03  0.2636000 0.303 0.151 1.000000e+00
## SNRNP200    1.345811e-03  0.2534552 0.265 0.127 1.000000e+00
## TMED9       1.348507e-03  0.2749583 0.295 0.146 1.000000e+00
## ZBTB37      1.352710e-03  0.2545735 0.303 0.151 1.000000e+00
## JAML        1.358989e-03 -0.4896275 0.114 0.254 1.000000e+00
## MAT2A       1.367087e-03  0.3991215 0.568 0.385 1.000000e+00
## BRD7        1.371212e-03  0.3062571 0.371 0.215 1.000000e+00
## SCAMP1      1.393814e-03  0.2662833 0.212 0.093 1.000000e+00
## TBCA        1.397421e-03  0.3042698 0.462 0.302 1.000000e+00
## USP36       1.399884e-03  0.4297943 0.311 0.176 1.000000e+00
## SCO1        1.419344e-03  0.2593675 0.174 0.063 1.000000e+00
## SPNS2       1.420139e-03 -0.3753368 0.030 0.132 1.000000e+00
## ELF4        1.425631e-03 -0.2713276 0.015 0.107 1.000000e+00
## TMEM123     1.444593e-03  0.2696749 0.432 0.259 1.000000e+00
## COMMD6      1.458393e-03  0.2901373 0.553 0.385 1.000000e+00
## ATL3        1.463118e-03  0.2969313 0.182 0.068 1.000000e+00
## VOPP1       1.463640e-03  0.2963402 0.280 0.141 1.000000e+00
## DIAPH1      1.485829e-03 -0.3672524 0.083 0.215 1.000000e+00
## ELAVL1      1.497494e-03  0.2709685 0.258 0.127 1.000000e+00
## AKR1B1      1.503574e-03  0.3171213 0.288 0.146 1.000000e+00
## MSMO1       1.508165e-03  0.3733124 0.212 0.093 1.000000e+00
## EMILIN2     1.512632e-03 -0.2603446 0.030 0.132 1.000000e+00
## UNC119      1.524576e-03 -0.3187328 0.030 0.132 1.000000e+00
## PDXDC1      1.547061e-03  0.3286682 0.295 0.151 1.000000e+00
## TRIM8       1.590576e-03  0.2706108 0.295 0.151 1.000000e+00
## LARP4       1.594064e-03  0.2746862 0.227 0.102 1.000000e+00
## NUBP1       1.603521e-03  0.2938141 0.167 0.063 1.000000e+00
## RALGDS      1.606496e-03 -0.5047709 0.061 0.176 1.000000e+00
## ACER3       1.614263e-03  0.2793849 0.227 0.102 1.000000e+00
## LRRFIP2     1.619403e-03 -0.4135834 0.083 0.210 1.000000e+00
## CRISPLD2    1.627506e-03 -0.3100548 0.023 0.117 1.000000e+00
## CISD3       1.636047e-03  0.2804724 0.152 0.054 1.000000e+00
## ABCA1       1.676319e-03  0.2651177 0.273 0.137 1.000000e+00
## HEXB        1.688324e-03  0.2853165 0.303 0.161 1.000000e+00
## HSPA6       1.724813e-03 -0.4271280 0.439 0.244 1.000000e+00
## ATP6V0D1    1.735335e-03 -0.4574897 0.182 0.327 1.000000e+00
## RYBP        1.742282e-03 -0.4857991 0.129 0.259 1.000000e+00
## PLEKHA1     1.751093e-03  0.2948568 0.220 0.098 1.000000e+00
## AC008105.3  1.784568e-03 -0.3248023 0.076 0.200 1.000000e+00
## SREBF1      1.815433e-03  0.2747877 0.182 0.073 1.000000e+00
## PDCD4       1.817425e-03  0.2702016 0.364 0.200 1.000000e+00
## PLP2        1.827838e-03 -0.3766301 0.023 0.117 1.000000e+00
## HMGN2       1.834305e-03 -0.5219914 0.439 0.551 1.000000e+00
## BACE1       1.846197e-03  0.2524354 0.114 0.029 1.000000e+00
## TAOK3       1.869843e-03  0.3435383 0.545 0.395 1.000000e+00
## APOBR       1.877877e-03 -0.3780169 0.045 0.156 1.000000e+00
## CCDC50      1.893145e-03  0.2600148 0.371 0.205 1.000000e+00
## EMC10       1.970325e-03  0.2939623 0.273 0.141 1.000000e+00
## RASGRP3     1.976416e-03  0.2803219 0.136 0.044 1.000000e+00
## ATP5F1C     1.992976e-03  0.2959001 0.227 0.107 1.000000e+00
## CTSO        2.016818e-03  0.3145640 0.174 0.068 1.000000e+00
## AC023157.3  2.017158e-03 -0.5375013 0.152 0.288 1.000000e+00
## CCDC82      2.017576e-03  0.2770577 0.227 0.107 1.000000e+00
## SCML1       2.018719e-03  0.2691715 0.197 0.083 1.000000e+00
## B4GALT5     2.033060e-03 -0.5262284 0.068 0.180 1.000000e+00
## CFL1        2.036226e-03 -0.4429117 0.659 0.717 1.000000e+00
## ZNF217      2.047489e-03  0.2818494 0.394 0.254 1.000000e+00
## ATP6V0A1    2.062558e-03  0.2659074 0.485 0.317 1.000000e+00
## PPP1R21     2.089603e-03  0.2857878 0.174 0.068 1.000000e+00
## KDM5A       2.106520e-03  0.2642640 0.394 0.229 1.000000e+00
## FAM133B     2.110957e-03  0.2728791 0.720 0.551 1.000000e+00
## EIF3M       2.114431e-03  0.3067644 0.265 0.141 1.000000e+00
## NDEL1       2.131464e-03 -0.3600127 0.083 0.205 1.000000e+00
## LINC01094   2.131756e-03  0.2537766 0.265 0.127 1.000000e+00
## C9orf72     2.143654e-03 -0.3895143 0.227 0.371 1.000000e+00
## LYRM2       2.164462e-03  0.3996123 0.197 0.088 1.000000e+00
## TTYH3       2.171402e-03  0.2619014 0.197 0.083 1.000000e+00
## INAFM2      2.172759e-03  0.3192373 0.288 0.156 1.000000e+00
## ATP5PD      2.208237e-03  0.2844209 0.288 0.156 1.000000e+00
## ABTB1       2.220312e-03 -0.4841983 0.091 0.210 1.000000e+00
## SOCS3       2.225891e-03 -0.4407244 0.053 0.166 1.000000e+00
## PSMB6       2.245927e-03  0.2971781 0.220 0.102 1.000000e+00
## OSM         2.250373e-03 -0.4211994 0.045 0.151 1.000000e+00
## EXOSC6      2.274454e-03  0.2536338 0.159 0.059 1.000000e+00
## RFLNB       2.323602e-03 -0.4777656 0.045 0.146 1.000000e+00
## MARCKSL1    2.330844e-03 -0.3444690 0.045 0.151 1.000000e+00
## SGK3        2.332222e-03  0.2987271 0.250 0.122 1.000000e+00
## GPR183      2.343441e-03  0.5943587 0.470 0.337 1.000000e+00
## ZNF106      2.343442e-03  0.2543175 0.492 0.327 1.000000e+00
## CCL2        2.374027e-03  0.6041761 0.220 0.102 1.000000e+00
## THAP9-AS1   2.419452e-03 -0.5684007 0.129 0.254 1.000000e+00
## APBB3       2.421969e-03  0.3116253 0.136 0.044 1.000000e+00
## CD276       2.460273e-03  0.2953690 0.197 0.083 1.000000e+00
## TMEM160     2.474325e-03  0.2745193 0.205 0.088 1.000000e+00
## FKBP1A      2.496252e-03 -0.4508536 0.341 0.468 1.000000e+00
## SH3BP2      2.510981e-03 -0.4595910 0.250 0.376 1.000000e+00
## AKAP9       2.523346e-03  0.2942254 0.568 0.410 1.000000e+00
## B4GALT1     2.587814e-03  0.2524812 0.386 0.234 1.000000e+00
## MRPS5       2.657139e-03  0.2970471 0.273 0.141 1.000000e+00
## UBN2        2.673458e-03  0.3155899 0.265 0.137 1.000000e+00
## LNPEP       2.697300e-03  0.2794179 0.424 0.268 1.000000e+00
## NFE2L1      2.722032e-03  0.2689732 0.242 0.117 1.000000e+00
## NME2        2.733522e-03  0.3104891 0.303 0.171 1.000000e+00
## KIF13A      2.750912e-03 -0.3258928 0.030 0.122 1.000000e+00
## TLE3        2.817552e-03 -0.5147879 0.258 0.380 1.000000e+00
## ATP5MC2     2.830086e-03  0.2862129 0.636 0.478 1.000000e+00
## UBE2D3      2.832720e-03 -0.3551165 0.689 0.732 1.000000e+00
## HLX         2.833900e-03 -0.3775567 0.038 0.137 1.000000e+00
## BHLHE40     2.850854e-03 -0.3579020 0.076 0.195 1.000000e+00
## STX11       2.875114e-03 -0.4548785 0.076 0.185 1.000000e+00
## TMEM87A     2.878820e-03  0.2923958 0.182 0.078 1.000000e+00
## NUFIP2      2.898018e-03  0.3133216 0.667 0.478 1.000000e+00
## MED6        2.906633e-03  0.2869197 0.189 0.083 1.000000e+00
## NBEAL2      2.906633e-03 -0.4439085 0.061 0.166 1.000000e+00
## TDP2        2.975263e-03 -0.5060683 0.098 0.220 1.000000e+00
## DEK         3.066358e-03  0.2738966 0.523 0.380 1.000000e+00
## SOCS6       3.099630e-03 -0.5896553 0.303 0.410 1.000000e+00
## GNLY        3.110477e-03 -0.3767643 0.030 0.122 1.000000e+00
## COX6A1      3.112163e-03  0.2593527 0.424 0.259 1.000000e+00
## GBP2        3.125630e-03 -0.3832492 0.167 0.302 1.000000e+00
## DENND1A     3.193685e-03  0.2533733 0.144 0.054 1.000000e+00
## HSPA2       3.282697e-03  0.4478128 0.174 0.073 1.000000e+00
## CD300LB     3.351783e-03 -0.3230166 0.053 0.161 1.000000e+00
## PAQR3       3.390795e-03  0.2704547 0.205 0.093 1.000000e+00
## ATP6AP2     3.394313e-03  0.3038865 0.432 0.268 1.000000e+00
## SUMO3       3.401872e-03  0.2660644 0.265 0.137 1.000000e+00
## HNRNPAB     3.419079e-03  0.2508502 0.311 0.180 1.000000e+00
## PGGT1B      3.454157e-03  0.2619072 0.205 0.093 1.000000e+00
## C5AR2       3.560647e-03  0.2533768 0.356 0.200 1.000000e+00
## HNRNPA2B1   3.659714e-03 -0.2731567 0.886 0.878 1.000000e+00
## RNASEH2B    3.673986e-03  0.2646558 0.258 0.132 1.000000e+00
## VENTX       3.682937e-03  0.2702273 0.318 0.180 1.000000e+00
## KRAS        3.722424e-03 -0.3877633 0.174 0.307 1.000000e+00
## SHOC2       3.734186e-03 -0.4604553 0.152 0.273 1.000000e+00
## IPCEF1      3.793593e-03 -0.5425755 0.159 0.273 1.000000e+00
## AC007032.1  3.803467e-03 -0.3348704 0.023 0.107 1.000000e+00
## USF3        3.806943e-03  0.2725793 0.167 0.068 1.000000e+00
## RAB8B       3.885242e-03 -0.4593937 0.235 0.366 1.000000e+00
## ISG20       3.889239e-03 -0.4987904 0.045 0.141 1.000000e+00
## VSIR        3.894803e-03 -0.4593845 0.409 0.512 1.000000e+00
## LTB         3.900577e-03 -0.2882312 0.038 0.132 1.000000e+00
## FAM111A     3.924951e-03  0.2751656 0.311 0.171 1.000000e+00
## CDK16       3.927411e-03 -0.2732346 0.023 0.107 1.000000e+00
## FGD4        3.954565e-03  0.2913032 0.492 0.337 1.000000e+00
## PNPLA8      3.955109e-03 -0.4859623 0.114 0.234 1.000000e+00
## GBP4        3.957330e-03  0.2610483 0.174 0.073 1.000000e+00
## FAM53C      3.999637e-03 -0.4219801 0.288 0.420 1.000000e+00
## RIT1        4.082604e-03 -0.4595481 0.129 0.249 1.000000e+00
## BLOC1S6     4.116226e-03  0.3089146 0.371 0.234 1.000000e+00
## SALL1       4.140597e-03 -0.4279849 0.038 0.127 1.000000e+00
## ZC3HAV1     4.171245e-03  0.2633194 0.447 0.288 1.000000e+00
## MGST2       4.171827e-03  0.2775844 0.341 0.205 1.000000e+00
## JMJD1C      4.203992e-03 -0.4248955 0.614 0.649 1.000000e+00
## RAB5IF      4.234589e-03 -0.4494726 0.121 0.239 1.000000e+00
## KDM3A       4.298333e-03 -0.3794933 0.106 0.229 1.000000e+00
## TFRC        4.302336e-03  0.6788475 0.280 0.166 1.000000e+00
## AGTRAP      4.407159e-03 -0.3114702 0.083 0.195 1.000000e+00
## TMEM52B     4.438841e-03 -0.2704345 0.023 0.102 1.000000e+00
## METTL7A     4.482327e-03 -0.4470874 0.227 0.346 1.000000e+00
## SPN         4.510804e-03 -0.3093385 0.106 0.229 1.000000e+00
## RAB5C       4.520962e-03  0.3071524 0.318 0.180 1.000000e+00
## DDX3X       4.526500e-03 -0.4100514 0.659 0.698 1.000000e+00
## SLC43A2     4.613428e-03 -0.3369341 0.182 0.307 1.000000e+00
## MGAT1       4.651195e-03  0.2678835 0.462 0.288 1.000000e+00
## THEMIS2     4.665445e-03  0.3410706 0.402 0.263 1.000000e+00
## HSP90B1     4.704369e-03  0.2661559 0.674 0.537 1.000000e+00
## STAT1       4.712811e-03  0.3009880 0.273 0.151 1.000000e+00
## LMNB1       4.733774e-03 -0.3784835 0.023 0.102 1.000000e+00
## NDUFA3      4.747553e-03  0.2959060 0.379 0.249 1.000000e+00
## HSPA1B      4.768151e-03 -0.7464085 0.856 0.800 1.000000e+00
## FAM49B      4.850139e-03 -0.4200047 0.553 0.624 1.000000e+00
## MOB3A       4.851312e-03 -0.4139247 0.136 0.254 1.000000e+00
## QSOX1       4.869011e-03 -0.4425731 0.076 0.180 1.000000e+00
## WDR11       4.911941e-03  0.2809489 0.341 0.205 1.000000e+00
## SIPA1L1     5.060697e-03 -0.4881033 0.098 0.210 1.000000e+00
## ZFAS1       5.103745e-03 -0.5712344 0.470 0.517 1.000000e+00
## TMEM183A    5.130763e-03 -0.2852505 0.038 0.127 1.000000e+00
## RNF19B      5.131569e-03 -0.4076718 0.068 0.171 1.000000e+00
## TENT5A      5.134167e-03  0.2804852 0.455 0.298 1.000000e+00
## LINC00996   5.209508e-03 -0.3441535 0.023 0.102 1.000000e+00
## HSPA9       5.324606e-03  0.3175825 0.364 0.239 1.000000e+00
## CACNB4      5.350038e-03 -0.3911512 0.083 0.190 1.000000e+00
## MKLN1       5.433111e-03  0.2911172 0.326 0.200 1.000000e+00
## DNAJB1      5.506515e-03 -0.5144219 0.636 0.727 1.000000e+00
## TMEM176B    5.554341e-03  0.3174501 0.303 0.176 1.000000e+00
## NLRP3       5.626387e-03 -0.3955493 0.205 0.341 1.000000e+00
## AZIN1       5.651777e-03 -0.3757838 0.227 0.356 1.000000e+00
## CUL3        5.735941e-03 -0.3328129 0.098 0.210 1.000000e+00
## COLGALT1    5.743751e-03  0.2541324 0.364 0.224 1.000000e+00
## UBE2S       5.861360e-03 -0.4361710 0.174 0.302 1.000000e+00
## MAVS        5.960193e-03  0.2501844 0.205 0.098 1.000000e+00
## RSRC1       5.969878e-03  0.3131417 0.326 0.205 1.000000e+00
## UQCRB       6.008862e-03  0.2641077 0.629 0.473 1.000000e+00
## CCNG2       6.082377e-03 -0.3493527 0.076 0.180 1.000000e+00
## PSMC4       6.165369e-03  0.2768752 0.212 0.107 1.000000e+00
## ETF1        6.353689e-03 -0.4486815 0.159 0.273 1.000000e+00
## TOM1        6.476686e-03 -0.3335540 0.053 0.146 1.000000e+00
## NOTCH2NL    6.511481e-03  0.3036034 0.227 0.122 1.000000e+00
## ZC3H18      6.532580e-03  0.2655193 0.152 0.063 1.000000e+00
## RTTN        6.678688e-03 -0.4269973 0.121 0.229 1.000000e+00
## SIN3B       6.715800e-03  0.2819989 0.189 0.093 1.000000e+00
## SENP7       6.762674e-03  0.2834755 0.280 0.151 1.000000e+00
## PTPN11      6.893257e-03  0.2777528 0.182 0.088 1.000000e+00
## QARS        6.917324e-03  0.3137845 0.129 0.049 1.000000e+00
## MTRNR2L8    7.035041e-03 -0.3112282 0.053 0.146 1.000000e+00
## CAND1       7.058921e-03  0.2633378 0.212 0.112 1.000000e+00
## CLEC4E      7.111235e-03 -0.5114982 0.106 0.210 1.000000e+00
## SLC25A45    7.119038e-03  0.2785659 0.159 0.068 1.000000e+00
## PELI1       7.256914e-03 -0.4784658 0.424 0.507 1.000000e+00
## JMJD6       7.304091e-03 -0.3743180 0.136 0.254 1.000000e+00
## SLC6A6      7.349405e-03 -0.4148603 0.129 0.244 1.000000e+00
## SGK1        7.428780e-03  0.2620075 0.765 0.649 1.000000e+00
## GPS2        7.736612e-03  0.2562082 0.121 0.044 1.000000e+00
## WASF2       7.773823e-03  0.2981314 0.417 0.273 1.000000e+00
## TPM3        8.058961e-03 -0.3030831 0.727 0.746 1.000000e+00
## PANK3       8.251560e-03  0.2665939 0.235 0.127 1.000000e+00
## UBE2J1      8.343366e-03 -0.3462210 0.114 0.224 1.000000e+00
## FRAT2       8.408066e-03 -0.3918611 0.053 0.141 1.000000e+00
## HSPA5       8.720718e-03 -0.8441687 0.303 0.405 1.000000e+00
## CCNH        8.727810e-03  0.2916851 0.341 0.210 1.000000e+00
## DUSP2       8.740314e-03 -0.4033833 0.045 0.132 1.000000e+00
## MTX1        8.778556e-03 -0.4048403 0.030 0.107 1.000000e+00
## MCFD2       8.845282e-03  0.3283910 0.121 0.044 1.000000e+00
## TFDP1       8.906868e-03 -0.3395183 0.030 0.107 1.000000e+00
## ATXN1       8.961054e-03  0.2577495 0.333 0.200 1.000000e+00
## MAPK6       9.065029e-03 -0.2919893 0.053 0.141 1.000000e+00
## ERO1A       9.168344e-03 -0.4419667 0.083 0.180 1.000000e+00
## LIMD2       9.182757e-03 -0.3414937 0.258 0.385 1.000000e+00
## IL10        9.323953e-03  0.2780242 0.144 0.059 1.000000e+00
## LRRC58      9.572480e-03  0.2846931 0.227 0.122 1.000000e+00
## UXT         9.636133e-03  0.2505589 0.220 0.117 1.000000e+00
## CERS2       9.809037e-03  0.2771358 0.197 0.098 1.000000e+00
## COQ7        9.854527e-03 -0.4880245 0.053 0.137 1.000000e+00
## LAMTOR5     9.882823e-03 -0.4164939 0.152 0.263 1.000000e+00
## SH2B3       9.883582e-03  0.3106249 0.280 0.176 1.000000e+00
## SERINC1     9.993106e-03  0.2718411 0.523 0.376 1.000000e+00
## MECP2       1.013069e-02 -0.3916272 0.227 0.351 1.000000e+00
## FAM208A     1.013989e-02  0.2806951 0.273 0.161 1.000000e+00
## LAMP2       1.081358e-02 -0.3029715 0.561 0.629 1.000000e+00
## PFN1        1.081909e-02 -0.5335334 0.674 0.683 1.000000e+00
## LSM12       1.082666e-02 -0.3255823 0.068 0.161 1.000000e+00
## ICAM1       1.094869e-02 -0.4464050 0.098 0.200 1.000000e+00
## PLXNC1      1.107703e-02 -0.3099897 0.068 0.161 1.000000e+00
## RASSF3      1.119262e-02 -0.3713868 0.159 0.268 1.000000e+00
## IMPAD1      1.126646e-02  0.2875327 0.189 0.098 1.000000e+00
## ITGA5       1.140447e-02 -0.2591284 0.091 0.195 1.000000e+00
## ITGAL       1.148465e-02 -0.4287615 0.121 0.215 1.000000e+00
## MTF2        1.187871e-02  0.3029560 0.235 0.132 1.000000e+00
## GK5         1.198108e-02 -0.2537083 0.038 0.117 1.000000e+00
## HIST1H2AC   1.210317e-02 -0.2786992 0.492 0.620 1.000000e+00
## KLF5        1.217413e-02 -0.3491925 0.030 0.102 1.000000e+00
## OXSR1       1.221682e-02 -0.3267775 0.083 0.180 1.000000e+00
## IFRD1       1.226985e-02 -0.3043236 0.515 0.629 1.000000e+00
## OSBPL3      1.231467e-02  0.2521055 0.212 0.112 1.000000e+00
## VPS26A      1.240541e-02  0.2555828 0.280 0.176 1.000000e+00
## ZBTB1       1.249038e-02  0.2661064 0.197 0.102 1.000000e+00
## SERTAD1     1.250811e-02 -0.3208910 0.106 0.205 1.000000e+00
## VPS28       1.257228e-02  0.2868869 0.326 0.205 1.000000e+00
## EVI2A       1.261021e-02 -0.4106251 0.417 0.488 1.000000e+00
## AGO3        1.263333e-02  0.2711491 0.227 0.127 1.000000e+00
## USP3        1.269535e-02 -0.3601244 0.220 0.332 1.000000e+00
## LHFPL2      1.272098e-02 -0.4495627 0.227 0.327 1.000000e+00
## FAM214B     1.281782e-02 -0.2955715 0.076 0.171 1.000000e+00
## USP15       1.307337e-02 -0.4438827 0.364 0.459 1.000000e+00
## SIPA1L2     1.314707e-02 -0.3414007 0.076 0.166 1.000000e+00
## NUTM2B-AS1  1.326356e-02  0.2760254 0.295 0.185 1.000000e+00
## LINC01736   1.341396e-02 -0.4559433 0.106 0.195 1.000000e+00
## PPDPF       1.344472e-02 -0.2720064 0.091 0.190 1.000000e+00
## NABP1       1.359474e-02 -0.3768656 0.424 0.498 1.000000e+00
## ARHGDIA     1.367697e-02 -0.4593839 0.318 0.405 1.000000e+00
## CCAR1       1.378240e-02  0.2819038 0.288 0.176 1.000000e+00
## ZNF10       1.390321e-02  0.2801215 0.106 0.039 1.000000e+00
## CWC25       1.393057e-02 -0.5317110 0.152 0.249 1.000000e+00
## LPAR2       1.395320e-02 -0.4290915 0.098 0.195 1.000000e+00
## NARF        1.426920e-02 -0.2590693 0.076 0.166 1.000000e+00
## FOSL2       1.452474e-02 -0.4685701 0.288 0.390 1.000000e+00
## CCL4        1.477014e-02 -0.5490221 0.083 0.180 1.000000e+00
## SEC22B      1.486724e-02 -0.4850766 0.189 0.293 1.000000e+00
## ANXA5       1.493597e-02  0.2517931 0.500 0.371 1.000000e+00
## SRSF5       1.494929e-02 -0.3185170 0.485 0.571 1.000000e+00
## PLCG2       1.499366e-02 -0.4577068 0.167 0.263 1.000000e+00
## CEBPZOS     1.517889e-02  0.3356903 0.167 0.083 1.000000e+00
## TUBB4B      1.525949e-02 -0.3262877 0.280 0.410 1.000000e+00
## ARMCX3      1.536198e-02  0.2533399 0.311 0.200 1.000000e+00
## RNF135      1.539190e-02  0.2802244 0.182 0.098 1.000000e+00
## RXRA        1.553194e-02 -0.4293963 0.189 0.288 1.000000e+00
## AC007952.4  1.559758e-02 -0.4219531 0.114 0.210 1.000000e+00
## OSBPL8      1.582399e-02 -0.4780210 0.258 0.351 1.000000e+00
## PMAIP1      1.639066e-02 -0.3394243 0.189 0.307 1.000000e+00
## SRPK2       1.643809e-02 -0.4677738 0.114 0.210 1.000000e+00
## ZNF562      1.691509e-02  0.2611905 0.174 0.088 1.000000e+00
## USP9X       1.702542e-02 -0.2743270 0.098 0.195 1.000000e+00
## CITED2      1.740950e-02  0.4334183 0.462 0.346 1.000000e+00
## FTL         1.749225e-02 -0.2921220 0.992 0.990 1.000000e+00
## ARID5B      1.778979e-02  0.2527825 0.152 0.073 1.000000e+00
## OAT         1.789180e-02 -0.3612223 0.121 0.215 1.000000e+00
## GOLGA7      1.803597e-02 -0.3611735 0.061 0.141 1.000000e+00
## RBM8A       1.806845e-02 -0.4283005 0.280 0.376 1.000000e+00
## IFI30       1.826323e-02 -0.4443772 0.152 0.249 1.000000e+00
## THADA       1.855257e-02  0.3829311 0.144 0.068 1.000000e+00
## TTF2        1.878492e-02  0.2823383 0.136 0.063 1.000000e+00
## CALM1       1.894261e-02 -0.3881800 0.500 0.566 1.000000e+00
## CCNDBP1     1.905276e-02 -0.3202647 0.076 0.161 1.000000e+00
## OLFML3      1.915474e-02 -0.4967629 0.083 0.166 1.000000e+00
## KIAA0232    1.918172e-02 -0.3222404 0.091 0.180 1.000000e+00
## XPO6        1.929559e-02 -0.5143099 0.220 0.312 1.000000e+00
## TTN-AS1     1.998447e-02 -0.3616816 0.098 0.185 1.000000e+00
## WDR1        2.018591e-02  0.2839205 0.288 0.180 1.000000e+00
## CBX6        2.027461e-02 -0.2616620 0.038 0.107 1.000000e+00
## IDI1        2.057129e-02  0.7415533 0.394 0.307 1.000000e+00
## ATP11B      2.080812e-02 -0.3571896 0.121 0.210 1.000000e+00
## TNFAIP2     2.094475e-02 -0.3830753 0.386 0.478 1.000000e+00
## CLEC12A     2.126043e-02 -0.3920210 0.068 0.146 1.000000e+00
## GABARAPL1   2.197961e-02 -0.3795536 0.295 0.385 1.000000e+00
## NRBF2       2.230799e-02 -0.5048265 0.152 0.239 1.000000e+00
## HSP90AA1    2.281059e-02 -0.4380864 0.970 0.961 1.000000e+00
## TOMM5       2.297374e-02  0.2974895 0.114 0.049 1.000000e+00
## WARS        2.317069e-02 -0.3208201 0.098 0.185 1.000000e+00
## MYL12B      2.336525e-02 -0.3559359 0.455 0.527 1.000000e+00
## TOB2        2.338289e-02  0.2594885 0.295 0.200 1.000000e+00
## ARPC4       2.341055e-02 -0.5047978 0.250 0.341 1.000000e+00
## LILRA2      2.424071e-02 -0.4010527 0.106 0.190 1.000000e+00
## MTMR3       2.435765e-02 -0.3581105 0.083 0.166 1.000000e+00
## ING3        2.465237e-02 -0.2769471 0.076 0.156 1.000000e+00
## CLK3        2.481639e-02 -0.2820151 0.098 0.180 1.000000e+00
## EXOC6B      2.498456e-02  0.2722962 0.106 0.044 1.000000e+00
## MRPL18      2.512945e-02 -0.3186186 0.106 0.195 1.000000e+00
## SLC36A4     2.513131e-02 -0.3149758 0.098 0.185 1.000000e+00
## AAMP        2.569933e-02  0.2960581 0.136 0.068 1.000000e+00
## MCF2L       2.604322e-02 -0.4210294 0.144 0.224 1.000000e+00
## LAMTOR4     2.609610e-02 -0.4312873 0.348 0.434 1.000000e+00
## C1orf43     2.694872e-02  0.2522190 0.250 0.156 1.000000e+00
## AL360012.1  2.773149e-02 -0.2605348 0.083 0.166 1.000000e+00
## WTAP        2.825732e-02 -0.3386294 0.258 0.346 1.000000e+00
## TRNAU1AP    2.834099e-02  0.2799732 0.205 0.117 1.000000e+00
## KLHL15      2.844160e-02 -0.3298143 0.068 0.141 1.000000e+00
## AC004687.1  2.849302e-02 -0.2908522 0.159 0.254 1.000000e+00
## CCL3L1      2.852152e-02 -0.4735816 0.053 0.127 1.000000e+00
## WDR45       2.872012e-02 -0.2803168 0.038 0.102 1.000000e+00
## NUP58       2.941139e-02 -0.4385816 0.212 0.302 1.000000e+00
## SLC44A2     2.945043e-02 -0.2911830 0.038 0.102 1.000000e+00
## NXF1        2.968125e-02 -0.3426103 0.189 0.293 1.000000e+00
## PAIP2       2.971274e-02 -0.3701930 0.265 0.366 1.000000e+00
## PSTPIP1     3.040477e-02 -0.4000109 0.083 0.161 1.000000e+00
## AC079015.1  3.046240e-02 -0.2948704 0.068 0.141 1.000000e+00
## IL10RA      3.046540e-02  0.2674435 0.424 0.327 1.000000e+00
## ELOB        3.100204e-02 -0.3915687 0.356 0.434 1.000000e+00
## GSTM2       3.122816e-02 -0.2549670 0.053 0.122 1.000000e+00
## IDS         3.125640e-02 -0.3402947 0.295 0.380 1.000000e+00
## TIAM1       3.126138e-02 -0.4081787 0.258 0.341 1.000000e+00
## DLEU2       3.143593e-02 -0.2996475 0.068 0.141 1.000000e+00
## TET3        3.225667e-02 -0.2636253 0.106 0.190 1.000000e+00
## HCLS1       3.229888e-02 -0.3658605 0.470 0.527 1.000000e+00
## NUS1        3.233310e-02  0.2843471 0.114 0.054 1.000000e+00
## ELOVL5      3.283605e-02 -0.3786444 0.144 0.224 1.000000e+00
## ADGRE2      3.300729e-02 -0.2729888 0.061 0.132 1.000000e+00
## ZC3H13      3.377744e-02  0.2924171 0.364 0.268 1.000000e+00
## OSER1       3.443203e-02 -0.4155888 0.098 0.176 1.000000e+00
## SAMSN1      3.677818e-02 -0.3758208 0.402 0.463 1.000000e+00
## ADAMTSL2    3.684458e-02 -0.2835950 0.061 0.127 1.000000e+00
## SRSF4       3.716405e-02 -0.3120862 0.205 0.302 1.000000e+00
## PHF21A      3.744216e-02 -0.3574609 0.114 0.195 1.000000e+00
## N4BP1       3.771157e-02 -0.3256376 0.114 0.190 1.000000e+00
## CDC42SE1    3.912773e-02 -0.4227392 0.371 0.424 1.000000e+00
## C3orf62     3.925272e-02 -0.3553080 0.045 0.107 1.000000e+00
## HSPA1A      4.028191e-02 -0.4990430 0.917 0.922 1.000000e+00
## TRMT112     4.082339e-02  0.2645043 0.280 0.200 1.000000e+00
## UBN1        4.091687e-02 -0.4437128 0.144 0.220 1.000000e+00
## FAM193B     4.091743e-02 -0.2886132 0.076 0.146 1.000000e+00
## TMEM167B    4.134463e-02 -0.3621542 0.106 0.180 1.000000e+00
## RAP1A       4.162299e-02 -0.3386315 0.205 0.288 1.000000e+00
## VPS37B      4.327532e-02 -0.3387933 0.136 0.210 1.000000e+00
## NSUN6       4.347095e-02 -0.2797835 0.129 0.210 1.000000e+00
## CHMP2A      4.378652e-02 -0.3216315 0.205 0.293 1.000000e+00
## PLBD1       4.384385e-02 -0.6158874 0.114 0.185 1.000000e+00
## SQOR        4.389740e-02 -0.2740276 0.053 0.117 1.000000e+00
## PPP2R2A     4.392590e-02 -0.2924182 0.159 0.244 1.000000e+00
## EIF4A3      4.547882e-02 -0.3387316 0.273 0.371 1.000000e+00
## IL1RAP      4.553896e-02 -0.3654076 0.106 0.180 1.000000e+00
## HIST1H2BF   4.557471e-02 -0.2777063 0.045 0.107 1.000000e+00
## IQSEC1      4.677805e-02 -0.4211375 0.144 0.220 1.000000e+00
## P2RY12      4.761999e-02 -0.6437484 0.295 0.337 1.000000e+00
## HTRA1       4.778120e-02 -0.4020905 0.250 0.327 1.000000e+00
## PPP1R18     4.877336e-02 -0.3880698 0.424 0.488 1.000000e+00
## OSBPL2      4.889271e-02 -0.3309557 0.091 0.161 1.000000e+00
## CHML        4.914535e-02  0.2547803 0.167 0.098 1.000000e+00
## SF1         4.965109e-02 -0.3853563 0.553 0.561 1.000000e+00
## ANKRD13A    5.004336e-02 -0.3908773 0.114 0.190 1.000000e+00
## AC025164.1  5.028186e-02 -0.2791323 0.167 0.249 1.000000e+00
## EZR         5.033330e-02 -0.5028523 0.197 0.268 1.000000e+00
## CIR1        5.120454e-02 -0.4999890 0.250 0.322 1.000000e+00
## SNHG9       5.253234e-02 -0.3246380 0.129 0.205 1.000000e+00
## LACTB       5.272314e-02 -0.2611874 0.106 0.180 1.000000e+00
## SHKBP1      5.283543e-02 -0.3043201 0.106 0.180 1.000000e+00
## RAB21       5.362700e-02 -0.3440555 0.242 0.327 1.000000e+00
## FLOT2       5.426902e-02 -0.4640448 0.136 0.205 1.000000e+00
## CSF2RB      5.450057e-02 -0.3039583 0.144 0.224 1.000000e+00
## BOD1L1      5.496900e-02 -0.4007887 0.523 0.541 1.000000e+00
## ID3         5.514638e-02  0.3801806 0.159 0.088 1.000000e+00
## SMIM27      5.535148e-02 -0.2666362 0.045 0.102 1.000000e+00
## CDK12       5.617411e-02 -0.3651508 0.235 0.312 1.000000e+00
## AREL1       5.638004e-02 -0.3705821 0.068 0.127 1.000000e+00
## CNPPD1      5.711235e-02 -0.2704767 0.083 0.151 1.000000e+00
## RPS6KB2     5.819732e-02 -0.2887664 0.045 0.102 1.000000e+00
## RBP7        5.924011e-02 -0.2548066 0.053 0.112 1.000000e+00
## ACTR3       6.077543e-02 -0.3689684 0.568 0.571 1.000000e+00
## H2AFY       6.162699e-02 -0.3348404 0.553 0.571 1.000000e+00
## CD53        6.306790e-02 -0.2541345 0.568 0.580 1.000000e+00
## PNPLA2      6.307175e-02 -0.4414136 0.197 0.263 1.000000e+00
## BID         6.336411e-02 -0.4093432 0.197 0.259 1.000000e+00
## IQGAP1      6.399575e-02 -0.2828339 0.591 0.580 1.000000e+00
## PXK         6.437868e-02 -0.2533527 0.083 0.151 1.000000e+00
## GPSM3       6.451250e-02 -0.3859351 0.417 0.468 1.000000e+00
## MPPE1       6.452671e-02 -0.2646001 0.098 0.166 1.000000e+00
## TMCO3       6.467293e-02 -0.3124027 0.068 0.127 1.000000e+00
## MAN2A2      6.530192e-02 -0.3142968 0.068 0.127 1.000000e+00
## PSTPIP2     6.593632e-02 -0.3210220 0.061 0.117 1.000000e+00
## PDXK        6.610520e-02 -0.2506573 0.159 0.239 1.000000e+00
## ERGIC1      6.648089e-02 -0.3333795 0.220 0.298 1.000000e+00
## MRC2        6.715230e-02 -0.4141906 0.136 0.200 1.000000e+00
## CSNK1D      6.770426e-02 -0.3187441 0.220 0.293 1.000000e+00
## TTC37       6.850484e-02  0.2574806 0.159 0.098 1.000000e+00
## ZNF394      6.869084e-02 -0.3686194 0.091 0.156 1.000000e+00
## GPR160      6.877924e-02 -0.3628293 0.076 0.137 1.000000e+00
## ASAP1       6.902836e-02 -0.2744530 0.182 0.263 1.000000e+00
## VEGFA       6.936989e-02 -0.2540904 0.348 0.429 1.000000e+00
## ZDHHC20     7.051537e-02 -0.2650482 0.098 0.166 1.000000e+00
## IST1        7.063067e-02 -0.2606650 0.258 0.346 1.000000e+00
## SIRT1       7.104071e-02 -0.2603235 0.106 0.176 1.000000e+00
## NADK        7.104991e-02 -0.3197813 0.129 0.200 1.000000e+00
## TFE3        7.130732e-02 -0.2657375 0.083 0.146 1.000000e+00
## SNX18       7.140860e-02 -0.2755129 0.235 0.322 1.000000e+00
## CPNE3       7.277295e-02 -0.5980930 0.144 0.210 1.000000e+00
## HIST2H2BE   7.287497e-02 -0.2858020 0.333 0.420 1.000000e+00
## SIGLEC8     7.321532e-02 -0.3470155 0.106 0.166 1.000000e+00
## RBM23       7.389337e-02 -0.3424680 0.152 0.220 1.000000e+00
## PRKCB       7.692469e-02 -0.3127668 0.265 0.332 1.000000e+00
## SRSF2       7.702653e-02 -0.2676901 0.492 0.556 1.000000e+00
## ALPK1       7.728609e-02 -0.3305036 0.136 0.205 1.000000e+00
## HIST3H2A    7.759440e-02  0.4326761 0.235 0.166 1.000000e+00
## PGAM1       7.780494e-02 -0.3862703 0.167 0.229 1.000000e+00
## EIF1AY      7.811285e-02 -0.4113751 0.114 0.176 1.000000e+00
## BACH1       8.112068e-02 -0.3776718 0.326 0.380 1.000000e+00
## CLIC1       8.274674e-02 -0.3223009 0.424 0.463 1.000000e+00
## SF3B1       8.393546e-02 -0.3061861 0.515 0.541 1.000000e+00
## AL034397.3  8.428777e-02 -0.3200068 0.197 0.268 1.000000e+00
## MAP4K4      8.500015e-02 -0.3206155 0.455 0.483 1.000000e+00
## MLF2        8.594973e-02 -0.3540537 0.197 0.268 1.000000e+00
## PPP4R1      8.803650e-02 -0.3112523 0.083 0.141 1.000000e+00
## EPB41       8.854987e-02 -0.3151616 0.121 0.185 1.000000e+00
## ATP5MPL     8.935014e-02 -0.3899369 0.417 0.439 1.000000e+00
## PPP2CA      8.992860e-02 -0.2953638 0.242 0.302 1.000000e+00
## MAP3K8      9.033882e-02 -0.2962566 0.159 0.229 1.000000e+00
## CCND3       9.072012e-02 -0.3119889 0.076 0.132 1.000000e+00
## NUP98       9.135292e-02 -0.2604045 0.152 0.220 1.000000e+00
## CFD         9.210421e-02 -0.2901249 0.258 0.332 1.000000e+00
## TSPO        9.231673e-02 -0.4940108 0.371 0.429 1.000000e+00
## CPVL        9.317010e-02  0.2515698 0.242 0.176 1.000000e+00
## SLC15A4     9.556434e-02 -0.3161091 0.174 0.239 1.000000e+00
## CITED4      9.557591e-02 -0.2699642 0.053 0.102 1.000000e+00
## SLC12A6     9.637017e-02 -0.3976296 0.121 0.176 1.000000e+00
## CAT         9.704904e-02 -0.2675877 0.129 0.195 1.000000e+00
## CHMP4B      9.813303e-02 -0.3297858 0.242 0.312 1.000000e+00
## ATF3        1.050029e-01  0.3678663 0.424 0.337 1.000000e+00
## RASSF5      1.055619e-01 -0.3025270 0.182 0.244 1.000000e+00
## RELT        1.071728e-01 -0.3294181 0.182 0.239 1.000000e+00
## KIAA1551    1.075588e-01 -0.4182502 0.379 0.415 1.000000e+00
## GYG1        1.077451e-01 -0.3593294 0.083 0.137 1.000000e+00
## RTF2        1.079670e-01 -0.2716056 0.182 0.244 1.000000e+00
## CBX4        1.089329e-01 -0.2816563 0.106 0.166 1.000000e+00
## MYH9        1.094271e-01 -0.2574924 0.273 0.332 1.000000e+00
## EHBP1L1     1.097843e-01 -0.3264916 0.220 0.278 1.000000e+00
## UBE2H       1.101869e-01 -0.3099955 0.106 0.166 1.000000e+00
## MAPK14      1.124078e-01 -0.3054808 0.152 0.215 1.000000e+00
## LSM8        1.133549e-01  0.3786068 0.341 0.263 1.000000e+00
## PTBP3       1.136166e-01 -0.2910767 0.470 0.527 1.000000e+00
## SH3KBP1     1.151796e-01 -0.3316896 0.333 0.371 1.000000e+00
## DDX21       1.153567e-01 -0.3912264 0.379 0.395 1.000000e+00
## ACAA1       1.155127e-01 -0.2707160 0.220 0.283 1.000000e+00
## KIF1B       1.164087e-01 -0.2646412 0.129 0.185 1.000000e+00
## ZYX         1.177358e-01 -0.3062965 0.273 0.341 1.000000e+00
## RRAGC       1.181595e-01 -0.2842765 0.114 0.171 1.000000e+00
## PRMT9       1.206989e-01  0.2731492 0.174 0.117 1.000000e+00
## SLU7        1.219615e-01 -0.2864882 0.205 0.268 1.000000e+00
## GPAT4       1.232226e-01 -0.2507789 0.106 0.161 1.000000e+00
## DENND5A     1.233759e-01 -0.3545729 0.242 0.302 1.000000e+00
## PPP2R5A     1.237399e-01 -0.2702069 0.076 0.127 1.000000e+00
## PHC2        1.273649e-01 -0.3525125 0.114 0.166 1.000000e+00
## H2AFZ       1.273828e-01 -0.4321661 0.258 0.312 1.000000e+00
## WIPF1       1.279873e-01 -0.3482009 0.462 0.498 1.000000e+00
## TSC22D2     1.280973e-01 -0.2599940 0.250 0.317 1.000000e+00
## PET100      1.290684e-01 -0.2544763 0.242 0.307 1.000000e+00
## CCNI        1.317546e-01 -0.3149381 0.492 0.498 1.000000e+00
## MAPKAPK2    1.342087e-01 -0.3182021 0.114 0.171 1.000000e+00
## SLC20A1     1.352502e-01 -0.2502475 0.348 0.405 1.000000e+00
## PGGHG       1.365035e-01 -0.3451706 0.121 0.171 1.000000e+00
## TAP1        1.369596e-01 -0.3033109 0.114 0.166 1.000000e+00
## ATP6V0B     1.408593e-01 -0.3069516 0.545 0.546 1.000000e+00
## CUTA        1.414551e-01 -0.2660171 0.114 0.166 1.000000e+00
## PPP2R2D     1.422766e-01 -0.2504929 0.083 0.132 1.000000e+00
## CHURC1      1.462707e-01 -0.3424881 0.250 0.302 1.000000e+00
## BUD31       1.468692e-01 -0.2976805 0.083 0.132 1.000000e+00
## FAM177A1    1.505919e-01 -0.2581394 0.136 0.190 1.000000e+00
## AKNA        1.517971e-01 -0.3043486 0.189 0.244 1.000000e+00
## AKIRIN1     1.537303e-01 -0.2687162 0.121 0.171 1.000000e+00
## VNN1        1.546781e-01 -0.5357588 0.106 0.156 1.000000e+00
## ZFAND2A     1.578377e-01 -0.6337076 0.295 0.346 1.000000e+00
## CCL3        1.583171e-01 -0.5198595 0.212 0.278 1.000000e+00
## STK38       1.584943e-01 -0.2615499 0.091 0.137 1.000000e+00
## TAX1BP1     1.602097e-01 -0.3638061 0.333 0.385 1.000000e+00
## PPP4R2      1.604946e-01 -0.3671640 0.189 0.239 1.000000e+00
## PIP4P1      1.610916e-01 -0.3542293 0.083 0.127 1.000000e+00
## ARHGEF2     1.625822e-01 -0.2792371 0.182 0.239 1.000000e+00
## DCAF11      1.631975e-01 -0.2916702 0.098 0.141 1.000000e+00
## RYR1        1.642858e-01 -0.3262271 0.197 0.244 1.000000e+00
## ATP6V1A     1.652599e-01 -0.2568903 0.189 0.239 1.000000e+00
## LPGAT1      1.653105e-01 -0.3381866 0.152 0.205 1.000000e+00
## MAPRE1      1.666349e-01 -0.2852835 0.182 0.234 1.000000e+00
## MAPK1       1.680942e-01 -0.2694768 0.258 0.317 1.000000e+00
## MBD6        1.682927e-01 -0.2857474 0.174 0.224 1.000000e+00
## TPM4        1.685221e-01 -0.2524581 0.220 0.273 1.000000e+00
## RAP1B       1.720654e-01 -0.2594622 0.356 0.390 1.000000e+00
## MT-ND4L     1.747990e-01 -0.2750901 0.576 0.566 1.000000e+00
## BAZ2B       1.753708e-01 -0.3595235 0.424 0.454 1.000000e+00
## C5AR1       1.762407e-01 -0.2878721 0.712 0.649 1.000000e+00
## MARF1       1.766988e-01 -0.3399787 0.265 0.307 1.000000e+00
## PSMB3       1.775604e-01 -0.3006273 0.227 0.278 1.000000e+00
## PACS2       1.796070e-01 -0.2632863 0.114 0.161 1.000000e+00
## UNC13D      1.821249e-01 -0.3176266 0.068 0.107 1.000000e+00
## MAP3K2      1.822935e-01 -0.3645137 0.432 0.449 1.000000e+00
## GUK1        1.824193e-01 -0.3208777 0.295 0.346 1.000000e+00
## SVIP        1.858961e-01 -0.3136422 0.068 0.107 1.000000e+00
## NLRP1       1.885500e-01 -0.2780061 0.167 0.210 1.000000e+00
## SP3         1.896984e-01 -0.3850688 0.250 0.283 1.000000e+00
## ATG16L2     2.004912e-01 -0.3487791 0.409 0.444 1.000000e+00
## USP32       2.034210e-01 -0.3181627 0.136 0.180 1.000000e+00
## NCF4        2.054959e-01 -0.3154213 0.333 0.361 1.000000e+00
## LIMS1       2.077319e-01 -0.2752912 0.235 0.283 1.000000e+00
## CCDC59      2.101081e-01 -0.2624824 0.114 0.156 1.000000e+00
## H2AFV       2.107227e-01 -0.3036824 0.235 0.278 1.000000e+00
## PGK1        2.116417e-01 -0.2950139 0.402 0.424 1.000000e+00
## CMIP        2.130183e-01 -0.2954443 0.356 0.410 1.000000e+00
## MIR29B2CHG  2.153302e-01 -0.2611277 0.258 0.302 1.000000e+00
## WDR74       2.189943e-01 -0.2527180 0.197 0.254 1.000000e+00
## TSC22D4     2.191817e-01 -0.2777791 0.098 0.141 1.000000e+00
## CX3CR1      2.198719e-01 -0.4246451 0.470 0.439 1.000000e+00
## PICALM      2.217263e-01 -0.2795041 0.402 0.420 1.000000e+00
## ATP6V1B2    2.279629e-01 -0.2939126 0.439 0.454 1.000000e+00
## SCLT1       2.307216e-01 -0.2652364 0.129 0.176 1.000000e+00
## POLB        2.307997e-01 -0.2767969 0.076 0.112 1.000000e+00
## VAPA        2.325191e-01 -0.2569405 0.386 0.420 1.000000e+00
## CAPG        2.356223e-01 -0.2713195 0.212 0.263 1.000000e+00
## GLRX        2.407759e-01 -0.3075459 0.212 0.259 1.000000e+00
## COX7B       2.414587e-01 -0.2911363 0.288 0.332 1.000000e+00
## TNF         2.458909e-01 -0.3371163 0.189 0.234 1.000000e+00
## HECA        2.485170e-01 -0.2969674 0.295 0.327 1.000000e+00
## NECAP1      2.497417e-01 -0.2858569 0.227 0.268 1.000000e+00
## GNG5        2.578895e-01 -0.2559228 0.379 0.400 1.000000e+00
## C3          2.585769e-01 -0.3465137 0.689 0.634 1.000000e+00
## TUBB2A      2.597426e-01  0.2544354 0.182 0.137 1.000000e+00
## SYF2        2.617161e-01 -0.3164683 0.318 0.346 1.000000e+00
## ZNF493      2.685011e-01 -0.3337624 0.091 0.127 1.000000e+00
## ZFAND2B     2.758226e-01 -0.2778415 0.106 0.141 1.000000e+00
## GAB2        2.768110e-01 -0.2760316 0.144 0.185 1.000000e+00
## CEP63       2.831026e-01 -0.2596836 0.144 0.180 1.000000e+00
## MSL1        2.846558e-01 -0.3368018 0.371 0.371 1.000000e+00
## ARHGAP26    2.886862e-01 -0.2576155 0.288 0.317 1.000000e+00
## HIST1H2BD   2.925408e-01 -0.2960730 0.220 0.249 1.000000e+00
## GMIP        2.954316e-01 -0.3054579 0.114 0.146 1.000000e+00
## WDR26       2.963466e-01 -0.3104688 0.311 0.332 1.000000e+00
## AC016745.2  2.981055e-01 -0.2740962 0.091 0.122 1.000000e+00
## ATP2B1      2.999347e-01 -0.3492514 0.417 0.410 1.000000e+00
## CHERP       3.028480e-01 -0.2536875 0.076 0.107 1.000000e+00
## ITGB2-AS1   3.098167e-01 -0.2969010 0.303 0.332 1.000000e+00
## BNIP2       3.158274e-01 -0.2653368 0.235 0.263 1.000000e+00
## TMEM119     3.160228e-01 -0.3407120 0.114 0.141 1.000000e+00
## AGTPBP1     3.191978e-01 -0.2537183 0.121 0.156 1.000000e+00
## AL031777.3  3.245235e-01 -0.3638461 0.258 0.288 1.000000e+00
## FAM160A2    3.245670e-01 -0.2620875 0.083 0.112 1.000000e+00
## STAT5B      3.275591e-01 -0.3046044 0.144 0.176 1.000000e+00
## LASP1       3.303183e-01 -0.4092227 0.242 0.254 1.000000e+00
## FCGR3A      3.375042e-01 -0.2825631 0.477 0.483 1.000000e+00
## GABARAPL2   3.513619e-01 -0.4164552 0.265 0.273 1.000000e+00
## ETS2        3.638864e-01 -0.2770568 0.462 0.439 1.000000e+00
## PTPRE       3.684718e-01 -0.2682221 0.621 0.585 1.000000e+00
## ARHGEF40    3.753856e-01 -0.3219838 0.189 0.210 1.000000e+00
## RAB31       3.773669e-01 -0.2752443 0.379 0.390 1.000000e+00
## TSPYL2      3.858464e-01  0.3457267 0.273 0.239 1.000000e+00
## MICAL1      3.903899e-01 -0.3073479 0.121 0.146 1.000000e+00
## HNRNPH2     4.008585e-01 -0.2541984 0.197 0.229 1.000000e+00
## CARF        4.228271e-01 -0.2652241 0.091 0.112 1.000000e+00
## SMNDC1      4.272940e-01 -0.2619864 0.159 0.185 1.000000e+00
## ZBTB10      4.292995e-01  0.2719995 0.144 0.117 1.000000e+00
## AL118516.1  4.385375e-01 -0.3154691 0.189 0.220 1.000000e+00
## ELL         4.582219e-01 -0.2961078 0.212 0.229 1.000000e+00
## ALOX5AP     4.714805e-01 -0.5066917 0.697 0.580 1.000000e+00
## METTL9      4.814222e-01 -0.2678711 0.280 0.288 1.000000e+00
## MLXIPL      5.039943e-01 -0.2704872 0.220 0.234 1.000000e+00
## ZCCHC2      5.235976e-01  0.4290574 0.136 0.117 1.000000e+00
## CCPG1       5.387469e-01 -0.3649669 0.318 0.307 1.000000e+00
## RICTOR      5.410036e-01 -0.2678644 0.432 0.400 1.000000e+00
## RCOR1       5.634421e-01 -0.2882770 0.235 0.244 1.000000e+00
## KATNBL1     5.822500e-01 -0.2548045 0.197 0.215 1.000000e+00
## PAK1        5.840369e-01 -0.2942294 0.273 0.278 1.000000e+00
## LRP10       5.933855e-01 -0.2539550 0.189 0.195 1.000000e+00
## AC108134.2  6.056097e-01 -0.3609390 0.242 0.220 1.000000e+00
## EIF4EBP2    6.127798e-01 -0.2599839 0.205 0.215 1.000000e+00
## GNAI3       6.182829e-01 -0.3121012 0.333 0.322 1.000000e+00
## HIST2H2BF   6.212744e-01 -0.2629042 0.333 0.341 1.000000e+00
## UGCG        8.581874e-01 -0.3902906 0.280 0.244 1.000000e+00
## H1FX        8.956886e-01 -0.6682320 0.394 0.322 1.000000e+00

DotPlot(MG22, feature = c("CD14", "MARCO", "F13A1", "RNASE1", "FCN1", "LYZ", "S100A9", "S100A8"))+ RotatedAxis()

FeaturePlot(MG22, feature = c("CD14", "MARCO", "F13A1", "RNASE1", "FCN1", "LYZ", "S100A9", "S100A8")) 

Annotation

I did literature mining of these different clusters markers and these were the best defining features that I could identify. exAM refers to ex vivo activated microglia which matched closely with the signature described in Daniele Mattei’s paper looking at the effects of enzymatic digestions of brain cells (https://www.mdpi.com/1422-0067/21/21/7944) and the preprint out of Beth Stevens lab (https://www.biorxiv.org/content/10.1101/2020.12.03.408542v1). It could be just activated microglia before post mortem processing as well. But the high Fos and Jun expression in this cluster are notable.

new.cluster.ids <- c("homeostatic", "'exAM'", "CD14- Mono", "CD14+ Monocytes", "T")
names(new.cluster.ids) <- levels(MG22)
MG22 <- RenameIdents(MG22, new.cluster.ids)
DimPlot(MG22, reduction = "umap", label = TRUE, pt.size = 0.5)

The number of cells in each cluster

table(Idents(MG22))

## 
##     homeostatic          'exAM'      CD14- Mono CD14+ Monocytes               T 
##            3120            2482             205             132              51

prop.table(table(Idents(MG22)))

## 
##     homeostatic          'exAM'      CD14- Mono CD14+ Monocytes               T 
##      0.52086811      0.41435726      0.03422371      0.02203673      0.00851419

Look at just the microglia clusters (without typical monocyte markers)

MG22.microglia <- subset(x = MG22, idents = c("homeostatic", "'exAM'"))

head(MG22.microglia$nCount_RNA)

## HIPP_AAACCCAAGGTGCCTC-1 HIPP_AAACCCATCTTCTTCC-1 HIPP_AAACGAAAGAGTTGTA-1 
##                    5124                    2397                    2707 
## HIPP_AAACGAAGTTCAGCGC-1 HIPP_AAACGCTCAAACTCTG-1 HIPP_AAAGAACCAACTCCAA-1 
##                    2615                    5436                    8555

median(MG22.microglia$nCount_RNA)

## [1] 2249

hist(MG22.microglia$nCount_RNA)

Feature Plots

Markers of interest to look at:

#General myeloid cell makers
myeloid_markers <- c("CD14", "AIF1")

#General microglial markers
mic_gen_markers <- c("P2RY12","C1QA" ,"C1QB", "C1QC", "GPR34")

#clusters based on Olah et al. 
mic1 <- c("EBI3", "CTSC", "PLTP", "CTSB", "TMEM176A", "MS4A6A", "VSIG4", "TREM2","MS4A4A", "TMEM176B")
mic2 <- c("BLNK","MEF2C", "CD81", "TLN2", "MAF", "SORL1", "ADGRG1", "PILRA")
mic3 <- c("FOS", "DUSP", "HSPB2", "DDIT4", "HSPA1A","HSPA1B", "JUN", "HSPE1", "RGS1")
mic4 <- c("MX1", "IFIT3", "BST2", "PARP9", "IFIT1", "IFITM3","ISG15","GMFB", "IRF7")
mic5 <- c("KDM6B", "LPL","NFKBID","CD83","ARID5A","TRIB1","C3AR1","ID2","MRC1","CCL2","CCL3","CCL4")
mic6 <- c("PLCG2","AP1B1","PLD4","SELPLG","SERPINE1","SRGAP2","TFRC")
mic7 <- c("PLAT","HLA−DMB","HTRA1","HLA−DMA","ITM2B","LIPA","CD68","HLA−DRB1","CD74")
mic8 <- c("PRDX1","GPX4","AIF1","BIN1","APOC1","FTL")
mic9 <- c("PCNA", "ORC6", "FEN1", "RFX2","MCM5","HMGB1","FAM111A","TK1","MKI67","UBE2C")

#other myeloid cells
mono <- c("FCN2", "VCAN", "LYZ")
Tcell <- c("CD3E","GZMB", "IL7R", "FGFBP2")
Bcell <- c("CCR7", "CD79A")
GFAPmyel <- c("AIF1","C1QA", "GFAP", "MBP", "SNAP25")
eryth <- c("HBA1")

#from katia 
homeostatic = c("CST3", "HEXB", "CTSD", "P2RY12", "C1QA", "C1QB", "CTSS", "CSF1R", "CX3CR1")
DAM_01 = c("APOE", "TYROBP", "CTSB", "LYZ2", "B2M", "FTH1")
DAM_02 = c("TREM2", "AXL", "CST7", "CTSL", "LPL", "CD9", "CSF1", "CCL6", "ITGAX", "CLEC7A", "LILRB4", "TIMP2")

#Patir core genes
patir.genes <- c("CX3CR1", "AIF1", "SCF1R", "TLR1", "C3AR1", "C1QA", "TREM2", "TYROBP", "CAPG", "WA5", "GPR34", "P2RY12", "P2RY13", "TMEM119")

Now we look at how the Olah microglial cluster defining markers distribute on our clusters. Olah mic1 are markers for Housekeeping cells for CNS parenchyma and the most abundant cell type in the Olah clusters

#compare to olah markers 
DotPlot(MG22, features = mic1)

FeaturePlot(MG22, features = mic1)

Olah mic2 markers are also housekeeping cell of the parenchyma. These markers are highly expressed in our clusters 0 thorugh 4 (MEF2C, CD81, MAF, SORL1)

DotPlot(MG22, features = mic2)

FeaturePlot(MG22, features = mic2)

Olah mic3 markers indicate distressed microglia cells. They mention that these are more of these cells in autopsy samples and less in surgery samples.

DotPlot(MG22, features = mic3)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: DUSP, HSPB2

FeaturePlot(MG22, features = mic3)

## Warning in FetchData(object = object, vars = c(dims, "ident", features), : The
## following requested variables were not found: DUSP, HSPB2

Olah mic4 markers are enriched in interferon response genes

DotPlot(MG22, features = mic4)

Olah mic 5 markers are labeled as “anti inflammatory response”

DotPlot(MG22, features = mic5)

FeaturePlot(MG22, features = mic5)

Olah mic 6 markers are also labeled as “anti-inflammatory response”

DotPlot(MG22, features = mic6)

Olah mic 7 markers are implicated perhaps in antigen presentaion.

DotPlot(MG22, features = mic7)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: HLA−DMB, HLA−DMA, HLA−DRB1

Olah mic 8 have upregulation of transcription factors but weren’t given a category title per se.

DotPlot(MG22, features = mic8)

Olah mic 9 for proliferating microglia

DotPlot(MG22, features = mic9)

Olah clusters 10 - 14 are monocytes, T cell, B cells, GFAP+ myeloid cells, erythrocytes.

DotPlot(MG22, features = mono)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: FCN2

DotPlot(MG22, features = Tcell )

DotPlot(MG22, features = Bcell )

DotPlot(MG22, features = GFAPmyel)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: SNAP25

DotPlot(MG22, features = eryth)

Look also at canonical homeostatic and DAM markers

DotPlot(MG22, features = homeostatic)

FeaturePlot(MG22,features = homeostatic)

DotPlot(MG22, features = DAM_01)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: LYZ2

FeaturePlot(MG22,features = DAM_01)

## Warning in FetchData(object = object, vars = c(dims, "ident", features), : The
## following requested variables were not found: LYZ2

DotPlot(MG22, features = DAM_02)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: CCL6

FeaturePlot(MG22,features = DAM_02)

## Warning in FetchData(object = object, vars = c(dims, "ident", features), : The
## following requested variables were not found: CCL6

Look also at Patir core genes (https://doi.org/10.1002/glia.23572)

DotPlot(MG22, features = patir.genes)

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: SCF1R, WA5

FeaturePlot(MG22,features = patir.genes)

## Warning in FetchData(object = object, vars = c(dims, "ident", features), : The
## following requested variables were not found: SCF1R, WA5

Nuclei

Cluster

NUC <- FindNeighbors(NUC, dims = 1:18)

## Computing nearest neighbor graph

## Computing SNN

NUC <- FindClusters(NUC, resolution = 0.5)

## Modularity Optimizer version 1.3.0 by Ludo Waltman and Nees Jan van Eck
## 
## Number of nodes: 3281
## Number of edges: 119430
## 
## Running Louvain algorithm...
## Maximum modularity in 10 random starts: 0.9157
## Number of communities: 15
## Elapsed time: 0 seconds

NUC <- RunUMAP(NUC, dims = 1:18)

## 23:34:43 UMAP embedding parameters a = 0.9922 b = 1.112

## 23:34:43 Read 3281 rows and found 18 numeric columns

## 23:34:43 Using Annoy for neighbor search, n_neighbors = 30

## 23:34:43 Building Annoy index with metric = cosine, n_trees = 50

## 0%   10   20   30   40   50   60   70   80   90   100%

## [----|----|----|----|----|----|----|----|----|----|

## **************************************************|
## 23:34:44 Writing NN index file to temp file /var/folders/b8/csskzk797yqdzbrxw6stvzhh0000gn/T//RtmprdwnBz/filefc952461cb69
## 23:34:44 Searching Annoy index using 1 thread, search_k = 3000
## 23:34:44 Annoy recall = 100%
## 23:34:45 Commencing smooth kNN distance calibration using 1 thread
## 23:34:45 Initializing from normalized Laplacian + noise
## 23:34:45 Commencing optimization for 500 epochs, with 133436 positive edges
## 23:34:49 Optimization finished

DimPlot(NUC, reduction = "umap")

NUC <- RunTSNE(object=NUC,  reduction.use = "pca", 
               dims.use = 1:18, 
               do.fast = T,
               tsne.method = "Rtsne", 
               seed.use = 2019,
               verbose=F) 
DimPlot(NUC, reduction = "tsne")

saveRDS(NUC, file="NUC_final.RDS")

Differential Expression with FindAllMarkers

NUC.markers <- FindAllMarkers(NUC, only.pos = TRUE, min.pct = 0.25, logfc.threshold = 0.25)

## Calculating cluster 0

## Calculating cluster 1

## Calculating cluster 2

## Calculating cluster 3

## Calculating cluster 4

## Calculating cluster 5

## Calculating cluster 6

## Calculating cluster 7

## Calculating cluster 8

## Calculating cluster 9

## Calculating cluster 10

## Calculating cluster 11

## Calculating cluster 12

## Calculating cluster 13

## Calculating cluster 14

top10markers <- NUC.markers %>% group_by(cluster) %>% top_n(n = 10, wt = avg_logFC)
top10markers_table <-data.frame(Cluster0 = top10markers[1:10,7], 
                                Cluster1 = top10markers[11:20,7],
                                Cluster2 = top10markers[21:30,7], 
                                Cluster3 = top10markers[31:40,7],
                                Cluster4 = top10markers[41:50,7], 
                                Cluster5 = top10markers[51:60,7],
                                Cluster6 = top10markers[61:70,7], 
                                Cluster7 = top10markers[71:80,7],
                                Cluster8 = top10markers[81:90,7],
                                Cluster9 = top10markers[91:100,7],
                                Cluster10 = top10markers[101:110,7],
                                Cluster11 = top10markers[111:120,7],
                                Cluster12 = top10markers[121:130,7],
                                Cluster13 = top10markers[131:140,7],
                                Cluster14 = top10markers[141:150,7]
                                )

top10markers_table

##          gene     gene.1     gene.2   gene.3   gene.4     gene.5 gene.6  gene.7
## 1      COL5A2       ENC1    TMEM144     GAD1 DLX6-AS1 AC109466.1  TSHZ2  SLC1A3
## 2  AJ009632.2     ACVR1C       MOBP      MAF    PROX1       RORB   RORB    AQP4
## 3     FAM19A1     ATP2B1       ST18     TAC1   CXCL14       PCP4 GABRG1    GJA1
## 4        MEG3    ONECUT2 AL359091.1     SOX6      VIP      LRRK1 COBLL1 SLC14A1
## 5        ANO3    EPB41L2        CNP  ZNF385D    CALB2     CLSTN2  FOXP2    APOE
## 6      BAIAP3 AC098850.3         TF     GAD2    ERBB4     CAMK2D    DCC  ATP1A2
## 7      R3HDM1      PDZD2       PLP1   SLC6A1    RGS12     CADPS2  TNNT2    GFAP
## 8        ENC1      PTK2B        MBP    ERBB4   THSD7A      STMN2  KCNH5  ADGRV1
## 9      CEP126     HS6ST3     CLDND1 CNTNAP3B     CNR1     TMSB10 COL6A1  SLC1A2
## 10 AC009041.2     TESPA1        SCD    SLIT2      CCK       CHGA   TBR1    CST3
##     gene.8  gene.9 gene.10 gene.11  gene.12 gene.13 gene.14
## 1    SSTR2    VCAN  SEMA3E     KIT    CSF1R   HTR2C  COL1A2
## 2  RASGRP1  PDGFRA   ZFHX3    SV2C     CD74   ITGA8     CFH
## 3   CTNNA2   SMOC1 B3GALT2    GAD2    PTPRC   TSHZ2   CLDN5
## 4      AK5   GPNMB  DIRAS2   NRIP3  HLA-DRA   ROBO3     VIM
## 5     ERC2    BCAN    TLE4  SLC6A1       C3 COL12A1     MGP
## 6      CCK  PTPRZ1    CHN1  CXCL14   ADAM28    TLE4     DCN
## 7   CLSTN2   OLIG1  DYNLL1    DNER     C1QB    CRYM  IFITM3
## 8   ARPP19     TNR     CLU    DGKD HLA-DRB1  TRMT9B    MYL9
## 9     LMO4    EPN2    NEFL   NTNG1      A2M    PCP4   TAGLN
## 10    SYT4 C1orf61    NEFM  PMEPA1     SPP1  COL6A1    APOD

top10.NUC <- NUC.markers %>% group_by(cluster) %>% top_n(n = 10, wt = avg_logFC)
DoHeatmap(NUC, features = top10.NUC$gene) + NoLegend() 

Gene Enrichment with grofiler

##   
## #### Cluster 0   
##  [1] "Biomarker: COL5A2"     "Biomarker: AJ009632.2" "Biomarker: FAM19A1"   
##  [4] "Biomarker: MEG3"       "Biomarker: ANO3"       "Biomarker: BAIAP3"    
##  [7] "Biomarker: R3HDM1"     "Biomarker: ENC1"       "Biomarker: CEP126"    
## [10] "Biomarker: AC009041.2"
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 1   
##  [1] "Biomarker: ENC1"       "Biomarker: ACVR1C"     "Biomarker: ATP2B1"    
##  [4] "Biomarker: ONECUT2"    "Biomarker: EPB41L2"    "Biomarker: AC098850.3"
##  [7] "Biomarker: PDZD2"      "Biomarker: PTK2B"      "Biomarker: HS6ST3"    
## [10] "Biomarker: TESPA1"    
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 2   
##  [1] "Biomarker: TMEM144"    "Biomarker: MOBP"       "Biomarker: ST18"      
##  [4] "Biomarker: AL359091.1" "Biomarker: CNP"        "Biomarker: TF"        
##  [7] "Biomarker: PLP1"       "Biomarker: MBP"        "Biomarker: CLDND1"    
## [10] "Biomarker: SCD"       
## <div id="htmlwidget-47778563dae988bc2b46" style="width:960px;height:500px;" class="datatables html-widget"></div>
## <script type="application/json" data-for="htmlwidget-47778563dae988bc2b46">{"x":{"filter":"none","extensions":["Buttons"],"caption":"<caption><\/caption>","data":[["539","668","5062"],["query_1","query_1","query_1"],[true,true,true],[5.6162015407843e-07,1.53040283841505e-06,5.9045753230014e-07],[50,10,8],[9,9,5],[4,3,3],[0.444444444444444,0.333333333333333,0.6],[0.08,0.3,0.375],["GO:0043209","GO:0019911","WP:WP2276"],["GO:CC","GO:MF","WP"],["myelin sheath","structural constituent of myelin sheath","Glial Cell Differentiation"],[18963,18694,7399],[2131,4645,90],["GO:0110165","GO:0005198","WP:000000"],[0.00285696172379697,0.00778515923901738,0.00300365746681081]],"container":"<table class=\"display\">\n  <thead>\n    <tr>\n      <th> <\/th>\n      <th>query<\/th>\n      <th>significant<\/th>\n      <th>p_value<\/th>\n      <th>term_size<\/th>\n      <th>query_size<\/th>\n      <th>intersection_size<\/th>\n      <th>precision<\/th>\n      <th>recall<\/th>\n      <th>term_id<\/th>\n      <th>source<\/th>\n      <th>term_name<\/th>\n      <th>effective_domain_size<\/th>\n      <th>source_order<\/th>\n      <th>parents<\/th>\n      <th>q_value<\/th>\n    <\/tr>\n  <\/thead>\n<\/table>","options":{"dom":"Bfrtip","buttons":["copy","csv","excel","pdf","print"],"scrollY":400,"scrollX":true,"scrollCollapse":true,"paging":false,"columnDefs":[{"className":"dt-center","targets":"_all"},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false}},"evals":[],"jsHooks":[]}</script>
##   
##   
## #### Cluster 3   
##  [1] "Biomarker: GAD1"     "Biomarker: MAF"      "Biomarker: TAC1"    
##  [4] "Biomarker: SOX6"     "Biomarker: ZNF385D"  "Biomarker: GAD2"    
##  [7] "Biomarker: SLC6A1"   "Biomarker: ERBB4"    "Biomarker: CNTNAP3B"
## [10] "Biomarker: SLIT2"   
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 4   
##  [1] "Biomarker: DLX6-AS1" "Biomarker: PROX1"    "Biomarker: CXCL14"  
##  [4] "Biomarker: VIP"      "Biomarker: CALB2"    "Biomarker: ERBB4"   
##  [7] "Biomarker: RGS12"    "Biomarker: THSD7A"   "Biomarker: CNR1"    
## [10] "Biomarker: CCK"     
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 5   
##  [1] "Biomarker: AC109466.1" "Biomarker: RORB"       "Biomarker: PCP4"      
##  [4] "Biomarker: LRRK1"      "Biomarker: CLSTN2"     "Biomarker: CAMK2D"    
##  [7] "Biomarker: CADPS2"     "Biomarker: STMN2"      "Biomarker: TMSB10"    
## [10] "Biomarker: CHGA"      
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 6   
##  [1] "Biomarker: TSHZ2"  "Biomarker: RORB"   "Biomarker: GABRG1"
##  [4] "Biomarker: COBLL1" "Biomarker: FOXP2"  "Biomarker: DCC"   
##  [7] "Biomarker: TNNT2"  "Biomarker: KCNH5"  "Biomarker: COL6A1"
## [10] "Biomarker: TBR1"  
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 7   
##  [1] "Biomarker: SLC1A3"  "Biomarker: AQP4"    "Biomarker: GJA1"   
##  [4] "Biomarker: SLC14A1" "Biomarker: APOE"    "Biomarker: ATP1A2" 
##  [7] "Biomarker: GFAP"    "Biomarker: ADGRV1"  "Biomarker: SLC1A2" 
## [10] "Biomarker: CST3"   
## <div id="htmlwidget-a1ecca66c43c6869c2bf" style="width:960px;height:500px;" class="datatables html-widget"></div>
## <script type="application/json" data-for="htmlwidget-a1ecca66c43c6869c2bf">{"x":{"filter":"none","extensions":["Buttons"],"caption":"<caption><\/caption>","data":[["6","7","8","9","10"],["query_1","query_1","query_1","query_1","query_1"],[true,true,true,true,true],[7.89357760766421e-07,1.59911262802268e-06,1.59911262802268e-06,1.84612255891947e-06,3.21543141919318e-06],[88,5,5,104,35],[10,10,10,10,10],[5,3,3,5,4],[0.5,0.3,0.3,0.5,0.4],[0.0568181818181818,0.6,0.6,0.0480769230769231,0.114285714285714],["GO:0010232","GO:0070779","GO:0070777","GO:0006835","GO:0140354"],["GO:BP","GO:BP","GO:BP","GO:BP","GO:BP"],["vascular transport","D-aspartate import across plasma membrane","D-aspartate transport","dicarboxylic acid transport","lipid import into cell"],[18092,18092,18092,18092,18092],[4509,19119,19117,2803,23148],[["GO:0003018","GO:0006810"],["GO:0070777","GO:0089718","GO:0140354"],["GO:0006869","GO:0015718","GO:0015740","GO:0042940"],"GO:0046942",["GO:0006869","GO:0098657"]],[0.00534710947143174,0.0108323889422257,0.0108323889422257,0.0125056342141205,0.0217813324336146]],"container":"<table class=\"display\">\n  <thead>\n    <tr>\n      <th> <\/th>\n      <th>query<\/th>\n      <th>significant<\/th>\n      <th>p_value<\/th>\n      <th>term_size<\/th>\n      <th>query_size<\/th>\n      <th>intersection_size<\/th>\n      <th>precision<\/th>\n      <th>recall<\/th>\n      <th>term_id<\/th>\n      <th>source<\/th>\n      <th>term_name<\/th>\n      <th>effective_domain_size<\/th>\n      <th>source_order<\/th>\n      <th>parents<\/th>\n      <th>q_value<\/th>\n    <\/tr>\n  <\/thead>\n<\/table>","options":{"dom":"Bfrtip","buttons":["copy","csv","excel","pdf","print"],"scrollY":400,"scrollX":true,"scrollCollapse":true,"paging":false,"columnDefs":[{"className":"dt-center","targets":"_all"},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false}},"evals":[],"jsHooks":[]}</script>
##   
##   
## #### Cluster 8   
##  [1] "Biomarker: SSTR2"   "Biomarker: RASGRP1" "Biomarker: CTNNA2" 
##  [4] "Biomarker: AK5"     "Biomarker: ERC2"    "Biomarker: CCK"    
##  [7] "Biomarker: CLSTN2"  "Biomarker: ARPP19"  "Biomarker: LMO4"   
## [10] "Biomarker: SYT4"   
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 9   
##  [1] "Biomarker: VCAN"    "Biomarker: PDGFRA"  "Biomarker: SMOC1"  
##  [4] "Biomarker: GPNMB"   "Biomarker: BCAN"    "Biomarker: PTPRZ1" 
##  [7] "Biomarker: OLIG1"   "Biomarker: TNR"     "Biomarker: EPN2"   
## [10] "Biomarker: C1orf61"
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 10   
##  [1] "Biomarker: SEMA3E"  "Biomarker: ZFHX3"   "Biomarker: B3GALT2"
##  [4] "Biomarker: DIRAS2"  "Biomarker: TLE4"    "Biomarker: CHN1"   
##  [7] "Biomarker: DYNLL1"  "Biomarker: CLU"     "Biomarker: NEFL"   
## [10] "Biomarker: NEFM"   
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 11   
##  [1] "Biomarker: KIT"    "Biomarker: SV2C"   "Biomarker: GAD2"  
##  [4] "Biomarker: NRIP3"  "Biomarker: SLC6A1" "Biomarker: CXCL14"
##  [7] "Biomarker: DNER"   "Biomarker: DGKD"   "Biomarker: NTNG1" 
## [10] "Biomarker: PMEPA1"
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 12   
##  [1] "Biomarker: CSF1R"    "Biomarker: CD74"     "Biomarker: PTPRC"   
##  [4] "Biomarker: HLA-DRA"  "Biomarker: C3"       "Biomarker: ADAM28"  
##  [7] "Biomarker: C1QB"     "Biomarker: HLA-DRB1" "Biomarker: A2M"     
## [10] "Biomarker: SPP1"    
## <div id="htmlwidget-0aba3cc1e88f50d7f734" style="width:960px;height:500px;" class="datatables html-widget"></div>
## <script type="application/json" data-for="htmlwidget-0aba3cc1e88f50d7f734">{"x":{"filter":"none","extensions":["Buttons"],"caption":"<caption><\/caption>","data":[["12","1261"],["query_1","query_1"],[true,true],[7.76792991561431e-06,6.12532679836001e-06],[1070,16],[10,10],[8,3],[0.8,0.3],[0.00747663551401869,0.1875],["GO:0002684","GO:0042613"],["GO:BP","GO:CC"],["positive regulation of immune system process","MHC class II protein complex"],[18092,18963],[1338,2067],[["GO:0002376","GO:0002682","GO:0048518"],"GO:0042611"],[0.0480213427383277,0.0378667702674616]],"container":"<table class=\"display\">\n  <thead>\n    <tr>\n      <th> <\/th>\n      <th>query<\/th>\n      <th>significant<\/th>\n      <th>p_value<\/th>\n      <th>term_size<\/th>\n      <th>query_size<\/th>\n      <th>intersection_size<\/th>\n      <th>precision<\/th>\n      <th>recall<\/th>\n      <th>term_id<\/th>\n      <th>source<\/th>\n      <th>term_name<\/th>\n      <th>effective_domain_size<\/th>\n      <th>source_order<\/th>\n      <th>parents<\/th>\n      <th>q_value<\/th>\n    <\/tr>\n  <\/thead>\n<\/table>","options":{"dom":"Bfrtip","buttons":["copy","csv","excel","pdf","print"],"scrollY":400,"scrollX":true,"scrollCollapse":true,"paging":false,"columnDefs":[{"className":"dt-center","targets":"_all"},{"orderable":false,"targets":0}],"order":[],"autoWidth":false,"orderClasses":false}},"evals":[],"jsHooks":[]}</script>
##   
##   
## #### Cluster 13   
##  [1] "Biomarker: HTR2C"   "Biomarker: ITGA8"   "Biomarker: TSHZ2"  
##  [4] "Biomarker: ROBO3"   "Biomarker: COL12A1" "Biomarker: TLE4"   
##  [7] "Biomarker: CRYM"    "Biomarker: TRMT9B"  "Biomarker: PCP4"   
## [10] "Biomarker: COL6A1" 
## [1] "No significant GO terms at p-value < 0.05."
##   
##   
## #### Cluster 14   
##  [1] "Biomarker: COL1A2" "Biomarker: CFH"    "Biomarker: CLDN5" 
##  [4] "Biomarker: VIM"    "Biomarker: MGP"    "Biomarker: DCN"   
##  [7] "Biomarker: IFITM3" "Biomarker: MYL9"   "Biomarker: TAGLN" 
## [10] "Biomarker: APOD"  
## [1] "No significant GO terms at p-value < 0.05."
## 

Feature Plots

Markers based on McKenzie et al 2018 - combined expression across multiple scRNA seq datasets in human and murine. also looked at pubmed searches and markers commonly used in the literature. nuclei markers https://www.nature.com/articles/s41598-018-27293-5#MOESM3

ast <- c("GFAP", "AQP4", "SLC1A2", "GJA1", "GJB6", "SLC4A4", "SLC39A12")
end <- c("APOLD1","FTL1","ABCB1", "RGS5", "PTPRB", "CD34", "VWF", "PALMD")
micro <- c("ITGAM", "CS3CR1", "CCL3", "CCL4", "CSF1R"," P2RY12", "C1QB", "PLEK", "GP183")
neuron <- c("VIP", "RELN", "GAD2", "SYNPR", "SCG2", "CNR1", "NPY")
oli <- c("MBP", "MOG", "MAG", "CNP", "PLP1", "MOBP", "CLDN22", "OPALIN", "UGT8", "ERMN", "BCAS1")
opc <- c("OLIG2", "OLIG1", "Sox10", "GPR17", "PDGFR", "VCAN", "SHC", "PCDH15", "TNR")

DoHeatmap(NUC, features = c(ast, end, micro, neuron, oli, opc))

## Warning in DoHeatmap(NUC, features = c(ast, end, micro, neuron, oli, opc)): The
## following features were omitted as they were not found in the scale.data slot
## for the RNA assay: SHC, PDGFR, Sox10, GP183, P2RY12, CS3CR1, FTL1

DoHeatmap(NUC, features = c(homeostatic, DAM_01, DAM_02))

## Warning in DoHeatmap(NUC, features = c(homeostatic, DAM_01, DAM_02)): The
## following features were omitted as they were not found in the scale.data slot
## for the RNA assay: CCL6, CST7, LYZ2

Neurons

#From Blue Lake 2016

Exc <- c( "SLC17A7", "SATB2")
In <- c("GAD1", "GAD2", "SLC32A1", "SLC6A1")

FeaturePlot(NUC, features = Exc)

FeaturePlot(NUC, features = In)

DotPlot(NUC, features = c(Exc, "RBFOX3"))

DotPlot(NUC, features = In)

These are the homeostatic microglial markers.

DotPlot(NUC, features =patir.genes ) + RotatedAxis()

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: SCF1R, WA5

DotPlot(MG22, features =patir.genes ) + RotatedAxis()

## Warning in FetchData(object = object, vars = features, cells = cells): The
## following requested variables were not found: SCF1R, WA5

Annotation

new.cluster.ids <- c("Ex1","Ex2", "Oligo","In1","In2", "Ex3", "Ex4","Astro", "Ex5", "OPCs", "Ex6", "In3", "MG", "Ex7", "Endo")
names(new.cluster.ids) <- levels(NUC)
NUC <- RenameIdents(NUC, new.cluster.ids)
DimPlot(NUC, reduction = "umap", label = TRUE, pt.size = 0.5)

table(Idents(NUC))

## 
##   Ex1   Ex2 Oligo   In1   In2   Ex3   Ex4 Astro   Ex5  OPCs   Ex6   In3    MG 
##   598   536   381   277   239   228   228   222   167   106    91    72    54 
##   Ex7  Endo 
##    44    38

prop.table(table(Idents(NUC)))

## 
##        Ex1        Ex2      Oligo        In1        In2        Ex3        Ex4 
## 0.18226151 0.16336483 0.11612313 0.08442548 0.07284365 0.06949101 0.06949101 
##      Astro        Ex5       OPCs        Ex6        In3         MG        Ex7 
## 0.06766230 0.05089912 0.03230722 0.02773545 0.02194453 0.01645840 0.01341055 
##       Endo 
## 0.01158183

Compare distribution of genes and transcripts of microglial cluster in the nuclei sample and microglia in the fresh microglia samples.

NUC.microglia <- subset(x = NUC, idents = c("MG"))

head(NUC.microglia$nCount_RNA)

## AAACGAATCTAGAACC-1 AACAACCCATGGTACT-1 AACCCAATCTCCTACG-1 AACGAAACAATCGTCA-1 
##               1553               1246               1019               2949 
## AACTTCTGTAGCACGA-1 AATGGCTCACCTGAAT-1 
##               4906               1934

median(NUC.microglia$nCount_RNA)

## [1] 1640

hist(NUC.microglia$nCount_RNA)

head(NUC.microglia$nFeature_RNA)

## AAACGAATCTAGAACC-1 AACAACCCATGGTACT-1 AACCCAATCTCCTACG-1 AACGAAACAATCGTCA-1 
##               1071                886                720               1857 
## AACTTCTGTAGCACGA-1 AATGGCTCACCTGAAT-1 
##               2648               1324

median(NUC.microglia$nFeature_RNA)

## [1] 1138.5

hist(NUC.microglia$nFeature_RNA)

head(MG22.microglia$nFeature_RNA)

## HIPP_AAACCCAAGGTGCCTC-1 HIPP_AAACCCATCTTCTTCC-1 HIPP_AAACGAAAGAGTTGTA-1 
##                    2060                    1137                    1125 
## HIPP_AAACGAAGTTCAGCGC-1 HIPP_AAACGCTCAAACTCTG-1 HIPP_AAAGAACCAACTCCAA-1 
##                    1151                    2278                    3050

median(MG22.microglia$nFeature_RNA)

## [1] 1228

hist(MG22.microglia$nFeature_RNA)

It appears that fresh microglia are able to pick up better resolution and more microglial cells that are able to subcluster. The micorglia from the nuclei isolation are limited and the isolation seems to bias for neurons and maybe oligos.