How to get plasmid or costruct sequence from fastq file in c.elegans Whole-genome Sequencing?
2
0
Entering edit mode
4.3 years ago
GiV17 ▴ 40

Hi,

I have a fastq file generated from a whole genome sequencing by C. elegans samples and I want to find in which site (or sites) a construct/plasmid has been integrated in the genome, comparing my sequence with a published genome in Ensembl (WB235).

I had a C. elegans strain with a construct integrated in the genome. From the construct design, I expect it to integrate into chromosome III of C. Elegans. So, by using a tipical SNPs and indels workflow from Varscan, I managed to identify Variants in all samples. Furthermore, I also analyzed the CNVs through BreakDance tool only in chr III.

But my question is: How can I find where the construct is integrated? I would like the position (Start and End) respect to a genome reference. It's possible?

Thanks for all.

sequencing sequence R • 2.2k views
ADD COMMENT
1
Entering edit mode

If I were you, I think I would do something like below.

  1. Extract flanking sequences of plasmid sequences from whole genome sequencing fastq files.
  2. Blat the flanking sequences.
ADD REPLY
0
Entering edit mode

Sorry, but I didn't understand... I have the construct sequence in fasta format, but I would like to get How this sequence is integrated in several sample fastq files. I don't have a flanking sequences... so How can I do?

ADD REPLY
1
Entering edit mode

If the plasmids were integrated into the genome and you ran whole genome sequencing, certain reads will have hybrid sequences (part of plasmid + part of genomic sequence).

Then, you go though your reads possessing the plasmid sequences, and some of them will have the hybrid sequences which are the "flanking sequence" I mentioned.

Extract those reads. Remove plasmid sequences. Then, align the flanking sequences on genome to obtain the genomic coordination.

ADD REPLY
0
Entering edit mode

Hi mbk0asis, so, to sum up, you advise me: 1- To map fastq files on Construct and to extract unmapped Paired end reads 2- To take the unmapped reads and to remap respect to genome of C. elegans and so I'll obtain the genomic coordination of the unmapped? Thanks

ADD REPLY
0
Entering edit mode

Map the fastq on construct.

Extract mapped reads.

Trim off the plasmid sequences.

Map the remaining part of genomic sequence to reference genome.

That's my rough idea.

ADD REPLY
1
Entering edit mode
4.3 years ago
Joe 20k

Here's an approach which might work:

  1. Align all your reads in your fastq to your known plasmid sequence (you might need to experiment with some stringency).

  2. de novo assemble your remaining reads to see if you regenerate the complete plasmid sequence.

  3. Hopefully, the reads which span the very edges of the plasmid in the genome will be retained in your new assembly (assuming they weren't thrown out by the mapping step if it was too stringent.

  4. Take the flanks of your new assembly, which with any luck will be a nice single contig containing a small amount of joining sequence.

  5. BLAST (or similar) your flanking sequences back against the reference/target genome which will give you the positiions of insertion.

The only problem I can forsee (other than not enough reads being retained after alignment as I mentioned), is if the flanking sequences are quite repetitive, in which case you might end up identifing multiple places within the genome.

Now here's what I would actually have done:

Design some primers internal to your plasmid pointing out along the genome (if you know how the plasmid integrates), then just send the DNA+Primer for Sanger sequencing for about $3. Basically end-sequence your joins, and you'd get more than enough sequence back that way to be certain of where the plasmid is.

ADD COMMENT
0
Entering edit mode
4.3 years ago
Lisa Ha ▴ 120

You can get the flanking sequences by extracting the reads that contain the beginning or end of the construct. Then you can search for the flanking sequences in the genome. To verify the position, you can do a PCR on the strain containing the construct.

ADD COMMENT
0
Entering edit mode

I Lisa Ha, Thank you for your answer, but I would like a bioinformatics tool to find the position of the construct respect to genome... I do not know if this is possible...

ADD REPLY
0
Entering edit mode

You are unlikely to find a tool 'ready made' for this. You're going to have to get your hands dirty.

ADD REPLY
0
Entering edit mode

Of course... it is clear! In fact I do not want a single tool but at least an idea for a strategy.

ADD REPLY
0
Entering edit mode

I doubt there is a single tool that does exactly what you want. You're going to have to put in a bit of effort. You can use grep on the command line to extract the reads that contain parts of your construct. Then map these to the genome and look at where the reads align (and stop aligning) with a genome viewer, something like IGV or the online Ensembl browser.

ADD REPLY

Login before adding your answer.

Traffic: 1363 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6