Package: sylsnet
Authors: Chunxuan Shao [aut, cre]
Compiled: 2023-10-24
Here is a step-by-step tutorial of the package synlet. This package processes data from cellBasedArrays. In this vignette, we will show a quick tour of using the synlet with dummy/stimulated data, including quality control, data visualization, and most importantly, hits selection.
First, Let’s have a look at the example data.
library(synlet)
data(example_dt)
head(example_dt)
## PLATE MASTER_PLATE WELL_CONTENT_NAME EXPERIMENT_TYPE
## 1: 1031121 P001 AAK1 si3 sample
## 2: 1031121 P001 PLK1 si1 control_positive
## 3: 1031121 P001 lipid only control_negative
## 4: 1031121 P001 scrambled control si1 control_negative
## 5: 1031121 P001 lipid only control_negative
## 6: 1031121 P001 AARSD1 si3 sample
## EXPERIMENT_MODIFICATION ROW_NAME COL_NAME READOUT
## 1: treatment C 10 455
## 2: treatment C 11 814
## 3: treatment C 12 537
## 4: treatment C 13 568
## 5: treatment C 14 566
## 6: treatment C 15 632
The example_dt is a data.table containing 8 columns, all of them are mandatory and please do NOT change the column names. Besides that, users are free to add new columns. However, please not that the new columns may not appear in the final results. Users need to generate example_dt like data.table from the screen results.
In the example_dt:
PLATE, the name of individual plate. Must be character value.
MASTER_PLATE, usually there are plate replicates designed in the RNAi screen experiments. Several individual PLATE having the identical layout are group to a MASTER_PLATE.
WELL_CONTENT_NAME, the name of siRNA used in each cell, in the format of “geneName si(number)”, as usually several siRNAs are available to control for the off-target effect.
EXPERIMENT_TYPE, indicates the cells are sample (interested knock-down), or negative control (control_negative), or positive control (control_positive).
EXPERIMENT_MODIFICATION, there are two conditions in synthetic lethal RNAi screen, named treatment and control.
ROW_NAME and COL_NAME, the position of siRNA reagent in the plate.
READOUT, the output of RNAi screen output for each cell in the plate. Here the data are randomized.
The example_dt contains three MASTER_PLATE, each MASTER_PLATE has three plates for treatment and control.
Z’ factor and Z factor are widely used quality metrics in RNAi experiments.
Z’ factor = 1- 3*(\(\delta_p\) + \(\delta_n\)) / |\(\mu_p\) - \(\mu_n\)|
Z factor = 1- 3*(\(\delta_s\) + \(\delta_n\)) / |\(\mu_s\) - \(\mu_n\)|
In Z’ factor calculation, \(\delta_p\), \(\delta_n\), \(\mu_p\), \(\mu_n\) are the standard deviation and mean of positive control and negative control signal, respectively; while in Z factor formular \(\delta_s\), \(\delta_n\), \(\mu_s\), \(\mu_n\) are standard deviation and mean of samples and negative control signal, respectively. From the definition, we could see that Z’ factor measures the quality of optimization of plates, and Z factor accesses the performance of screen on actual samples. Generally, the plates with value >= 0.5 are considered as excellent assay. See [1] and [2] for more information and discussion.
By default \(\mu\) in the denominator are mean value of signals, zFactor function offers the option to use median instead, which could be more robust in certain conditions.
Here is an example, we specify negative control is “scrambled control si1”, and positive control is “PLK1 si1”. These informations are stored in EXPERIMENT_TYPE column.
res <- zFactor(example_dt, negativeCon = "scrambled control si1", positiveCon = "PLK1 si1")
## (II) Number of plates to calculate Z/Z' factor: 18
head(res)
## zFactor zPrimeFactor
## 1031121 -10.63422 -41.36325
## 1031122 -33.61931 -605.05780
## 1031123 -23.36894 -30.10812
## 1031124 -32.31388 -14.53770
## 1031125 -115.06929 -137.87408
## 1031126 -23.08053 -51.94333
As the READOUT are shuffled data, not surprising the Z and Z’ factor are negative value.
Synlet could plot the screen results in the format of heatmap, in which a dot represents a single cell. All plates are organized together in single figure. This provides a general view of RNAi screen quality.
plateHeatmap(example_dt)
Usually in each plate there are negative and positive control siRNAs, which set the bound of siRNA effect. The scatter plot provide a direct way to examine the bound of controls. It is possible to specify multiple positive / negative controls in the function scatterPlot, and the output of all plates are arranged in a single plot.
scatterPlot(example_dt, controlOnly = FALSE, control_name = c("PLK1 si1", "scrambled control si1", "lipid only"))
## Warning: Removed 2 rows containing missing values (`geom_point()`).
## Knock-down effect It is often intuitive to look at the knock-down effect of a single gene in RNAi screen experiments. The function siRNAPlot plot the normalized and raw signals, the positive and negative control signals, and Z’ factor of plates into a single graph. siRNAPlot provides the option to specify control siRNAs and normalization methods. The following codes show the knock-down effect of “AAK1”.
We need to calculate the Z’ factor based on mean and median of signals. By default, a pdf file named “AAK1.pdf” will be generated in the working directory, the return value contains all subplot and could be plotted separately.
zF_mean <- zFactor(example_dt, negativeCon = "scrambled control si1", positiveCon = "PLK1 si1")
## (II) Number of plates to calculate Z/Z' factor: 18
zF_med <- zFactor(example_dt, negativeCon = "scrambled control si1", positiveCon = "PLK1 si1", useMean = FALSE)
## (II) Number of plates to calculate Z/Z' factor: 18
siRNAPlot("AAK1", example_dt,
controlsiRNA = c("lipid only", "scrambled control si1"),
FILEPATH = ".", zPrimeMed = zF_med, zPrimeMean = zF_mean,
treatment = "treatment", control = "control",
normMethod = c("PLATE", "lipid only", "scrambled control si1"))
## (==) Processing: AAK1
## $full_plot
##
## $p_norm_boxpl
##
## $p_norm_barpl
##
## $p_raw_barpl
##
## $p_cont_boxplot
##
## $p_zprimer_barpl
The main goal of synthetic lethal RNAi screen experiments is to identify interesting genes led to reliable difference in mortality between treatment and control plates, which is a difficult task because of cell heterogeneity, reagent efficiency and intrinsic character of genes. Synlet tries to improve the results of hits selection by employing several algorithms that explore data from different directions, including student’s t-test, median ± k median absolute deviation, rank products and redundant siRNA activity (RSA).
Student’s t-test is commonly used to test whether the mean from two samples are identical, thus it could be a helpful strategy in identifying synthetic lethal genes [3]. Synlet applies t-test to robust B-score value calculated from raw data of each plates, and the BH method is used to correct for the multiple comparisons.
We start with B-score calculation.
bscore_res <- sapply(unique(example_dt$MASTER_PLATE)[1], bScore, example_dt, treatment = "treatment", control = "control", simplify = FALSE)
## (II) Processing PLATE:1031613
## 1: 33731
## Final: 33589
## (II) Processing PLATE:1031614
## 1: 32981
## 2: 32617.25
## Final: 32546.69
## (II) Processing PLATE:1031615
## 1: 37003
## 2: 36509
## Final: 36360.38
## (II) Processing PLATE:1031121
## 1: 34257
## Final: 33965.5
## (II) Processing PLATE:1031122
## 1: 35666
## Final: 35464.75
## (II) Processing PLATE:1031123
## 1: 26858
## Final: 26757
head(bscore_res$P001)
## 1031121 1031122 1031123 1031613 1031614 1031615
## AAK1 si3 -0.379059117 1.4344181 -0.8484694 -0.1088711 -1.14224946 -0.0854374
## AAK1 si2 0.703408672 -0.4497696 -0.2156487 -0.7355607 0.01035855 0.2937760
## AAK1 si1 -2.395497309 -1.2288465 -0.8126849 1.0424412 -0.59587991 -0.1999736
## AANAT si2 -0.004689391 -0.4785758 3.0162568 -1.7154010 2.40476344 1.8867009
## AANAT si1 0.740923801 -0.1256998 -0.5395926 0.6960949 0.15450039 1.8628214
## AANAT si3 -0.353267466 1.8881158 -1.4963572 -0.3769663 0.24491823 -0.1728052
bscore_res is a list containing B-score of plates belonging to the same master plate. The first three plates are treament, and the following three plates are control.
Now let’s apply the student’s t-test to the B-score and combine the results together.
t_res <- tTest(bscore_res$P001, 3, 3)
head(t_res)
## pvalue Treat_Cont p_adj
## AAK1 si3 0.55670002 0.51448255 0.8761181
## AAK1 si2 0.75470866 0.15647216 0.9592663
## AAK1 si1 0.08470045 -1.56120548 0.8761181
## AANAT si2 0.99366232 -0.01435729 0.9936623
## AANAT si1 0.24037426 -0.87926175 0.8761181
## AANAT si3 0.91960090 0.11444817 0.9789864
The columns of bscore.combined are self-explanatory.
We now apply the madSelect to our example data. madSelection is a list containing results for each master plate, and we combine the results together to madSelection.c. In the setting of synthetic lethal screen we are looking for the genes have a stronger knock-down effect, thus “Yes” means the ratio < Median - k*MAD.
mad_res <- sapply((unique(example_dt$MASTER_PLATE)),
madSelect,
example_dt,
control = "control",
treatment = "treatment",
simplify = FALSE)
head(mad_res$P001)
## MASTER_PLATE treat_cont_ratio treat_median control_median Hits
## A2M si1 P001 0.9427558 1.0201578 1.0821018 No
## A2M si2 P001 0.5500522 0.6324181 1.1497418 No
## A2M si3 P001 1.0504126 0.8768822 0.8347979 No
## A4GALT si1 P001 1.3721300 0.9530558 0.6945813 No
## A4GALT si2 P001 2.0811647 1.6950044 0.8144499 No
## A4GALT si3 P001 1.3238316 1.5751037 1.1898067 No
The rank products is a non-parametric statistic method proposed to find consistently up/down-regulated genes between treatment and controls replicates in microarray data, and has been successfully used in the RNAi screen data [5]. It has several advantages over the parametric Student’s t-test, including clear biological meaning, fewer assumptions of the data and improved performance. Detailed information about rank products could be found in [6] and [7].
It is straight forward to use rank product in synlet, as showed in the following codes.
rp_res <- sapply(unique(example_dt$MASTER_PLATE),
rankProdHits,
example_dt,
control = "control",
treatment = "treatment",
simplify = FALSE)
## Rank Product analysis for unpaired case
##
##
## done Rank Product analysis for unpaired case
##
##
## done Rank Product analysis for unpaired case
##
##
## done
head(rp_res)
## $P001
## MASTER_PLATE pvalue_treat_lowerthan_cont FDR_treat_lowerthan_cont
## A2M si1 P001 0.045970011 1.2608917
## A2M si2 P001 0.118764788 1.4251775
## A2M si3 P001 0.402154304 1.0159688
## A4GALT si1 P001 0.203272445 1.1150946
## A4GALT si2 P001 0.931748963 1.0164534
## A4GALT si3 P001 0.447919402 0.9555614
## A4GNT si1 P001 0.247157546 1.0785057
## A4GNT si2 P001 0.973494556 0.9837419
## A4GNT si3 P001 0.889614535 1.0106863
## AACS si1 P001 0.344574528 0.9874375
## AACS si2 P001 0.669885112 0.9818163
## AACS si3 P001 0.957828269 0.9940704
## AADAC si1 P001 0.277001956 1.0428309
## AADAC si2 P001 0.327148757 0.9663471
## AADAC si3 P001 0.321490942 0.9644728
## AADAT si1 P001 0.414486957 0.9947687
## AADAT si2 P001 0.438913232 0.9686361
## AADAT si3 P001 0.916569700 1.0113873
## AAK1 si1 P001 0.072501742 1.1600279
## AAK1 si2 P001 0.712908578 0.9707691
## AAK1 si3 P001 0.538217861 1.0032799
## AANAT si1 P001 0.431862912 0.9871152
## AANAT si2 P001 0.470602246 0.9412045
## AANAT si3 P001 0.692320107 0.9773931
## AARS si1 P001 0.819242483 0.9830910
## AARS si2 P001 0.636278246 1.0013559
## AARS si3 P001 0.987333412 0.9925027
## AARS2 si1 P001 0.840416914 0.9839027
## AARS2 si2 P001 0.317575091 0.9678479
## AARS2 si3 P001 0.582856386 0.9816529
## AARSD1 si1 P001 0.288699996 1.0078254
## AARSD1 si2 P001 0.313000009 0.9692904
## AARSD1 si3 P001 0.532013721 1.0113528
## AASDH si1 P001 0.227961502 1.0675270
## AASDH si2 P001 0.568521390 0.9746081
## AASDH si3 P001 0.415652702 0.9852508
## AASDHPPT si1 P001 0.989858805 0.9898588
## AASDHPPT si2 P001 0.302151055 0.9668834
## AASDHPPT si3 P001 0.444165178 0.9581990
## AASS si1 P001 0.648708512 0.9964163
## AASS si2 P001 0.298891394 0.9726635
## AASS si3 P001 0.719118097 0.9392563
## AATK si1 P001 0.630682575 1.0007525
## AATK si2 P001 0.845375755 0.9837100
## AATK si3 P001 0.644742720 0.9983113
## ABAT si1 P001 0.208590097 1.0824135
## ABAT si2 P001 0.463362880 0.9464433
## ABAT si3 P001 0.151725170 1.3241469
## ABCA1 si1 P001 0.883857036 1.0101223
## ABCA1 si2 P001 0.377680345 1.0213328
## ABCA1 si3 P001 0.406422742 0.9877616
## ABCA10 si1 P001 0.134304921 1.2893272
## ABCA10 si2 P001 0.166806276 1.1438145
## ABCA10 si3 P001 0.473738955 0.9281416
## ABCA12 si1 P001 0.046761181 1.1222683
## ABCA12 si2 P001 0.435387638 0.9834638
## ABCA12 si3 P001 0.251473812 1.0272973
## ABCA13 si1 P001 0.168269897 1.1140628
## ABCA13 si2 P001 0.466153610 0.9421210
## ABCA13 si3 P001 0.544843782 0.9962858
## ABCA2 si1 P001 0.304558067 0.9586090
## ABCA2 si2 P001 0.010640373 1.0214758
## ABCA2 si3 P001 0.718185951 0.9509773
## ABCA3 si1 P001 0.188853691 1.1696745
## ABCA3 si2 P001 0.952226670 0.9990575
## ABCA3 si3 P001 0.209611079 1.0590876
## ABCA4 si1 P001 0.940022519 1.0026907
## ABCA4 si2 P001 0.210157039 1.0346193
## ABCA4 si3 P001 0.362091277 0.9931646
## ABCA5 si1 P001 0.759903677 0.9536046
## ABCA5 si2 P001 0.623362559 1.0057614
## ABCA5 si3 P001 0.713305822 0.9644698
## ABCA6 si1 P001 0.935135992 1.0086860
## ABCA6 si2 P001 0.443836353 0.9683702
## ABCA6 si3 P001 0.539142307 0.9953396
## ABCA7 si1 P001 0.346200055 0.9775060
## ABCA7 si2 P001 0.755114375 0.9538287
## ABCA7 si3 P001 0.691233295 0.9830874
## ABCA8 si1 P001 0.820364328 0.9783227
## ABCA8 si2 P001 0.547657068 0.9919826
## ABCA8 si3 P001 0.459365541 0.9483676
## ABCA9 si1 P001 0.437860200 0.9775484
## ABCA9 si2 P001 0.155704108 1.2997908
## ABCA9 si3 P001 0.964540314 0.9850624
## ABCB1 si1 P001 0.064653157 1.1284915
## ABCB1 si2 P001 0.898967289 1.0093668
## ABCB1 si3 P001 0.250875929 1.0471343
## ABCB10 si1 P001 0.502152130 0.9641321
## ABCB10 si2 P001 0.276605079 1.0621635
## ABCB10 si3 P001 0.384098071 1.0102305
## ABCB11 si1 P001 0.688979907 0.9871951
## ABCB11 si2 P001 0.166459298 1.1837106
## ABCB11 si3 P001 0.959196728 0.9848437
## ABCB4 si1 P001 0.080974938 1.1105134
## ABCB4 si2 P001 0.280557238 1.0359036
## ABCB4 si3 P001 0.932937104 1.0119996
## ABCB5 si1 P001 0.403629868 0.9935504
## ABCB5 si2 P001 0.951370363 1.0036435
## ABCB5 si3 P001 0.338147323 0.9837013
## ABCB6 si1 P001 0.558847380 0.9843917
## ABCB6 si2 P001 0.564931877 0.9860629
## ABCB6 si3 P001 0.699675400 0.9734614
## ABCB7 si1 P001 0.125933797 1.4223111
## ABCB7 si2 P001 0.863044929 0.9922433
## ABCB7 si3 P001 0.001115915 0.2142557
## ABCB8 si1 P001 0.696499397 0.9761159
## ABCB8 si2 P001 0.063355406 1.2164238
## ABCB8 si3 P001 0.937379053 1.0054569
## ABCB9 si1 P001 0.283686630 1.0276950
## ABCB9 si2 P001 0.796529067 0.9866683
## ABCB9 si3 P001 0.476041245 0.9232315
## ABCC1 si1 P001 0.143677413 1.3136221
## ABCC1 si2 P001 0.161137069 1.1899353
## ABCC1 si3 P001 0.449310910 0.9376923
## ABCC10 si1 P001 0.203065877 1.1467250
## ABCC10 si2 P001 0.569217741 0.9671664
## ABCC10 si3 P001 0.226832599 1.0887965
## ABCC11 si1 P001 0.891536804 1.0069122
## ABCC11 si2 P001 0.719036038 0.9455816
## ABCC11 si3 P001 0.549932895 0.9776585
## ABCC12 si1 P001 0.019142805 0.9188547
## ABCC12 si2 P001 0.607788648 0.9973968
## ABCC12 si3 P001 0.418527223 0.9799662
## ABCC2 si1 P001 0.717776960 0.9570359
## ABCC2 si2 P001 0.536361895 1.0096224
## ABCC2 si3 P001 0.156233858 1.2498709
## MARCH1 si1 P001 0.716326800 0.9617814
## MARCH1 si2 P001 0.940710380 0.9978806
## MARCH1 si3 P001 0.033188720 1.0620390
## MARCH2 si1 P001 0.968566323 0.9839404
## MARCH2 si2 P001 0.799767643 0.9780598
## MARCH2 si3 P001 0.657336281 0.9783610
## MARCH3 si1 P001 0.790618746 0.9857065
## MARCH3 si2 P001 0.426647791 0.9869443
## MARCH3 si3 P001 0.906216632 1.0115907
## MARCH4 si1 P001 0.076150864 1.1246897
## MARCH4 si2 P001 0.705623787 0.9677126
## MARCH4 si3 P001 0.684567922 0.9882484
## MARCH5 si1 P001 0.806394482 0.9799224
## MARCH5 si2 P001 0.704549895 0.9731912
## MARCH5 si3 P001 0.847275819 0.9799817
## MARCH6 si1 P001 0.738898471 0.9457900
## MARCH6 si2 P001 0.394227983 1.0092236
## MARCH6 si3 P001 0.831935759 0.9859979
## MARCH8 si1 P001 0.283780801 1.0089984
## MARCH8 si2 P001 0.205174695 1.0942650
## MARCH8 si3 P001 0.680235450 0.9894334
## MARCH9 si1 P001 0.297215890 0.9838871
## MARCH9 si2 P001 0.294293128 0.9913032
## MARCH9 si3 P001 0.642538856 1.0029875
## PDE12 si1 P001 0.392511228 1.0184075
## PDE12 si2 P001 0.229986244 1.0513657
## PDE12 si3 P001 0.289021460 0.9909307
## RBFOX1 si1 P001 0.179928694 1.1515436
## RBFOX1 si2 P001 0.907663699 1.0073493
## RBFOX1 si3 P001 0.926733080 1.0167586
## SEPT1 si1 P001 0.192524102 1.1201402
## SEPT1 si2 P001 0.250781869 1.0700026
## SEPT1 si3 P001 0.954882636 0.9963993
## SEPT10 si1 P001 0.082478052 1.0557191
## SEPT10 si2 P001 0.268236632 1.0510497
## SEPT10 si3 P001 0.725432383 0.9411015
## SEPT11 si1 P001 0.262368964 1.0494759
## SEPT11 si2 P001 0.604653220 1.0008053
## SEPT11 si3 P001 0.382509357 1.0200250
## SEPT12 si1 P001 0.654499690 0.9894798
## SEPT12 si2 P001 0.649801086 0.9901731
## SEPT12 si3 P001 0.817685432 0.9873937
## SEPT14 si1 P001 0.629304845 1.0068878
## SEPT14 si2 P001 0.190819047 1.1449143
## SEPT14 si3 P001 0.448822988 0.9469672
## SEPT2 si1 P001 0.655385220 0.9830778
## SEPT2 si2 P001 0.664175222 0.9809357
## SEPT2 si3 P001 0.472103708 0.9344733
## SEPT3 si1 P001 0.059053972 1.2598181
## SEPT3 si2 P001 0.032289437 1.2399144
## SEPT3 si3 P001 0.126021461 1.3442289
## SEPT4 si1 P001 0.013357283 0.8548661
## SEPT4 si2 P001 0.403575777 1.0063188
## SEPT4 si3 P001 0.728527270 0.9387734
## SEPT5 si1 P001 0.156663153 1.2031730
## SEPT5 si2 P001 0.565122144 0.9775086
## SEPT5 si3 P001 0.797338690 0.9813399
## SEPT6 si1 P001 0.597778929 0.9980309
## SEPT6 si2 P001 0.750862464 0.9547390
## SEPT6 si3 P001 0.834368157 0.9828140
## SEPT7 si1 P001 0.548293026 0.9838529
## SEPT7 si2 P001 0.359253363 0.9996615
## SEPT7 si3 P001 0.958344473 0.9892588
## SEPT9 si1 P001 0.134278434 1.3569189
## SEPT9 si2 P001 0.622568464 1.0129928
## SEPT9 si3 P001 0.242893681 1.0845485
## treat_cont_log2FC
## A2M si1 -0.642752923
## A2M si2 -0.471593565
## A2M si3 -0.230002531
## A4GALT si1 0.210760723
## A4GALT si2 0.806341462
## A4GALT si3 -0.098293743
## A4GNT si1 -0.493250613
## A4GNT si2 1.005688215
## A4GNT si3 0.538373735
## AACS si1 -0.326538033
## AACS si2 0.044541120
## AACS si3 1.101071893
## AADAC si1 -0.386775496
## AADAC si2 -0.301894192
## AADAC si3 -0.047546282
## AADAT si1 -0.281953405
## AADAT si2 -0.237892597
## AADAT si3 0.678590774
## AAK1 si1 -0.704448724
## AAK1 si2 0.155395115
## AAK1 si3 0.018494773
## AANAT si1 -0.210320819
## AANAT si2 -0.092241794
## AANAT si3 0.310013278
## AARS si1 0.399325479
## AARS si2 0.330109831
## AARS si3 1.222128665
## AARS2 si1 0.439350505
## AARS2 si2 -0.220401132
## AARS2 si3 0.366154231
## AARSD1 si1 0.080687495
## AARSD1 si2 -0.292866666
## AARSD1 si3 -0.099566839
## AASDH si1 -0.488771511
## AASDH si2 -0.035187708
## AASDH si3 -0.192514564
## AASDHPPT si1 1.107468954
## AASDHPPT si2 -0.340073035
## AASDHPPT si3 -0.145538893
## AASS si1 0.018303544
## AASS si2 -0.379009240
## AASS si3 0.115531491
## AATK si1 0.230937090
## AATK si2 0.332582178
## AATK si3 0.071720559
## ABAT si1 -0.526149889
## ABAT si2 -0.148450169
## ABAT si3 -0.615228726
## ABCA1 si1 0.445900038
## ABCA1 si2 -0.318467263
## ABCA1 si3 -0.166233923
## ABCA10 si1 -0.482137418
## ABCA10 si2 -0.620878314
## ABCA10 si3 0.270638991
## ABCA12 si1 -0.817903254
## ABCA12 si2 -0.064621811
## ABCA12 si3 -0.461765543
## ABCA13 si1 -0.172482198
## ABCA13 si2 -0.146709800
## ABCA13 si3 0.244694446
## ABCA2 si1 -0.367364021
## ABCA2 si2 -1.115620107
## ABCA2 si3 0.221568591
## ABCA3 si1 -0.531452458
## ABCA3 si2 0.768639095
## ABCA3 si3 -0.044547897
## ABCA4 si1 0.654654294
## ABCA4 si2 -0.396901028
## ABCA4 si3 0.172958015
## ABCA5 si1 0.272237146
## ABCA5 si2 0.117514691
## ABCA5 si3 0.192767370
## ABCA6 si1 0.547037914
## ABCA6 si2 -0.199620500
## ABCA6 si3 -0.048623911
## ABCA7 si1 -0.328832254
## ABCA7 si2 0.296397958
## ABCA7 si3 0.262053855
## ABCA8 si1 0.386022884
## ABCA8 si2 0.073919573
## ABCA8 si3 0.053689888
## ABCA9 si1 -0.134418255
## ABCA9 si2 -0.333186803
## ABCA9 si3 0.996799028
## ABCB1 si1 -0.702966064
## ABCB1 si2 0.451240547
## ABCB1 si3 -0.375187690
## ABCB10 si1 -0.089700755
## ABCB10 si2 0.090081742
## ABCB10 si3 -0.039085283
## ABCB11 si1 0.069130778
## ABCB11 si2 -0.553401541
## ABCB11 si3 1.111462655
## ABCB4 si1 -0.232055905
## ABCB4 si2 -0.430213990
## ABCB4 si3 0.727083667
## ABCB5 si1 -0.209274257
## ABCB5 si2 0.860034059
## ABCB5 si3 -0.253454182
## ABCB6 si1 -0.115116689
## ABCB6 si2 -0.069700738
## ABCB6 si3 0.354165952
## ABCB7 si1 -0.635154070
## ABCB7 si2 0.332469490
## ABCB7 si3 -1.746733374
## ABCB8 si1 0.795241401
## ABCB8 si2 -0.776562822
## ABCB8 si3 0.553172527
## ABCB9 si1 -0.125837597
## ABCB9 si2 0.393606060
## ABCB9 si3 -0.079843701
## ABCC1 si1 -0.655892262
## ABCC1 si2 -0.618615145
## ABCC1 si3 0.111667615
## ABCC10 si1 -0.105748945
## ABCC10 si2 0.405956629
## ABCC10 si3 -0.516073326
## ABCC11 si1 0.666982925
## ABCC11 si2 0.118372055
## ABCC11 si3 -0.069326278
## ABCC12 si1 -0.705370770
## ABCC12 si2 0.226446070
## ABCC12 si3 -0.251067085
## ABCC2 si1 0.269424998
## ABCC2 si2 -0.149498040
## ABCC2 si3 -0.431769121
## MARCH1 si1 0.145177413
## MARCH1 si2 0.692886101
## MARCH1 si3 -0.290463445
## MARCH2 si1 0.774885211
## MARCH2 si2 0.445513970
## MARCH2 si3 0.072672116
## MARCH3 si1 0.565269868
## MARCH3 si2 -0.211169787
## MARCH3 si3 0.743504779
## MARCH4 si1 -0.636347499
## MARCH4 si2 0.196685408
## MARCH4 si3 0.046064546
## MARCH5 si1 0.405883950
## MARCH5 si2 0.100751490
## MARCH5 si3 0.551943221
## MARCH6 si1 0.102292402
## MARCH6 si2 -0.142691220
## MARCH6 si3 0.578722602
## MARCH8 si1 -0.167719423
## MARCH8 si2 -0.120903732
## MARCH8 si3 0.064732870
## MARCH9 si1 0.308909616
## MARCH9 si2 -0.421638514
## MARCH9 si3 0.140926982
## PDE12 si1 -0.297765639
## PDE12 si2 -0.319725388
## PDE12 si3 -0.382383482
## RBFOX1 si1 -0.400956406
## RBFOX1 si2 0.496498011
## RBFOX1 si3 0.974177127
## SEPT1 si1 -0.512433396
## SEPT1 si2 -0.170563030
## SEPT1 si3 0.615587790
## SEPT10 si1 -0.375013554
## SEPT10 si2 -0.438648032
## SEPT10 si3 0.291851777
## SEPT11 si1 -0.464925918
## SEPT11 si2 0.058689055
## SEPT11 si3 -0.205555480
## SEPT12 si1 -0.006655006
## SEPT12 si2 0.004863198
## SEPT12 si3 0.353635539
## SEPT14 si1 -0.001675321
## SEPT14 si2 -0.567674976
## SEPT14 si3 -0.148223434
## SEPT2 si1 0.172286382
## SEPT2 si2 0.318410231
## SEPT2 si3 -0.167537814
## SEPT3 si1 -0.860700922
## SEPT3 si2 -1.085586408
## SEPT3 si3 -0.549254745
## SEPT4 si1 -1.187162297
## SEPT4 si2 -0.130367211
## SEPT4 si3 0.385862747
## SEPT5 si1 -0.124258087
## SEPT5 si2 -0.034968284
## SEPT5 si3 0.192394438
## SEPT6 si1 0.291836886
## SEPT6 si2 0.205058935
## SEPT6 si3 0.334184351
## SEPT7 si1 -0.132248579
## SEPT7 si2 -0.050601315
## SEPT7 si3 0.730347563
## SEPT9 si1 -0.660659591
## SEPT9 si2 0.073371845
## SEPT9 si3 -0.497496944
##
## $P002
## MASTER_PLATE pvalue_treat_lowerthan_cont FDR_treat_lowerthan_cont
## ABCC3 si1 P002 0.494979511 1.0003796
## ABCC3 si2 P002 0.094560208 1.0086422
## ABCC3 si3 P002 0.825365816 0.9966682
## ABCC4 si1 P002 0.386811262 1.0609680
## ABCC4 si2 P002 0.330831075 1.1143784
## ABCC4 si3 P002 0.914996698 0.9925388
## ABCC5 si1 P002 0.184501621 1.1070097
## ABCC5 si2 P002 0.380398403 1.0900969
## ABCC5 si3 P002 0.698096004 0.9855473
## ABCC6 si1 P002 0.611637943 1.0123663
## ABCC6 si2 P002 0.272145993 1.1611562
## ABCC6 si3 P002 0.922681636 0.9841937
## ABCC8 si1 P002 0.454031468 0.9794836
## ABCC8 si2 P002 0.426574182 1.0367373
## ABCC8 si3 P002 0.377775532 1.1158908
## ABCC9 si1 P002 0.327576723 1.1231202
## ABCC9 si2 P002 0.963578215 0.9893423
## ABCC9 si3 P002 0.958305232 0.9892183
## ABCD1 si1 P002 0.011162317 0.5357912
## ABCD1 si2 P002 0.596474210 1.0225272
## ABCD1 si3 P002 0.747338378 1.0034194
## ABCD2 si1 P002 0.003449163 0.3311196
## ABCD2 si2 P002 0.278201315 1.1364820
## ABCD2 si3 P002 0.195524454 1.1375968
## ABCD3 si1 P002 0.291866365 1.0573272
## ABCD3 si2 P002 0.600939824 1.0121092
## ABCD3 si3 P002 0.119233646 0.9953417
## ABCD4 si1 P002 0.426813663 1.0243528
## ABCD4 si2 P002 0.010945096 0.7004862
## ABCD4 si3 P002 0.943910379 0.9849500
## ABCE1 si1 P002 0.462152498 0.9750910
## ABCE1 si2 P002 0.354984650 1.0993073
## ABCE1 si3 P002 0.881927115 0.9960589
## ABCF1 si1 P002 0.061153868 0.7338464
## ABCF1 si2 P002 0.279290862 1.1171634
## ABCF1 si3 P002 0.815699384 1.0104147
## ABCF2 si1 P002 0.103164356 0.9903778
## ABCF2 si2 P002 0.899930202 0.9817420
## ABCF2 si3 P002 0.563327695 1.0500866
## ABCF3 si1 P002 0.347034346 1.0923048
## ABCF3 si2 P002 0.230364326 1.1954041
## ABCF3 si3 P002 0.862997225 1.0103382
## ABCG1 si1 P002 0.839041064 1.0068493
## ABCG1 si2 P002 0.672071875 0.9775591
## ABCG1 si3 P002 0.016615641 0.6380406
## ABCG2 si1 P002 0.313191127 1.0933218
## ABCG2 si2 P002 0.885947999 0.9889652
## ABCG2 si3 P002 0.870678954 0.9950617
## ABCG4 si1 P002 0.918122984 0.9848023
## ABCG4 si2 P002 0.935789653 0.9926608
## ABCG4 si3 P002 0.749055620 0.9987408
## ABCG5 si1 P002 0.390378954 1.0410105
## ABCG5 si2 P002 0.539566239 1.0359672
## ABCG5 si3 P002 0.769496006 0.9656420
## ABCG8 si1 P002 0.134774551 0.9583968
## ABCG8 si2 P002 0.595614673 1.0302524
## ABCG8 si3 P002 0.868120395 1.0101765
## ABHD1 si1 P002 0.530473645 1.0287974
## ABHD1 si2 P002 0.440736801 1.0073984
## ABHD1 si3 P002 0.998192954 0.9981930
## ABHD11 si1 P002 0.668319105 0.9795211
## ABHD11 si2 P002 0.113259135 0.9884434
## ABHD11 si3 P002 0.587849696 1.0354784
## ABHD14B si1 P002 0.159907391 1.0234073
## ABHD14B si2 P002 0.757047486 0.9887967
## ABHD14B si3 P002 0.119425151 0.9171852
## ABHD2 si1 P002 0.605712048 1.0112758
## ABHD2 si2 P002 0.197717252 1.0846204
## ABHD2 si3 P002 0.518429117 1.0368582
## ABHD3 si1 P002 0.435819506 1.0204554
## ABHD3 si2 P002 0.025306383 0.6941179
## ABHD3 si3 P002 0.474534117 0.9903321
## ABHD4 si1 P002 0.759082132 0.9847552
## ABHD4 si2 P002 0.963909006 0.9844177
## ABHD4 si3 P002 0.696948392 0.9912155
## ABHD5 si1 P002 0.868962646 1.0050652
## ABHD5 si2 P002 0.936627941 0.9880910
## ABHD5 si3 P002 0.389320274 1.0528098
## ABHD6 si1 P002 0.970309983 0.9857117
## ABHD6 si2 P002 0.405517590 1.0381250
## ABHD6 si3 P002 0.459582833 0.9804434
## ABHD8 si1 P002 0.299972502 1.0665689
## ABHD8 si2 P002 0.666925791 0.9849981
## ABHD8 si3 P002 0.738905944 1.0061698
## ABI2 si1 P002 0.280225090 1.0980248
## ABI2 si2 P002 0.413702367 1.0451428
## ABI2 si3 P002 0.686567286 0.9911347
## ABI3 si1 P002 0.584109679 1.0481220
## ABI3 si2 P002 0.713469814 0.9998993
## ABI3 si3 P002 0.618590539 1.0065202
## ABL1 si1 P002 0.938818928 0.9849904
## ABL1 si2 P002 0.884698135 0.9933453
## ABL1 si3 P002 0.640986156 0.9924947
## ABL2 si1 P002 0.763349136 0.9836445
## ABL2 si2 P002 0.197547616 1.1155630
## ABL2 si3 P002 0.585241859 1.0404300
## ABLIM1 si1 P002 0.067976048 0.7677295
## ABLIM1 si2 P002 0.664075567 0.9961134
## ABLIM1 si3 P002 0.793376467 0.9891447
## ABLIM2 si1 P002 0.822366439 0.9993314
## ABLIM2 si2 P002 0.816440312 1.0048496
## ABLIM2 si3 P002 0.356128896 1.0853452
## ABLIM3 si1 P002 0.056101077 0.7693862
## ABLIM3 si2 P002 0.849580807 1.0131647
## ABLIM3 si3 P002 0.266412655 1.1895635
## ABO si1 P002 0.630682575 1.0007525
## ABO si2 P002 0.282479126 1.0634508
## ABO si3 P002 0.119375518 0.9550041
## ABP1 si1 P002 0.577460355 1.0559275
## ABP1 si2 P002 0.270754864 1.1814758
## ABP1 si3 P002 0.873694001 0.9925991
## ABR si1 P002 0.252732971 1.2131183
## ABR si2 P002 0.447488689 0.9990445
## ABR si3 P002 0.851025262 1.0086225
## ABTB2 si1 P002 0.718638721 0.9998452
## ABTB2 si2 P002 0.378302899 1.1005175
## ABTB2 si3 P002 0.730465934 1.0017819
## ACAA1 si1 P002 0.659687267 0.9973225
## ACAA1 si2 P002 0.029707946 0.5703926
## ACAA1 si3 P002 0.893832187 0.9806616
## ACAA2 si1 P002 0.855900419 1.0081772
## ACAA2 si2 P002 0.025957681 0.6229843
## ACAA2 si3 P002 0.819768072 1.0025189
## ACACA si1 P002 0.356382640 1.0691479
## ACACA si2 P002 0.447427113 1.0106589
## ACACA si3 P002 0.581338093 1.0529898
## ACACB si1 P002 0.666436463 0.9919054
## ACACB si2 P002 0.483466055 0.9981235
## ACACB si3 P002 0.453019502 0.9884062
## ACAD10 si1 P002 0.432824571 1.0259545
## ACAD10 si2 P002 0.033093517 0.5776323
## ACAD10 si3 P002 0.426205272 1.0491207
## ACAD11 si1 P002 0.600650432 1.0205742
## ACAD11 si2 P002 0.053731755 0.7935767
## ACAD11 si3 P002 0.650869134 0.9918006
## ACAD8 si1 P002 0.768727298 0.9710240
## ACAD8 si2 P002 0.485043610 0.9907274
## ACAD8 si3 P002 0.028389680 0.6056465
## ACAD9 si1 P002 0.403339281 1.0465019
## ACAD9 si2 P002 0.591642087 1.0326844
## ACAD9 si3 P002 0.634407977 0.9984126
## ACADL si1 P002 0.751388002 0.9949414
## ACADL si2 P002 0.890218948 0.9823106
## ACADL si3 P002 0.398893297 1.0491440
## ACADM si1 P002 0.645601082 0.9916433
## ACADM si2 P002 0.437807704 1.0127600
## ACADM si3 P002 0.380437748 1.0741772
## ACADS si1 P002 0.125319155 0.9254338
## ACADS si2 P002 0.384891740 1.0710031
## ACADS si3 P002 0.149127666 0.9873280
## ACADSB si1 P002 0.045772969 0.7323675
## ACADSB si2 P002 0.744343320 1.0064360
## ACADSB si3 P002 0.346696608 1.1094291
## ACADVL si1 P002 0.722781915 0.9983750
## ACADVL si2 P002 0.627779753 1.0044476
## ACADVL si3 P002 0.097093224 0.9811526
## ACAT1 si1 P002 0.451020637 0.9953559
## ACAT1 si2 P002 0.340999975 1.1288275
## ACAT1 si3 P002 0.549151512 1.0439316
## ACAT2 si1 P002 0.766903764 0.9816368
## ACAT2 si2 P002 0.146943844 1.0076149
## ACAT2 si3 P002 0.345382727 1.1239573
## ACBD4 si1 P002 0.627325357 1.0121552
## ACBD4 si2 P002 0.001072556 0.2059307
## ACBD4 si3 P002 0.638949270 0.9973842
## ACBD5 si1 P002 0.752098473 0.9890610
## ACBD5 si2 P002 0.254530141 1.1919460
## ACBD5 si3 P002 0.282876360 1.0444666
## ACBD7 si1 P002 0.058775215 0.7523228
## ACBD7 si2 P002 0.281834000 1.0822426
## ACBD7 si3 P002 0.998156137 1.0033821
## ASIC1 si1 P002 0.613687474 1.0070769
## ASIC1 si2 P002 0.980532424 0.9908538
## ASIC1 si3 P002 0.957501397 0.9937312
## ASIC2 si1 P002 0.424862942 1.0593985
## ASIC2 si2 P002 0.109324842 0.9995414
## ASIC2 si3 P002 0.272786496 1.1385871
## ASIC3 si1 P002 0.869425082 0.9995785
## ASIC3 si2 P002 0.551040348 1.0372524
## ASIC3 si3 P002 0.232042602 1.1724258
## ASIC4 si1 P002 0.247652676 1.2192132
## ASIC4 si2 P002 0.916278255 0.9883451
## ASIC4 si3 P002 0.696115943 0.9974199
## ASIC5 si1 P002 0.521271523 1.0317952
## ASIC5 si2 P002 0.887836480 0.9853445
## ASIC5 si3 P002 0.022798699 0.7295584
## EPHX3 si1 P002 0.563727755 1.0407282
## EPHX3 si2 P002 0.214656423 1.1448343
## EPHX3 si3 P002 0.767645625 0.9760792
## EPHX4 si1 P002 0.523706346 1.0260369
## EPHX4 si2 P002 0.265481841 1.2136313
## EPHX4 si3 P002 0.181925496 1.1267644
## treat_cont_log2FC
## ABCC3 si1 0.044456432
## ABCC3 si2 -0.760802637
## ABCC3 si3 0.358396601
## ABCC4 si1 0.267931801
## ABCC4 si2 -0.327492585
## ABCC4 si3 0.779763082
## ABCC5 si1 -0.076314991
## ABCC5 si2 -0.282309124
## ABCC5 si3 0.204441978
## ABCC6 si1 0.013722404
## ABCC6 si2 0.393094853
## ABCC6 si3 0.578421782
## ABCC8 si1 -0.005542042
## ABCC8 si2 -0.174731490
## ABCC8 si3 -0.218958943
## ABCC9 si1 -0.219310783
## ABCC9 si2 0.654085148
## ABCC9 si3 0.677203667
## ABCD1 si1 -1.125179351
## ABCD1 si2 0.193299533
## ABCD1 si3 0.981209240
## ABCD2 si1 -1.431438610
## ABCD2 si2 -0.217317723
## ABCD2 si3 -0.456788966
## ABCD3 si1 0.819880560
## ABCD3 si2 -0.024680790
## ABCD3 si3 -0.629989525
## ABCD4 si1 -0.090501041
## ABCD4 si2 -1.381591360
## ABCD4 si3 0.639788489
## ABCE1 si1 -0.112360388
## ABCE1 si2 -0.080918449
## ABCE1 si3 0.437538631
## ABCF1 si1 -0.921272129
## ABCF1 si2 -0.102947288
## ABCF1 si3 0.585318877
## ABCF2 si1 -0.312888370
## ABCF2 si2 0.774845946
## ABCF2 si3 0.093205920
## ABCF3 si1 -0.314432073
## ABCF3 si2 -0.391734685
## ABCF3 si3 0.854337230
## ABCG1 si1 0.585084242
## ABCG1 si2 0.090394283
## ABCG1 si3 -1.251182307
## ABCG2 si1 0.017459657
## ABCG2 si2 0.379954871
## ABCG2 si3 0.422075318
## ABCG4 si1 0.587763766
## ABCG4 si2 0.773145755
## ABCG4 si3 0.216871300
## ABCG5 si1 -0.280006137
## ABCG5 si2 0.129182363
## ABCG5 si3 0.472440767
## ABCG8 si1 -0.680042026
## ABCG8 si2 0.017963946
## ABCG8 si3 0.347747876
## ABHD1 si1 0.162969077
## ABHD1 si2 0.085302532
## ABHD1 si3 1.393706791
## ABHD11 si1 0.083918972
## ABHD11 si2 -0.717691619
## ABHD11 si3 -0.036293016
## ABHD14B si1 -0.439170996
## ABHD14B si2 0.248597742
## ABHD14B si3 -0.597586959
## ABHD2 si1 0.065038842
## ABHD2 si2 -0.312503983
## ABHD2 si3 -0.081481447
## ABHD3 si1 -0.205419896
## ABHD3 si2 -1.179902584
## ABHD3 si3 -0.076286182
## ABHD4 si1 0.263866452
## ABHD4 si2 0.667375117
## ABHD4 si3 0.114754756
## ABHD5 si1 0.703921794
## ABHD5 si2 0.591713192
## ABHD5 si3 -0.057177808
## ABHD6 si1 0.670304416
## ABHD6 si2 -0.193788714
## ABHD6 si3 -0.168639069
## ABHD8 si1 -0.071761444
## ABHD8 si2 0.066793513
## ABHD8 si3 0.503993578
## ABI2 si1 -0.350737178
## ABI2 si2 -0.170799350
## ABI2 si3 0.161161001
## ABI3 si1 -0.032458777
## ABI3 si2 0.114109652
## ABI3 si3 0.169608088
## ABL1 si1 0.810998941
## ABL1 si2 0.756045242
## ABL1 si3 0.097997447
## ABL2 si1 0.258869064
## ABL2 si2 -0.523514706
## ABL2 si3 0.163684520
## ABLIM1 si1 -0.683457325
## ABLIM1 si2 0.125183999
## ABLIM1 si3 0.534979480
## ABLIM2 si1 0.363442411
## ABLIM2 si2 0.315786286
## ABLIM2 si3 0.001937280
## ABLIM3 si1 -0.880693390
## ABLIM3 si2 0.898054996
## ABLIM3 si3 -0.364434288
## ABO si1 0.206667577
## ABO si2 -0.177199044
## ABO si3 -0.281218574
## ABP1 si1 0.014887152
## ABP1 si2 0.036565994
## ABP1 si3 0.455269051
## ABR si1 -0.443655545
## ABR si2 -0.118274133
## ABR si3 0.285198327
## ABTB2 si1 1.162764808
## ABTB2 si2 -0.282304116
## ABTB2 si3 0.149517027
## ACAA1 si1 0.052764006
## ACAA1 si2 -1.012411291
## ACAA1 si3 0.486341376
## ACAA2 si1 0.465151903
## ACAA2 si2 -0.734906859
## ACAA2 si3 0.342300507
## ACACA si1 -0.211815849
## ACACA si2 -0.100766584
## ACACA si3 -0.005799791
## ACACB si1 0.089491198
## ACACB si2 0.241067890
## ACACB si3 -0.051409880
## ACAD10 si1 -0.036458110
## ACAD10 si2 -0.957827735
## ACAD10 si3 -0.071581846
## ACAD11 si1 0.027718084
## ACAD11 si2 -0.885675167
## ACAD11 si3 0.051988851
## ACAD8 si1 0.200749107
## ACAD8 si2 -0.103544464
## ACAD8 si3 -0.770541899
## ACAD9 si1 -0.243170224
## ACAD9 si2 0.043175028
## ACAD9 si3 0.026779967
## ACADL si1 0.243643945
## ACADL si2 0.449085351
## ACADL si3 -0.180493820
## ACADM si1 0.025026745
## ACADM si2 -0.192121106
## ACADM si3 -0.257786830
## ACADS si1 -0.640906681
## ACADS si2 -0.272716052
## ACADS si3 0.056707168
## ACADSB si1 -0.473784934
## ACADSB si2 0.507064969
## ACADSB si3 -0.097113174
## ACADVL si1 0.365873720
## ACADVL si2 0.036025181
## ACADVL si3 -0.489483387
## ACAT1 si1 -0.072012857
## ACAT1 si2 -0.285856336
## ACAT1 si3 0.032240911
## ACAT2 si1 0.565624932
## ACAT2 si2 -0.676201521
## ACAT2 si3 -0.182614353
## ACBD4 si1 0.249093646
## ACBD4 si2 -1.743848199
## ACBD4 si3 0.104378536
## ACBD5 si1 0.306742737
## ACBD5 si2 -0.396079974
## ACBD5 si3 -0.292262982
## ACBD7 si1 -0.575856568
## ACBD7 si2 0.189933993
## ACBD7 si3 1.400125973
## ASIC1 si1 0.312425544
## ASIC1 si2 0.869660062
## ASIC1 si3 0.712949246
## ASIC2 si1 -0.145883683
## ASIC2 si2 -0.374238897
## ASIC2 si3 -0.282852897
## ASIC3 si1 0.523846473
## ASIC3 si2 -0.023835670
## ASIC3 si3 -0.492121964
## ASIC4 si1 -0.438901408
## ASIC4 si2 0.491036204
## ASIC4 si3 0.209365687
## ASIC5 si1 0.179450683
## ASIC5 si2 0.881292671
## ASIC5 si3 -0.927144040
## EPHX3 si1 0.138202543
## EPHX3 si2 -0.251764851
## EPHX3 si3 0.275231494
## EPHX4 si1 -0.104575185
## EPHX4 si2 -0.403204759
## EPHX4 si3 -0.548263677
##
## $P003
## MASTER_PLATE pvalue_treat_lowerthan_cont FDR_treat_lowerthan_cont
## ACE si1 P003 0.791227985 0.9676164
## ACE si2 P003 0.530531188 1.0394080
## ACE si3 P003 0.287733643 1.1048972
## ACE2 si1 P003 0.127036473 0.9756401
## ACE2 si2 P003 0.502903147 1.0849147
## ACE2 si3 P003 0.551831243 1.0090628
## ACHE si1 P003 0.285646898 1.1192695
## ACHE si2 P003 0.112615417 0.9400939
## ACHE si3 P003 0.527457183 1.0549144
## ACIN1 si1 P003 0.229469889 1.0490052
## ACIN1 si2 P003 0.435716297 1.1305071
## ACIN1 si3 P003 0.021331636 0.5119593
## ACLY si5 P003 0.992742476 0.9927425
## ACLY si6 P003 0.196641111 0.9935551
## ACLY si7 P003 0.527489450 1.0441028
## ACMSD si1 P003 0.404390222 1.1418077
## ACMSD si2 P003 0.675118826 0.9819910
## ACMSD si3 P003 0.184168140 1.0715237
## ACO1 si1 P003 0.605847915 0.9942120
## ACO1 si2 P003 0.536830245 1.0411253
## ACO1 si3 P003 0.982635349 0.9877800
## ACO2 si1 P003 0.247868635 1.0816086
## ACO2 si2 P003 0.649948610 0.9825995
## ACO2 si3 P003 0.696955678 0.9767554
## ACOT11 si1 P003 0.913371236 0.9907756
## ACOT11 si2 P003 0.626226561 0.9775244
## ACOT11 si3 P003 0.526206448 1.0748047
## ACOT12 si1 P003 0.887076008 0.9845005
## ACOT12 si2 P003 0.826866904 0.9739782
## ACOT12 si3 P003 0.819399634 0.9771722
## ACOT2 si1 P003 0.287851279 1.0836754
## ACOT2 si2 P003 0.288805926 1.0663603
## ACOT2 si3 P003 0.751879079 0.9820461
## ACOT4 si1 P003 0.002001684 0.3843234
## ACOT4 si2 P003 0.335603856 1.1109645
## ACOT4 si3 P003 0.010509558 1.0089176
## ACOT7 si1 P003 0.382807814 1.1307554
## ACOT7 si2 P003 0.550309344 1.0159557
## ACOT7 si3 P003 0.673660946 0.9949454
## ACOT8 si1 P003 0.549815056 1.0248980
## ACOT8 si2 P003 0.556818816 0.9991515
## ACOT8 si3 P003 0.492086979 1.0736443
## ACOT9 si1 P003 0.566862021 1.0077547
## ACOT9 si2 P003 0.265726260 1.1091183
## ACOT9 si3 P003 0.190354795 1.0152256
## ACOX1 si1 P003 0.867845284 0.9801547
## ACOX1 si2 P003 0.696779198 0.9836883
## ACOX1 si3 P003 0.755252654 0.9732115
## ACOX2 si1 P003 0.902737702 0.9904322
## ACOX2 si2 P003 0.526897071 1.0648867
## ACOX2 si3 P003 0.752510099 0.9762293
## ACOX3 si1 P003 0.222319632 1.0411066
## ACOX3 si2 P003 0.945063051 1.0024978
## ACOX3 si3 P003 0.851766290 0.9676871
## ACOXL si1 P003 0.345903510 1.1068912
## ACOXL si2 P003 0.847377725 0.9684317
## ACOXL si3 P003 0.423241292 1.1286434
## ACP1 si1 P003 0.410914416 1.1270795
## ACP1 si2 P003 0.555144217 1.0055442
## ACP1 si3 P003 0.821923524 0.9741316
## ACP2 si1 P003 0.657867712 0.9868016
## ACP2 si2 P003 0.055326000 0.8852160
## ACP2 si3 P003 0.093613786 0.8558975
## ACP5 si1 P003 0.077908379 0.9349005
## ACP5 si2 P003 0.407863319 1.1349240
## ACP5 si3 P003 0.970390475 0.9963367
## ACP6 si1 P003 0.424411171 1.1162595
## ACP6 si2 P003 0.976183719 0.9916787
## ACP6 si3 P003 0.629724019 0.9750565
## ACPL2 si1 P003 0.482258676 1.0893372
## ACPL2 si2 P003 0.445512673 1.0827650
## ACPL2 si3 P003 0.609463555 0.9916695
## ACPP si1 P003 0.511765169 1.0680317
## ACPP si2 P003 0.490256725 1.0819459
## ACPP si3 P003 0.489974807 1.0938972
## ACPT si1 P003 0.829861790 0.9715455
## ACPT si2 P003 0.699313879 0.9729584
## ACPT si3 P003 0.437584955 1.1202175
## ACR si1 P003 0.413223728 1.1174501
## ACR si2 P003 0.186405357 1.0526420
## ACR si3 P003 0.708705309 0.9719387
## ACSBG1 si1 P003 0.590494641 0.9945173
## ACSBG1 si2 P003 0.800837899 0.9670495
## ACSBG1 si3 P003 0.580335299 0.9860564
## ACSBG2 si1 P003 0.019563819 0.5366076
## ACSBG2 si2 P003 0.792868839 0.9634862
## ACSBG2 si3 P003 0.620516376 0.9928262
## ACSL1 si1 P003 0.816003905 0.9792047
## ACSL1 si2 P003 0.482004825 1.1017253
## ACSL1 si3 P003 0.913311091 0.9963394
## ACSL3 si1 P003 0.118949615 0.9515969
## ACSL3 si2 P003 0.868865419 0.9755682
## ACSL3 si3 P003 0.763519713 0.9644460
## ACSL4 si1 P003 0.252368428 1.0767720
## ACSL4 si2 P003 0.041862978 0.7306993
## ACSL4 si3 P003 0.313774304 1.0569240
## ACSL5 si1 P003 0.216804043 1.0406594
## ACSL5 si2 P003 0.843669912 0.9699678
## ACSL5 si3 P003 0.202202675 0.9954593
## ACSL6 si1 P003 0.173887121 1.0769783
## ACSL6 si2 P003 0.440452517 1.0982712
## ACSL6 si3 P003 0.383446472 1.1154806
## ACSM1 si1 P003 0.450916173 1.0821988
## ACSM1 si2 P003 0.696349213 0.9903633
## ACSM1 si3 P003 0.294061607 1.0455524
## ACSM2A si1 P003 0.746644061 0.9818881
## ACSM2A si2 P003 0.459820336 1.0899445
## ACSM2A si3 P003 0.166661760 1.1428235
## ACSM2B si1 P003 0.728038791 0.9707184
## ACSM2B si2 P003 0.595149382 0.9936407
## ACSM2B si3 P003 0.076150864 0.9747311
## ACSM3 si1 P003 0.310720722 1.0653282
## ACSM3 si2 P003 0.167268174 1.1074307
## ACSM3 si3 P003 0.709409105 0.9592010
## ACSS1 si1 P003 0.078249232 0.8837560
## ACSS1 si2 P003 0.959620050 1.0013427
## ACSS1 si3 P003 0.064015656 0.9454620
## ACSS2 si1 P003 0.895644173 0.9882970
## ACSS2 si2 P003 0.310225787 1.0829700
## ACSS2 si3 P003 0.371892710 1.1156781
## ACTA1 si1 P003 0.574669251 0.9940225
## ACTA1 si2 P003 0.580067288 0.9944011
## ACTA1 si3 P003 0.647039728 0.9938530
## ACTA2 si1 P003 0.089455145 0.8587694
## ACTA2 si2 P003 0.625006388 0.9836166
## ACTA2 si3 P003 0.884707466 0.9875804
## ACTC1 si1 P003 0.352077507 1.1081784
## ACTC1 si2 P003 0.784027000 0.9711818
## ACTC1 si3 P003 0.106576646 0.9301235
## ACTG1 si1 P003 0.604456172 1.0004792
## ACTG1 si2 P003 0.951587267 0.9983866
## ACTG1 si3 P003 0.833142652 0.9694751
## ACTG2 si1 P003 0.015284348 0.5869189
## ACTG2 si2 P003 0.540460268 1.0274096
## ACTG2 si3 P003 0.080519657 0.8588763
## ACTL6A si1 P003 0.787120984 0.9687643
## ACTL6A si2 P003 0.647141230 0.9861200
## ACTL6A si3 P003 0.463597727 1.0854971
## ACTL6B si1 P003 0.018132317 0.5802341
## ACTL6B si2 P003 0.469168800 1.0853061
## ACTL6B si3 P003 0.544809705 1.0255241
## ACTL7A si1 P003 0.947311756 0.9993619
## ACTL7A si2 P003 0.401499782 1.1505665
## ACTL7A si3 P003 0.038173025 0.7329221
## ACTL7B si1 P003 0.358285698 1.1095299
## ACTL7B si2 P003 0.915288372 0.9872773
## ACTL7B si3 P003 0.013357283 0.8548661
## ACTL8 si1 P003 0.438377488 1.1074800
## ACTL8 si2 P003 0.944358103 1.0073153
## ACTL8 si3 P003 0.838918687 0.9703156
## ACTN1 si1 P003 0.961145823 0.9975135
## ACTN1 si2 P003 0.245319090 1.0953783
## ACTN1 si3 P003 0.194353782 1.0085385
## ACTN2 si1 P003 0.755710725 0.9673097
## ACTN2 si2 P003 0.674801182 0.9890216
## ACTN2 si3 P003 0.088171160 0.8909928
## ACTN3 si1 P003 0.369387083 1.1257511
## ACTN3 si2 P003 0.526089814 1.0861209
## ACTN3 si3 P003 0.932866161 1.0006162
## ACTN4 si1 P003 0.567240936 0.9991767
## ACTN4 si2 P003 0.142702890 1.0538060
## ACTN4 si3 P003 0.540077416 1.0369486
## ACTR10 si1 P003 0.278005655 1.1356827
## ACTR10 si2 P003 0.145570220 1.0351660
## ACTR10 si3 P003 0.732175371 0.9695012
## ACTR1A si1 P003 0.609774083 0.9838372
## ACTR1A si2 P003 0.280007352 1.1200294
## ACTR1A si3 P003 0.620789306 0.9850541
## ACTR1B si1 P003 0.662883993 0.9866180
## ACTR1B si2 P003 0.442633858 1.0895603
## ACTR1B si3 P003 0.014405579 0.6914678
## ACTR2 si1 P003 0.506110745 1.0797029
## ACTR2 si2 P003 0.677888911 0.9713035
## ACTR2 si3 P003 0.189712491 1.0407085
## ACTR3 si1 P003 0.510098261 1.0762513
## ACTR3 si2 P003 0.773843379 0.9647917
## ACTR3 si3 P003 0.070253520 0.9634768
## ACTR3B si1 P003 0.980389960 0.9907099
## ACTR3B si2 P003 0.168074470 1.0756766
## ACTR3B si3 P003 0.709157104 0.9656607
## ACTR5 si1 P003 0.722432937 0.9699799
## ACTR5 si2 P003 0.677665724 0.9782844
## ACTR5 si3 P003 0.702420388 0.9702497
## ACTR6 si1 P003 0.962533129 0.9935826
## ACTR6 si2 P003 0.766990807 0.9624983
## ACTR6 si3 P003 0.574432943 1.0026466
## ACTR8 si1 P003 0.973666550 0.9943829
## ACTR8 si2 P003 0.292223001 1.0586192
## ACTR8 si3 P003 0.760809727 0.9673872
## ACTRT1 si1 P003 0.177039322 1.0622359
## ACTRT1 si2 P003 0.030367156 0.6478327
## ACTRT1 si3 P003 0.338215530 1.1006336
## treat_cont_log2FC
## ACE si1 0.4333006257
## ACE si2 0.0844739612
## ACE si3 -0.3573622807
## ACE2 si1 -0.7174841691
## ACE2 si2 -0.0827676448
## ACE2 si3 0.0008850614
## ACHE si1 -0.3303339557
## ACHE si2 -0.4400466650
## ACHE si3 -0.0033495591
## ACIN1 si1 -0.3650085502
## ACIN1 si2 -0.1594153394
## ACIN1 si3 -0.7978322059
## ACLY si5 1.2271537311
## ACLY si6 0.0529250362
## ACLY si7 0.1591417376
## ACMSD si1 0.1084406175
## ACMSD si2 0.3926564565
## ACMSD si3 -0.4486922126
## ACO1 si1 0.1483454150
## ACO1 si2 0.1696506295
## ACO1 si3 1.2268067911
## ACO2 si1 -0.0656122562
## ACO2 si2 0.1183854575
## ACO2 si3 0.5361326215
## ACOT11 si1 0.6240093555
## ACOT11 si2 0.1079269297
## ACOT11 si3 0.0140340074
## ACOT12 si1 0.6149884511
## ACOT12 si2 0.3440757295
## ACOT12 si3 0.6261636653
## ACOT2 si1 -0.4232798535
## ACOT2 si2 -0.2857313234
## ACOT2 si3 0.1997114304
## ACOT4 si1 -1.5487654112
## ACOT4 si2 -0.2491704073
## ACOT4 si3 -1.1836534492
## ACOT7 si1 -0.2417836535
## ACOT7 si2 0.8413503641
## ACOT7 si3 0.1604829400
## ACOT8 si1 -0.0678154423
## ACOT8 si2 0.0567421156
## ACOT8 si3 -0.1254724884
## ACOT9 si1 0.2121646700
## ACOT9 si2 -0.2872966756
## ACOT9 si3 -0.4846061632
## ACOX1 si1 0.9414341896
## ACOX1 si2 0.1464588201
## ACOX1 si3 0.3368252241
## ACOX2 si1 0.8584981777
## ACOX2 si2 -0.0682056897
## ACOX2 si3 0.1560271195
## ACOX3 si1 -0.5248625579
## ACOX3 si2 0.6773259294
## ACOX3 si3 0.3489700311
## ACOXL si1 -0.2339952356
## ACOXL si2 0.3975575592
## ACOXL si3 -0.1604547870
## ACP1 si1 -0.2064055333
## ACP1 si2 0.0470607804
## ACP1 si3 0.3260261754
## ACP2 si1 0.3656135646
## ACP2 si2 -0.5006882318
## ACP2 si3 -0.1988024211
## ACP5 si1 -0.4493272952
## ACP5 si2 -0.1454494014
## ACP5 si3 0.9595259606
## ACP6 si1 -0.0375874581
## ACP6 si2 0.8092553053
## ACP6 si3 0.5770228485
## ACPL2 si1 -0.0474140972
## ACPL2 si2 0.1639714839
## ACPL2 si3 0.3786313293
## ACPP si1 -0.0657335781
## ACPP si2 -0.0243645932
## ACPP si3 0.0766522398
## ACPT si1 0.7915266506
## ACPT si2 0.2843823942
## ACPT si3 -0.1379115652
## ACR si1 0.1818829434
## ACR si2 -0.5071285833
## ACR si3 0.2793989609
## ACSBG1 si1 0.0923281457
## ACSBG1 si2 0.2400015930
## ACSBG1 si3 0.0169762389
## ACSBG2 si1 -0.9772830927
## ACSBG2 si2 0.3602033442
## ACSBG2 si3 0.0745930243
## ACSL1 si1 0.3072016332
## ACSL1 si2 -0.1057354937
## ACSL1 si3 0.8890896789
## ACSL3 si1 -0.2896470546
## ACSL3 si2 0.6198467815
## ACSL3 si3 0.5110131501
## ACSL4 si1 -0.3532183393
## ACSL4 si2 -0.4947397017
## ACSL4 si3 -0.1765487391
## ACSL5 si1 -0.3427417837
## ACSL5 si2 0.5874810075
## ACSL5 si3 -0.1756598468
## ACSL6 si1 -0.4833815365
## ACSL6 si2 0.2483484671
## ACSL6 si3 -0.2105453677
## ACSM1 si1 -0.0870632038
## ACSM1 si2 0.6480416783
## ACSM1 si3 0.3509841762
## ACSM2A si1 0.1913379991
## ACSM2A si2 -0.0224704843
## ACSM2A si3 -0.5338234584
## ACSM2B si1 0.3206642701
## ACSM2B si2 0.0993041555
## ACSM2B si3 -0.7110579035
## ACSM3 si1 -0.0247322491
## ACSM3 si2 -0.5618993265
## ACSM3 si3 0.1726793531
## ACSS1 si1 -0.4391968757
## ACSS1 si2 1.0505062086
## ACSS1 si3 -0.3754022324
## ACSS2 si1 0.6750008057
## ACSS2 si2 -0.2563903946
## ACSS2 si3 0.0844092083
## ACTA1 si1 -0.0219258643
## ACTA1 si2 0.0913887628
## ACTA1 si3 0.0789910417
## ACTA2 si1 -0.4651611780
## ACTA2 si2 0.1379289136
## ACTA2 si3 0.5736643529
## ACTC1 si1 -0.2134078378
## ACTC1 si2 0.2787356982
## ACTC1 si3 -0.1841597411
## ACTG1 si1 0.0095072351
## ACTG1 si2 0.7992892016
## ACTG1 si3 0.3714543008
## ACTG2 si1 -1.0421048718
## ACTG2 si2 0.0049684099
## ACTG2 si3 -0.7818945453
## ACTL6A si1 0.4001044742
## ACTL6A si2 0.3142103417
## ACTL6A si3 -0.1485049454
## ACTL6B si1 -1.0919136548
## ACTL6B si2 0.2068464107
## ACTL6B si3 -0.0302153582
## ACTL7A si1 0.6629860859
## ACTL7A si2 -0.2245243478
## ACTL7A si3 -0.9229206161
## ACTL7B si1 -0.3142015814
## ACTL7B si2 0.5603251529
## ACTL7B si3 -1.1022927224
## ACTL8 si1 -0.1804733626
## ACTL8 si2 0.7223029128
## ACTL8 si3 0.4421572016
## ACTN1 si1 0.8786988787
## ACTN1 si2 -0.4561402610
## ACTN1 si3 -0.4618440214
## ACTN2 si1 0.3055256903
## ACTN2 si2 0.1665056189
## ACTN2 si3 -0.6317479432
## ACTN3 si1 -0.0831323549
## ACTN3 si2 -0.0705066496
## ACTN3 si3 0.8414429371
## ACTN4 si1 0.0645177829
## ACTN4 si2 -0.2799849839
## ACTN4 si3 -0.0582432500
## ACTR10 si1 -0.4171049987
## ACTR10 si2 -0.3741925025
## ACTR10 si3 0.1296686032
## ACTR1A si1 0.5244760888
## ACTR1A si2 -0.0601996130
## ACTR1A si3 0.4207186758
## ACTR1B si1 0.0792681460
## ACTR1B si2 -0.2059447676
## ACTR1B si3 -0.9761028687
## ACTR2 si1 0.0308812841
## ACTR2 si2 0.0872939033
## ACTR2 si3 -0.3558643671
## ACTR3 si1 0.0224556455
## ACTR3 si2 0.3834879848
## ACTR3 si3 -0.4760454396
## ACTR3B si1 1.1666001129
## ACTR3B si2 -0.5666646068
## ACTR3B si3 0.2149889541
## ACTR5 si1 0.2231762421
## ACTR5 si2 0.1570354783
## ACTR5 si3 0.2561972114
## ACTR6 si1 0.8895109323
## ACTR6 si2 0.4258319547
## ACTR6 si3 0.0124899564
## ACTR8 si1 0.6958281997
## ACTR8 si2 -0.2480679708
## ACTR8 si3 0.3057647497
## ACTRT1 si1 -0.5235631450
## ACTRT1 si2 -1.0014429414
## ACTRT1 si3 -0.2965465714
We could select interesting siRNA hits by FDR or pvalue in the rankp.c,
Redundant siRNA activity (RSA) is a novel method proposed for RNAi screen data, which systemically employs the information provided by multiple siRNAs targeting a single gene to reduce the off-target and improve confirmation rate. Briefly, RSA calculates a P-value for the rank distribution of multiple siRNAs silenced the same gene under the background of all siRNA signals in the experiment by iterative hypergeometric distribution formula. Compared to the methods mentioned above, siRNAs targeted the same genes have identical P-value, and genes with several moderately effect siRNAs may have smaller P-value than genes with fewer strong effect siRNAs. Synlet provides a wrapper function to use the RSA R codes [8].
rsaHits(example_dt, treatment = "treatment", control = "control",
normMethod = "PLATE", LB = 0.2, UB = 0.8, revHits = FALSE,
Bonferroni = FALSE, outputFile = "RSAhits.csv")
## File written to:RSAhits.csv
## Summary of RSA:
## Total #Genes = 192
## Total #Wells = 576
## Total #Hit Genes = 104
## Total #Hit Wells = 133
The meaning of column names:
Gene_ID,Well_ID,Score: columns from input spreadsheet
LogP: OPI p-value in log10, i.e., -2 means 0.0
OPI_Hit: whether the well is a hit, 1 means yes, 0 means no
#hitWell: number of hit wells for the gene
#totalWell: total number of wells for the gene
OPI_Rank: ranking column to sort all wells for hit picking
Cutoff_Rank: ranking column to sort all wells based on Score in the simple activity-based method
Note: a rank value of 999999 means the well is not a hit
Users could pick up hits based on rsa.res, or calculate the FDR from LogP to set a stringent cut-off.
We have went through the process of hits selection in synlet. However, the siRNA hits picked by each algorithms may be different. How to get a reasonable gene lists for further verification? Here are some simple suggestions:
Rule of thumb, an interesting genes should show significant knock-down effect in at least two siRNA to avoid off-target effect.
High confidence: genes are picked up by at least two methods.
Low confidence: Get a union of siRNA identified as hits, then select genes targeted by at least 2 siRNA.
It is always helpful to check the interested genes by siRNAPlot.
[1] Zhang J.H., Chung T.D. & Oldenburg K.R. A simple statistical parameter for use in evaluation and validation of high throughput screening assays. J. Biomol. Screen. B, 4 67-73 (1999).
[2] Birmingham,A. et al. (2009) Statistical methods for analysis of high-throughput RNA interference screens. Nat Methods, 6, 569–575.
[3] Whitehurst,A.W. et al. (2007) Synthetic lethal screen identification of chemosensitizer loci in cancer cells. Nature, 446, 815–819.
[4] Chung,N. et al. (2008) Median Absolute Deviation to Improve Hit Selection for Genome-Scale RNAi Screens. Journal of Biomolecular Screening, 13, 149–158.
[5] Rieber,N. et al. (2009) RNAither, an automated pipeline for the statistical analysis of high-throughput RNAi screens. Bioinformatics, 25, 678–679.
[6] Breitling,R. et al. (2004) Rank products: a simple, yet powerful, new method to detect differentially regulated genes in replicated microarray experiments. FEBS Lett, 573, 83–92.
[7] Hong,F. et al. (2006) RankProd: a bioconductor package for detecting differentially expressed genes in meta-analysis. Bioinformatics, 22, 2825–2827.
[8] Koenig, R. et al. A probability-based approach for the analysis of large-scale RNAi screens. Nat Methods 4, 847-849 (2007).
sessionInfo()
## R version 4.3.1 (2023-06-16)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 22.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.18-bioc/R/lib/libRblas.so
## LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.10.0
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## time zone: America/New_York
## tzcode source: system (glibc)
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] synlet_2.2.0 BiocStyle_2.30.0
##
## loaded via a namespace (and not attached):
## [1] gtable_0.3.4 jsonlite_1.8.7 Rmpfr_0.9-3
## [4] dplyr_1.1.3 compiler_4.3.1 BiocManager_1.30.22
## [7] Rcpp_1.0.11 tidyselect_1.2.0 magick_2.8.1
## [10] RankProd_3.28.0 jquerylib_0.1.4 scales_1.2.1
## [13] yaml_2.3.7 fastmap_1.1.1 ggplot2_3.4.4
## [16] R6_2.5.1 labeling_0.4.3 generics_0.1.3
## [19] patchwork_1.1.3 knitr_1.44 tibble_3.2.1
## [22] bookdown_0.36 munsell_0.5.0 bslib_0.5.1
## [25] pillar_1.9.0 RColorBrewer_1.1-3 rlang_1.1.1
## [28] utf8_1.2.4 cachem_1.0.8 xfun_0.40
## [31] sass_0.4.7 cli_3.6.1 withr_2.5.1
## [34] magrittr_2.0.3 digest_0.6.33 grid_4.3.1
## [37] gmp_0.7-2 lifecycle_1.0.3 vctrs_0.6.4
## [40] evaluate_0.22 glue_1.6.2 data.table_1.14.8
## [43] farver_2.1.1 fansi_1.0.5 colorspace_2.1-0
## [46] rmarkdown_2.25 tools_4.3.1 pkgconfig_2.0.3
## [49] htmltools_0.5.6.1