Part of Learn. For the methods behind each step, read the free book Epigenomics Data Analysis with EpiNexus.
Welcome to EpiNexus! This guide will walk you through using the web interface to analyse your ChIP-seq, histone modification, or other epigenomic data — no command-line or coding experience needed.
Drag and drop your FASTQ.gz files — or BAM files with their peak files — or import public data from GEO, SRA or ENA, and choose your reference genome.
A full epigenomic pipeline runs automatically — from raw QC to publication-ready results.
Browse volcano plots, heatmaps, super-enhancer rankings, a genome browser and more — interactive, with tables and the main figures ready to download.
See mapping rates, duplication, FRiP scores and more for every sample.
Find out which genes and genomic features your peaks overlap.
Identify peaks that are gained or lost between treatment and control.
Compare H3K27ac and H3K27me3 changes gene by gene: activated, repressed, poised, or stable.
Detect large regulatory regions that drive cell identity genes.
Connect enhancer peaks to their likely target genes.
Bring in GEO, SRA or ENA samples by accession, straight from the archive.
See which transcription-factor motifs and biological processes your regions point to.
Find differentially methylated regions in WGBS or RRBS data, with bisulfite QC.
You will need your data files ready to upload. EpiNexus accepts:
Raw sequencing reads straight from the sequencer. If your data hasn't been aligned yet, start here — EpiNexus handles everything from raw reads onward.
Aligned sequencing reads (already mapped to a genome), with the peaks called from them (.narrowPeak or .broadPeak). EpiNexus calls peaks only from FASTQ files, so a BAM-based analysis needs its peak files too.
Samples published in GEO, SRA or ENA can be imported by accession instead of uploaded — the reads go straight from the archive into your project.
Bisulfite FASTQ files (WGBS or RRBS), or Bismark coverage (.cov) files if you have already aligned them.
Once you have your files, follow the steps in the sidebar — starting with Create a Project.
A project is a container for one experiment — it groups your samples, analysis runs, and results together in one place.
After signing in at app.epinexus.io, click All Projects in the sidebar. You will see your existing projects, or an empty page if this is your first time.
The New Project button is at the top of the page.
Project Name — Give it a descriptive name, e.g. “H3K27ac Treatment vs Control - March 2026”.
Description (optional) — Add notes about the experiment for your records.
Reference Genome — The list shows the genomes installed on the service, for example:
All 17 configured genomes are listed on the Pipelines page; if the one you need isn't offered, ask us. Not sure which genome? Check with whoever aligned your data, or look at the BAM file header. Most recent human data uses hg38.
Your project is ready. Its page shows where you are — Setup, Pipeline, Results — and what to do next.
Keep things organised by creating separate projects for different experiments. For example, one project for “H3K27ac in liver cancer” and another for “H3K4me3 in normal liver”.
After creating a project, bring in your sequencing data — uploaded from your computer or imported from a public archive — and describe your samples. EpiNexus handles all the bioinformatics processing from there.
In your project, click Upload Data in the sidebar (under Setup). It has two tabs: From my computer and From a public repository (GEO / SRA / ENA).
Drag and drop files, or click Browse Files. Depending on what you have:
FASTQ files (.fastq.gz) — Raw reads from the sequencer. Upload both read 1 and read 2 files for paired-end data.
BAM files — Aligned reads, one file per sample (e.g. “Treatment_Rep1.bam”), together with their peak files.
Files upload in the background in 50 MB pieces and pick up where they stopped if the connection drops. When a file is complete it shows Transfer verified, plus a format check for FASTQ, BAM and peak files.
EpiNexus calls peaks only from FASTQ files, so a BAM-based analysis uses the peaks you provide. Name each peak file after its BAM — e.g. “Treatment_Rep1.narrowPeak” (or .broadPeak) beside “Treatment_Rep1.bam” — and EpiNexus pairs them automatically.
On the From a public repository tab, type one or more accessions (GSE, GSM, SRR/SRX/SRP, PRJNA…) and click Look up. Check the table — target, condition and replicate are filled in from the archive's records; correct anything that's wrong — and click Import. The files go straight from the archive to your project and are checked against the archive's MD5 checksums; nothing is downloaded to your computer.
Open Samples in the sidebar. Each sample needs a name, its target (the histone mark or factor, or input/IgG for a control), its condition, its replicate number, and its files. Use Add Sample, or Import CSV with the columns SampleID, fastq_r1, fastq_r2, Factor, Condition and Replicate. Imported public samples are registered for you. Note that importing a CSV replaces all the samples in the project.
| File type | Typical size | What it contains |
|---|---|---|
| FASTQ | 5 – 30 GB | Raw sequencing reads straight from the sequencer |
| BAM | 1 – 10 GB | Aligned sequencing reads mapped to the genome |
| narrowPeak | 0.5 – 5 MB | Regions where your mark or factor is enriched (sharp peaks) |
| broadPeak | 1 – 20 MB | Regions where your mark is enriched (broad domains) |
Each account has 100 GB of storage; deleting files or projects frees it again.
BAM and FASTQ files can be very large, so a stable connection helps. If an upload is interrupted, it resumes from where it stopped — you can even close the browser and resume later.
Once your samples are registered and their files uploaded, you set up a run in four short steps. EpiNexus chooses the analysis steps that fit your assay and data.
Click Run Pipeline in the sidebar (or New Run on the project page).
Choose your assay — ChIP-seq, CUT&Tag, CUT&RUN, ATAC-seq or DNA Methylation — and the pipeline:
For FASTQ files: alignment, peak calling, QC and every downstream analysis — everything from raw reads to results.
For BAM files with their peak files: QC, annotation, differential analysis and the other downstream steps.
Tick the samples to include and check the reference genome. For DNA methylation, put each sample in a group (e.g. treated, control) and choose the reference group.
Every pair of conditions is compared by default — untick the comparisons you don't need. The other defaults suit most experiments; see Settings & Options for what you can change.
EpiNexus checks your design before anything runs — controls, replicates, the genome, and for methylation the library type — and lists anything to fix. Give the run a name if you like, then click Launch Pipeline. The run executes on the server: you can close the browser and come back later.
The run page shows each step as it starts and finishes, with its duration and any warnings; the Logs tab has the details, and you can cancel a run while it is going. Runs that start from FASTQ files take several hours, mostly for alignment; runs that start from BAM files are shorter. If alignment fails, the run stops; if a later step fails, the other steps still finish and the run page says what went wrong.
When a run completes, click View analysis on the run page, or open any page under Analysis in the sidebar. Each page has a Viewing results from picker to switch between runs. The run page's Results tab lists every output file, and links to the MultiQC report.
Metrics for each sample — reads, mapping rate, duplication, FRiP, library complexity, fragment sizes, TSS enrichment — graded Good, Marginal or Poor against published thresholds, with a plain-language interpretation. DNA-methylation runs show bisulfite QC here.
Where your peaks fall — promoter, UTR, exon, intron, distal intergenic — their distance to the nearest TSS and their nearest gene, with a sortable table.
Peaks that change significantly between conditions: volcano and MA plots, PCA, heatmaps, and a table with fold change, FDR and direction (gained or lost). Download every peak as a CSV, or open any peak in the genome browser.
Tornado heatmaps and average profiles of each sample's signal around peaks, and a sample-correlation matrix (runs from FASTQ files).
Signal tracks, peaks and reads in an IGV browser, with a link to share, an IGV session for the desktop app, and a button to open the tracks in UCSC.
Gene Ontology terms enriched near gained or lost regions, or near each sample's peaks (GREAT-style), by biological process, molecular function and cellular component.
For H3K27ac, MED1 or BRD4 samples: the ROSE ranking (hockey-stick plot) and a table of super-enhancers with their location, signal and rank.
When the project has differential results for both H3K27ac and H3K27me3: each gene classed as activated, repressed, both gained, poised or stable, a diagram of how the two marks changed, and a gene table.
Transcription-factor motifs from JASPAR 2024 enriched in your peaks (MEME AME), with sequence logos to compare them.
Each enhancer's predicted target genes from the Activity-by-Contact (ABC) model, and the genes with the most enhancers.
Fragment-size distributions and regions that become more open or more closed between conditions.
Differentially methylated regions (dmrseq) with their nearest gene: the CpGs tested, mean methylation per group, a volcano plot, a DMR table, and bisulfite QC with M-bias plots.
The sidebar shows the pages that fit your project's assays — for example, a DNA-methylation project shows Quality Control and Methylation.
Tables download as tab-separated files (they open in Excel); the main figures download as PNG or SVG. Under Output in the sidebar, Reports builds a report with figures and a methods section (HTML or Markdown, or PDF from your browser), and Export downloads a zip of the run's peaks, tables, signal tracks and alignments, a GEO submission package, and the files to reproduce the run (parameters, commands and methods).
Not sure what these results mean? Here is a plain-language guide to the key concepts.
A peak is a region of the genome where your histone mark is enriched compared to the background. Think of it as a “hotspot” of epigenetic activity. Each peak has a location (chromosome, start, end) and a signal strength.
Fraction of Reads in Peaks — what percentage of your sequencing reads land inside called peaks. Higher is better, and what counts as good depends on the assay. EpiNexus grades ChIP-seq as good from 5% and failing below 1% (ENCODE), ATAC-seq as good from 30% and failing below 20% (ENCODE), and CUT&Tag and CUT&RUN as good from 30% and failing below 10% (EpiNexus's own guide; no published standard exists). Low FRiP may mean low enrichment or too much background.
A differential peak is one that is significantly more (gained) or less (lost) enriched in one condition versus another. For example, an H3K27ac peak that appears only in treated cells but not in control cells is a “gained” peak — suggesting that genomic region became activated by the treatment.
False Discovery Rate — the expected share of false positives among the results called significant. With an FDR cut-off of 0.1, about 10% of the peaks called differential are expected to be false positives. Lower FDR = more confident. EpiNexus calls a peak differential when its FDR is below 0.1 and its fold change is more than √2 (|log2 fold change| > 0.5).
Super-enhancers are clusters of nearby enhancers that have unusually high levels of H3K27ac (or of the co-activators MED1 and BRD4). They control important cell-identity genes and are often disrupted in cancer. They are identified by ranking all stitched enhancer regions by signal and finding the “inflection point” above which regions are exceptionally strong.
TSS = Transcription Start Site, the beginning of a gene. The distance from a peak to the nearest TSS helps you understand whether a peak is at a promoter (close to TSS, within ~3 kb) or at a distal enhancer (far from TSS, often 10–1000 kb away).
The defaults suit most experiments and are set for your assay. These are the options on the Parameters step when you launch a run:
| Setting | Default | What it does |
|---|---|---|
| Differential comparisons | Every pair of conditions | Which conditions to compare, and which side is the reference (control-like names such as “control”, “WT” or “DMSO” are taken as the reference). Untick all to skip differential analysis. |
| Paired design | Off | Tick when samples with the same replicate number were processed together (e.g. littermates); the model then accounts for the pairing. |
| Peak caller | MACS3 (ChIP-seq, ATAC-seq); SEACR (CUT&Tag, CUT&RUN) | GoPeaks is also available for paired-end CUT&Tag and CUT&RUN. |
| Normalization | TMM (ChIP-seq), library size (ATAC-seq), RLE (CUT&Tag, CUT&RUN) | How samples are scaled before comparing them. These methods assume most regions don't change; when a mark changes globally, use spike-in calibration. |
| Q-value threshold | 0.05 (ChIP-seq, ATAC-seq), 0.01 (CUT&Tag, CUT&RUN) | The peak-calling significance cut-off. |
| Spike-in normalization | On with E. coli (CUT&Tag, CUT&RUN); off otherwise | For libraries containing DNA or chromatin from another species. Each read is assigned to the genome it matches best, and samples are scaled by their spike-in reads. |
| Remove duplicates from targets | Off | CUT&Tag and CUT&RUN only: identical fragments are expected signal there, so duplicates are kept in target samples (IgG controls are always deduplicated). |
| Control depth scaling | MACS default | ChIP-seq only: scale a shallow input up instead of the deeper library down. |
| Samples that fail QC | Keep and flag | Keep failing samples and flag the results, or exclude them from the statistics. |
| Consensus peaks for differential | Each sample's peaks | Or only the peaks that reproduce between replicates. |
| Library type | Detect from the reads | DNA methylation only: RRBS or WGBS. RRBS reads start at an MspI site, which EpiNexus checks at the Review step; RRBS is never deduplicated. |
These follow published practice and are the same for every run:
| What | Value |
|---|---|
| Differential peaks | FDR < 0.1 and |log2 fold change| > 0.5 |
| Narrow or broad peaks | Chosen per mark (e.g. H3K27ac and H3K4me3 narrow; H3K27me3, H3K9me3, H3K36me3 and H3K4me1 broad) |
| Promoter window | ± 3 kb around a TSS |
| Super-enhancer stitching | 12.5 kb, with peaks within ± 2.5 kb of a TSS left out (ROSE) |
| Peak–gene links | Up to 1 Mb |
| DMRs (methylation) | CpGs with at least 4 reads in every sample; region FDR ≤ 0.05 (dmrseq) |
Most experiments don't need to. Untick comparisons you aren't interested in; tick the paired design when replicates were processed in pairs; turn on spike-in calibration if your libraries carry a spike-in, especially when the mark may change globally.
Open the run page: the steps table shows which step failed and why, and the Logs tab has the details. The most common causes are: (1) wrong genome selected — make sure your files match the chosen reference, (2) corrupted or truncated upload — try re-uploading, (3) very low quality data with no peaks detected. If none of these apply, reach out via the Contact page.
Large FASTQ.gz files can take a while to upload, depending on your connection. Uploads resume from where they stopped: re-open the Upload Data page and add the same file again. If the page reports that your storage is full, delete files or projects you no longer need (each account has 100 GB).
This usually means the data has very low signal-to-noise. Verify that (1) the correct assay type is selected, (2) you have enough sequencing depth (at least 10–20M reads for ChIP-seq), and (3) the antibody or enrichment step worked as expected. Check the QC report for library complexity and FRiP scores.
Make sure you have at least two conditions with replicates. If you only have one sample per condition, the statistical test cannot estimate variance. Also check that each sample's condition is right on the Samples page and that the comparison you want was ticked when you launched the run.
Minor differences in peak counts or differential calls are expected between tools, because of different thresholds, filters and statistics. EpiNexus uses MACS3 for ChIP-seq and ATAC-seq, SEACR for CUT&Tag and CUT&RUN, and DESeq2 statistics (csaw for broad marks); its fixed thresholds are listed under Settings & Options, and every run's methods text and recorded commands show exactly what was done.
Interactive plots with very large datasets (>100k peaks) can be demanding. Try using Chrome or Firefox with hardware acceleration enabled. Closing other tabs can also free up memory. If a specific plot is slow, use the built-in filters to reduce the data shown.
Open the run (or the project’s Runs page) and click Delete run. Any finished run — completed, failed or cancelled — can be deleted; cancel a run that is still going first. A deleted run disappears at once and can be restored from Recently deleted on the Runs page for 7 days; after that its result files are removed for good, and backup copies expire within a further 30 days.
Open All Projects, find the project, and click its trash icon. You will be asked to confirm. Deletion is permanent — all uploaded files, runs and results are removed, including recently deleted runs, and backup copies expire within 30 days. A project with a run still going can’t be deleted: cancel the run (or wait for it to finish) first.