Question: RNA-seq from Co-culture, can I filter out the reads that map to one species and assemble what is left?
0
gravatar for Rößti
2.2 years ago by
Rößti30
United Kingdom
Rößti30 wrote:

Hi,

I have recently performed RNA-seq on two bacterial species both on their own and in co-culture. Unfortunately for me, one of the species (A) doesn't yet have a published genome that it maps nicely to... I have used rnaSPAdes to assemble a reference for A from it's monoculture growth. What I wanted to know is whether there is a way that I could filter out the reads of the other species (B) in the co-culture condition, and assemble what is left. A grows much better when co-cultured, so I want to be sure that I am not missing anything, that isn't being expressed on its own!

I hope that makes sense

sequencing rna-seq • 996 views
ADD COMMENTlink modified 2.2 years ago by Jake Warner750 • written 2.2 years ago by Rößti30
1
gravatar for Devon Ryan
2.2 years ago by
Devon Ryan91k
Freiburg, Germany
Devon Ryan91k wrote:

If you have a decent genome or transcriptome for B and the genomes/transcriptomes are different enough then sure, align everything against B and then take the unaligned reads and assemble them. Many mappers have an option to write unaligned reads to a separate file, which would be convenient for this.

ADD COMMENTlink written 2.2 years ago by Devon Ryan91k
1
gravatar for genomax
2.2 years ago by
genomax71k
United States
genomax71k wrote:

You could give BBSplit from BBMap suite a try with the genome you have at hand. You are going to have trouble with genes that are shared/homologous with both species. BBSplit allows you multiple options. You can choose to assign the reads to both or collect them in an ambiguous pool that you can look at separately.

Seal.sh from BBMap suite is also able to separate reads based on k-mer profiles. If it helps you may need to be tried out.

ADD COMMENTlink modified 2.2 years ago • written 2.2 years ago by genomax71k
1
gravatar for Jake Warner
2.2 years ago by
Jake Warner750
Jake Warner750 wrote:

You could map the mixed AB reads to B's genome and keep the unaligned reads. This approach is obviously dependent on the stringency of the mapping and the similarity of their genomes...

Here's an example:

bowtie2 --phred33 --very-sensitive --threads 8 --time -a -x B_genome.fa \
-1 AB_R1.fastq -2 AB_R2.fastq \
--un-conc A.fastq \
--al-conc B.fastq \
> /dev/null
ADD COMMENTlink written 2.2 years ago by Jake Warner750
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1178 users visited in the last hour