def test_simple_genes_by_gene_annotation(self): """Test web api backend call /gene/ """ # http://www.broadinstitute.org/oncotator/gene/MAPK1/ datasource_list = DatasourceFactory.createDatasources(self._determine_db_dir(), "hg19", isMulticore=False) annotator = Annotator() for ds in datasource_list: annotator.addDatasource(ds) txs = annotator.retrieve_transcripts_by_genes(["MAPK1"]) self.assertTranscriptsFound(txs) mut_dict = annotator.annotate_genes_given_txs(txs) self.assertTrue(len(mut_dict.keys()) == 1)
def test_simple_genes_by_region_annotation(self): """Test web api backend call /genes/ """ # http://www.broadinstitute.org/oncotator/genes/chr22_22112223_22312558/ # Two genes: chr22:22,112,223-22,312,558 datasource_list = DatasourceFactory.createDatasources(self._determine_db_dir(), "hg19", isMulticore=False) annotator = Annotator() for ds in datasource_list: annotator.addDatasource(ds) # Here is what the API would call.... txs = annotator.retrieve_transcripts_by_region("22", 22112223, 22312558) self.assertTranscriptsFound(txs) mut_dict = annotator.annotate_genes_given_txs(txs) # Each mut will be for a separate gene for gene in mut_dict.keys(): mut = mut_dict[gene] alt_accessions = mut["UniProt_alt_uniprot_accessions"].split("|") tcgascape_amp_peaks = mut["TCGAScape_Amplification_Peaks"].split("|") tcgascape_del_peaks = mut["TCGAScape_Deletion_Peaks"].split("|") tumorscape_amp_peaks = mut["TUMORScape_Amplification_Peaks"].split("|") tumorscape_del_peaks = mut["TUMORScape_Deletion_Peaks"].split("|") full_name = mut["HGNC_Approved Name"] cosmic = { "tissue_types_affected": mut["COSMIC_Tissue_tissue_types_affected"], "total_alterations_in_gene": mut["COSMIC_Tissue_tissue_types_affected"], } alt_aliases = list( itertools.chain([mut["HGNC_Previous Symbols"].split(", "), mut["HGNC_Synonyms"].split(", ")]) ) location = mut["HGNC_Chromosome"] uniprot_accession = mut["UniProt_uniprot_accession"] transcripts = mut["transcripts"] self.assertTrue(transcripts is not None) self.assertTrue(len(transcripts) > 0) self.assertTrue(transcripts.startswith("ENST")) strand = mut["strand"] klass = mut["class"] uniprot_experimentals = mut["UniProt_AA_experimental_info"].split("|") self.assertTrue(uniprot_experimentals is not None) uniprot_natural_variations = mut["UniProt_AA_natural_variation"].split("|") uniprot_regions = mut["UniProt_AA_region"].split("|") uniprot_sites = mut["UniProt_AA_site"].split("|") uniprot_go_biological_processes = mut["UniProt_GO_Biological_Process"].split("|") uniprot_go_cellular_components = mut["UniProt_GO_Cellular_Component"].split("|") self.assertTrue(uniprot_go_cellular_components is not None) uniprot_go_molecular_functions = mut["UniProt_GO_Molecular_Function"].split("|") pass
def test_simple_genes_by_region_annotation(self): """Test web api backend call /genes/ """ # http://www.broadinstitute.org/oncotator/genes/chr22_22112223_22312558/ # Two genes: chr22:22,112,223-22,312,558 datasource_list = DatasourceFactory.createDatasources(self._determine_db_dir(), "hg19", isMulticore=False) annotator = Annotator() for ds in datasource_list: annotator.addDatasource(ds) # Here is what the API would call.... txs = annotator.retrieve_transcripts_by_region("22", 22112223, 22312558) self.assertTranscriptsFound(txs) mut_dict = annotator.annotate_genes_given_txs(txs) # Each mut will be for a separate gene for gene in mut_dict.keys(): mut = mut_dict[gene] alt_accessions = mut['UniProt_alt_uniprot_accessions'].split("|") tcgascape_amp_peaks = mut['TCGAScape_Amplification_Peaks'].split("|") tcgascape_del_peaks = mut['TCGAScape_Deletion_Peaks'].split("|") tumorscape_amp_peaks = mut['TUMORScape_Amplification_Peaks'].split("|") tumorscape_del_peaks = mut['TUMORScape_Deletion_Peaks'].split("|") full_name = mut['HGNC_Approved Name'] cosmic = {"tissue_types_affected": mut['COSMIC_Tissue_tissue_types_affected'], "total_alterations_in_gene": mut["COSMIC_Tissue_tissue_types_affected"]} alt_aliases = list(itertools.chain([mut["HGNC_Previous Symbols"].split(", "), mut["HGNC_Synonyms"].split(", ")])) location = mut["HGNC_Chromosome"] uniprot_accession = mut["UniProt_uniprot_accession"] transcripts = mut['transcripts'] self.assertTrue(transcripts is not None) self.assertTrue(len(transcripts) > 0) self.assertTrue(transcripts.startswith('ENST')) strand = mut['strand'] klass = mut['class'] uniprot_experimentals = mut['UniProt_AA_experimental_info'].split("|") self.assertTrue(uniprot_experimentals is not None) uniprot_natural_variations = mut['UniProt_AA_natural_variation'].split("|") uniprot_regions = mut['UniProt_AA_region'].split("|") uniprot_sites = mut['UniProt_AA_site'].split("|") uniprot_go_biological_processes = mut["UniProt_GO_Biological_Process"].split("|") uniprot_go_cellular_components = mut["UniProt_GO_Cellular_Component"].split("|") self.assertTrue(uniprot_go_cellular_components is not None) uniprot_go_molecular_functions = mut["UniProt_GO_Molecular_Function"].split("|") pass