Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • format uniref90.xml to database for BLAST

    we want to get a formated database for BLAST from Uniref90.
    we referred an article which used an early release Uniref90, say version 10.0
    ftp.uniprot.org only provide the XML format of this version.

    How can we convert XML uniref90.xml to uniref90.fasta ? or, can formatdb take XML file as input file?
    Thanks!

    XML file:
    ftp://ftp.uniprot.org/pub/databases/...ref10.0.tar.gz

  • #2
    You should be able to convert the UniProt XML to FASTA using Biopython,

    Code:
    from Bio import SeqIO
    count = SeqIO.convert("uniref90.xml", "uniprot-xml", "converted.fasta", "fasta")
    print("Converted %i records" % count)

    Comment


    • #3
      Thank you for your quick reply. I'll try it out.

      Comment


      • #4
        Code:
        >>> from Bio import SeqIO
        >>> count = SeqIO.convert("uniref90.xml", "uniprot-xml", "uniref90converted.fasta", "fasta")
        >>> print("Converted %i records" % count)
        Converted 0 records

        We checked uniref90.xml:
        more uniref90.xml
        Code:
        <?xml version="1.0" encoding="ISO-8859-1" ?>
        <UniRef90 xmlns="http://uniprot.org/uniref" 
        xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" 
        xsi:schemaLocation="http://uniprot.org/uniref http://www.uniprot.org/support/docs/uniref.xsd" 
         releaseDate="2007-03-06" version="10.0"> 
        <entry id="UniRef90_Q3ASY8" updated="2007-03-06">
        <name>Cluster: Parallel beta-helix repeat</name>
        <property type="member count" value="1"/>
        <property type="common taxon" value="Chlorobium chlorochromatii CaD3"/>
        <property type="common taxon ID" value="340177"/>
        <representativeMember>
        <dbReference type="UniProtKB ID" id="Q3ASY8_CHLCH">
        <property type="UniProtKB accession" value="Q3ASY8"/>
        <property type="UniParc ID" value="UPI00005D5563"/>
        <property type="UniRef100 ID" value="UniRef100_Q3ASY8"/>
        <property type="UniRef50 ID" value="UniRef50_Q3ASY8"/>
        <property type="protein name" value="Parallel beta-helix repeat"/>
        <property type="source organism" value="Chlorobium chlorochromatii (strain CaD3)"/>
        <property type="NCBI taxonomy" value="340177"/>
        <property type="length" value="36805"/>
        <property type="isSeed" value="true"/>
        </dbReference>
        <sequence length="36805" checksum="A7A8EA21B9345FF9">
        MKPRFYIEQLEPRILLSGDILSELVPLLSSREASQMQSDYLLEHPEARRVAPLSAVEAAR
        ....
        Could you help us, thanks.
        Last edited by emanlee; 09-22-2013, 04:55 PM.

        Comment


        • #5
          Originally posted by emanlee View Post
          Could you help us, thanks.
          Wouldn't it just be much easier to download the UniRef90 FASTA file directly?

          Comment


          • #6
            Originally posted by kmcarr View Post
            Wouldn't it just be much easier to download the UniRef90 FASTA file directly?
            Indeed, I should have doubled checked that really didn't exist.

            As to the Biopython conversion failing, that is probably a bug - I'd have replied earlier but missed the thread reply alert - sorry.

            Comment


            • #7
              The file linked by kmcarr does not refer to a "version 10.0" that emanlee was asking for. Perhaps that is not important.

              Comment


              • #8
                OK then... first this is how I just extracted the uniref90.xml file from the FTP site (multiple levels of bundling!):
                Code:
                $ wget ftp://ftp.uniprot.org/pub/databases/uniprot/previous_releases/release10.0/uniref/uniref10.0.tar.gz
                ...
                $ tar -zxvf uniref10.0.tar.gz 
                uniref100.tar
                uniref50.tar
                $ tar -xvf uniref90.tar 
                README
                uniref90.dtd
                uniref90.xml.gz
                $ gunzip uniref90.xml.gz
                And here is what the start of the file looks like for me too (same as emanlee reported):
                Code:
                $ head -n 25 uniref90.xml 
                <?xml version="1.0" encoding="ISO-8859-1" ?>
                <UniRef90 xmlns="http://uniprot.org/uniref" 
                xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" 
                xsi:schemaLocation="http://uniprot.org/uniref http://www.uniprot.org/support/docs/uniref.xsd" 
                 releaseDate="2007-03-06" version="10.0"> 
                <entry id="UniRef90_Q3ASY8" updated="2007-03-06">
                <name>Cluster: Parallel beta-helix repeat</name>
                <property type="member count" value="1"/>
                <property type="common taxon" value="Chlorobium chlorochromatii CaD3"/>
                <property type="common taxon ID" value="340177"/>
                <representativeMember>
                <dbReference type="UniProtKB ID" id="Q3ASY8_CHLCH">
                <property type="UniProtKB accession" value="Q3ASY8"/>
                <property type="UniParc ID" value="UPI00005D5563"/>
                <property type="UniRef100 ID" value="UniRef100_Q3ASY8"/>
                <property type="UniRef50 ID" value="UniRef50_Q3ASY8"/>
                <property type="protein name" value="Parallel beta-helix repeat"/>
                <property type="source organism" value="Chlorobium chlorochromatii (strain CaD3)"/>
                <property type="NCBI taxonomy" value="340177"/>
                <property type="length" value="36805"/>
                <property type="isSeed" value="true"/>
                </dbReference>
                <sequence length="36805" checksum="A7A8EA21B9345FF9">
                MKPRFYIEQLEPRILLSGDILSELVPLLSSREASQMQSDYLLEHPEARRVAPLSAVEAAR
                ACMVVVQSEAPSLLTEDGLMYPFEVGVGEERSSEANAEPTLAADFSADYTFSKSEWDALE
                And here's how many records there seem to be according to grep:
                Code:
                $ grep -c "^<entry id" uniref90.xml 
                2781437
                Biopython 1.61 and 1.62 do appear to have a problem parsing this - I suspect the XML is different in some way to what we expect.

                Update: Raised here: http://lists.open-bio.org/pipermail/...er/010909.html
                Last edited by maubp; 09-27-2013, 07:52 AM.

                Comment

                Latest Articles

                Collapse

                • seqadmin
                  Strategies for Sequencing Challenging Samples
                  by seqadmin


                  Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
                  03-22-2024, 06:39 AM
                • seqadmin
                  Techniques and Challenges in Conservation Genomics
                  by seqadmin



                  The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

                  Avian Conservation
                  Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
                  03-08-2024, 10:41 AM

                ad_right_rmr

                Collapse

                News

                Collapse

                Topics Statistics Last Post
                Started by seqadmin, Yesterday, 06:37 PM
                0 responses
                10 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, Yesterday, 06:07 PM
                0 responses
                9 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-22-2024, 10:03 AM
                0 responses
                49 views
                0 likes
                Last Post seqadmin  
                Started by seqadmin, 03-21-2024, 07:32 AM
                0 responses
                67 views
                0 likes
                Last Post seqadmin  
                Working...
                X