# -*-Perl-*-
## Bioperl Test Harness Script for Modules
##
# $Id: UniGene.t,v 1.13 2003/09/11 23:52:59 andrew Exp $

# Before `make install' is performed this script should be runnable with
# `make test'. After `make install' it should work as `perl test.t'

use strict;
use vars qw($NUMTESTS);
#use lib '.','./blib/lib'; # make test should take care of this

my $error;

BEGIN { 
    # to handle systems with no installed Test module
    # we include the t dir (where a copy of Test.pm is located)
    # as a fallback
    eval { require Test; };
    $error = 0;
    if( $@ ) {
	use lib 't';
    }
    use Test;

    $NUMTESTS = 61;
    plan tests => $NUMTESTS;
    
}

use Bio::Cluster::UniGene;
use Bio::ClusterIO;

my ($str, $unigene); # predeclare variables for strict


$str = Bio::ClusterIO->new('-file' => Bio::Root::IO->catfile("t","data","unigene.data"), '-format' => "unigene");

ok $str;

ok ( defined ($unigene = $str->next_cluster()));

# check interface implementations to be sure
ok $unigene->isa("Bio::Cluster::UniGeneI");
ok $unigene->isa("Bio::ClusterI");
ok $unigene->isa("Bio::IdentifiableI");
ok $unigene->isa("Bio::DescribableI");


# test specific instance of unigene record provided in the unigene.data file
ok($unigene->unigene_id, 'Hs.2');
ok($unigene->title, 'N-acetyltransferase 2 (arylamine N-acetyltransferase)');
ok($unigene->gene, 'NAT2');
ok($unigene->cytoband,'8p22');
ok($unigene->gnm_terminus,'S');
ok($unigene->scount,25);
ok(scalar @{ $unigene->locuslink }, 1);
ok(scalar @{ $unigene->chromosome }, 1);
ok(scalar @{ $unigene->express }, 7);
ok(scalar @{ $unigene->sts }, 10);
ok(scalar @{ $unigene->txmap }, 0);
ok(scalar @{ $unigene->protsim } , 4);
ok(scalar @{ $unigene->sequences },25);

ok($unigene->locuslink->[0], '10');
ok($unigene->chromosome->[0], '8');
ok($unigene->express->[0], 'liver');
ok($unigene->sts->[0], 'ACC=GDB:386004 UNISTS=157141');
ok($unigene->protsim->[0], 'ORG=Escherischia coli; PROTGI=16129422; PROTID=ref:NP_415980.1; PCT=24; ALN=255');

my ($seq1) = $unigene->next_seq;
ok($seq1->display_id, 'D90040');
#ok($seq1->desc, 'ACC=D90042 NID=g219415 PID=g219416');

# test recognition of species
ok $unigene->species;
skip (! $unigene->species, $unigene->species->binomial, "Homo sapiens");

# test accessors of interfaces
ok ($seq1->namespace, "GenBank");
ok ($seq1->authority, "NCBI");
ok ($seq1->alphabet, "dna");
my $n = 1; # we've seen already one seq
while($seq1 = $unigene->next_seq()) {
    $n++;
}
ok ($n, 25);
ok ($unigene->size(), 25);
ok (scalar($unigene->get_members()), 25);
ok ($unigene->description, 'N-acetyltransferase 2 (arylamine N-acetyltransferase)');
ok ($unigene->display_id, "Hs.2");
ok ($unigene->namespace, "UniGene");
ok ($unigene->authority, "NCBI");

$unigene->unigene_id('Hs.50');
ok($unigene->unigene_id, 'Hs.50', 'unigene_id was ' . $unigene->unigene_id);

$unigene->title('title_test');
ok($unigene->title, 'title_test', 'title was ' . $unigene->title);

$unigene->gene('gene_test');
ok($unigene->gene, 'gene_test', 'gene was ' . $unigene->gene);

$unigene->cytoband('cytoband_test');
ok($unigene->cytoband, 'cytoband_test', 'cytoband was ' . $unigene->cytoband);

$unigene->gnm_terminus('gnm_terminus_test');
ok($unigene->gnm_terminus, 'gnm_terminus_test', 'gnm_terminus was ' . $unigene->gnm_terminus);

$unigene->scount('scount_test');
ok($unigene->scount, 'scount_test', 'scount was ' . $unigene->scount);

my $seq = $unigene->next_seq;
$seq = $unigene->next_seq;
ok($seq->isa('Bio::PrimarySeqI'), 1,'expected a Bio::PrimarySeq object but got a ' . ref($seq));
my $accession = $seq->accession_number;
ok($accession, 'BC015878');
my $version = $seq->seq_version();
ok($version, "");


# test the sequence parsing is working
my $ac = $seq->annotation();
my $simple_ann_object;
($simple_ann_object) = $ac->get_Annotations('seqtype');
ok $simple_ann_object;
ok($simple_ann_object->value(), 'mRNA', 'seqtype was ' . $simple_ann_object->value);	



# tests not specific to unigene record provided in the unigene.data file
my @locuslink_test = ( "58473", "5354" );
$unigene->locuslink(\@locuslink_test);
my @locuslink_results;
while (my $locuslink = $unigene->next_locuslink) {
	push @locuslink_results, $locuslink;
}
ok scalar(@locuslink_results), 2, 'expected locuslink to have 2 entries but it had ' . scalar(@locuslink_results);
my $locuslink = shift @locuslink_results;
ok $locuslink, '58473', 'expected 58473 but got ' . $locuslink;


my @express_test = qw( kidney heart liver spleen );
$unigene->express(\@express_test);
my @express_results;
while (my $tissue = $unigene->next_express) {
	push @express_results, $tissue;
}
ok scalar(@express_results), 4, 'expected express to have 4 entries but it had ' . scalar(@express_results);

my @chromosome_test = ( "7", "11" );
$unigene->chromosome(\@chromosome_test);
my @chromosome_results;
while (my $chromosome = $unigene->next_chromosome) {
	push @chromosome_results, $chromosome;
}
ok scalar(@chromosome_results), 2, 'expected chromosome to have 2 entries but it had ' . scalar(@chromosome_results);
my $chromosome = shift @chromosome_results;
ok $chromosome, '7', 'expected 7 but got ' . $chromosome;

my @sts_test = ( "ACC=- NAME=sts-D90276 UNISTS=37687", "ACC=G29786 NAME=SHGC-35230 UNISTS=58455" );
$unigene->sts(\@sts_test);
my @sts_results;
while (my $sts = $unigene->next_sts) {
	push @sts_results, $sts;
}
ok scalar(@sts_results), 2, 'expected sts to have 2 entries but it had ' . scalar(@sts_results);
my $sts = shift @sts_results;
ok $sts, 'ACC=- NAME=sts-D90276 UNISTS=37687', 'expected ACC=- NAME=sts-D90276 UNISTS=37687 but got ' . $sts;

my @txmap_test = ("D19S425-D19S418; MARKER=sts-D90276; RHPANEL=GB4" , "D19S425-D19S418; MARKER=stSG41396; RHPANEL=GB4");
$unigene->txmap(\@txmap_test);
my @txmap_results;
while (my $txmap = $unigene->next_txmap) {
	push @txmap_results, $txmap;
}
ok scalar(@txmap_results), 2, 'expected txmap to have 2 entries but it had ' . scalar(@txmap_results);
my $txmap = shift @txmap_results;
ok $txmap, 'D19S425-D19S418; MARKER=sts-D90276; RHPANEL=GB4', 'expected D19S425-D19S418; MARKER=sts-D90276; RHPANEL=GB4 but got ' . $txmap;

my @protsim_test = ("ORG=Homo sapiens; PROTGI=107211; PROTID=pir:A40428; PCT=100; ALN=243" , "ORG=Mus musculus; PROTGI=2497288; PROTID=sp:Q61400; PCT=42; ALN=143");
$unigene->protsim(\@protsim_test);
my @protsim_results;
while (my $protsim = $unigene->next_protsim) {
    push @protsim_results, $protsim;
}
ok scalar(@protsim_results), 2, 'expected protsim to have 2 entries but it had ' . scalar(@protsim_results);
my $protsim = shift @protsim_results;
ok $protsim, 'ORG=Homo sapiens; PROTGI=107211; PROTID=pir:A40428; PCT=100; ALN=243', 'expected ORG=Homo sapiens; PROTGI=107211; PROTID=pir:A40428; PCT=100; ALN=243 but got ' . $protsim;



# do a quick test on Rn record included as the next cluster in the test data file because it has version numbers tacked on the end of the accession numbers in each seq line - NCBI has started doing this now (Sept 2003).

$unigene = $str->next_cluster();
$seq = $unigene->next_seq;
ok($seq->isa('Bio::PrimarySeqI'), 1,'expected a Bio::PrimarySeq object but got a ' . ref($seq));
$version = $seq->seq_version();
ok($version, '1');