package jalview.ext.ensembl;
import static org.testng.AssertJUnit.assertEquals;
+import static org.testng.AssertJUnit.assertFalse;
+import static org.testng.AssertJUnit.assertTrue;
+import static org.testng.internal.junit.ArrayAsserts.assertArrayEquals;
import jalview.datamodel.Alignment;
import jalview.datamodel.AlignmentI;
+import jalview.datamodel.SequenceFeature;
import jalview.datamodel.SequenceI;
import jalview.io.AppletFormatAdapter;
import jalview.io.FastaFile;
import jalview.io.FileParse;
+import jalview.io.gff.SequenceOntologyFactory;
+import jalview.io.gff.SequenceOntologyLite;
import java.lang.reflect.Method;
import java.net.MalformedURLException;
import java.util.List;
import org.testng.Assert;
+import org.testng.annotations.AfterClass;
+import org.testng.annotations.BeforeClass;
import org.testng.annotations.DataProvider;
import org.testng.annotations.Test;
+ "NRDQIIFMVGRGYLSPDLSKVRSNCPKAMKRLMAECLKKKRDERPLFPQILASIELLARS\n"
+ "LPKIHRSASEPSLNRAGFQTEDFSLYACASPKTPIQAGGYGAFPVH" } };
- @DataProvider(name = "queries")
- public Object[][] createQueryData(Method m)
+ @BeforeClass
+ public void setUp()
{
- return new Object[][] { { "CCDS5863.1" }, { "ENSP00000288602" } };
+ SequenceOntologyFactory.setInstance(new SequenceOntologyLite());
}
- @Test(dataProvider = "queries")
- public void testIsValidReference(String query) throws Exception
+ @AfterClass
+ public void tearDown()
{
- EnsemblSequenceFetcher esq = new EnsemblProtein();
- Assert.assertTrue(esq.isValidReference(query),
- "Expected reference string " + query
- + " to be valid for regex "
- + esq.getAccessionValidator().toString());
+ SequenceOntologyFactory.setInstance(null);
}
@DataProvider(name = "ens_seqs")
: "DOWN or unreachable ******************* BAD!"));
}
- /**
- * Tests for the method that computes all peptide variants given codon
- * variants
- */
@Test(groups = "Functional")
- public void testComputePeptideVariants()
+ public void getGenomicRangesFromFeatures()
{
- String[][] codonVariants = new String[][] { { "A" }, { "G" }, { "T" } };
- /*
- * AGT codes for S - this is not included in the variants returned
- */
- List<String> variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "S");
- assertEquals("[]", variants.toString());
+ }
- // S is reported if it differs from the current value (A):
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "A");
- assertEquals("[S]", variants.toString());
+ @Test(groups = "Functional")
+ public void testIsTranscriptIdentifier()
+ {
+ EnsemblSeqProxy testee = new EnsemblGene();
+ assertFalse(testee.isTranscriptIdentifier(null));
+ assertFalse(testee.isTranscriptIdentifier(""));
+ assertFalse(testee.isTranscriptIdentifier("ENSG00000012345"));
+ assertTrue(testee.isTranscriptIdentifier("ENST00000012345"));
+ assertTrue(testee.isTranscriptIdentifier("ENSMUST00000012345"));
+ assertFalse(testee.isTranscriptIdentifier("enst00000012345"));
+ assertFalse(testee.isTranscriptIdentifier("ENST000000123456"));
+ assertFalse(testee.isTranscriptIdentifier("ENST0000001234"));
+ }
- /*
- * synonymous variant is not reported
- */
- codonVariants = new String[][] { { "A" }, { "G" }, { "C", "T" } };
- // AGC and AGT both code for S
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "s");
- assertEquals("[]", variants.toString());
+ @Test(groups = "Functional")
+ public void testIsGeneIdentifier()
+ {
+ EnsemblSeqProxy testee = new EnsemblGene();
+ assertFalse(testee.isGeneIdentifier(null));
+ assertFalse(testee.isGeneIdentifier(""));
+ assertFalse(testee.isGeneIdentifier("ENST00000012345"));
+ assertTrue(testee.isGeneIdentifier("ENSG00000012345"));
+ assertTrue(testee.isGeneIdentifier("ENSMUSG00000012345"));
+ assertFalse(testee.isGeneIdentifier("ensg00000012345"));
+ assertFalse(testee.isGeneIdentifier("ENSG000000123456"));
+ assertFalse(testee.isGeneIdentifier("ENSG0000001234"));
+ }
- /*
- * equivalent variants are only reported once
- */
- codonVariants = new String[][] { { "C" }, { "T" },
- { "A", "C", "G", "T" } };
- // CTA CTC CTG CTT all code for L
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "S");
- assertEquals("[L]", variants.toString());
+ /**
+ * Test the method that appends a single allele's reverse complement to a
+ * string buffer
+ */
+ @Test(groups = "Functional")
+ public void testReverseComplementAllele()
+ {
+ StringBuilder sb = new StringBuilder();
+ EnsemblSeqProxy.reverseComplementAllele(sb, "G"); // comp=C
+ EnsemblSeqProxy.reverseComplementAllele(sb, "g"); // comp=c
+ EnsemblSeqProxy.reverseComplementAllele(sb, "C"); // comp=G
+ EnsemblSeqProxy.reverseComplementAllele(sb, "T"); // comp=A
+ EnsemblSeqProxy.reverseComplementAllele(sb, "A"); // comp=T
+ assertEquals("C,c,G,A,T", sb.toString());
- /*
- * vary codons 1 and 2; variant products are sorted and non-redundant
- */
- codonVariants = new String[][] { { "a", "C" }, { "g", "T" }, { "A" } };
- // aga ata cga cta code for R, I, R, L
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "S");
- assertEquals("[I, L, R]", variants.toString());
+ sb = new StringBuilder();
+ EnsemblSeqProxy.reverseComplementAllele(sb, "-GATt"); // revcomp=aATC-
+ EnsemblSeqProxy.reverseComplementAllele(sb, "hgmd_mutation");
+ EnsemblSeqProxy.reverseComplementAllele(sb, "PhenCode_variation");
+ assertEquals("aATC-,hgmd_mutation,PhenCode_variation", sb.toString());
+ }
- /*
- * vary codons 2 and 3
- */
- codonVariants = new String[][] { { "a" }, { "g", "T" }, { "A", "c" } };
- // aga agc ata atc code for R, S, I, I
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "S");
- assertEquals("[I, R]", variants.toString());
+ /**
+ * Test the method that computes the reverse complement of the alleles in a
+ * sequence_variant feature
+ */
+ @Test(groups = "Functional")
+ public void testReverseComplementAlleles()
+ {
+ String alleles = "C,G,-TAC,HGMD_MUTATION,gac";
+ SequenceFeature sf = new SequenceFeature("sequence_variant", alleles,
+ 1, 2, 0f, null);
+ sf.setValue("alleles", alleles);
+ sf.setAttributes("x=y,z;alleles=" + alleles + ";a=b,c");
+
+ EnsemblSeqProxy.reverseComplementAlleles(sf);
+ String revcomp = "G,C,GTA-,HGMD_MUTATION,gtc";
+ // verify description is updated with reverse complement
+ assertEquals(revcomp, sf.getDescription());
+ // verify alleles attribute is updated with reverse complement
+ assertEquals(revcomp, sf.getValue("alleles"));
+ // verify attributes string is updated with reverse complement
+ assertEquals("x=y,z;alleles=" + revcomp + ";a=b,c", sf.getAttributes());
+ }
+
+ @Test(groups = "Functional")
+ public void testSortFeatures()
+ {
+ SequenceFeature sf1 = new SequenceFeature("", "", 10, 15, 0f, null);
+ SequenceFeature sf2 = new SequenceFeature("", "", 8, 12, 0f, null);
+ SequenceFeature sf3 = new SequenceFeature("", "", 8, 13, 0f, null);
+ SequenceFeature sf4 = new SequenceFeature("", "", 11, 11, 0f, null);
+ SequenceFeature[] sfs = new SequenceFeature[] { sf1, sf2, sf3, sf4 };
- /*
- * vary codons 1 and 3
- */
- codonVariants = new String[][] { { "a", "t" }, { "a" }, { "t", "g" } };
- // aat aag tat tag code for N, K, Y, STOP - STOP sorted to end
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "S");
- assertEquals("[K, N, Y, STOP]", variants.toString());
+ // sort by start position ascending (forward strand)
+ // sf2 and sf3 tie and should not be reordered by sorting
+ EnsemblSeqProxy.sortFeatures(sfs, true);
+ assertArrayEquals(new SequenceFeature[] { sf2, sf3, sf1, sf4 }, sfs);
- /*
- * vary codons 1, 2 and 3
- */
- codonVariants = new String[][] { { "a", "t" }, { "G", "C" },
- { "t", "g" } };
- // agt agg act acg tgt tgg tct tcg code for S, R, T, T, C, W, S, S
- variants = EnsemblSeqProxy.computePeptideVariants(codonVariants, "S");
- assertEquals("[C, R, T, W]", variants.toString());
+ // sort by end position descending (reverse strand)
+ EnsemblSeqProxy.sortFeatures(sfs, false);
+ assertArrayEquals(new SequenceFeature[] { sf1, sf3, sf2, sf4 }, sfs);
}
}
\ No newline at end of file