phc-ingestion 0.8.26__tar.gz → 0.8.27__tar.gz

{phc-ingestion-0.8.26 → phc-ingestion-0.8.27}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: phc-ingestion
-Version: 0.8.26
+Version: 0.8.27
 Summary: Functions for LifeOmic PHC genomic ingestions
 License: MIT
 Author-email: LifeOmic Development <development@lifeomic.com>

phc-ingestion-0.8.27/ingestion/nextgen/util/nextgen_specific_genes.py

@@ -0,0 +1,28 @@
+from typing import TypedDict
+
+
+class GeneWithLocation(TypedDict):
+    gene: str
+    chr: str
+    start: int
+    end: int
+
+
+nextgen_specific_genes: list[GeneWithLocation] = [
+    {"gene": "IGK", "chr": "chr2", "start": 88852034, "end": 90258119},
+    {"gene": "NSD2", "chr": "chr4", "start": 1792518, "end": 1940193},
+    {"gene": "CCND3", "chr": "chr6", "start": 41920534, "end": 42562008},
+    {"gene": "MYC", "chr": "chr8", "start": 125309416, "end": 129673293},
+    {"gene": "CCND1", "chr": "chr11", "start": 69090733, "end": 69656860},
+    {"gene": "IGH", "chr": "chr14", "start": 105578834, "end": 109902208},
+    {"gene": "MAF", "chr": "chr16", "start": 78428398, "end": 79615096},
+    {"gene": "MAFB", "chr": "chr20", "start": 39039005, "end": 40688948},
+    {"gene": "IGL", "chr": "chr22", "start": 22012552, "end": 22965858},
+]
+
+
+def maybe_get_matching_gene_for_location(chr: str, position: int) -> str | None:
+    for gene in nextgen_specific_genes:
+        if gene["chr"] == chr and gene["start"] <= position <= gene["end"]:
+            return gene["gene"]
+    return None

{phc-ingestion-0.8.26 → phc-ingestion-0.8.27}/ingestion/nextgen/util/process_structural.py

@@ -5,6 +5,7 @@ from typing import TypedDict
 from ingestion.shared_util.coords_to_genes import coords_to_genes
 from ingestion.nextgen.util.alteration_table import extract_variant_table
 from ingestion.nextgen.util.interpretation import map_interpretation
+from ingestion.nextgen.util.nextgen_specific_genes import maybe_get_matching_gene_for_location
 from ingestion.shared_util.open_maybe_gzipped import open_maybe_gzipped
 
 
@@ -42,6 +43,23 @@ def is_del_dup_or_ins(variant: list[str]) -> bool:
     return any([x in variant[2] for x in ["MantaDEL", "MantaDUP", "MantaINS"]])
 
 
+def get_gene_from_coords(
+    chromosome: str, start_position: str, end_position: str, log: Logger
+) -> str:
+    """
+    A number of genes of interest with specific start and end positions have been provided.
+    If a variant falls within the start and end positions of one of those genes of interest, that gene will be used.
+    Otherwise, we fall back to the standard gene lookup.
+    """
+    center_position = int((int(start_position) + int(end_position)) / 2)
+
+    gene = maybe_get_matching_gene_for_location(chromosome, center_position)
+    if gene:
+        return gene
+
+    return coords_to_genes("GRCh38", chromosome, center_position, log)
+
+
 def process_structural(
     sv_in_file: str, xml_in_file, root_path: str, prefix: str, log: Logger
 ) -> str | None:
@@ -74,9 +92,7 @@ def process_structural(
                 effect = "insertion"
 
             # Get genes from coordinates using center point of start and end positions
-            gene1 = coords_to_genes(
-                "GRCh38", chromosome1, int((int(start_position1) + int(end_position1)) / 2), log
-            )
+            gene1 = get_gene_from_coords(chromosome1, start_position1, end_position1, log)
             gene2 = "N/A"
 
         else:
@@ -89,12 +105,8 @@ def process_structural(
             effect = "translocation"
 
             # Get genes from coordinates using center point of start and end positions
-            gene1 = coords_to_genes(
-                "GRCh38", chromosome1, int((int(start_position1) + int(end_position1)) / 2), log
-            )
-            gene2 = coords_to_genes(
-                "GRCh38", chromosome2, int((int(start_position2) + int(end_position2)) / 2), log
-            )
+            gene1 = get_gene_from_coords(chromosome1, start_position1, end_position1, log)
+            gene2 = get_gene_from_coords(chromosome2, start_position2, end_position2, log)
 
         # Scrape interpretation
         interpretation = "unknown"

{phc-ingestion-0.8.26 → phc-ingestion-0.8.27}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "phc-ingestion"
-version = "0.8.26"
+version = "0.8.27"
 description = "Functions for LifeOmic PHC genomic ingestions"
 authors = [
     { name = "LifeOmic Development", email = "development@lifeomic.com" },