https://www.mdpi.com/2073-4425/11/7/756/htm

https://shadowcaster.readthedocs.io/en/latest/

https://github.com/dani2s/ShadowCaster_testData

Address of the bookmark: https://github.com/dani2s/ShadowCaster

https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-020-6568-2

Address of the bookmark: https://github.com/HuffordLab/GenomeQC

More at https://www.youtube.com/watch?v=vdGuy1GAJrQ

Address of the bookmark: https://sourceforge.net/projects/accnet/

Address of the bookmark: https://github.com/mkirsche/Jasmine

• recovering signals in networks after a stimulation (cascade network reverse engineering),
• analysing periodic signals.

Address of the bookmark: https://github.com/fbertran/Patterns

DrugBank:

a bioinformatics- cheminformatics resource combining detailed drug data with comprehensive drug target information with >4900 drug (~3500 experimental) and >1500 non-redundant protein entries http://www.drugbank.ca/

Drug-Target Network:

network data of 890 drugs and 394 target human proteins http://www.nature.com/nbt/journal/v25/ n10/suppinfo/nbt1338_S1.html

Drug-Therapy Network:

three layers of drug-therapy networks according to the ATC classification http://www.biomedcentral.com/1471-2210/8/5/additional/

FDA Orange Book:

approved drug products with therapeutic equivalence evaluations http://www.fda.gov/cder/ob/HIDdb: Thomson Investigational drugs database including information on 107000 patents, 25000 investigational drugs and 80000 chemical structures http://scientific.thomson.com/products/iddb/HOMIM: a knowledgebase of human genes and genetic disorders http://www.ncbi.nlm.nih.gov/ sites/entrez?db=omim

PDTD:

3D drug target structure database with a target identification option http://www.dddc.ac.cn/pdtd/

Predicted drug targets:

a set of 1383 predicted drug targets http://www.biomedcentral.com/1471-2105/8/353/additional/ [25] Protein ligand network: a network of 4208 ligands and ~15000 binding sites http://pbil.kaist.ac.kr/~parkkw/Lnet/

TDR Targets Database:

identification and ranking targets against neglected tropical diseases http://tdrtargets.org/

Therapeutic Target Database:

lists >1500 therapeutic targets, disease conditions and corresponding drugs http://xin.cz3.nus.edu.sg/group/cjttd/ttd.asp

More at https://github.com/AntonBankevich/jumboDB

Address of the bookmark: https://github.com/AntonBankevich/jumboDB

What Does DESeq2 Do?
DESeq2 analyzes count data—the number of sequencing reads that map to each gene. It:

Normalizes the data to account for sequencing depth and library size.

Estimates variance (dispersion) for each gene.

Fits a model to compare groups (e.g., control vs treated).

Calculates fold-changes and p-values to determine significance.

Installing DESeq2

You can install DESeq2 via Bioconductor in R:

if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("DESeq2")

Inputs Needed

A count matrix: genes as rows, samples as columns (raw counts, not normalized).

A sample metadata table (also called colData): defines the condition/group for each sample.

Example:
# Count matrix (rows = genes, columns = samples)
counts <- read.csv("counts.csv", row.names = 1)

# Sample metadata
colData <- data.frame(
row.names = colnames(counts),
condition = c("control", "control", "treated", "treated")
)

DESeq2 Workflow

1. Load the package
library(DESeq2)
2. Create a DESeqDataSet object
dds <- DESeqDataSetFromMatrix(countData = counts,
colData = colData,
design = ~ condition)
3. Run the differential expression analysis
dds <- DESeq(dds)
4. Get the results
res <- results(dds)
head(res)
This gives a table with:

log2FoldChange: how much expression changed

pvalue: statistical significance

padj: adjusted p-value (FDR corrected)

Visualization (Optional but Powerful)

MA Plot
plotMA(res, ylim = c(-2, 2))

Volcano Plot (custom)
library(ggplot2)
res$significant <- res$padj < 0.05
ggplot(res, aes(x=log2FoldChange, y=-log10(padj), color=significant)) +
geom_point() +
theme_minimal()

Heatmap of Top Genes
library(pheatmap)
topgenes <- head(order(res$padj), 20)
vsd <- vst(dds, blind=FALSE)
pheatmap(assay(vsd)[topgenes, ])

Tips for Best Results
Use raw counts (not normalized or TPM/RPKM values).

Have replicates: DESeq2 relies on variance estimates, so at least 3 per group is ideal.

Watch out for batch effects—include them in your design if needed (e.g., ~ batch + condition).

Summary

Step Purpose
DESeqDataSetFromMatrix() Load your data into DESeq2
DESeq() Run the differential expression analysis
results() Extract the output (log fold change, p-values, etc.)
plotMA() / ggplot2 / pheatmap Visualize the results

Final Thoughts
DESeq2 is an essential tool for RNA-seq data analysis. It abstracts away much of the complexity of statistical modeling, while still giving you control when needed. Whether you're a bioinformatician or a wet-lab biologist, DESeq2 offers both ease of use and analytical power.

ashr - https://github.com/stephens999/ashr, https://cran.r-project.org/web/packages/ashr/index.html

consensusDE - https://bioconductor.org/packages/release/bioc/html/consensusDE.html

DESeq2 - https://bioconductor.org/packages/release/bioc/html/DESeq2.html

edgeR - https://bioconductor.org/packages/release/bioc/html/edgeR.html

limma - https://kasperdanielhansen.github.io/genbioconductor/html/limma.html  https://bioconductor.org/packages/release/bioc/html/limma.html

MetaCycle - https://cran.r-project.org/web/packages/MetaCycle/index.html, https://github.com/gangwug/MetaCycle

RUVSeq - https://bioconductor.org/packages/release/bioc/html/RUVSeq.html

SARTools - https://github.com/PF2-pasteur-fr/SARTools

tximport - https://github.com/mikelove/tximport