[Bioperl-l] Install help errors

snoze pa snoze.pa at gmail.com
Tue Jul 8 19:00:05 UTC 2008


Hilmar,

 I really had tough time dealing with bioperl-db and then loading the
sequences. As we discussed earlier(long back) that DBSOURCE is
crashing every time so I removed the gen bank installation. Currently
I am loading uniprot_sprot.dat file which I downloaded from swissprot
site. Based on the guideline in bio-perl website I am able to manage
the bioperl/biosql/bioperl-db/ installation and then try to load the
swissprot file,  uniprot_sprot.dat,  officially uniprot_sprot.dat.gz
latest release july 7th 2008. It always give me a warning

"insert in Bio::DB::BioSQL::SeqFeatureAdaptor (driver) failed,"

which I think okie as long the loading seq is fine. I can see it from
SELECT * FROM biosequence

but now it also crashes after loading 214 sequence from uniprot_sprot.dat file.

Here is the error message

---------------------------------------------------
Could not store P21215:
------------- EXCEPTION: Bio::Root::Exception -------------
MSG: The supplied lineage does not start near 'Clostridium sp.' (I was
supplied 'sp. ATCC29733 | Clostridium | Clostridiaceae | Clostridiales
| Clostridia |Firmicutes | Bacteria')
STACK: Error::throw
STACK: Bio::Root::Root::throw
/usr/local/lib/perl5/site_perl/5.10.0/Bio/Root/Root.pm:359
STACK: Bio::Species::classification
/usr/local/lib/perl5/site_perl/5.10.0/Bio/Species.pm:174
STACK: Bio::DB::Persistent::PersistentObject::AUTOLOAD
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/Persistent/PersistentObject.pm:552
STACK: Bio::DB::BioSQL::SpeciesAdaptor::populate_from_row
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/SpeciesAdaptor.pm:281
STACK: Bio::DB::BioSQL::BasePersistenceAdaptor::_build_object
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/BasePersistenceAdaptor.pm:1305
STACK: Bio::DB::BioSQL::BasePersistenceAdaptor::_find_by_unique_key
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/BasePersistenceAdaptor.pm:973
STACK: Bio::DB::BioSQL::BasePersistenceAdaptor::find_by_unique_key
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/BasePersistenceAdaptor.pm:852
STACK: Bio::DB::BioSQL::BasePersistenceAdaptor::create
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/BasePersistenceAdaptor.pm:182
STACK: Bio::DB::Persistent::PersistentObject::create
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/Persistent/PersistentObject.pm:244
STACK: Bio::DB::BioSQL::BasePersistenceAdaptor::create
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/BasePersistenceAdaptor.pm:169
STACK: Bio::DB::BioSQL::BasePersistenceAdaptor::store
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/BioSQL/BasePersistenceAdaptor.pm:251
STACK: Bio::DB::Persistent::PersistentObject::store
/usr/local/lib/perl5/site_perl/5.10.0/Bio/DB/Persistent/PersistentObject.pm:271
STACK: /usr/local/bioperl-db/scripts/biosql/load_seqdatabase.pl:630
-----------------------------------------------------------

 at /usr/local/bioperl-db/scripts/biosql/load_seqdatabase.pl line 643


If I remove this entry from uniprot_sprot.dat file then it gives
similar error in other IDs

Any Idea,

 Thanks a lot for your reply and help

On Tue, Jul 8, 2008 at 1:35 PM, Hilmar Lapp <hlapp at gmx.net> wrote:
>
> On Jul 8, 2008, at 11:29 AM, snoze pa wrote:
>
>> I agree with you. There is no clean installation help of the
>> bioperl/biosql/bioperl-db available. If you read the biosql INSTALL
>> file and follow what is written then you will find database name
>> biosql. However if you read the mysql entry in same page then you will
>> find the name of the data base is bioseqdb. For new bee in the area it
>> is not only difficult but also  confusing.
>
> That's useful feedback. In fact it sounds like there is a documentation bug
> (inconsistent documentation is a bug, IMHO). Would you mind filing a bug
> report (BioSQL has its own queue now)?
>
>>
>>
>> I am trying a clean installation of bioperl/biosql/bioseqdb but still
>> there is no success. If I assume the installation is done then it
>> always crash on DBSOURCE entry.  The load_sequencedb.pl crashes every
>> time when I encounter with sequences like
>
> Could you post the error message? The entry you quote in reality is a
> Swissprot (or Uniprot) entry sucked up into GenPept, so it's possible that
> it's the BioPerl genbank parser that's failing, not bioperl-db. I recall
> that we've had issues with these records before, though I don't recall what
> the actual cause was, nor the recommended fix. Is there a requirement that
> prevents you from loading these entries in their native Uniprot format?
>
>        -hilmar
>
>>
>>
>> LOCUS       P27912                   792 aa            linear   VRL
>> 10-JUN-2008
>> DEFINITION  Genome polyprotein [Contains: Protein C (Core protein) (Capsid
>>           protein); prM; Peptide pr; Small envelope protein M (Matrix
>>           protein); Envelope protein E; Non-structural protein 1 (NS1)].
>> ACCESSION   P27912
>> VERSION     P27912.1  GI:130422
>> DBSOURCE    swissprot: locus POLG_DEN1A, accession P27912;
>>           class: standard.
>>           created: Aug 1, 1992.
>>           sequence updated: Aug 1, 1992.
>>           annotation updated: Jun 10, 2008.
>>           xrefs: D00502.1, BAA00394.1, B32401
>>           xrefs (non-sequence databases): HSSP:Q88653, SMR:P27912,
>>           GO:0005789, InterPro:IPR011999, InterPro:IPR013754,
>>           InterPro:IPR001122, InterPro:IPR000069, InterPro:IPR001157,
>>           InterPro:IPR002535, InterPro:IPR000336, Gene3D:G3DSA:2.60.98.10,
>>           Gene3D:G3DSA:2.60.40.350, Pfam:PF01003, Pfam:PF02832,
>> Pfam:PF00869,
>>           Pfam:PF01004, Pfam:PF00948, Pfam:PF01570
>> KEYWORDS    Capsid protein; Cleavage on pair of basic residues;
>> Endoplasmic
>>           reticulum; Envelope protein; Glycoprotein; Membrane; Secreted;
>>           Transmembrane; Viral nucleoprotein; Virion.
>>
>>
>> It can not parse the DBSOURCE entry.  May be I am doing something
>> wrong but It will great if someone can try to load
>> P27912, entry of genbank to bioseqdb database.  Please let me know if
>> you can load the this genbank file.
>>
>> Thank you all for your Help.
>> Best S
>>
>>
>>
>> On Mon, Jul 7, 2008 at 4:20 PM, Sendu Bala <bix at sendu.me.uk> wrote:
>>>
>>> Hilmar Lapp wrote:
>>>>
>>>> You may use any name for the database - is there some piece of
>>>> documentation that prompted you to conclude that this parameter isn't
>>>> under
>>>> your control? There should be a command-line parameter (--namespace) in
>>>> load_seqdatabase.pl
>>>
>>> For what its worth, I was caught out by this as well. It would be a
>>> little
>>> more friendly if the documentation for biosql and bioperl-db were
>>> consistent
>>> in their suggestions for what to call the database.
>>>
>>> When you don't really have any idea what's going on, you at least want to
>>> be
>>> able to copy/paste the suggested commands/code ;)
>>>
>>>
>>>> On Jul 7, 2008, at 12:23 PM, snoze pa wrote:
>>>>
>>>>> Dear Users,
>>>>> Thank you very much for you timely support. I highly appreciate your
>>>>> help in installing the software.  Meanwhile during the  installation
>>>>> of  bioperl-db and  loading load_seqdatabase.pl, it ask me to use the
>>>>> dbname "biosql" while the database created in during  bioperl
>>>>> installation is "bioseqdb". Shall I use bioseqdb instead of biosql as
>>>>> dbname.
>>>
>> _______________________________________________
>> Bioperl-l mailing list
>> Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
> --
> ===========================================================
> : Hilmar Lapp  -:-  Durham, NC  -:-  hlapp at gmx dot net :
> ===========================================================
>
>
>
>



More information about the Bioperl-l mailing list