ARACNE (Algorithm for the Reconstruction of Accurate Cellular Networks), a novel algorithm, using microarray expression profiles, specifically designed to scale up to the complexity of regulatory networks in mammalian cells, yet general enough to address a wider range of network deconvolution problems. This method uses an information theoretic approach to eliminate the vast majority of indirect interactions typically inferred by pairwise analysis.
On synthetic datasets ARACNE achieves extremely low error rates and significantly outperforms established methods, such as Relevance Networks and Bayesian Networks. Application to the deconvolution of genetic networks in human B cells demonstrates ARACNE’s ability to infer validated transcriptional targets of the c-MYC proto-oncogene.
Usage.txt: This file is used by the native aracne2 binaries compiled from C++ source to provide ARACNE2 usage summary. Please copy this file to the same directory as the binary (This can be ignored if the entire source distribution is downloaded)
aracne2.exe: PE32 executable for MS Windows (console) Intel 80386 32-bit
aracne2: ELF 64-bit LSB executable, AMD x86-64, version 1 (SYSV), for GNU/Linux 2.6.9, dynamically linked (uses shared libs), for GNU/Linux 2.6.9, not stripped
aracne2.macosx: Mach-O 64-bit executable x86_64
aracne2.jar: Java executable jarfile
ARACNE.src.tar.gz: C++ source only with Makefiles
ARACNE-java.src.tar.gz: Java source only with ANT build file
aracne.zip: Java Graphic User Interface (GUI) for loading adjacency matrices and drawing network diagrams using a built-in Cytoscape plugin. Please set a 'JAVA_HOME' environment variable pointing to your JDK and use the launch_aracne scripts in the distribution to start the application
See http://www.nature.com/nprot/journal/v1/n2/abs/nprot.2006.106.html or Download PDF and Supplemental Documents for detailed usage instructions.
aracne [OPTIONS] ... or
java -jar ARACNE-java.jar [OPTIONS]
-i <file> Input gene expression profile dataset
-o <file> Output file name (optional) [*]
-j <file> Existing adjacency matrix (.adj) file
-a <fixed_bandwidth|variable_bandwidth|adaptive_partitioning> Algorithm (fixed_bandwidth | variable_bandwidth | adaptive_partitioning), default: adaptive_partitioning
-k <kernel width> Kernel width (accurate method only), default: determined by program
-b <# bins> No. of bins (fast method only), default: 6
-t <threshold> MI threshold, default: 0-p
P-value for MI threshold (e.g. 1e-7), default: 1 [**] -e <tolerance> DPI tolerance, default: 1 -h <probeId> Hub gene (only MI w/ hub gene will be computed), default: NONE -r <sample number> Use resampling arrays -s <file> A file containing a list of probes for which a subnetwork will be constructed, default: NONE -l <file> A file containing a list of probes annotated as transcription factors in the input dataset, default: NONE [***] -c <+/-probeId %> Conditional network reconstruction, default: NONE [****] [format: "+24 0.35", "-1973_s_at 0.4"] -f <mean> <cv> Gene filter by the mean and coefficient of variance (cv) of the expression values, default: mean=0, cv=0 -H <ARACNE_HOME> To specify where the ARACNE configuration files locates, default: current working directory --help Display this help and exit [*] If no output file is specified by the user, an output will be automatically generated in the same directory as the input file by appending some of the parameter values, such as kernel width, MI threshold, tolerance and so on, at the end of the input file name, and changing the file extension to ".adj". [**] If the "-t" option is supplied, it will enforce the program to use the specified MI threshold, therefore the "-p" option will be ignored. Otherwise, the program will automatically determines the MI threshold given the p-value. The default, p-value=1, will preserve all pairwise MI. [***] This option is ideal for transcriptional network reconstruction. If provided, DPI will not remove any connection of a transcription factor (TF) by connections between two probes not annotated as TFs. This option is often used in conjunction with '-s', which specifies a list of probes that are either the same or a subset of the probes specified by '-l'. [****] Conditional network reconstructs the network given a specified probe being most expressed or least expressed. In the format that follows "-c", "probeId" indicate the probe to be conditioned on; "+" or "-" specify whether the upper or lower tail of the probe's expression should be used as the condition, and "%" is a percentage between (0, 1) specifying the proportion of samples used as the conditioning subset. Example useage:"-c +24 0.35", "-c -1973_s_at 0.4
1) Reverse engineering cellular networks. Nature Protocols 1, 662 - 671 (2006). (Download PDF) (Supplemental Documents)
2) Reverse engineering of regulatory networks in human B cells. Nature Genetics. 2005 Apr;37(4):382-90. (Download PDF)
3) ARACNE: An Algorithm for the Reconstruction of Gene Regulatory Networks in a Mammalian Cellular Context. In press in BMC Bioinformatic. (Download PDF)
4) On The Reconstruction of Interaction Networks with Applications to Transcriptional Regulation.
http://arxiv.org/abs/q-bio.MN/0410036. Accepted in NIPS 2005
5) Conditional Network Analysis Identifies Candidate Regulator Genes in Human B Cells.http://arxiv.org/abs/q-bio/?0411003. Submitted to RECOMB 2005