cool-seq-tool 0.5.1__py3-none-any.whl → 0.7.0__py3-none-any.whl

cool_seq_tool/mappers/mane_transcript.py

@@ -25,7 +25,9 @@ from cool_seq_tool.mappers.liftover import LiftOver
 from cool_seq_tool.schemas import (
     AnnotationLayer,
     Assembly,
-    ResidueMode,
+    CoordinateType,
+    GenomicTxMetadata,
+    ManeGeneData,
     Strand,
     TranscriptPriority,
 )
@@ -71,10 +73,10 @@ class CdnaRepresentation(DataRepresentation):
 class GenomicRepresentation(BaseModel):
     """Define object model for genomic representation"""
 
-    refseq: str
     pos: tuple[int, int]
-    status: TranscriptPriority
-    alt_ac: str
+    mane_genes: list[ManeGeneData] = []
+    status: Literal["grch38"] = TranscriptPriority.GRCH38.value
+    ac: str
 
 
 class ProteinAndCdnaRepresentation(BaseModel):
@@ -100,7 +102,7 @@ class ManeTranscript:
         A handful of resources are required for initialization, so when defaults are
         enough, it's easiest to let the core CoolSeqTool class handle it for you:
 
-        >>> from cool_seq_tool.app import CoolSeqTool
+        >>> from cool_seq_tool import CoolSeqTool
         >>> mane_mapper = CoolSeqTool().mane_transcript
 
         Note that most methods are defined as Python coroutines, so they must be called
@@ -108,7 +110,7 @@ class ManeTranscript:
 
         >>> import asyncio
         >>> result = asyncio.run(mane_mapper.g_to_grch38("NC_000001.11", 100, 200))
-        >>> result["ac"]
+        >>> result.ac
         'NC_000001.11'
 
         See the :ref:`Usage section <async_note>` for more information.
@@ -128,7 +130,7 @@ class ManeTranscript:
         self.liftover = liftover
 
     @staticmethod
-    def _get_reading_frame(pos: int) -> int:
+    def get_reading_frame(pos: int) -> int:
         """Return reading frame number. Only used on c. coordinate.
 
         :param pos: cDNA position
@@ -181,13 +183,12 @@ class ManeTranscript:
         pos = self._p_to_c_pos(start_pos, end_pos)
         return ac, pos
 
-    async def _c_to_g(self, ac: str, pos: tuple[int, int]) -> dict | None:
+    async def _c_to_g(self, ac: str, pos: tuple[int, int]) -> GenomicTxMetadata | None:
         """Get g. annotation from c. annotation.
 
         :param ac: cDNA accession
         :param pos: [cDNA pos start, cDNA pos end]
-        :return: Gene, Transcript accession and position change,
-            Altered transcript accession and position change, Strand
+        :return: Metadata for genomic and transcript accessions
         """
         # UTA does not store ENST versions
         # So we want to make sure version is valid
@@ -219,13 +220,13 @@ class ManeTranscript:
             ac, pos, AnnotationLayer.CDNA, coding_start_site=coding_start_site
         )
 
-    async def _liftover_to_38(self, genomic_tx_data: dict) -> None:
+    async def _liftover_to_38(self, genomic_tx_data: GenomicTxMetadata) -> None:
         """Liftover genomic_tx_data to hg38 assembly.
 
-        :param genomic_tx_data: Dictionary containing gene, nc_accession, alt_pos, and
-            strand. This will be mutated in-place if not GRCh38 assembly.
+        :param genomic_tx_data: Metadata for genomic and transcript accessions. This
+            will be mutated in-place if not GRCh38 assembly.
         """
-        descr = await self.uta_db.get_chr_assembly(genomic_tx_data["alt_ac"])
+        descr = await self.uta_db.get_chr_assembly(genomic_tx_data.alt_ac)
         if descr is None:
             # already grch38
             return
@@ -234,14 +235,14 @@ class ManeTranscript:
         query = f"""
             SELECT DISTINCT alt_ac
             FROM {self.uta_db.schema}.tx_exon_aln_v
-            WHERE tx_ac = '{genomic_tx_data['tx_ac']}';
+            WHERE tx_ac = '{genomic_tx_data.tx_ac}';
             """  # noqa: S608
         nc_acs = await self.uta_db.execute_query(query)
         nc_acs = [nc_ac[0] for nc_ac in nc_acs]
-        if nc_acs == [genomic_tx_data["alt_ac"]]:
+        if nc_acs == [genomic_tx_data.alt_ac]:
             _logger.warning(
                 "UTA does not have GRCh38 assembly for %s",
-                genomic_tx_data["alt_ac"].split(".")[0],
+                genomic_tx_data.alt_ac.split(".")[0],
             )
             return
 
@@ -257,7 +258,7 @@ class ManeTranscript:
         )
 
         # Change alt_ac to most recent
-        if genomic_tx_data["alt_ac"].startswith("EN"):
+        if genomic_tx_data.alt_ac.startswith("EN"):
             order_by_cond = "ORDER BY alt_ac DESC;"
         else:
             order_by_cond = """
@@ -267,50 +268,49 @@ class ManeTranscript:
         query = f"""
             SELECT alt_ac
             FROM {self.uta_db.schema}.genomic
-            WHERE alt_ac LIKE '{genomic_tx_data['alt_ac'].split('.')[0]}%'
+            WHERE alt_ac LIKE '{genomic_tx_data.alt_ac.split('.')[0]}%'
             {order_by_cond}
             """  # noqa: S608
         nc_acs = await self.uta_db.execute_query(query)
-        genomic_tx_data["alt_ac"] = nc_acs[0][0]
+        genomic_tx_data.alt_ac = nc_acs[0][0]
 
     def _set_liftover(
         self,
-        genomic_tx_data: dict,
+        genomic_tx_data: GenomicTxMetadata,
         key: str,
         chromosome: str,
         liftover_to_assembly: Assembly,
     ) -> None:
         """Update genomic_tx_data to have coordinates for given assembly.
 
-        :param genomic_tx_data: Dictionary containing gene, nc_accession, alt_pos, and
-            strand
+        :param genomic_tx_data: Metadata for genomic and transcript accessions
         :param key: Key to access coordinate positions
         :param chromosome: Chromosome, must be prefixed with ``chr``
         :param liftover_to_assembly: Assembly to liftover to
         """
+        coords = getattr(genomic_tx_data, key)
         liftover_start_i = self.liftover.get_liftover(
-            chromosome, genomic_tx_data[key][0], liftover_to_assembly
+            chromosome, coords[0], liftover_to_assembly
         )
         if liftover_start_i is None:
             _logger.warning(
                 "Unable to liftover position %s on %s",
-                genomic_tx_data[key][0],
+                coords[0],
                 chromosome,
             )
             return
 
         liftover_end_i = self.liftover.get_liftover(
-            chromosome, genomic_tx_data[key][1], liftover_to_assembly
+            chromosome, coords[1], liftover_to_assembly
         )
         if liftover_end_i is None:
             _logger.warning(
                 "Unable to liftover position %s on %s",
-                genomic_tx_data[key][1],
+                coords[1],
                 chromosome,
             )
             return
-
-        genomic_tx_data[key] = liftover_start_i[1], liftover_end_i[1]
+        setattr(genomic_tx_data, key, (liftover_start_i[1], liftover_end_i[1]))
 
     async def _get_and_validate_genomic_tx_data(
         self,
@@ -320,7 +320,7 @@ class ManeTranscript:
         | Literal[AnnotationLayer.GENOMIC] = AnnotationLayer.CDNA,
         coding_start_site: int | None = None,
         alt_ac: str | None = None,
-    ) -> dict | None:
+    ) -> GenomicTxMetadata | None:
         """Get and validate genomic_tx_data
 
         :param tx_ac: Accession on c. coordinate
@@ -328,7 +328,8 @@ class ManeTranscript:
         :param annotation_layer: Annotation layer for ``ac`` and ``pos``
         :param coding_start_site: Coding start site
         :param alt_ac: Accession on g. coordinate
-        :return: genomic_tx_data if found and validated, else None
+        :return: Metadata for genomic and transcript accessions if found and validated,
+            else None
         """
         genomic_tx_data = await self.uta_db.get_genomic_tx_data(
             tx_ac, pos, annotation_layer, alt_ac=alt_ac
@@ -341,14 +342,14 @@ class ManeTranscript:
                 annotation_layer,
             )
             return None
-        genomic_tx_data["coding_start_site"] = coding_start_site
+        genomic_tx_data.coding_start_site = coding_start_site
 
         if not alt_ac:
             # Only want to liftover if alt_ac not provided. If alt_ac is provided,
             # it means user wants to stick with the queried assembly
-            og_alt_exon_id = genomic_tx_data["alt_exon_id"]
+            og_alt_exon_id = genomic_tx_data.alt_exon_id
             await self._liftover_to_38(genomic_tx_data)
-            liftover_alt_exon_id = genomic_tx_data["alt_exon_id"]
+            liftover_alt_exon_id = genomic_tx_data.alt_exon_id
 
             # Validation check: Exon structure
             if og_alt_exon_id != liftover_alt_exon_id:
@@ -466,14 +467,14 @@ class ManeTranscript:
         :return: Transcript data
         """
         if found_result:
-            tx_g_pos = g["alt_pos_range"]
-            tx_pos_range = g["tx_pos_range"]
+            tx_g_pos = g.alt_pos_range
+            tx_pos_range = g.tx_pos_range
         else:
             result = await self.uta_db.get_tx_exon_aln_v_data(
                 refseq_c_ac,
-                g["alt_pos_change_range"][0],
-                g["alt_pos_change_range"][1],
-                alt_ac=alt_ac if alt_ac else g["alt_ac"],
+                g.alt_pos_change_range[0],
+                g.alt_pos_change_range[1],
+                alt_ac=alt_ac if alt_ac else g.alt_ac,
                 use_tx_pos=False,
             )
 
@@ -483,18 +484,18 @@ class ManeTranscript:
                 )
                 return None
             result = result[-1]
-            tx_g_pos = result[5], result[6]  # alt_start_i, alt_end_i
-            tx_pos_range = result[2], result[3]  # tx_start_i, tx_end_i
+            tx_g_pos = result.alt_start_i, result.alt_end_i
+            tx_pos_range = result.tx_start_i, result.tx_end_i
 
         cds_start_end = await self.uta_db.get_cds_start_end(refseq_c_ac)
         if not cds_start_end:
             return None
         coding_start_site = cds_start_end[0]
 
-        g_pos = g["alt_pos_change_range"]  # start/end genomic change
+        g_pos = g.alt_pos_change_range  # start/end genomic change
         g_pos_change = g_pos[0] - tx_g_pos[0], tx_g_pos[1] - g_pos[1]
 
-        if g["strand"] == Strand.NEGATIVE:
+        if g.strand == Strand.NEGATIVE:
             g_pos_change = (tx_g_pos[1] - g_pos[0], g_pos[1] - tx_g_pos[0])
 
         c_pos_change = (
@@ -506,10 +507,10 @@ class ManeTranscript:
             c_pos_change = c_pos_change[1], c_pos_change[0]
 
         return self._get_c_data(
-            gene=g["gene"],
+            gene=g.gene,
             cds_start_end=cds_start_end,
             c_pos_change=c_pos_change,
-            strand=g["strand"],
+            strand=g.strand,
             alt_ac=alt_ac,
             status=status,
             refseq_c_ac=refseq_c_ac,
@@ -531,8 +532,8 @@ class ManeTranscript:
         """
         for pos, pos_index in [(start_pos, 0), (end_pos, 1)]:
             if pos is not None:
-                og_rf = self._get_reading_frame(pos)
-                new_rf = self._get_reading_frame(transcript_data.pos[pos_index])
+                og_rf = self.get_reading_frame(pos)
+                new_rf = self.get_reading_frame(transcript_data.pos[pos_index])
 
                 if og_rf != new_rf:
                     _logger.warning(
@@ -561,7 +562,7 @@ class ManeTranscript:
         | GenomicRepresentation,
         expected_ref: str,
         anno: AnnotationLayer,
-        residue_mode: ResidueMode,
+        coordinate_type: CoordinateType,
     ) -> bool:
         """Return whether or not reference changes are the same.
 
@@ -573,7 +574,7 @@ class ManeTranscript:
             position change
         :param expected_ref: Reference at position given during input
         :param anno: Annotation layer we are starting from
-        :param residue_mode: Residue mode for ``start_pos`` and ``end_pos``
+        :param coordinate_type: Coordinate type for ``start_pos`` and ``end_pos``
         :return: ``True`` if reference check passes. ``False`` otherwise.
         """
         if anno == AnnotationLayer.CDNA:
@@ -581,7 +582,7 @@ class ManeTranscript:
             end_pos += coding_start_site
 
         ref, _ = self.seqrepo_access.get_reference_sequence(
-            ac, start=start_pos, end=end_pos, residue_mode=residue_mode
+            ac, start=start_pos, end=end_pos, coordinate_type=coordinate_type
         )
         if ref is None:
             return False
@@ -597,7 +598,7 @@ class ManeTranscript:
                 mane_transcript.refseq,
                 start=mane_start_pos,
                 end=mane_end_pos if mane_start_pos != mane_end_pos else None,
-                residue_mode=residue_mode,
+                coordinate_type=coordinate_type,
             )
             if not mane_ref:
                 _logger.info("Unable to validate reference for MANE Transcript")
@@ -618,7 +619,7 @@ class ManeTranscript:
 
         return True
 
-    def _validate_index(
+    def validate_index(
         self, ac: str, pos: tuple[int, int], coding_start_site: int
     ) -> bool:
         """Validate that positions actually exist on accession
@@ -632,7 +633,10 @@ class ManeTranscript:
         end_pos = pos[1] + coding_start_site
         return bool(
             self.seqrepo_access.get_reference_sequence(
-                ac, start=start_pos, end=end_pos, residue_mode=ResidueMode.INTER_RESIDUE
+                ac,
+                start=start_pos,
+                end=end_pos,
+                coordinate_type=CoordinateType.INTER_RESIDUE,
             )[0]
         )
 
@@ -689,7 +693,7 @@ class ManeTranscript:
         start_annotation_layer: AnnotationLayer,
         gene: str | None = None,
         ref: str | None = None,
-        residue_mode: ResidueMode = ResidueMode.RESIDUE,
+        coordinate_type: CoordinateType = CoordinateType.RESIDUE,
         mane_transcripts: set | None = None,
         alt_ac: str | None = None,
         end_annotation_layer: EndAnnotationLayer | None = None,
@@ -699,8 +703,8 @@ class ManeTranscript:
         information.
 
         >>> import asyncio
-        >>> from cool_seq_tool.app import CoolSeqTool
-        >>> from cool_seq_tool.schemas import AnnotationLayer, ResidueMode
+        >>> from cool_seq_tool import CoolSeqTool
+        >>> from cool_seq_tool.schemas import AnnotationLayer, CoordinateType
         >>> mane_mapper = CoolSeqTool().mane_transcript
         >>> mane_transcripts = {
         ...     "ENST00000646891.2",
@@ -714,7 +718,7 @@ class ManeTranscript:
         ...         599,
         ...         gene="BRAF",
         ...         start_annotation_layer=AnnotationLayer.PROTEIN,
-        ...         residue_mode=ResidueMode.INTER_RESIDUE,
+        ...         coordinate_type=CoordinateType.INTER_RESIDUE,
         ...         mane_transcripts=mane_transcripts,
         ...     )
         ... )
@@ -731,7 +735,7 @@ class ManeTranscript:
         :param start_annotation_layer: Starting annotation layer
         :param gene: HGNC gene symbol
         :param ref: Reference at position given during input
-        :param residue_mode: Residue mode for ``start_pos`` and ``end_pos``
+        :param coordinate_type: Coordinate type for ``start_pos`` and ``end_pos``
         :param mane_transcripts: Attempted mane transcripts that were not compatible
         :param alt_ac: Genomic accession
         :param end_annotation_layer: The end annotation layer. If not provided, will be
@@ -767,8 +771,8 @@ class ManeTranscript:
             )
 
         lcr_result = None
-        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, residue_mode)
-        residue_mode = ResidueMode.INTER_RESIDUE
+        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, coordinate_type)
+        coordinate_type = CoordinateType.INTER_RESIDUE
 
         is_p_or_c_start_anno = True
         if start_annotation_layer == AnnotationLayer.PROTEIN:
@@ -856,7 +860,7 @@ class ManeTranscript:
                         {},
                         ref,
                         AnnotationLayer.PROTEIN,
-                        residue_mode,
+                        coordinate_type,
                     )
                 elif start_annotation_layer == AnnotationLayer.CDNA:
                     valid_references = self._validate_references(
@@ -867,7 +871,7 @@ class ManeTranscript:
                         {},
                         ref,
                         AnnotationLayer.CDNA,
-                        residue_mode,
+                        coordinate_type,
                     )
                 else:
                     valid_references = self._validate_references(
@@ -878,7 +882,7 @@ class ManeTranscript:
                         {},
                         ref,
                         AnnotationLayer.GENOMIC,
-                        residue_mode,
+                        coordinate_type,
                     )
 
                 if not valid_references:
@@ -902,7 +906,7 @@ class ManeTranscript:
                         gene,
                         row["pro_ac"],
                         lcr_c_data.pos,
-                        g["strand"],
+                        g.strand,
                         lcr_c_data.status,
                     )
                     coding_start_site = 0
@@ -910,7 +914,7 @@ class ManeTranscript:
                 ac = lcr_result.refseq or lcr_result.ensembl
                 pos = lcr_result.pos
 
-                if not self._validate_index(ac, pos, coding_start_site):
+                if not self.validate_index(ac, pos, coding_start_site):
                     _logger.warning(
                         "%s are not valid positions on %s with coding start site %s",
                         pos,
@@ -924,7 +928,7 @@ class ManeTranscript:
                     gene,
                     row["pro_ac"],
                     lcr_c_data.pos,
-                    g["strand"],
+                    g.strand,
                     lcr_c_data.status,
                 ),
                 cdna=lcr_c_data,
@@ -936,7 +940,7 @@ class ManeTranscript:
                 cds = lcr_result_dict[k].get("coding_start_site", 0)
                 ac = lcr_result_dict[k]["refseq"] or lcr_result_dict[k]["ensembl"]
                 pos = lcr_result_dict[k]["pos"]
-                if not self._validate_index(ac, pos, cds):
+                if not self.validate_index(ac, pos, cds):
                     valid = False
                     _logger.warning(
                         "%s are not valid positions on %s with coding start site %s",
@@ -959,13 +963,22 @@ class ManeTranscript:
         gene: str | None = None,
         ref: str | None = None,
         try_longest_compatible: bool = False,
-        residue_mode: Literal[ResidueMode.RESIDUE]
-        | Literal[ResidueMode.INTER_RESIDUE] = ResidueMode.RESIDUE,
+        coordinate_type: Literal[CoordinateType.RESIDUE]
+        | Literal[CoordinateType.INTER_RESIDUE] = CoordinateType.RESIDUE,
     ) -> DataRepresentation | CdnaRepresentation | None:
-        """Return MANE transcript.
-
-        >>> from cool_seq_tool.app import CoolSeqTool
-        >>> from cool_seq_tool.schemas import AnnotationLayer, ResidueMode
+        """Return MANE representation
+
+        If ``start_annotation_layer`` is ``AnnotationLayer.PROTEIN``, will return
+            ``AnnotationLayer.PROTEIN`` representation.
+        If ``start_annotation_layer`` is ``AnnotationLayer.CDNA``, will return
+            ``AnnotationLayer.CDNA`` representation.
+        If ``start_annotation_layer`` is ``AnnotationLayer.GENOMIC`` will return
+            ``AnnotationLayer.CDNA`` representation if ``gene`` is provided and
+            ``AnnotationLayer.GENOMIC`` GRCh38 representation if ``gene`` is NOT
+            provided.
+
+        >>> from cool_seq_tool import CoolSeqTool
+        >>> from cool_seq_tool.schemas import AnnotationLayer, CoordinateType
         >>> import asyncio
         >>> mane_mapper = CoolSeqTool().mane_transcript
         >>> result = asyncio.run(
@@ -973,7 +986,7 @@ class ManeTranscript:
         ...         "NP_004324.2",
         ...         599,
         ...         AnnotationLayer.PROTEIN,
-        ...         residue_mode=ResidueMode.INTER_RESIDUE,
+        ...         coordinate_type=CoordinateType.INTER_RESIDUE,
         ...     )
         ... )
         >>> result.gene, result.refseq, result.status
@@ -983,17 +996,21 @@ class ManeTranscript:
         :param start_pos: Start position change
         :param end_pos: End position change
         :param start_annotation_layer: Starting annotation layer.
-        :param gene: HGNC gene symbol
+        :param gene: HGNC gene symbol.
+            If ``gene`` is not provided and ``start_annotation_layer`` is
+            ``AnnotationLayer.GENOMIC``, will return GRCh38 representation.
+            If ``gene`` is provided and ``start_annotation_layer`` is
+            ``AnnotationLayer.GENOMIC``, will return cDNA representation.
         :param ref: Reference at position given during input
         :param try_longest_compatible: ``True`` if should try longest compatible remaining
             if mane transcript was not compatible. ``False`` otherwise.
-        :param ResidueMode residue_mode: Starting residue mode for ``start_pos`` and
-            ``end_pos``. Will always return coordinates in inter-residue
+        :param CoordinateType coordinate_type: Starting Coordinate type for
+            ``start_pos`` and ``end_pos``. Will always return inter-residue coordinates
         :return: MANE data or longest transcript compatible data if validation
             checks are correct. Will return inter-residue coordinates. Else, ``None``.
         """
-        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, residue_mode)
-        residue_mode = ResidueMode.INTER_RESIDUE
+        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, coordinate_type)
+        coordinate_type = CoordinateType.INTER_RESIDUE
         if ref:
             ref = ref[: end_pos - start_pos]
 
@@ -1012,7 +1029,7 @@ class ManeTranscript:
             if g is None:
                 return None
             # Get mane data for gene
-            mane_data = self.mane_transcript_mappings.get_gene_mane_data(g["gene"])
+            mane_data = self.mane_transcript_mappings.get_gene_mane_data(g.gene)
             if not mane_data:
                 return None
 
@@ -1039,10 +1056,8 @@ class ManeTranscript:
                 if not mane:
                     continue
 
-                if not mane.alt_ac:
-                    g_alt_ac = g.get("alt_ac")
-                    if g_alt_ac:
-                        mane.alt_ac = g_alt_ac
+                if not mane.alt_ac and g.alt_ac:
+                    mane.alt_ac = g.alt_ac
 
                 valid_reading_frame = self._validate_reading_frames(
                     c_ac, c_pos[0], c_pos[1], mane
@@ -1058,13 +1073,13 @@ class ManeTranscript:
                 if ref:
                     valid_references = self._validate_references(
                         ac,
-                        g["coding_start_site"],
+                        g.coding_start_site,
                         start_pos,
                         end_pos,
                         mane,
                         ref,
                         start_annotation_layer,
-                        residue_mode,
+                        coordinate_type,
                     )
                     if not valid_references:
                         continue
@@ -1078,8 +1093,8 @@ class ManeTranscript:
                         end_pos,
                         AnnotationLayer.PROTEIN,
                         ref=ref,
-                        gene=g["gene"],
-                        residue_mode=residue_mode,
+                        gene=g.gene,
+                        coordinate_type=coordinate_type,
                         mane_transcripts=mane_transcripts,
                     )
                 return await self.get_longest_compatible_transcript(
@@ -1088,34 +1103,61 @@ class ManeTranscript:
                     AnnotationLayer.CDNA,
                     ref=ref,
                     gene=g["gene"],
-                    residue_mode=residue_mode,
+                    coordinate_type=coordinate_type,
                     mane_transcripts=mane_transcripts,
                 )
             return None
         if start_annotation_layer == AnnotationLayer.GENOMIC:
+            if not gene:
+                return await self.g_to_grch38(
+                    ac,
+                    start_pos,
+                    end_pos,
+                    get_mane_genes=True,
+                    coordinate_type=coordinate_type,
+                )
+
             return await self.g_to_mane_c(
-                ac, start_pos, end_pos, gene=gene, residue_mode=residue_mode
+                ac, start_pos, end_pos, gene, coordinate_type=coordinate_type
             )
         _logger.warning("Annotation layer not supported: %s", start_annotation_layer)
         return None
 
-    async def g_to_grch38(self, ac: str, start_pos: int, end_pos: int) -> dict | None:
+    async def g_to_grch38(
+        self,
+        ac: str,
+        start_pos: int,
+        end_pos: int,
+        get_mane_genes: bool = False,
+        coordinate_type: CoordinateType = CoordinateType.RESIDUE,
+    ) -> GenomicRepresentation | None:
         """Return genomic coordinate on GRCh38 when not given gene context.
 
         :param ac: Genomic accession
         :param start_pos: Genomic start position
         :param end_pos: Genomic end position
-        :return: NC accession, start and end pos on GRCh38 assembly
+        :param get_mane_genes: ``True`` if mane genes for genomic position should be
+            included in response. ``False``, otherwise.
+        :param coordinate_type: Coordinate type for ``start_pos`` and ``end_pos``
+        :return: GRCh38 genomic representation (accession and start/end inter-residue
+            position)
         """
-        if end_pos is None:
-            end_pos = start_pos
+        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, coordinate_type)
 
         # Checking to see what chromosome and assembly we're on
         descr = await self.uta_db.get_chr_assembly(ac)
         if not descr:
             # Already GRCh38 assembly
-            if self._validate_index(ac, (start_pos, end_pos), 0):
-                return {"ac": ac, "pos": (start_pos, end_pos)}
+            if self.validate_index(ac, (start_pos, end_pos), 0):
+                return GenomicRepresentation(
+                    ac=ac,
+                    pos=(start_pos, end_pos),
+                    mane_genes=self.mane_transcript_mappings.get_genomic_mane_genes(
+                        ac, start_pos + 1, end_pos
+                    )
+                    if get_mane_genes
+                    else [],
+                )
             return None
         chromosome, assembly = descr
         is_same_pos = start_pos == end_pos
@@ -1145,13 +1187,21 @@ class ManeTranscript:
         newest_ac = await self.uta_db.get_newest_assembly_ac(ac)
         if newest_ac:
             ac = newest_ac[0]
-            if self._validate_index(ac, (start_pos, end_pos), 0):
-                return {"ac": ac, "pos": (start_pos, end_pos)}
+            if self.validate_index(ac, (start_pos, end_pos), 0):
+                return GenomicRepresentation(
+                    ac=ac,
+                    pos=(start_pos, end_pos),
+                    mane_genes=self.mane_transcript_mappings.get_genomic_mane_genes(
+                        ac, start_pos + 1, end_pos
+                    )
+                    if get_mane_genes
+                    else [],
+                )
         return None
 
     @staticmethod
     def get_mane_c_pos_change(
-        mane_tx_genomic_data: dict, coding_start_site: int
+        mane_tx_genomic_data: GenomicTxMetadata, coding_start_site: int
     ) -> tuple[int, int]:
         """Get mane c position change
 
@@ -1159,12 +1209,12 @@ class ManeTranscript:
         :param coding_start_site: Coding start site
         :return: cDNA pos start, cDNA pos end
         """
-        tx_pos_range = mane_tx_genomic_data["tx_pos_range"]
-        alt_pos_change = mane_tx_genomic_data["alt_pos_change"]
+        tx_pos_range = mane_tx_genomic_data.tx_pos_range
+        pos_change = mane_tx_genomic_data.pos_change
 
         mane_c_pos_change = (
-            tx_pos_range[0] + alt_pos_change[0] - coding_start_site,
-            tx_pos_range[1] - alt_pos_change[1] - coding_start_site,
+            tx_pos_range[0] + pos_change[0] - coding_start_site,
+            tx_pos_range[1] - pos_change[1] - coding_start_site,
         )
 
         if mane_c_pos_change[0] > mane_c_pos_change[1]:
@@ -1176,16 +1226,13 @@ class ManeTranscript:
         ac: str,
         start_pos: int,
         end_pos: int,
-        gene: str | None = None,
-        residue_mode: ResidueMode = ResidueMode.RESIDUE,
-    ) -> GenomicRepresentation | CdnaRepresentation | None:
+        gene: str,
+        coordinate_type: CoordinateType = CoordinateType.RESIDUE,
+    ) -> CdnaRepresentation | None:
         """Return MANE Transcript on the c. coordinate.
 
-        If an arg for ``gene`` is provided, lifts to GRCh38, then gets MANE cDNA
-        representation.
-
         >>> import asyncio
-        >>> from cool_seq_tool.app import CoolSeqTool
+        >>> from cool_seq_tool import CoolSeqTool
         >>> cst = CoolSeqTool()
         >>> result = asyncio.run(
         ...     cst.mane_transcript.g_to_mane_c(
@@ -1198,33 +1245,16 @@ class ManeTranscript:
         <TranscriptPriority.MANE_SELECT: 'mane_select'>
         >>> del cst
 
-        Locating a MANE transcript requires a ``gene`` symbol argument -- if none is
-        given, this method will only lift over to genomic coordinates on GRCh38.
-
         :param ac: Transcript accession on g. coordinate
         :param start_pos: genomic start position
         :param end_pos: genomic end position
         :param gene: HGNC gene symbol
-        :param residue_mode: Starting residue mode for ``start_pos`` and ``end_pos``.
-            Will always return coordinates in inter-residue.
-        :return: MANE Transcripts with cDNA change on c. coordinate if gene
-            is provided. Else, GRCh38 data
+        :param coordinate_type: Starting Coordinate type for ``start_pos`` and
+            ``end_pos``. Will always return inter-residue coordinates.
+        :return: MANE Transcripts with cDNA change on c. coordinate
         """
-        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, residue_mode)
-        residue_mode = ResidueMode.INTER_RESIDUE
-
-        # If gene not provided, return GRCh38
-        if not gene:
-            grch38 = await self.g_to_grch38(ac, start_pos, end_pos)
-            if not grch38:
-                return None
-
-            return GenomicRepresentation(
-                refseq=grch38["ac"],
-                pos=grch38["pos"],
-                status=TranscriptPriority.GRCH38,
-                alt_ac=grch38["ac"],
-            )
+        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, coordinate_type)
+        coordinate_type = CoordinateType.INTER_RESIDUE
 
         if not await self.uta_db.validate_genomic_ac(ac):
             _logger.warning("Genomic accession does not exist: %s", ac)
@@ -1238,12 +1268,18 @@ class ManeTranscript:
             mane_c_ac = current_mane_data["RefSeq_nuc"]
 
             # Liftover to GRCh38
-            grch38 = await self.g_to_grch38(ac, start_pos, end_pos)
+            grch38 = await self.g_to_grch38(
+                ac,
+                start_pos,
+                end_pos,
+                get_mane_genes=False,
+                coordinate_type=coordinate_type,
+            )
             mane_tx_genomic_data = None
             if grch38:
                 # GRCh38 -> MANE C
                 mane_tx_genomic_data = await self.uta_db.get_mane_c_genomic_data(
-                    mane_c_ac, grch38["ac"], grch38["pos"][0], grch38["pos"][1]
+                    mane_c_ac, grch38.ac, grch38.pos[0], grch38.pos[1]
                 )
 
             if not grch38 or not mane_tx_genomic_data:
@@ -1255,15 +1291,13 @@ class ManeTranscript:
                     continue
                 _logger.info("Not using most recent assembly")
 
-            coding_start_site = mane_tx_genomic_data["coding_start_site"]
-            coding_end_site = mane_tx_genomic_data["coding_end_site"]
+            coding_start_site = mane_tx_genomic_data.coding_start_site
+            coding_end_site = mane_tx_genomic_data.coding_end_site
             mane_c_pos_change = self.get_mane_c_pos_change(
                 mane_tx_genomic_data, coding_start_site
             )
 
-            if not self._validate_index(
-                mane_c_ac, mane_c_pos_change, coding_start_site
-            ):
+            if not self.validate_index(mane_c_ac, mane_c_pos_change, coding_start_site):
                 _logger.warning(
                     "%s are not valid positions on %s with coding start site %s",
                     mane_c_pos_change,
@@ -1284,7 +1318,7 @@ class ManeTranscript:
                 ),
                 refseq_c_ac=current_mane_data["RefSeq_nuc"],
                 ensembl_c_ac=current_mane_data["Ensembl_nuc"],
-                alt_ac=grch38["ac"] if grch38 else None,
+                alt_ac=grch38.ac if grch38 else None,
             )
         return None
 
@@ -1294,7 +1328,7 @@ class ManeTranscript:
         start_pos: int,
         end_pos: int,
         gene: str | None = None,
-        residue_mode: ResidueMode = ResidueMode.RESIDUE,
+        coordinate_type: CoordinateType = CoordinateType.RESIDUE,
         try_longest_compatible: bool = False,
     ) -> dict | None:
         """Given GRCh38 genomic representation, return protein representation.
@@ -1307,8 +1341,8 @@ class ManeTranscript:
         :param start_pos: Start position
         :param end_pos: End position
         :param gene: HGNC gene symbol
-        :param residue_mode: Starting residue mode for ``start_pos`` and ``end_pos``. Will
-            always return coordinates as inter-residue.
+        :param coordinate_type: Starting Coordinate type for ``start_pos`` and
+            ``end_pos``. Will always return inter-residue coordinates.
         :param try_longest_compatible: ``True`` if should try longest compatible remaining
             if mane transcript(s) not compatible. ``False`` otherwise.
         :return: If successful, return MANE data or longest compatible remaining (if
@@ -1327,8 +1361,8 @@ class ManeTranscript:
             return None
 
         # Step 2: Get inter-residue position
-        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, residue_mode)
-        residue_mode = ResidueMode.INTER_RESIDUE
+        start_pos, end_pos = get_inter_residue_pos(start_pos, end_pos, coordinate_type)
+        coordinate_type = CoordinateType.INTER_RESIDUE
 
         # Step 3: Try getting MANE protein representation
         mane_transcripts = set()  # Used if getting longest compatible remaining
@@ -1344,16 +1378,14 @@ class ManeTranscript:
                 continue
 
             # Get MANE C positions
-            coding_start_site = mane_tx_genomic_data["coding_start_site"]
-            coding_end_site = mane_tx_genomic_data["coding_end_site"]
+            coding_start_site = mane_tx_genomic_data.coding_start_site
+            coding_end_site = mane_tx_genomic_data.coding_end_site
             mane_c_pos_change = self.get_mane_c_pos_change(
                 mane_tx_genomic_data, coding_start_site
             )
 
             # Validate MANE C positions
-            if not self._validate_index(
-                mane_c_ac, mane_c_pos_change, coding_start_site
-            ):
+            if not self.validate_index(mane_c_ac, mane_c_pos_change, coding_start_site):
                 _logger.warning(
                     "%s are not valid positions on %s with coding start site %s",
                     mane_c_pos_change,
@@ -1367,7 +1399,7 @@ class ManeTranscript:
                 cdna=self._get_c_data(
                     (coding_start_site, coding_end_site),
                     mane_c_pos_change,
-                    mane_tx_genomic_data["strand"],
+                    mane_tx_genomic_data.strand,
                     TranscriptPriority(
                         "_".join(current_mane_data["MANE_status"].split()).lower()
                     ),
@@ -1383,7 +1415,7 @@ class ManeTranscript:
                 start_pos,
                 end_pos,
                 AnnotationLayer.GENOMIC,
-                residue_mode=residue_mode,
+                coordinate_type=coordinate_type,
                 alt_ac=alt_ac,
                 end_annotation_layer=EndAnnotationLayer.PROTEIN_AND_CDNA,
                 mane_transcripts=mane_transcripts,