X-Git-Url: http://source.jalview.org/gitweb/?a=blobdiff_plain;f=test%2Fjalview%2Fext%2Fensembl%2FEnsemblGeneTest.java;h=446b4f766e463fad94af2b56ef1e932948c8e25a;hb=b9bde68b87a581a1f21e47cb19839791a192134a;hp=edecc23258f2da7225e4ceaae665b7f41553207b;hpb=87a85be7fc7678455c298287349b03fdd12fd3ad;p=jalview.git diff --git a/test/jalview/ext/ensembl/EnsemblGeneTest.java b/test/jalview/ext/ensembl/EnsemblGeneTest.java index edecc23..446b4f7 100644 --- a/test/jalview/ext/ensembl/EnsemblGeneTest.java +++ b/test/jalview/ext/ensembl/EnsemblGeneTest.java @@ -25,6 +25,8 @@ import static org.testng.AssertJUnit.assertFalse; import static org.testng.AssertJUnit.assertTrue; import jalview.api.FeatureSettingsModelI; +import jalview.bin.Cache; +import jalview.datamodel.Sequence; import jalview.datamodel.SequenceDummy; import jalview.datamodel.SequenceFeature; import jalview.datamodel.SequenceI; @@ -53,6 +55,7 @@ public class EnsemblGeneTest @BeforeClass(alwaysRun = true) public void setUp() { + Cache.loadProperties("test/jalview/io/testProps.jvprops"); SequenceOntologyFactory.setInstance(new SequenceOntologyLite()); } @@ -75,17 +78,9 @@ public class EnsemblGeneTest genomic.setEnd(50000); String geneId = "ABC123"; - // gene at (start+20000) length 501 - // should be ignored - the first 'gene' found defines the whole range - // (note features are found in position order, not addition order) - SequenceFeature sf = new SequenceFeature("gene", "", 20000, 20500, 0f, - null); - sf.setValue("ID", "gene:" + geneId); - sf.setStrand("+"); - genomic.addSequenceFeature(sf); - // gene at (start + 10500) length 101 - sf = new SequenceFeature("gene", "", 10500, 10600, 0f, null); + SequenceFeature sf = new SequenceFeature("gene", "", 10500, 10600, 0f, + null); sf.setValue("ID", "gene:" + geneId); sf.setStrand("+"); genomic.addSequenceFeature(sf); @@ -115,17 +110,9 @@ public class EnsemblGeneTest genomic.setEnd(50000); String geneId = "ABC123"; - // gene at (start+20000) length 501 - // should be ignored - the first 'gene' found defines the whole range - // (real data would only have one such feature) - SequenceFeature sf = new SequenceFeature("ncRNA_gene", "", 20000, - 20500, 0f, null); - sf.setValue("ID", "gene:" + geneId); - sf.setStrand("-"); - genomic.addSequenceFeature(sf); - // gene at (start + 10500) length 101 - sf = new SequenceFeature("gene", "", 10500, 10600, 0f, null); + SequenceFeature sf = new SequenceFeature("gene", "", 10500, 10600, 0f, + null); sf.setValue("ID", "gene:" + geneId); sf.setStrand("+"); genomic.addSequenceFeature(sf); @@ -173,7 +160,8 @@ public class EnsemblGeneTest // NMD_transcript_variant treated like transcript in Ensembl SequenceFeature sf3 = new SequenceFeature("NMD_transcript_variant", "", 22000, 22500, 0f, null); - sf3.setValue("Parent", "gene:" + geneId); + // id matching should not be case-sensitive + sf3.setValue("Parent", "gene:" + geneId.toLowerCase()); sf3.setValue("transcript_id", "transcript3"); genomic.addSequenceFeature(sf3); @@ -211,22 +199,24 @@ public class EnsemblGeneTest sf.setValue("ID", "gene:" + geneId); assertFalse(testee.retainFeature(sf, geneId)); - sf.setType("transcript"); + sf = new SequenceFeature("transcript", "", 20000, 20500, 0f, null); sf.setValue("Parent", "gene:" + geneId); assertTrue(testee.retainFeature(sf, geneId)); - sf.setType("mature_transcript"); + sf = new SequenceFeature("mature_transcript", "", 20000, 20500, 0f, + null); sf.setValue("Parent", "gene:" + geneId); assertTrue(testee.retainFeature(sf, geneId)); - sf.setType("NMD_transcript_variant"); + sf = new SequenceFeature("NMD_transcript_variant", "", 20000, 20500, + 0f, null); sf.setValue("Parent", "gene:" + geneId); assertTrue(testee.retainFeature(sf, geneId)); sf.setValue("Parent", "gene:XYZ"); assertFalse(testee.retainFeature(sf, geneId)); - sf.setType("anything"); + sf = new SequenceFeature("anything", "", 20000, 20500, 0f, null); assertTrue(testee.retainFeature(sf, geneId)); } @@ -235,34 +225,48 @@ public class EnsemblGeneTest * accession id as ID */ @Test(groups = "Functional") - public void testIdentifiesSequence() + public void testGetIdentifyingFeatures() { String accId = "ABC123"; - EnsemblGene testee = new EnsemblGene(); + SequenceI seq = new Sequence(accId, "HIBEES"); // gene with no ID not valid - SequenceFeature sf = new SequenceFeature("gene", "", 1, 2, 0f, null); - assertFalse(testee.identifiesSequence(sf, accId)); + SequenceFeature sf1 = new SequenceFeature("gene", "", 1, 2, 0f, null); + seq.addSequenceFeature(sf1); // gene with wrong ID not valid - sf.setValue("ID", "gene:XYZ"); - assertFalse(testee.identifiesSequence(sf, accId)); + SequenceFeature sf2 = new SequenceFeature("gene", "", 1, 2, 0f, null); + sf2.setValue("ID", "gene:XYZ"); + seq.addSequenceFeature(sf2); // gene with right ID is valid - sf.setValue("ID", "gene:" + accId); - assertTrue(testee.identifiesSequence(sf, accId)); + SequenceFeature sf3 = new SequenceFeature("gene", "", 1, 2, 0f, null); + sf3.setValue("ID", "gene:" + accId); + seq.addSequenceFeature(sf3); // gene sub-type with right ID is valid - sf.setType("snRNA_gene"); - assertTrue(testee.identifiesSequence(sf, accId)); + SequenceFeature sf4 = new SequenceFeature("snRNA_gene", "", 1, 2, 0f, null); + sf4.setValue("ID", "gene:" + accId); + seq.addSequenceFeature(sf4); // transcript not valid: - sf.setType("transcript"); - assertFalse(testee.identifiesSequence(sf, accId)); + SequenceFeature sf5 = new SequenceFeature("transcript", "", 1, 2, 0f, null); + sf5.setValue("ID", "gene:" + accId); + seq.addSequenceFeature(sf5); // exon not valid: - sf.setType("exon"); - assertFalse(testee.identifiesSequence(sf, accId)); + SequenceFeature sf6 = new SequenceFeature("exon", "", 1, 2, 0f, null); + sf6.setValue("ID", "gene:" + accId); + seq.addSequenceFeature(sf6); + + List sfs = new EnsemblGene() + .getIdentifyingFeatures(seq, accId); + assertFalse(sfs.contains(sf1)); + assertFalse(sfs.contains(sf2)); + assertTrue(sfs.contains(sf3)); + assertTrue(sfs.contains(sf4)); + assertFalse(sfs.contains(sf5)); + assertFalse(sfs.contains(sf6)); } /** @@ -291,4 +295,28 @@ public class EnsemblGeneTest assertEquals(-1, fc.compare("coding_exon", "feature_variant")); assertEquals(1f, fc.getTransparency()); } + + @Test(groups = "Network") + public void testGetGeneIds() + { + /* + * ENSG00000158828 gene id PINK1 human + * ENST00000321556 transcript for the same gene - should not be duplicated + * P30419 Uniprot identifier for ENSG00000136448 + * ENST00000592782 transcript for Uniprot gene - should not be duplicated + * BRAF - gene name resolvabe (at time of writing) for 6 model species + */ + String ids = "ENSG00000158828 ENST00000321556 P30419 ENST00000592782 BRAF"; + EnsemblGene testee = new EnsemblGene(); + List geneIds = testee.getGeneIds(ids); + assertEquals(8, geneIds.size()); + assertTrue(geneIds.contains("ENSG00000158828")); + assertTrue(geneIds.contains("ENSG00000136448")); + assertTrue(geneIds.contains("ENSG00000157764")); // BRAF human + assertTrue(geneIds.contains("ENSMUSG00000002413")); // mouse + assertTrue(geneIds.contains("ENSRNOG00000010957")); // rat + assertTrue(geneIds.contains("ENSXETG00000004845")); // xenopus + assertTrue(geneIds.contains("ENSDARG00000017661")); // zebrafish + assertTrue(geneIds.contains("ENSGALG00000012865")); // chicken + } }