Method and apparatus for the access to bioinformatics data structured in access units

Inventors

Baluch, Mohamed KhosoAlberti, ClaudioZoia, GiorgioRenzi, Daniele

Assignees

Koninklijke Philips NV

Interested in licensing this patent?

MTEC can help explore whether this patent might be available for licensing for your application.

Publication Number

US-11763918-B2

Patent

Publication Date

2023-09-19

Expiration Date


Abstract

Method and apparatus for the coding and selective access of compressed genomic sequence data produced by genomic sequencing machines. The coding process is based on aligning sequence reads with respect to pre-existing or constructed reference sequences, on classifying and coding the sequence reads by means of sets of descriptors, and further partitioning the descriptor sets into access units of different types. Efficient selective access to specific genomic regions with the guarantee of retrieving all sequence reads mapped to those regions, is provided by: signaling the type of data mapping configuration used to store or transmit the descriptor sets, determining the minimum number of access units that need to be retrieved and decoded to access a genomic region, providing a master index table that contain all information for optimizing the data access process.

Core Innovation

The invention relates to selectively accessing compressed genomic data structured into entropy coded Access Units, where each AU includes genomic records associated with a genomic region. The system uses an indexing table disjointly coded from the Access Units, from which it extracts values of start coordinate s_i and end coordinate e_i for the Access Units, enabling selective access by identifying which AUs cover a queried genomic interval.

The invention additionally extracts a threshold value MaxD from the dataset header of the file, where MaxD specifies the maximum value that the pairing distance of read pairs can assume. This MaxD governs whether read pairs are coded as single entity or separately coded, and the same threshold value MaxD is parsed for decoding.

For efficient selective access and file organization, the invention uses a Master Index Table including AU start and end coordinates and byte offsets to descriptor blocks, along with descriptor stream sizes or AU sizes depending on the selected storage mode. The disclosed file layout supports selective access using modes in which blocks of data for different descriptors of the same AU are stored contiguously, or blocks belonging to the same descriptors across different AUs are stored in contiguous regions.

Claims Coverage

The independent claims cover three core aspects of the disclosed approach: selective access to entropy coded Access Units by parsing AU start and end coordinates from an indexing table, coding and decoding using a dataset-header threshold MaxD that controls coding of read pairs as single entities versus separately, and AU/range structure supported by mode-dependent organization. Across the independent claims, the inventive features focus on indexing-table-driven AU identification and MaxD-controlled handling of read pairs.

Indexing-table extraction of Access Unit coordinates for selective access

extracting values of start coordinate s_i and end coordinate e_i of said Access Units AU_i from an indexing table contained in the file and disjointly coded from said Access Units AU_i

Dataset-header threshold MaxD limiting pairing distance of read pairs

extracting a threshold value MaxD from the dataset header of said file, wherein said threshold value MaxD specifies the maximum value that the pairing distance of read pairs can assume, when read pairs are present in said file, such that said read pairs are coded as single entity or are separately coded.

Disjointly coded indexing table for coding Access Unit coordinate values

coding values of start coordinate s_i and end coordinate e_i of said Access Units AU_i in an indexing table disjointly coded from said Access Units for the case of storing information in a file

Dataset-header threshold MaxD for coding read pairs as single entity or separately

coding a threshold value MaxD in the genomic dataset header of said file, wherein said threshold value MaxD specifies the maximum value that the pairing distance of read pairs can assume, when read pairs are present in said genomic data to be coded, such that said read pairs are coded as single entity or are separately coded.

Indexing-table extraction of Access Unit coordinates for decoding

parsing values of start coordinate s_i and end coordinate e_i of said Access Units AU_i from an indexing table disjointly coded from said Access Units for the case of information stored in a file

Dataset-header threshold MaxD for decoding read pairs as single entity or separately

parsing a threshold value MaxD from the genomic dataset header of said file, wherein said threshold value MaxD specifies the maximum value that the pairing distance of read pairs can assume, when read pairs are present in said genomic encoded data, such that said read pairs are coded as single entity or are separately coded.

Across the independent claims, the inventive coverage centers on AU-level selective identification using an indexing table containing AU start and end coordinates disjointly coded from the entropy coded AUs, and on using a dataset-header threshold MaxD that specifies a maximum pairing distance for read pairs to determine whether read pairs are coded together as a single entity or separately. The decoding claim mirrors the coding claim by parsing the same AU coordinate values and MaxD from the genomic dataset header and disjoint indexing table.

Stated Advantages

Minimizing data volume and decoding/seek operations by identifying the minimum set of AUs to decode to retrieve all reads overlapping a queried genomic interval.

Selective access enabled by using an indexing table disjointly coded from entropy coded Access Units.

Decoder-side selection enabled by class-based mapping and pointer-based AU identification without entropy decoding the Access Units.

Documented Applications

Selectively accessing compressed genomic data structured into entropy coded Access Units for retrieval of reads overlapping a queried genomic interval.

Coding genomic data arranged into entropy coded Access Units using a dataset-header threshold MaxD to determine coding of read pairs as single entity or separately.

Decoding compressed coded genomic data arranged into entropy coded Access Units by parsing AU coordinate values and a dataset-header threshold MaxD to interpret read pairs as single entity or separately.

JOIN OUR MAILING LIST

Stay Connected with MTEC

Keep up with active and upcoming solicitations, MTEC news and other valuable information.