In this HiCExplorer tutorial we will generate and plot a Hi-C contact matrix. For this the following steps are necessary to be performed:

1. Map the Hi-C reads to the reference genome
2. Creation of a Hi-C matrix
3. Plotting the Hi-C matrix
4. Correction of Hi-C matrix
5. TAD Calling
6. A/B compartments computation
7. pyGenomeTracks visualization
8. Loop detection

After a corrected Hi-C matrix is created other tools can be used to visualize it, call TADS or compare it with other matrices.

Agenda

In this tutorial, we will deal with:

1. Data upload
2. Reads mapping
3. Creation of a Hi-C matrix
4. Plotting the Hi-C matrix
5. Correction of Hi-C matrix
   1. Plotting the corrected Hi-C matrix
   2. Load new data
6. TAD calling
7. A/B compartments computation
8. Integrating Hi-C and other data
9. Loop detection
10. Conclusion

Data upload

Hands-on: Data upload

1. Create a new history

   Tip: Creating a new history

   To create a new history simply click the new-history icon at the top of the history panel:

   

2. Import from Zenodo.

   Tip: Importing data via links

   1. Copy the link location
   2. Open the Galaxy Upload Manager

   3. Select Paste/Fetch Data

      Below are the links to the read files that can be copied and pasted in the upload manager.

      https://zenodo.org/record/1183661/files/HiC_S2_1p_10min_lowU_R1.fastq.gz
      https://zenodo.org/record/1183661/files/HiC_S2_1p_10min_lowU_R2.fastq.gz
      

   4. Paste the links into the text field
   5. Press Start

3. Rename the data set to something meaningful, e.g. HiC_S2_1p_10min_lowU_R1 and HiC_S2_1p_10min_lowU_R2. By default, when data is imported via its link, Galaxy names it with its URL.

Comment: Get data from public sources

HiCExplorer needs as input the forward and reverse strand of a pair end read which are mapped independently. A usual start point for a typical analysis is the given GSE number of a publication, e.g. GSE63525 for Rao 2014. To get the actual data, go to NCBI and search for the GSE number. In the section ‘Samples’ the GSM numbers of all samples are given. Select the correct one for you, and go to the European Nucleotide Archive and enter the GSM number. Select a matching result e.g. SRX764936 and download the data given in the row ‘FASTQ files (FTP)’ the forward and reverse strand. It is important to have the forward and reverse strand individual as a FASTQ file and to map it individually, HiCExplorer can not work with interleaved files.

Reads mapping

Mates have to be mapped individually to avoid mapper specific heuristics designed for standard paired-end libraries.

We have used the HiCExplorer successfully with bwa, bowtie2 and hisat2. In this tutorial we will be using Map with BWA-MEM tool. It is important to remember to:

• use local mapping, in contrast to end-to-end. A fraction of Hi-C reads are chimeric and will not map end-to-end thus, local mapping is important to increase the number of mapped reads
• tune the aligner parameters to penalize deletions and insertions. This is important to avoid aligned reads with gaps if they happen to be chimeric.
• If bowtie2 or hisat2 are used, --reorder option and as a file format bam_native needs to be used. Regular bam files are sorted by Galaxy and can not be used as an input for HiCExplorer.

Hands-on: Mapping reads

1. Map with Bowtie tool: Run Bowtie on both strands HiC_S2_1p_10min_lowU_R1 and HiC_S2_1p_10min_lowU_R2 with:
   • “Is this single or paired library” to Single-end
   • Set multiple data sets
   • “FASTQ file” to HiC_S2_1p_10min_lowU_R1and HiC_S2_1p_10min_lowU_R2
   • “Will you select a reference genome from your history or use a built-in index?” to Use a built-in index
   • “Select a reference genome” to dm3
   • “Do you want to tweak SAM/BAM Options?” to Yes
   • “Reorder output to reflect order of the input file” to Yes
2. Rename the output of the tool according to the corresponding files: R1.bam and R2.bam

Creation of a Hi-C matrix

Once the reads have been mapped the Hi-C matrix can be built.

For this step we will use hicBuildMatrix tool, which builds the matrix of read counts over the bins in the genome, considering the sites around the given restriction site.

Hands-on: hicBuildMatrix

1. hicBuildMatrix tool: Run hicBuildMatrix on the R1.bam and R2.bam from previous step with modifying the following parameters:
   • “1: Sam/Bam files to process” to R1.bam
   • “2: Sam/Bam files to process” to R2.bam
   • “Choose to use a restriction cut file or a bin size” to Bin size
   • “Bin size in bp” to 10000

   • “Sequence of the restriction site” to GATC

     Comment

     hicBuildMatrix creates two files, a bam file containing only the valid Hi-C read pairs and a matrix containing the Hi-C contacts at the given resolution. The bam file is useful to check the quality of the Hi-C library on the genome browser. A good Hi-C library should contain piles of reads near the restriction fragment sites. In the QC folder a html file is saved with plots containing useful information for the quality control of the Hi-C sample like the number of valid pairs, duplicated pairs, self-ligations etc. Usually, only 25%-40% of the reads are valid and used to build the Hi-C matrix mostly because of the reads that are on repetitive regions that need to be discarded.

     Comment

     Normally 25% of the total reads are selected. The output matrices have counts for the genomic regions. The extension of output matrix files is .h5. A quality report is created in e.g. hicMatrix/R1_10kb_QC, have a look at the report hicQC.html.

2. Rename the output to 10 kb contact matrix.

Plotting the Hi-C matrix

A 10kb bin matrix is too large to plot, it’s better to reduce the resolution. We usually run out of memory for a 1 kb or a 10 kb bin matrix and the time to plot it is very long (minutes instead of seconds). In order to reduce the resolution we use the tool hicMergeMatrixBins.

hicMergeMatrixBins merges the bins into larger bins of given number (specified by –numBins). We will merge 100 bins in the original (uncorrected) matrix and then correct it. The new bin size is going to be 10.000 bp * 100 = 1.000.000 bp = 1 Mb

Hands-on: hicMergeMatrixBins

1. hicMergeMatrixBins tool: Run hicMergeMatrixBins on the output from previous step setting the following parameter:
   • “Number of bins to merge” to 100

2. Rename the output to 1 MB contact matrix.

3. hicPlotMatrix tool: Run hicPlotMatrix on the output from hicMergeMatrixBins 1 MB contact matrix adjusting the parameters:
   • “Plot title” to Hi-C matrix for dm3
   • “Remove masked bins from the matrix” to True
   • “Plot the log1p of the matrix values: True
   • “Chromosomes to include (and order to plot in)” to chr2L
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr2R
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr3L
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr3R
   • ”+ Insert Chromosomes to include (and order to plot in):” to chrX

   Tip: log1p

   Because of the large differences in counts found in the matrix, it is better to plot the counts using the –log1p option.

The resulting plot of the 1 Mb contact matrix should look like:

Correction of Hi-C matrix

hicCorrectMatrix corrects the matrix counts in an iterative manner. For correcting the matrix, it’s important to remove the unassembled scaffolds (e.g. NT_) and keep only chromosomes, as scaffolds create problems with matrix correction. Therefore we use the chromosome names (chr2R, chr2L, chr3R, chr3L, chrX) here.

Matrix correction works in two steps: first a histogram containing the sum of contact per bin (row sum) is produced. This plot needs to be inspected to decide the best threshold for removing bins with lower number of reads. The second steps removes the low scoring bins and does the correction.

Hands-on: Matrix diagnostic

1. hicCorrectMatrix tool: Run hicCorrectMatrix on the output from hicBuildMatrix 10 kb contact matrix adjusting the parameters:
   • “Range restriction (in bp)” to Diagnostic plot
   • “Chromosomes to include (and order to plot in)” to chr2L
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr2R
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr3L
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr3R
   • ”+ Insert Chromosomes to include (and order to plot in):” to chrX

The output of the program prints a threshold suggestion that is usually accurate but is better to revise the histogram plot. The threshold is visualized in the plot as a black vertical line.

In our case the distribution describes the counts per bin of a genomic distance. To remove all bins with a z-score threshold less / more than X means to remove all bins which have less / more counts than X of mean of their specific distribution in units of the standard deviation. Looking at the distribution, we can select the value of -1.6 (lower end) and 1.8 (upper end) to remove. This is given by the –filterThreshold option in hicCorrectMatrix set to ‘correct matrix’ mode.

Hands-on: Matrix correction

1. hicCorrectMatrix tool: Run hicCorrectMatrix on the original matrix 10 kb contact matrix adjusting the parameters:
   • “Range restriction (in bp)” to Correct matrix
   • “Normalize each chromosome separately” to True
   • “Remove bins of low coverage” to -1.6
   • “Remove bins of large coverage” to 1.8
   • “Chromosomes to include (and order to plot in)” to chr2L
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr2R
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr3L
   • ”+ Insert Chromosomes to include (and order to plot in):” to chr3R
   • ”+ Insert Chromosomes to include (and order to plot in):” to chrX
2. Rename the corrected matrix to 10 kb corrected contact matrix.

It can happen that the correction stops with:

ERROR:iterative correction:*Error* matrix correction produced extremely large values.
This is often caused by bins of low counts. Use a more stringent filtering of bins.

This can be solved by a more stringent z-score values for the filter threshold or by a look at the plotted matrix. For example, chromosomes with 0 reads in its bins can be excluded from the correction by not defining it for the set of chromosomes that should be corrected (parameter ‘Include chromosomes’).

Plotting the corrected Hi-C matrix

We can now plot chromosome 2L with the corrected matrix.

Hands-on: Plotting the corrected Hi-C matrix

1. hicPlotMatrix tool: Run hicPlotMatrix on 10 kb corrected contact matrix adjusting the parameters:
   • “Plot title” to Hi-C matrix for dm3
   • “Plot per chromosome” to False
   • “Plot only this region” to chr2L
   • “Plot the log1p of the matrix values” to True

Load new data

The steps so far would have led to long run times if real data would have been used. We therefore prepared a new matrix for you, corrected contact matrix dm3 large. Please load it into your history.

TAD calling

“The partitioning of chromosomes into topologically associating domains (TADs) is an emerging concept that is reshaping our understanding of gene regulation in the context of physical organization of the genome” (Ramírez et al. 2017).

TAD calling works in two steps: First HiCExplorer computes a TAD-separation score based on a z-score matrix for all bins. Then those bins having a local minimum of the TAD-separation score are evaluated with respect to the surrounding bins to assign a p-value. Then a cutoff is applied to select the bins more likely to be TAD boundaries.

hicFindTADs tries to identify sensible parameters but those can be change to identify more stringent set of boundaries.

Hands-on: Finding TADs

1. hicFindTADs tool: Run hicFindTADs on corrected contact matrix dm3 large adjusting the parameters:
   • “Minimum window length (in bp) to be considered to the left and to the right of each Hi-C bin.” to 30000
   • “Maximum window length (in bp) to be considered to the left and to the right of each Hi-C bin.” to 100000
   • “Step size when moving from minDepth to maxDepth” to 10000
   • “Multiple Testing Corrections” to False discovery rate
   • “q-value” to 0.05
   • “Minimum threshold of the difference between the TAD-separation score of a putative boundary and the mean of the TAD-sep. score of surrounding bins.” to 0.001

As an output we get the boundaries, domains and scores separated files. We will use in the plot later only the TAD-score file.

A/B compartments computation

Hands-on: Computing A / B compartments

1. hicPCA tool: Run hicPCA adjusting the parameters:
   • “Matrix to compute on” to corrected contact matrix dm3 large
   • “Output file format” to bigwig
   • “Return internally used Pearson matrix” to Yes

Hands-on: Plotting the pearson matrix and PCA track

1. hicPlotMatrix tool: Run hicPlotMatrix on pearson_matrix from PCA computation adjusting the parameters:
   • “Plot title” to Pearson matrix and PC1
   • “Chromosomes to include” to chr2L
   • “Color map to use for the heatmap” to gist_heat
   • “Datatype of eigenvector file” to bigwig
   • “Eigenvector file” to hicPCA on [...] PC1
2. hicPlotMatrix tool: Run hicPlotMatrix on pearson_matrix from PCA computation adjusting the parameters:
   • “Plot title” to Pearson matrix and PC2
   • “Chromosomes to include” to chr2L
   • “Color map to use for the heatmap” to gist_heat
   • “Datatype of eigenvector file” to bigwig
   • “Eigenvector file” to hicPCA on [...] PC2

The first principal component correlates with the chromosome arms, while the second component correlates with A/B compartments.

Integrating Hi-C and other data

We can plot the TADs for a given chromosomal region. For this we will use pyGenomeTracks.

For the next step we need additional data tracks. Please load dm3_genes.bed, H3K27me3.bw, H3K36me3.bw and H4K16ac.bw to your history.

Hands-on: Plotting TADs

1. pyGenomeTracks tool: Run pyGenomeTracks adjusting the parameters:
   • “Region of the genome to limit the operation” to chr2L:14500000-16500000
   • “Choose style of the track” to TAD visualization
     • “Plot title” to HiC dm3 chr2L:14500000-16500000
     • “Matrix to compute on.” to the corrected matrix from hicCorrectMatrix step
     • “Depth” to 750000
     • “Height” to 4
     • “Boundaries file” to hicFindTads on data XX: TAD domains
     • “Show x labels” to Yes
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Bedgraph matrix track
     • “Plot title” to TAD separation score
     • “Track file bedgraph format” to hicFindTads on data XX: TAD information in bm file
     • “Height” to 4
     • “Set type to lines” to ‘True’
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Bigwig track
     • “Plot title” to PC1
     • “Track file bigwig format” the first computed hicPCA result
     • “Height” to 1.5
     • “Color of track” to a color of your choice
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Bigwig track
     • “Plot title” to PC2
     • “Track file bigwig format” the second computed hicPCA result
     • “Height” to 1.5
     • “Color of track” to a color of your choice
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Bigwig track
     • “Plot title” to H3K36me3
     • “Track file bigwig format” to H3K36me3
     • “Height” to 1.5
     • “Color of track” to a color of your choice
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Bigwig track
     • “Plot title” to H3K27me3
     • “Track file bigwig format” to H3K27me3
     • “Height” to 1.5
     • “Color of track” to a color of your choice
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Bigwig track
     • “Plot title” to H4K16ac
     • “Track file bigwig format” to H4K16ac
     • “Height” to 1.5
     • “Color of track” to a color of your choice
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Gene track
     • “Plot title” to dm3 genes
     • “Track file bedgraph format” the imported .bed file
     • “Height” to 3
     • “Type” to genes
     • “Gene rows” to 15
     • “Color of track” to a color of your choice
   • “+Insert Include tracks in your plot”
     • “Choose style of the track” to Vlines track
     • “Track file bed format” to hicFindTads on data XX: TAD domains

The resulting image should look like this one:

Loop detection

In Hi-C data, the term loop refers to a 3D structure which represents enhancer-promoter, gene, architectural or polycomb-mediated interactions. These interactions have the characteristics to be enriched in a single region compared to the local background. These loops are also called long-range interactions with an expected maximum distance of 2 MB (see Rao et al. 2014).

To compute loops, we have to import a new data set from the shared library to our history: GM12878-MboI-allreps-filtered.10kb.cool or download it via the FTP cool files server. (FTP links seems to be blocked, therefore the full URL: ftp://cooler.csail.mit.edu/coolers/hg19/Rao2014-GM12878-MboI-allreps-filtered.10kb.cool)

This dataset is from the human cell GM12878, mapped to hg19 and of 10 kb resolution. We use a new file because to detect loop structures the read coverage is required to be in the hunderts of million; this was not the case for the previous used drosophila dataset.

Hands-on: Matrix information

1. hicInfo tool: Run hicInfo adjusting the parameters:
   • “Select” Multiple datasets
   • “Matrix to compute on” to corrected contact matrix dm3 large and GM12878-MboI-allreps-filtered.10kb.cool

We now investigate the result of hicInfo and see that the new imported file is having 1.2 billion non-zero elements, while the drosophila Hi-C interaction matrix has around 12 million non-zero elements.

Hands-on: Computing loops

1. hicDetectLoops tool: Run hicDetectLoops adjusting the parameters:
   • “Matrix to compute on” to GM12878-MboI-allreps-filtered.10kb.cool
   • “Peak width” to 6
   • “Window size” to 10
   • “P-value preselection” to 0.01
   • “P-value” to 0.01
   • “Chromosomes to include” to chr1

The detection of the loops is based on a pre-selection of interactions, a p-value given a continuous negative binomial distribution over all interactions of a relative distance is computed. All interactions are filtered with a threshold (p-value preselection) to retrieve loop candidates. In a second step, the selected peak candidate is compared against its background using a Wilcoxon rank-sum test.

As an output we get a loop file containing the positions of both anchor points of the loop and the p-value of the used statistical test.

Hands-on: Plotting of loops

1. hicPlotMatrix tool: Run hicPlotMatrix adjusting the parameters:
   • “Matrix to compute on” to GM12878-MboI-allreps-filtered.10kb.cool
   • “Plot title” to Loops
   • “Plot only this region” to chr1:18000000-22000000
   • “Plot the log1p of the matrix values” to Yes
   • “Add detected loops” to Computed loops of the previous step
   • “DPI for image” to 300

Conclusion

In this tutorial we used HiCExplorer to analyze drosophila melanogaster cells. We mapped chimeric reads and created a contact matrix, to reduce noise this contact matrix was normalized. We showed how to visualize a contact matrix and how we can investigate topological associating domains and relate them to additional data like gene tracks. Moreover, we used a human Hi-C interaction matrix to compute loop structures.

To improve your learned skills we offer an additional tutorial based on mouse stem cells: following work.

You've Finished the Tutorial

Please also consider filling out the Feedback Form as well!

Key points

• Hi-C helps to investigate the 3D structure of the genome and to uncover folding principles of chromatin.

• In order to build a Hi-C contact matrix the reads have to be mapped to the reference genome.

• Based on a contact matrix open and closed chromatin (A/B compartments) and topological associating domains (TADs) can be computed.

Frequently Asked Questions

References

1. Rao, S. S. P., M. H. Huntley, N. C. Durand, E. K. Stamenova, I. D. Bochkov et al., 2014 A 3D Map of the Human Genome at Kilobase Resolution Reveals Principles of Chromatin Looping. Cell 159: 1665–1680. 10.1016/j.cell.2014.11.021
2. Ramírez, F., V. Bhardwaj, J. Villaveces, L. Arrigoni, B. A. Grüning et al., 2017 High-resolution TADs reveal DNA sequences underlying genome organization in flies. 10.1101/115063

Feedback

Did you use this material as an instructor? Feel free to give us feedback on how it went.
Did you use this material as a learner or student? Click the form below to leave feedback.

Citing this Tutorial

1. Joachim Wolff, Fidel Ramirez, Vivek Bhardwaj, Ekaterina Polkh, Hi-C analysis of Drosophila melanogaster cells using HiCExplorer (Galaxy Training Materials). https://training.galaxyproject.org/training-material/topics/epigenetics/tutorials/hicexplorer/tutorial.html Online; accessed TODAY
2. Hiltemann, Saskia, Rasche, Helena et al., 2023 Galaxy Training: A Powerful Framework for Teaching! PLOS Computational Biology 10.1371/journal.pcbi.1010752
3. Batut et al., 2018 Community-Driven Data Analysis Training for Biology Cell Systems 10.1016/j.cels.2018.05.012

BibTeX

@misc{epigenetics-hicexplorer,
author = "Joachim Wolff and Fidel Ramirez and Vivek Bhardwaj and Ekaterina Polkh",
	title = "Hi-C analysis of Drosophila melanogaster cells using HiCExplorer (Galaxy Training Materials)",
	year = "",
	month = "",
	day = "",
	url = "\url{https://training.galaxyproject.org/training-material/topics/epigenetics/tutorials/hicexplorer/tutorial.html}",
	note = "[Online; accessed TODAY]"
}
@article{Hiltemann_2023,
	doi = {10.1371/journal.pcbi.1010752},
	url = {https://doi.org/10.1371%2Fjournal.pcbi.1010752},
	year = 2023,
	month = {jan},
	publisher = {Public Library of Science ({PLoS})},
	volume = {19},
	number = {1},
	pages = {e1010752},
	author = {Saskia Hiltemann and Helena Rasche and Simon Gladman and Hans-Rudolf Hotz and Delphine Larivi{\`{e}}re and Daniel Blankenberg and Pratik D. Jagtap and Thomas Wollmann and Anthony Bretaudeau and Nadia Gou{\'{e}} and Timothy J. Griffin and Coline Royaux and Yvan Le Bras and Subina Mehta and Anna Syme and Frederik Coppens and Bert Droesbeke and Nicola Soranzo and Wendi Bacon and Fotis Psomopoulos and Crist{\'{o}}bal Gallardo-Alba and John Davis and Melanie Christine Föll and Matthias Fahrner and Maria A. Doyle and Beatriz Serrano-Solano and Anne Claire Fouilloux and Peter van Heusden and Wolfgang Maier and Dave Clements and Florian Heyl and Björn Grüning and B{\'{e}}r{\'{e}}nice Batut and},
	editor = {Francis Ouellette},
	title = {Galaxy Training: A powerful framework for teaching!},
	journal = {PLoS Comput Biol}
}

                   

Funding

These individuals or organisations provided funding support for the development of this resource

ELIXIR Europe

de.NBI

UFR

CRC992

Galaxy Administrators: Install the missing tools

You can use Ephemeris's shed-tools install command to install the tools used in this tutorial.

shed-tools install [-g GALAXY] [-a API_KEY] -t <(curl https://training.galaxyproject.org/training-material/api/topics/epigenetics/tutorials/hicexplorer/tutorial.json | jq .admin_install_yaml -r)

Alternatively you can copy and paste the following YAML

---
install_tool_dependencies: true
install_repository_dependencies: true
install_resolver_dependencies: true
tools:
- name: hicexplorer_hicbuildmatrix
  owner: bgruening
  revisions: 707f691c974c
  tool_panel_section_label: HiCExplorer
  tool_shed_url: https://toolshed.g2.bx.psu.edu/
- name: hicexplorer_hiccorrectmatrix
  owner: bgruening
  revisions: f7d344dacfeb
  tool_panel_section_label: HiCExplorer
  tool_shed_url: https://toolshed.g2.bx.psu.edu/
- name: hicexplorer_hicfindtads
  owner: bgruening
  revisions: 43332e51a712
  tool_panel_section_label: HiCExplorer
  tool_shed_url: https://toolshed.g2.bx.psu.edu/
- name: hicexplorer_hicmergematrixbins
  owner: bgruening
  revisions: 8926550196ce
  tool_panel_section_label: HiCExplorer
  tool_shed_url: https://toolshed.g2.bx.psu.edu/
- name: hicexplorer_hicpca
  owner: bgruening
  revisions: ea50f7bba657
  tool_panel_section_label: HiCExplorer
  tool_shed_url: https://toolshed.g2.bx.psu.edu/
- name: hicexplorer_hicplotmatrix
  owner: bgruening
  revisions: d398f5373c70
  tool_panel_section_label: HiCExplorer
  tool_shed_url: https://toolshed.g2.bx.psu.edu/
- name: hicexplorer_hicplottads
  owner: bgruening
  revisions: 94fff1a6a899
  tool_panel_section_label: Graph/Display Data
  tool_shed_url: https://toolshed.g2.bx.psu.edu/

Feedback

5 stars	2
4 stars	6
0 stars	1

February 2024

• 4 stars: Liked: Details Disliked: For someone trying without support, there are some mandatory fields in hicbuildmatrix that are not explained: 1. Need to add the hicfindrestsite step to generate the compulsory contact bed file 2. Need to add a dangling sequence for the restriction enzyme for hicbuildmatrix

January 2023

• 4 stars: Liked: The tutorial was simple to follow and was a great balance between guidance and letting do things myself Disliked: Some of the steps seem outdated with some parameter having changed names or being not available anymore. But it doesn't prevent the completion of the tutorial.
• 5 stars: Liked: I like the detailed instructions of this tool. Disliked: When I tred to follow the steps of hicBuildMatrix and ran it on usegalaxy. It showed that I need to provide BED file with all restriction cut places. Could you tell me how to choose to use the bin size or where I can find the BED file?

February 2021

• 4 stars: Liked: Covered all the steps Disliked: Parameters could have been explained

July 2020

• 4 stars: Liked: step wise explanation was very clear Disliked: comparison of two HiC datasets

December 2019

• 0 stars: Disliked: Note from @jennaj: Noticed mismatched tools across tuto components. The "Reads mapping" step description states "We have used the HiCExplorer successfully with bwa, bowtie2 and hisat2. In this tutorial we will be using Map with BWA-MEM tool." *However* the "Hands-on: Mapping reads" box has the mapping tool specified as "Map with Bowtie". The tool name doesn't fully match a Galaxy wrapped tool but looks as if it was intended to match "Map with Bowtie for Illumina" tool from some earlier tutorial revision, but the tool options/settings are actually for "Bowtie2" (tweak SAM/BAM output). The tuto workflow uses "Map with BWA-MEM (Galaxy Version 0.8.0)" which isn't available at EU or ORG (or that version is hidden in the tool panel + tool versions menu). --------- Punchline ... three different tools are mixed up, at the first step of the tuto after loading the initial fastq inputs. Probably should adjust to make all for either Bowtie2 or BWA-MEM using a version available at EU (so it can be run there). Be nice to have this work at (at least) one of the usegalaxy.* servers :) ORG doesn't include HiC tools. Will ticket this and whatever else is found after reviewing the remainder of steps.

February 2019

• 4 stars: Liked: Nothing bad, just I do not have sufficient background knowledge to comprehend everything. Nevertheless, very well-structured for a beginner to learn.

September 2018

• 5 stars: Liked: perfect step by step !!!! Disliked: maybe use human data ??