I am using my university's HPC cluster to de novo assemble paired-end HiSeq 400 reads
I've used Trimmomatic for adapter trimming and quality-based trimming
I used BBmap's kmercountexact.sh to determine my preferred k-mer value for assembly
Now I want to use ABySS for my de novo asssembly - the HPC uses SLURM scheduler

    module load openmpi
    Module openmpi/2.0.1 loaded 
    module load abyss
    Module abyss/1.9.0 loaded 

    srun --partition=high --mem=24000 --time=12:00:00 --nodes=1 abyss-pe np=8 name=EthFoc-11_trimmomatic7_corr k=41 in='EthFoc-11_S285_L007_trimmomatic7_1P.00.0_0.cor.fastq EthFoc-11_S285_L007_trimmomatic7_2P.00.0_0.cor.fastq'

There is a lot of text in STDOUT - available at [http://txt.do/d61vl][1]
Briefly, in STDOUT the last lines that indicate the the run was still working read as follows:

    Assembling...
    0: Assembled 51879929 k-mer in 81891 contigs.
    Assembled 51879929 k-mer in 81891 contigs.
    Concatenating fasta files to EthFoc-11_trimmomatic7_corr-1.fa
    Concatenating fasta files to EthFoc-11_trimmomatic7_corr-bubbles.fa
    Done.

But soon after this in the STDOUT, the run terminates with the following error message:

    Concatenating fasta files to EthFoc-11_trimmomatic7_corr-1.fa
    error: `contigs-0.fa': No such file or directory
    make: *** [EthFoc-11_trimmomatic7_corr-1.fa] Error 1
    srun: error: c11-96: task 0: Exited with exit code 2

I went through several biostars posts on ABySS run errors, but I don't think I have a direct solution to my problem...
https://www.biostars.org/p/140610/#140656, https://www.biostars.org/p/243511/, https://www.biostars.org/p/168100/, https://www.biostars.org/p/231310/

Could it be a share access mis-configuration on the HPCC? - https://www.biostars.org/p/178246/

The files generated from this run are listed below:

-rw-rw-r--  1 aksrao aksrao   42 Aug 18 00:01 EthFoc-11_trimmomatic7_corr-1.dot

-rw-rw-r--  1 aksrao aksrao    0 Aug 17 23:16 EthFoc-11_trimmomatic7_corr-1.fa

-rw-rw-r--  1 aksrao aksrao 1.3M Aug 17 23:15 EthFoc-11_trimmomatic7_corr-bubbles.fa


I wanted to test whether removing MPI from the equation will allow the run to completion, I tried

    srun --partition=high --mem=24000 --time=12:00:00 --nodes=1 abyss-pe name=EthFoc-11_trimmomatic7_corr k=41 in='EthFoc-11_S285_L007_trimmomatic7_1P.00.0_0.cor.fastq EthFoc-11_S285_L007_trimmomatic7_2P.00.0_0.cor.fastq'

    srun: job 13800791 queued and waiting for resources
    srun: job 13800791 has been allocated resources

And the error I see happens almost right away, with the STDOUT looking as follows:

    abyss-filtergraph  --dot   -k41 -g EthFoc-11_trimmomatic7_corr-2.dot1 EthFoc-11_trimmomatic7_corr-1.dot EthFoc-11_trimmomatic7_corr-1.fa >EthFoc-11_trimmomatic7_corr-1.path
    abyss-filtergraph: ../Graph/DotIO.h:302: std::istream& read_dot(std::istream&, Graph&, BetterEP) [with Graph = DirectedGraph<ContigProperties, Distance>; BetterEP = DisallowParallelEdges; std::istream = std::basic_istream<char>]: Assertion `num_vertices(g) > 0' failed.
    /bin/bash: line 1: 27510 Aborted                 abyss-filtergraph --dot -k41 -g EthFoc-11_trimmomatic7_corr-2.dot1 EthFoc-11_trimmomatic7_corr-1.dot EthFoc-11_trimmomatic7_corr-1.fa > EthFoc-11_trimmomatic7_corr-1.path
    make: *** [EthFoc-11_trimmomatic7_corr-1.path] Error 134
    make: *** Deleting file `EthFoc-11_trimmomatic7_corr-1.path'
    srun: error: c11-91: task 0: Exited with exit code 2

What am I doing wrong? And how can I fix it? Since I am only a couple of days new to genome assembly and an hour into ABySS use, the more detailed your reply, the more useful it might be for me. Thanks!


  [1]: http://txt.do/d61vl
Hi Anand,

I'm not sure what the problem is, but I can hopefully provide some hints.

It looks like your `abyss-pe` command is correct.  I suspect your problems are related your cluster job submission parameters.   Learning to run MPI jobs on a cluster usually requires a bit of experimenting with job submission flags.  If you have an IT department, you should ask them if they have any example scripts showing how to run MPI jobs on your cluster. Also, I would recommend first testing that you can successfully run a simple MPI program before trying ABySS.  For example, this page provides a MPI "Hello, World!" program: https://hpcc.usc.edu/support/documentation/examples-of-mpi-programs/.  You would have to paste the code into a file and compile it yourself, but it probably worth the effort.

In the log of a successful ABySS run, you should see multiple `ABYSS-P` processes ("ranks" in MPI terminology) running in parallel and communicating with each other.  The processes can be running on different cluster nodes. Each MPI process (rank) writes its own temporary `contig-<rank>.fa` file, so if you are running a job with 8 processes (np=8), you would expect to see the following in your assembly directory:

<pre>
contig-0.fa
contig-1.fa
contig-2.fa
contig-3.fa
contig-4.fa
contig-5.fa
contig-6.fa
contig-7.fa
</pre>  

When the parallel runs of ABYSS-P finish, these files are concatenated together into a single FASTA file and then removed.

At the beginning of a ABYSS-P log for np=8, you should see something like:
<pre>
0: Running on host c11-96
1: Running on host c11-96
2: Running on host c11-96
3: Running on host c11-96
4: Running on host c11-96
5: Running on host c11-96
6: Running on host c11-96
7: Running on host c11-96
</pre>

whereas in http://textuploader.com/d61vl, you are just seeing:
<pre>
0: Running on host c11-96
</pre>

in multiple runs of `ABYSS-P`.

It is appears that multiple independent ABySS jobs are being started rather than a single job with 8 processes ("ranks").   If they are all running in the same directory, they will overwrite each other's `contig-0.fa` files.  Also, each independent job will delete its `contig-0.fa` file once the concatenating step is finished, which is likely why you are seeing:

<pre>
error: `contigs-0.fa': No such file or directory
</pre>