Question: Extracting and joining exons from multiple sequence alignments using PERL
1
gravatar for Eskimo
3.9 years ago by
Eskimo90
Eskimo90 wrote:

Hi all,

I would really appreciate any help you can provide with the following coding problem (using PERL). I have an aligned multiple sequence alignment (fasta format) as my input, along with a file containing the co-ordinates of all open reading frames (ORFs) and the exons that make them up, I am wanting to slice out the ORFs according to those co-ordinates. So far I have code that works well for individual ORFs (and have no problem with the reverse complementing etc) but the problem lies in extracting and concatenating multi-exon ORFs.

Thus far my code reads in the multiple sequence alignment as follows

use Bio::SimpleAlign;
use Bio::AlignIO;
$str = Bio::AlignIO->new(-file => $inputfilename, -format => 'fasta');
$aln = $str->next_aln();

and deals with the splicing as follows

$mini = $aln->slice($array[0], $array[1]);  
$out = Bio::AlignIO->new(-file => $array[3],
                                     -format => 'fasta'); 
$out->write_aln($mini);

an example of the input file containing co-ordinates looks like this

Start Stop Strand Name Note
24 89 + ORF1 exon1
165 560 - ORF2 exon1
680 1004 + ORF3 exon1
1240 1760 + ORF3 exon2
1790 2360 + ORF3 exon3
2600 2900 - ORF4 exon1
2850 3100 + ORF 5 exon1

Would anyone know of a clever way to extract the individual exons for ORF3 and then concatenate the files (side by side to ensure that the multiple sequence alignment is not comprised)? My initial thought was to change the co-ordinate file to a GFF type file and use Bio::Tools::GFF but I don't think this is compatible with a multiple sequence alignment as the input.

Any help would be hugely appreciated!

ADD COMMENTlink modified 3.9 years ago • written 3.9 years ago by Eskimo90
0
gravatar for novice
3.9 years ago by
novice950
United States
novice950 wrote:

I'm not sure what you mean by "side by side" concatenation, but here's an idea to include exons from the same ORF in one file:

my @orfs;
unless ($array[3] ~~ @orfs) {
    $out = Bio::AlignIO->new(-file => $array[3], -format => 'fasta');
    push @orfs, $array[3];
}

$mini = $aln->slice($array[0], $array[1]); 
$out->write_aln($mini);
ADD COMMENTlink written 3.9 years ago by novice950
0
gravatar for anp375
3.9 years ago by
anp375170
anp375170 wrote:

Can you make a 2D array with blank spaces?

ADD COMMENTlink written 3.9 years ago by anp375170
0
gravatar for Eskimo
3.9 years ago by
Eskimo90
Eskimo90 wrote:

Hi folks, thanks for the answers so far. I just want to clarify my question a little further. Assume my multiple sequence alignment looks like this

genomeA

ACTCGAGCTATCGATCGATCATGCGAGCGCTACTAAATTTCATCGAGCGTATTCTATCTGAGCTAGCATCTTCA

genomeB

ACTCGAGCTATCGATCGATCATGCG---GCTACTATCTCATTCATCGAGCGTATTCGTGCTGAGCTAGCATCTTCA

genomeC

ACTCGAGCTATCTGCGATCATGCGAGCGCTACTATCTCATTCATCGAGCGTAATGTATCTGCTCTAGCATCTTCA

If I have a single exon between the aligned position 5 and 25 I can splice out the MSA in that region with

$mini = $aln->slice(5, 25);
$out = Bio::AlignIO->new(-file => $array[3], -format => 'fasta'); $out->write_aln($mini);

which would give me

genomeA

GAGCTATCGATCGATCATGC

genomeB

GAGCTATCGATCGATCATGC

genomeC

GAGCTATCTGCGATCATGCG

If I then have a second exon between position 30 and 35 I would splice that out with

$mini = $aln->slice(30, 35);
$out = Bio::AlignIO->new(-file => $array[3], -format => 'fasta'); $out->write_aln($mini)

which would give me

genomeA

GCTAC

genomeB

GCTAC

genomeC

GCTAC

What I then want to do use join the two MSA outputs into one so i get

genomeA

GAGCTATCGATCGATCATGCGCTAC

genomeB

GAGCTATCGATCGATCATGCGCTAC

genomeC

GAGCTATCTGCGATCATGCGGCTAC

(Note that I don't want the asterisks, I am just indicating where the sequences are being joined)

What I need is a piece of code that will concatenate the spliced regions from the same ORF (i.e. ORF3 exon1,2&3) but will then write this to the outfile and start again when it reaches a new ORF (i.e. ORF4).

Hopefully that is a little clearer.

Any further comments would be most appreciated.

ADD COMMENTlink written 3.9 years ago by Eskimo90
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: 1043 users visited in the last hour