逐行解析 XML 文件

Question

所以这就是问题所在。我正在尝试解析来自 GenBank 的 XML 信息文件。该文件包含多个 DNA 序列的信息。我已经对 genbacnk 的另外两种 xml 格式（TINY xml 和 INSD xml）完成了此操作，但纯 xml 让我头疼。这是我的程序应该如何工作。下载 xml 格式的文件，其中包含有关 GenBank 中 X 个序列的信息。运行我的 perl 脚本，逐行搜索该 xml 文件并将我想要的信息以 fasta 格式打印到新文件中。这是这样的： >Sequence_name_and_information\n strings\n >sequence_name.... 如此反复，直到获得 xml 文件中的所有序列。但我的问题是，在纯 xml 中，序列本身位于基因或序列基因座的标识符之前。基因或序列基因座应与“>”位于同一行。这是我从打开文件并解析它开始得到的代码：

open( New_File, "+>$PWD_file/$new_file" ) or die "\n\nCouldn't create file. Check permissions on location.\n\n";

    while ( my $lines = <INSD> ) {
        foreach ($lines) {
            if (m/<INSDSeq_locus>.*<\/INSDSeq_locus>/) {
                $lines =~ s/<INSDSeq_locus>//g and $lines =~ s/<\/INSDSeq_locus>//g and $lines =~ s/[a-z, |]//g; #this last bit may cause a bug of removing the letters in the genbank accession number
                $lines =~ s/ //g;
                chomp($lines);
                print New_File ">$lines\_";
            } elsif (m/<INSDSeq_organism>.*<\/INSDSeq_organism>/) {
                $lines =~ s/<INSDSeq_organism>//g and $lines =~ s/<\/INSDSeq_organism>//g;
                $lines =~ s/(\.|\?|\-| )/_/g;
                $lines =~ s/_{2,}/_/g;
                $lines =~ s/_{1,}$//;
                $lines =~ s/^>*_{1,}//; 
                $lines =~ s/\s{2}//g;
                chomp($lines);
                print New_File "$lines\n";
            } elsif (m/<INSDSeq_sequence>.*<\/INSDSeq_sequence>/) {
                $lines =~ s/<INSDSeq_sequence>//g and $lines =~ s/<\/INSDSeq_sequence>//g;
                $lines =~ s/ //g;
                chomp($lines);
                print New_File "$lines\n";
            }
        }
    }
    close INSD;
    close New_File;
}

有两个地方可以找到基因/位点信息。该信息位于以下两个标签之间：LOCUS_NAME 或 GENE_NAME。会有一个，或者另一个。如果一个有信息，另一个将是空的。无论哪种情况，都需要添加到 >....... 行的末尾。

谢谢，

AlphaA

PS--我尝试通过打开“$NA”，“>”将该信息打印到“文件”中序列，然后继续程序，找到基因信息，将其打印到 >行，然后读取 $NA 文件并将其打印到 > 之后的行；线。我希望这一点很清楚。

Answer 1

在我看来，您应该将 XSLT 与 XPath 导航到您需要的数据。

正如 @Brian 所建议的，使用已建立的 XML 解析技术和库会更容易。

甚至还有一个 用于 XSLT 的 Perl 库

Answer 2

使用 XML 解析器。我不是生物学家，我不确定你想要的最终格式，但以此为起点应该很简单。匿名子中的 $_[1] 包含一个哈希引用，从上面我可以看出，我认为您希望通过解析所需标签的父标签来保存所有内容。以您想要的格式打印 $_[1] 的元素应该很容易：

use strict;
use warnings;

use XML::Rules;
use Data::Dumper;

my @rules = (
  _default => '',
  'INSDSeq_locus,INSDSeq_organism,INSDSeq_sequence' => 'content',
  INSDSeq  => sub { delete $_[1]{_content}; print Dumper $_[1]; return },
);

my $p = XML::Rules->new(rules => \@rules);
$p->parsefile('sequence.gbc.xml');

这样就可以很容易地只打印您想要的标签。或者，如果您想要一些其他标签，我真正可能做的是这样的（如果您只是逐个打印元素，则根本不需要 @tags 变量）

my @tags = qw(
  INSDSeq_locus
  INSDSeq_organism
  INSDSeq_sequence
);

my @rules = (
  _default => 'content',
  # Elements are, e.g. $_[1]{INSDSeq_locus}
  INSDSeq  => sub { print "$_: $_[1]{$_}\n" for @tags; return; },
);

：

my $p = XML::Rules->new(rules => \@rules, stripspaces => 4);