Summary

This page covers information on how to use the new 'Pawsey Bio - Ubuntu 22.04 - 2023-03' image for Nimbus. Instructions on how to choose this image when creating your instance can be found here. This Bio-image is created to cater to bioinformatics users who prefer to have their instances pre-installed with software, tools and datasets commonly used in the bioinformatics domain, including over 8000 Biocontainer tools.

Commonly used software

Instructions

On this page, we will only cover instructions for how to use access and use Biocontainers and reference genome data sets. For instructions for the software listed above, please see the software's original documentation page.

1. Biocontainers and Reference Genome data

CernVM-FS is a read-only file system that was developed by another supercomputing centre (Cern). It allows files such as container tools, reference datasets and other shared resources that are commonly used by many researchers to be accessed, added to, and updated in the one location. At Pawsey, we currently cache the Biocontainer tools and reference genome datasets that are on Galaxy Project's repositories. 

To use the Biocontainer tools, you can skip this step and proceed to the next section, 2. Using the Biocontainer tools.

To view these repositories, you can do the following:

1. List the Biocontainer tools repository:

   Note
   Note: It may take a minute or two to load the folders. When you have done it once, it will not take as long to show again.

   
   

   Tip
   The /cvmfs/singularity.galaxyproject.org/all subdirectory is where the entire list of 8000+ Biocontainers can be found, with the alphabetical subdirectories being symlinks to them.

   
   

   Code Block

   $ ls -la /cvmfs/singularity.galaxyproject.org total 140 drwxr-xr-x 33 cvmfs cvmfs 4096 Mar 25 2020 . drwxr-xr-x 3 cvmfs cvmfs 4096 May 12 2022 1 drwxr-xr-x 3 cvmfs cvmfs 4096 Jun 29 2020 2 drwxr-xr-x 4 cvmfs cvmfs 4096 Feb 10 07:09 3 drwxr-xr-x 24 cvmfs cvmfs 4096 Feb 23 2021 a drwxr-xr-x 4 cvmfs cvmfs 4096 Mar 7 19:06 all drwxr-xr-x 23 cvmfs cvmfs 4096 Sep 2 2022 b drwxr-xr-x 26 cvmfs cvmfs 4096 Feb 10 15:05 c drwxr-xr-x 22 cvmfs cvmfs 4096 Feb 10 15:05 d drwxr-xr-x 23 cvmfs cvmfs 4096 May 1 2020 e drwxr-xr-x 21 cvmfs cvmfs 4096 Feb 10 19:06 f drwxr-xr-x 25 cvmfs cvmfs 4096 Feb 5 2022 g drwxr-xr-x 19 cvmfs cvmfs 4096 Feb 10 21:21 h drwxr-xr-x 18 cvmfs cvmfs 4096 Feb 10 21:21 i drwxr-xr-x 14 cvmfs cvmfs 4096 Jul 10 2020 j drwxr-xr-x 21 cvmfs cvmfs 4096 Feb 25 2021 k drwxr-xr-x 16 cvmfs cvmfs 4096 Feb 10 22:54 l drwxr-xr-x 26 cvmfs cvmfs 4096 Feb 10 22:54 m drwxr-xr-x 20 cvmfs cvmfs 4096 Feb 11 05:13 n drwxr-xr-x 14 cvmfs cvmfs 4096 May 7 2020 o drwxr-xr-x 24 cvmfs cvmfs 4096 Aug 24 2021 p drwxr-xr-x 12 cvmfs cvmfs 4096 Feb 27 2021 q drwxr-xr-x 27 cvmfs cvmfs 4096 Jun 28 2022 r drwxr-xr-x 26 cvmfs cvmfs 4096 Apr 7 2020 s drwxr-xr-x 24 cvmfs cvmfs 4096 Feb 11 23:17 t drwxr-xr-x 13 cvmfs cvmfs 4096 Feb 11 23:17 u drwxr-xr-x 18 cvmfs cvmfs 4096 Feb 11 23:17 v drwxr-xr-x 17 cvmfs cvmfs 4096 Jul 28 2021 w drwxr-xr-x 16 cvmfs cvmfs 4096 Apr 7 2020 x drwxr-xr-x 4 cvmfs cvmfs 4096 Feb 28 2021 y drwxr-xr-x 10 cvmfs cvmfs 4096 Feb 11 23:17 z

   
   

2. List the reference genome sets and other data files:

   Code Block
   $ ls -la /cvmfs/data.galaxyproject.org/ total 14 drwxr-xr-x 4 cvmfs cvmfs 4096 Mar 31 2016 . -rw-r--r-- 1 cvmfs cvmfs 21 Oct 24 2018 .cvmfsdirtab drwxr-xr-x 210 cvmfs cvmfs 4096 Apr 21 2022 byhand drwxr-xr-x 18 cvmfs cvmfs 4096 Nov 24 2020 managed

   Please note that the data sets may not be comprehensive, and this service is not meant to replace your current methods for accessing public datasets.

3. To use these data files for your analyses, copy the absolute file path in your workflow/pipeline. For example, with the reference genome Hg38, the file can be found in the following location, and specifically under the seq sub directory:

   Code Block

   $ ls -la /cvmfs/data.galaxyproject.org/byhand/hg38 total 46 drwxrwxr-x 10 cvmfs cvmfs 4096 Apr 22 2016 . drwxr-xr-x 210 cvmfs cvmfs 4096 Apr 21 2022 .. -rw-r--r-- 1 cvmfs cvmfs 0 Apr 22 2016 .cvmfscatalog drwxrwxr-x 3 cvmfs cvmfs 4096 Jan 21 2015 download drwxrwxr-x 6 cvmfs cvmfs 4096 Jan 20 2015 hg38canon drwxrwxr-x 6 cvmfs cvmfs 4096 Jan 20 2015 hg38female drwxrwxr-x 6 cvmfs cvmfs 4096 Jan 20 2015 hg38full drwxrwxr-x 2 cvmfs cvmfs 4096 Mar 18 2014 liftOver drwxrwxr-x 2 cvmfs cvmfs 4096 Mar 18 2014 picard_index drwxrwxr-x 2 cvmfs cvmfs 4096 Mar 18 2014 sam_index drwxrwxr-x 2 cvmfs cvmfs 4096 Apr 1 2016 seq $ ls -la /cvmfs/data.galaxyproject.org/byhand/hg38/seq total 10108046 drwxrwxr-x 2 cvmfs cvmfs 4096 Apr 1 2016 . drwxrwxr-x 10 cvmfs cvmfs 4096 Apr 22 2016 .. -rw-rw-r-- 1 cvmfs cvmfs 136 Mar 18 2014 README -rw-rw-r-- 1 cvmfs cvmfs 835393456 Mar 18 2014 hg38.2bit lrwxrwxrwx 1 cvmfs cvmfs 11 May 17 2014 hg38.fa -> hg38full.fa -rw-r--r-- 1 cvmfs cvmfs 19327 Aug 24 2015 hg38.fa.fai -rw-rw-r-- 1 cvmfs cvmfs 3150052305 Mar 17 2014 hg38canon.fa -rw-rw-r-- 1 cvmfs cvmfs 3091680335 Mar 17 2014 hg38female.fa -rw-r--r-- 1 cvmfs cvmfs 757 Apr 1 2016 hg38female.fa.fai -rw-rw-r-- 1 cvmfs cvmfs 3273481150 Mar 18 2014 hg38full.fa

   
   

4. So the full absolute path for the Hg38 sequence file would be:

   Code Block
   /cvmfs/data.galaxyproject.org/byhand/hg38/seq/hg38full.fa

   
   

sudo apt-get autoremove cvmfs
sudo apt-get purge cvmfs
sudo rm -rf /etc/cvmfs/ 
cd /home/ubuntu
git clone https://github.com/PawseySC/Pawsey-CernVM-FS.git
cd Pawsey-CernVM-FS
sudo ./install-cvmfs.sh install

2. Using the Biocontainer tools 

Singularity-HPC (SHPC) is a software for container modules. In this Pawsey Bio - Ubuntu 22.04 - 2023-03 image, we have integrated the use of SHPC seamlessly with CernVM-FS. This means that you can easily access and use over 8000 Biocontainers (and up to the latest versions) without needing to understand container syntax.

1.  To search for versions and information on a particular tool, e.g. for cuttlefish, use shpc show:

   Code Block

   $ shpc show quay.io/biocontainers/cuttlefish url: https://biocontainers.pro/tools/cuttlefish maintainer: '@vsoch' description: shpc-registry automated BioContainers addition for cuttlefish latest: 2.2.0--hf1761c0_0: sha256:63cdd7778b144a37684ae53b8e760ed00852f3010aa79292b3f1a6a6470f0992 tags: 2.1.0--hf1761c0_0: sha256:aa009abd48c372125e060d39f49f1690be74b6dac276d451bf1cc4c847a914d6 2.1.1--hf1761c0_0: sha256:8bccced83dd6bbf87843cf08851563c812ef7c36afda2efbcf0d54f9102b913f 2.2.0--hf1761c0_0: sha256:63cdd7778b144a37684ae53b8e760ed00852f3010aa79292b3f1a6a6470f0992 docker: quay.io/biocontainers/cuttlefish aliases: cuttlefish: /usr/local/bin/cuttlefish

   
   

2. To check for availability and to load the tool, use the module avail and module load commands:

   Note
   Note that on first use, the tool might take a 30 seconds or so to run the command as the container is being accessed from the filesystem for the first time

   
   

   Code Block
   $ module avail cuttlefish ------------------------------- /home/ubuntu/singularity-hpc/modules ------------------------------- quay.io/biocontainers/cuttlefish/2.2.0--hf1761c0_0/module $ module load quay.io/biocontainers/cuttlefish/2.2.0--hf1761c0_0/module $ $ cuttlefish --version cuttlefish 2.2.0 Supported commands: `build`, `help`, `version`. Usage: cuttlefish build [options]

   
   

3. To check for the list of modules loaded:

   Note
   This list will be cleared whenever you log out of your instance. After logging back in, you will need to reload the module for it to be on the list for use.

   
   

   Code Block
   $ module list Currently Loaded Modules: 1) quay.io/biocontainers/cuttlefish/2.2.0--hf1761c0_0/module

   
   

4. To use another version of the tool available from the above shpc show list, use shpc to install the module, ensuring to use and keep the cvmfs path of the container:

   Code Block
   $ sudo shpc install quay.io/biocontainers/cuttlefish:2.1.1--hf1761c0_0 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.1.1--hf1761c0_0 --keep-path Module quay.io/biocontainers/cuttlefish:2.1.1--hf1761c0_0 was created.

   
   

5. Update the modules cache, then, run the same module commands to load the tool in this version:

   Code Block

   $ /usr/local/lmod/lmod/libexec/update_lmod_system_cache_files -d /opt/mData/cacheDir -t /opt/mData/cacheTS.txt /home/ubuntu/singularity-hpc/modules $ module avail cuttlefish ------------------------------- /home/ubuntu/singularity-hpc/modules ------------------------------- quay.io/biocontainers/cuttlefish/2.1.1--hf1761c0_0/module quay.io/biocontainers/cuttlefish/2.2.0--hf1761c0_0/module (L,D) Where: L: Module is loaded D: Default Module   $ module load quay.io/biocontainers/cuttlefish/2.1.1--hf1761c0_0/module The following have been reloaded with a version change: 1) quay.io/biocontainers/cuttlefish/2.2.0--hf1761c0_0/module => quay.io/biocontainers/cuttlefish/2.1.1--hf1761c0_0/module

    You will notice that the previous version of the tool is now swapped out for the version you just loaded

   Code Block
   $ module list   Currently Loaded Modules: 1) quay.io/biocontainers/cuttlefish/2.1.1--hf1761c0_0/module

   
   

3. Adding a local SHPC registry

If there are versions (usually older ones) of a Biocontainer tool that is present in the cvmfs repository but not on the shpc show list (i.e. the default recipe), you can create a local SHPC registry and add/update a recipe file for the Biocontainer tool.

1. Clone the remote SHPC-registry and add it as a local registry:

   Code Block
   $ cd /home/ubuntu $ git clone https://github.com/singularityhub/shpc-registry.git $ sudo shpc config add registry /home/ubuntu/shpc-registry/ Warning: Check with shpc config edit - ordering of list can change. Added registry to /home/ubuntu/shpc-registry/

   
   

2. Look up for available versions of the tool, e.g. cuttlefish:

   Code Block

   $ ls /cvmfs/singularity.galaxyproject.org/all/cuttlefish* /cvmfs/singularity.galaxyproject.org/all/cuttlefish:1.0.0--h2e03b76_0 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:1.0.0--h2e03b76_1 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.0.0--h95f258a_0 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.0.0--hf1761c0_1 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.1.0--hf1761c0_0 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.1.1--hf1761c0_0 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.2.0--hf1761c0_0

   
   

3. To add a different version to the recipe file for the tool, e.g. 2.0.0--hf1761c0_1:

   Code Block

   $ shpc add /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.0.0--hf1761c0_1 quay.io/biocontainers/cuttlefish:2.0.0--hf1761c0_1 --registry /home/ubuntu/shpc-registry/ quay.io/biocontainers/cuttlefish:2.0.0--hf1761c0_1 already exists and will be updated! Registry entry quay.io/biocontainers/cuttlefish was added! Before shpc install, edit: /home/ubuntu/shpc-registry/quay.io/biocontainers/cuttlefish/container.yaml

   
   

4. To install the tool:

   Warning

   If you are using the now deprecated 'Pawsey Bio - Ubuntu 20.04 - 2021-11' image, run the following steps first: Code Block cd /home/ubuntu git clone https://github.com/singularityhub/singularity-hpc.git mkdir /home/ubuntu/singularity-hpc/modules sudo shpc config set module_base /home/ubuntu/singularity-hpc/modules cat >> ~/.bashrc <<'EOF' module use /home/ubuntu/singularity-hpc/modules EOF source ~/.bashrc By setting your module_base to this new location, all new container modules will be installed to this path.

   
   

   Code Block
   $ sudo shpc install quay.io/biocontainers/cuttlefish:2.0.0--hf1761c0_1 /cvmfs/singularity.galaxyproject.org/all/cuttlefish:2.0.0--hf1761c0_1 --keep-path Module biocontainers/cuttlefish:2.0.0--hf1761c0_1 was created.

   
   

5. Now when you do a module avail, the newly installed 2.0.0--hf1761c0_1 version will be available:

   Code Block
   $ module avail cuttlefish   ------------------------------- /home/ubuntu/singularity-hpc/modules ------------------------------- quay.io/biocontainers/cuttlefish/2.0.0--hf1761c0_1/module

   
   

6. Since you have created your own local registry, shpc will default to your local registry whenever you do a look up with shpc show. To look up the full list of Biocontainer tools with the latest versions, you will need to add a flag to point to the remote (Github) shpc-registry in your search:

   Tip
   The shpc-registry is kept up-to-date with the latest versions of all Biocontainers on a nightly update.

   
   

   Code Block
   $ shpc show quay.io/biocontainers/cuttlefish --registry https://github.com/singularityhub/shpc-registry

   
   

4. Using Nextflow

Nextflow makes use of containers to run your workflows sequentially. Each step of your workflow is called a process. For each process, Nextflow pulls the appropriate container required for that step to run it. You can prevent Nextflow from pulling the container and using what is present on your instance to save time and space.

To do so, you would create an additional config file to point Nextflow to either 1) the paths of the containers on /cvmfs/singularity.galaxyproject.org, or 2) the module paths on your instance. Nextflow prioritises this custom config file above the default nextflow.config file(s), if present, in other directories for your workflow.

Nextflow pipelines for Bioinformatics

Nextflow has a repository of pipelines that are available through https://nf-co.re. These are becoming increasingly popular, as more peer-reviewed pipelines are added by the community. A couple of popular ones include:

• https://nf-co.re/rnaseq/3.10.1
• https://nf-co.re/sarek

Configuring Nextflow to use Biocontainers

Please note this is only an example of how you can configure your Nextflow workflow to use the containers available from your instance.

1. 1. Suppose you are using the nfcore/rnaseq pipeline. Note that in the main.nf for the FASTQC process, there are a few parameters for the tool. These will be over-written by the config file that you will create in the next step.

      Code Block

      $ cat rnaseq/modules/nf-core/fastqc/main.nf process FASTQC { tag "$meta.id" label 'process_medium' conda "bioconda::fastqc=0.11.9" container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ? 'https://depot.galaxyproject.org/singularity/fastqc:0.11.9--0' : 'quay.io/biocontainers/fastqc:0.11.9--0' }" input: tuple val(meta), path(reads) output: tuple val(meta), path("*.html"), emit: html tuple val(meta), path("*.zip") , emit: zip path "versions.yml" , emit: versions when: task.ext.when == null || task.ext.when script: def args = task.ext.args ?: '' def prefix = task.ext.prefix ?: "${meta.id}" // Make list of old name and new name pairs to use for renaming in the bash while loop def old_new_pairs = reads instanceof Path || reads.size() == 1 ? [[ reads, "${prefix}.${reads.extension}" ]] : reads.withIndex().collect { entry, index -> [ entry, "${prefix}_${index + 1}.${entry.extension}" ] } def rename_to = old_new_pairs*.join(' ').join(' ') def renamed_files = old_new_pairs.collect{ old_name, new_name -> new_name }.join(' ') """ printf "%s %s\\n" $rename_to | while read old_name new_name; do [ -f "\${new_name}" ] || ln -s \$old_name \$new_name done fastqc $args --threads $task.cpus $renamed_files cat <<-END_VERSIONS > versions.yml "${task.process}": fastqc: \$( fastqc --version | sed -e "s/FastQC v//g" ) END_VERSIONS """ stub: def prefix = task.ext.prefix ?: "${meta.id}" """ touch ${prefix}.html touch ${prefix}.zip cat <<-END_VERSIONS > versions.yml "${task.process}": fastqc: \$( fastqc --version | sed -e "s/FastQC v//g" ) END_VERSIONS """ }

      
      

   2. Nextflow prioritises a custom config file over any other config files or values defined in the workflow files. To ensure that Nextflow uses the existing container for fastqc, you would create and use a custom config file, choosing either of the two ways:

      Ui tabs
      Ui tab title CVMFS path Create a cvmfs_path.config file with the following: Then run the workflow with this config: Ui tab title Modules path Create a modules_path.config file with the following: Then run the workflow with this config:

      For the other processes in your workflow, you will need to also add to the corresponding *_path.config file to ensure that every process has the path for the existing container to run that part of the workflow from. More info on config files for processes can be found here: https://www.nextflow.io/docs/latest/config.html#scope-process