diff --git a/bin/ontobio-parse-assocs.py b/bin/ontobio-parse-assocs.py index f74174da..b465af85 100755 --- a/bin/ontobio-parse-assocs.py +++ b/bin/ontobio-parse-assocs.py @@ -80,7 +80,9 @@ def main(): parser.add_argument("-g", "--gpi", type=str, required=False, default=None, help="GPI file") parser.add_argument("-m", "--metadata_dir", type=dir_path, required=False, - help="Path to metadata directory") + help="Path to metadata directory") + parser.add_argument("--retracted_pub_set", type=argparse.FileType('r'), required=False, + help="Path to retracted publications file") parser.add_argument("-l", "--rule", action="append", required=None, default=[], dest="rule_set", help="Set of rules to be run. Default is no rules to be run, with the exception \ of gorule-0000027 and gorule-0000020. See command line documentation in the \ @@ -143,11 +145,17 @@ def main(): rule_set = assocparser.RuleSet.ALL goref_metadata = None - ref_species_metadata = None + ref_species_metadata = None if args.metadata_dir: absolute_metadata = os.path.abspath(args.metadata_dir) goref_metadata = metadata.yamldown_lookup(os.path.join(absolute_metadata, "gorefs")) ref_species_metadata = metadata.yaml_set(absolute_metadata, "go-reference-species.yaml", "taxon_id") + + retracted_pub_set = None + if args.retracted_pub_set: + retracted_pub_set = metadata.retracted_pub_set(args.retracted_pub_set.name) + elif args.metadata_dir: + retracted_pub_set = metadata.retracted_pub_set_from_meta(absolute_metadata) # set configuration filtered_evidence_file = open(args.filtered_file, "w") if args.filtered_file else None @@ -164,6 +172,7 @@ def main(): gpi_authority_path=args.gpi, goref_metadata=goref_metadata, ref_species_metadata=ref_species_metadata, + retracted_pub_set=retracted_pub_set, rule_set=rule_set ) p = None diff --git a/bin/validate.py b/bin/validate.py index ed024015..25cd0f83 100755 --- a/bin/validate.py +++ b/bin/validate.py @@ -34,6 +34,7 @@ from ontobio.validation import tools from ontobio.validation import rules + from typing import Dict, Set # logging.basicConfig(format="%(asctime)s - %(name)s - %(levelname)s: %(message)s", level=logging.WARNING) @@ -210,7 +211,7 @@ def create_parser(config, group, dataset, format="gaf"): """ @tools.gzips -def produce_gaf(dataset, source_gaf, ontology_graph, gpipaths=None, paint=False, group="unknown", rule_metadata=None, goref_metadata=None, ref_species_metadata=None, db_entities=None, group_idspace=None, format="gaf", suppress_rule_reporting_tags=[], annotation_inferences=None, group_metadata=None, extensions_constraints=None, rule_contexts=[], gaf_output_version="2.2", rule_set=assocparser.RuleSet.ALL): +def produce_gaf(dataset, source_gaf, ontology_graph, gpipaths=None, paint=False, group="unknown", rule_metadata=None, goref_metadata=None, ref_species_metadata=None, retracted_pub_set=None,db_entities=None, group_idspace=None, format="gaf", suppress_rule_reporting_tags=[], annotation_inferences=None, group_metadata=None, extensions_constraints=None, rule_contexts=[], gaf_output_version="2.2", rule_set=assocparser.RuleSet.ALL): filtered_associations = open(os.path.join(os.path.split(source_gaf)[0], "{}_noiea.gaf".format(dataset)), "w") config = assocparser.AssocParserConfig( ontology=ontology_graph, @@ -221,6 +222,7 @@ def produce_gaf(dataset, source_gaf, ontology_graph, gpipaths=None, paint=False, rule_metadata=rule_metadata, goref_metadata=goref_metadata, ref_species_metadata=ref_species_metadata, + retracted_pub_set=retracted_pub_set, entity_idspaces=db_entities, group_idspace=group_idspace, suppress_rule_reporting_tags=suppress_rule_reporting_tags, @@ -493,7 +495,8 @@ def cli(ctx, verbose): @click.option("--only-dataset", default=None) @click.option("--gaf-output-version", default="2.2", type=click.Choice(["2.1", "2.2"])) @click.option("--rule-set", "-l", "rule_set", default=[assocparser.RuleSet.ALL], multiple=True) -def produce(ctx, group, metadata_dir, gpad, ttl, target, ontology, exclude, base_download_url, suppress_rule_reporting_tag, skip_existing_files, gaferencer_file, only_dataset, gaf_output_version, rule_set): +@click.option("--retracted_pub_set", type=click.Path(exists=True), default=None, required=False, help="Path to retracted publications file") +def produce(ctx, group, metadata_dir, gpad, ttl, target, ontology, exclude, base_download_url, suppress_rule_reporting_tag, skip_existing_files, gaferencer_file, only_dataset, gaf_output_version, rule_set, retracted_pub_set): logger.info("Logging is verbose") products = { @@ -529,7 +532,7 @@ def produce(ctx, group, metadata_dir, gpad, ttl, target, ontology, exclude, base db_entities = metadata.database_entities(absolute_metadata) group_ids = metadata.groups(absolute_metadata) - extensions_constraints = metadata.extensions_constraints_file(absolute_metadata) + extensions_constraints = metadata.extensions_constraints_file(absolute_metadata) gaferences = None if gaferencer_file: @@ -539,6 +542,12 @@ def produce(ctx, group, metadata_dir, gpad, ttl, target, ontology, exclude, base if rule_set == (assocparser.RuleSet.ALL,): rule_set = assocparser.RuleSet.ALL + retracted_pubs = None + if retracted_pub_set: + retracted_pubs = metadata.retracted_pub_set(retracted_pub_set) + else: + retracted_pubs = metadata.retracted_pub_set_from_meta(absolute_metadata) + for dataset_metadata, source_gaf in downloaded_gaf_sources: dataset = dataset_metadata["dataset"] # Set paint to True when the group is "paint". @@ -550,6 +559,7 @@ def produce(ctx, group, metadata_dir, gpad, ttl, target, ontology, exclude, base rule_metadata=rule_metadata, goref_metadata=goref_metadata, ref_species_metadata=ref_species_metadata, + retracted_pub_set=retracted_pubs, db_entities=db_entities, group_idspace=group_ids, suppress_rule_reporting_tags=suppress_rule_reporting_tag, @@ -634,13 +644,14 @@ def paint(group, dataset, metadata, target, ontology): absolute_target = os.path.abspath(target) os.makedirs(os.path.join(absolute_target, "groups"), exist_ok=True) paint_metadata = metadata.dataset_metadata_file(absolute_metadata, "paint") + paint_src_gaf = check_and_download_mixin_source(paint_metadata, dataset, absolute_target) click.echo("Loading ontology: {}...".format(ontology)) ontology_graph = OntologyFactory().create(ontology) gpi_path = os.path.join(absolute_target, "groups", dataset, "{}.gpi".format(dataset)) - click.echo("Using GPI at {}".format(gpi_path)) + click.echo("Using GPI at {}".format(gpi_path)) paint_gaf = produce_gaf("paint_{}".format(dataset), paint_src_gaf, ontology_graph, gpipath=gpi_path) @@ -650,7 +661,8 @@ def paint(group, dataset, metadata, target, ontology): @click.option("--ontology", type=click.Path(), required=True) @click.option("--gaferencer-file", "-I", type=click.Path(exists=True), default=None, required=False, help="Path to Gaferencer output to be used for inferences") -def rule(metadata_dir, out, ontology, gaferencer_file): +@click.option("--retracted_pub_set", type=click.Path(exists=True), default=None, required=False, help="Path to retracted publications file") +def rule(metadata_dir, out, ontology, gaferencer_file, retracted_pub_set): absolute_metadata = os.path.abspath(metadata_dir) click.echo("Loading ontology: {}...".format(ontology)) @@ -659,6 +671,12 @@ def rule(metadata_dir, out, ontology, gaferencer_file): goref_metadata = metadata.yamldown_lookup(os.path.join(absolute_metadata, "gorefs")) gorule_metadata = metadata.yamldown_lookup(os.path.join(absolute_metadata, "rules")) ref_species_metadata = metadata.yaml_set(absolute_metadata, "go-reference-species.yaml", "taxon_id") + retracted_pubs = None + if retracted_pub_set: + retracted_pubs = metadata.retracted_pub_set(retracted_pub_set) + else: + retracted_pubs = metadata.retracted_pub_set_from_meta(absolute_metadata) + click.echo("Found {} GO Rules".format(len(gorule_metadata.keys()))) @@ -673,6 +691,7 @@ def rule(metadata_dir, out, ontology, gaferencer_file): ontology=ontology_graph, goref_metadata=goref_metadata, ref_species_metadata=ref_species_metadata, + retracted_pub_set=retracted_pubs, entity_idspaces=db_entities, group_idspace=group_ids, annotation_inferences=gaferences, diff --git a/ontobio/io/assocparser.py b/ontobio/io/assocparser.py index 8f00c72e..bf0a490e 100644 --- a/ontobio/io/assocparser.py +++ b/ontobio/io/assocparser.py @@ -234,6 +234,7 @@ def __init__(self, ref_species_metadata=None, group_metadata=None, dbxrefs=None, + retracted_pub_set=None, suppress_rule_reporting_tags=[], annotation_inferences=None, extensions_constraints=None, @@ -258,6 +259,7 @@ def __init__(self, self.goref_metadata = goref_metadata self.ref_species_metadata = ref_species_metadata self.group_metadata = group_metadata + self.retracted_pub_set = retracted_pub_set self.suppress_rule_reporting_tags = suppress_rule_reporting_tags self.annotation_inferences = annotation_inferences self.entity_idspaces = entity_idspaces diff --git a/ontobio/io/qc.py b/ontobio/io/qc.py index f15e73f0..f08668a2 100644 --- a/ontobio/io/qc.py +++ b/ontobio/io/qc.py @@ -421,6 +421,20 @@ def test(self, annotation: association.GoAssociation, config: assocparser.AssocP return self._result(bool(withfrom)) else: return self._result(True) + +class GoRule22(GoRule): + + def __init__(self): + super().__init__("GORULE:0000022", "Check for, and filter, annotations made to retracted publications", FailMode.HARD) + + def test(self, annotation: association.GoAssociation, config: assocparser.AssocParserConfig, group=None) -> TestResult: + if config.retracted_pub_set is not None: + references = annotation.evidence.has_supporting_reference + for ref in references: + ref = str(ref) + if ref in config.retracted_pub_set: + return self._result(False) + return self._result(True) class GoRule26(GoRule): @@ -952,6 +966,7 @@ def test(self, annotation: association.GoAssociation, config: assocparser.AssocP "GoRule16": GoRule16(), "GoRule17": GoRule17(), "GoRule18": GoRule18(), + "GoRule22": GoRule22(), "GoRule26": GoRule26(), "GoRule28": GoRule28(), "GoRule29": GoRule29(), diff --git a/ontobio/validation/metadata.py b/ontobio/validation/metadata.py index 8b370635..b5be0c1d 100644 --- a/ontobio/validation/metadata.py +++ b/ontobio/validation/metadata.py @@ -182,7 +182,32 @@ def yaml_set(metadata, yaml_file_name, field) -> Set[str]: except Exception as e: raise click.ClickException("Could not find or read {}: {}".format(yaml_path, str(e))) - return set([yaml[field] for yaml in yaml_list]) + return set([yaml[field] for yaml in yaml_list]) + +def retracted_pub_set_from_meta(metadata) -> Set: + retracted_path = os.path.join(metadata, "retracted-publications.txt") + if os.access(retracted_path, os.R_OK): + return retracted_pub_set_use_abspath(retracted_path) + else: + return set() + +def retracted_pub_set(abspath_retracted_file) -> Set: + return retracted_pub_set_use_abspath(os.path.abspath(abspath_retracted_file)) + +def retracted_pub_set_use_abspath(abspath_retracted_file) -> Set: + try: + retracted_pubs = None + with open(abspath_retracted_file, "r") as f: + retracted_pubs = set() + for line in f: + li=line.strip() + if not li.startswith("!"): + if "," in li: + li = li.partition(',')[0] + retracted_pubs.add(li) + return retracted_pubs + except Exception as e: + raise click.ClickException("Could not find or read {}: {}".format(abspath_retracted_file, str(e))) diff --git a/tests/test_qc.py b/tests/test_qc.py index 90d92368..ef4bd800 100644 --- a/tests/test_qc.py +++ b/tests/test_qc.py @@ -354,6 +354,20 @@ def test_go_rule_18(): test_result = qc.GoRule18().test(assoc, all_rules_config()) assert test_result.result_type == qc.ResultType.PASS +def test_go_rule22(): + config = assocparser.AssocParserConfig( + ontology=ontology, + retracted_pub_set={"RETRACTED:1234","PMID:37772366"}, + rule_set=assocparser.RuleSet.ALL + ) + assoc = make_annotation(goid="GO:1234567", evidence="IBA", references="PMID:12345").associations[0] + test_result = qc.GoRule22().test(assoc, config) + assert test_result.result_type == qc.ResultType.PASS + + assoc = make_annotation(goid="GO:1234567", evidence="IBA", references="PMID:37772366").associations[0] + test_result = qc.GoRule22().test(assoc, config) + assert test_result.result_type == qc.ResultType.ERROR + def test_go_rule26(): config = assocparser.AssocParserConfig( @@ -819,7 +833,7 @@ def test_all_rules(): assoc = gafparser.to_association(a).associations[0] test_results = qc.test_go_rules(assoc, config).all_results - assert len(test_results.keys()) == 26 + assert len(test_results.keys()) == 27 assert test_results[qc.GoRules.GoRule26.value].result_type == qc.ResultType.PASS assert test_results[qc.GoRules.GoRule29.value].result_type == qc.ResultType.PASS