Question: How can I use BLAST to extract chloroplast sequences from DNA reads?
0
gravatar for AcademicDialysis
3.5 years ago by
United States
AcademicDialysis60 wrote:

I'm trying to extract the chloroplast sequences from my reads, as Whole Genome Sequencing was used to produce them.

This paper: http://www.sciencedirect.com/science/article/pii/S2214540013000169 mentions that to do this, they BLASTed their reads against all of the known genomes in the same family. For me, this family would be Fabacaea.

Does anyone know of a quicker way to do this besides manually downloading every FASTA file containing Fabacaea chloroplast sequences from NCBI? Or of a better way to extract chloroplast sequences from my reads? I do know that chloroplast DNA should be more abundant than other DNA because it is more highly repeated than nuclear or mitochondrial DNA.

Info about reads: 300bp average, paired-end reads from Illumina MiSeq

 

Thanks in advance!

ADD COMMENTlink modified 3.5 years ago by 5heikki7.8k • written 3.5 years ago by AcademicDialysis60
4
gravatar for Pierre Lindenbaum
3.5 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum114k wrote:

search NCBI for chloroplast + Fabacaea

http://www.ncbi.nlm.nih.gov/nuccore?term=%28%22chloroplast%22[Filter]%29%20AND%20%22fabaceae%22[Organism]

and download the sequences as fasta.

Index the fasta with `bwa index` and map with `bwa mem`

 

ADD COMMENTlink written 3.5 years ago by Pierre Lindenbaum114k

Since the closest reference is just in the same family, I don't think my consensus sequences would be very large. Should I just blast and then use those reads to do de novo assembly? Or should I still use bwa and just use all of those sequences as reference?

 

Thanks!

ADD REPLYlink written 3.5 years ago by AcademicDialysis60
3
gravatar for 5heikki
3.5 years ago by
5heikki7.8k
Finland
5heikki7.8k wrote:

Assuming the chloroplast genome differs from the host in GC% and codon usage, the quickest way would be to bin the reads based on tetramer frequencies.
 

ADD COMMENTlink written 3.5 years ago by 5heikki7.8k
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: 769 users visited in the last hour