International Business Machines Corporation

Filippo Utro

Laxmi Priya Parida

Kahn Rhrissorrakrai

Chaya Levovitz

Aldo Guzman Saenz

Ajay K Royyuru

Boaz Carmeli

Niina S Haiminen

Zeev Waks

Takahiko Koyama

Bilal Erhan

Enrico Siragusa

Myson Burch

Aritra Bose

A computer-implemented method includes receiving at a processor an input set of multiway data. The method determines a set of persistent homology barcodes based on the multiway data and identifies at least a first significant persistent homology barcode in the set of persistent homology barcodes. A representative cycle of the first significant persistent homology is returned and an orthonormal basis for a boundary of the multiway data is computed. The method obtains a harmonic representative by computing a projection of the representative cycle to an orthogonal complement and generates a set of feature vectors using the harmonic representative. Each feature vector has a magnitude corresponding to an impact that feature has on a given condition. The method then ranks the feature vectors by feature vector magnitude.

Ranking cross-disorder features in multiway data interactions

A computer-implemented method includes to determine a cell, tissue or a lesion representation in cell-free DNA comprises inputting, to a processor, cell-free DNA (cfDNA) genomic profiles from one or more fluid biopsy samples from a patient and one or more genomic profiles from one or more cells, tissues or lesions from the patient; constructing, by the processor, a plurality of synthetic fluid hypotheses (SFs); comparing, by the processor, each of the plurality of SFs to the cfDNA genomic profiles to determine goodness of fit, of each of the plurality of SFs; selecting, by the processor, a subset of the plurality of SFs, wherein each SF of the subset of SFs has a minimum distance in goodness of fit compared to the cfDNA genomic profile; and outputting, by the processor, based on the subset of SFs, a cell, tissue or a lesion representation in the cfDNA of the patient.

Determining cell, tissue, or lesion representations in cell-free DNA

Methods and systems for training a machine learning model are described. A processor can transform single cell data in a first space into projection data in a second space having a dimensionality lower than or equal to the first space. The processor can produce a cover having a plurality of sets of the projection data. The processor can determine a plurality of transition paths among the plurality of sets. A transition path can represent a transition from one cell state to another cell state. The processor can translate the transition paths from the second dimensional space to the first dimensional space. The processor can extract features from the transition paths in the first dimensional space. The processor can generate training data using the features, and use the training data to train a machine learning model for classifying cell state transitions.

Cell state transition features from single cell data

A computer-implemented method is disclosed which includes receiving biological sample information from one or more subjects at a first time period. The method further includes receiving biological sample information from the one or more subjects at a second time period. The method further includes comparing the biological sample information at the second time period with the biological sample information at the first time period. The method further includes generating a precedence graph based on results of the comparison. The method further includes determining one or more actions based on the precedence graph.

Methods and systems for determining drug resistance using a precedence graph

A computer-implemented method includes generating, by a processor, a set of training data for each phenotype in a database including a set of subjects. The set of training data is generated by dividing genomic information of N subjects selected with or without repetition into windows, computing a distribution of genomic events in the windows for each of N subjects, and extracting, for each window, a tensor that represents the distribution of genomic events for each of N subjects. A set of test data is generated for each phenotype in the database, a distribution of genomic events in windows for each phenotype is computed, and a tensor is extracted for each window that represents a distribution of genomic events for each phenotype. The method includes classifying each phenotype of the test data with a classifier, and assigning a phenotype to a patient.

Single sample genetic classification via tensor motifs

A computer-implemented method incudes calculating, by a processor, based on sequence data for a tumor from a subject at a plurality of time points, a mutation frequency for each of a plurality of SSVs at each of the time points to provide a plurality of time-resolved mutation frequencies (between 0 and 1) for each of the plurality of SSVs, the sequence data including a plurality of simple somatic variations (SSVs) at each of the time points; binning, by the processor, the plurality of time-resolved mutation frequencies for each SSV at each of the time points to provide a matrix of SSVs and time points; converting, by the processor, the matrix cells to pseudo-clones; and constructing, by the processor, a time-series tumor evolution tree from the pseudo-clones, wherein each time point in the time-series evolution tree represents an event in the subject's cancer treatment.

Time-series phylogenetic tumor evolution trees

A computer-implemented method includes determining, by a processor, from a time-series evolution tree comprising one or more clones at each of the plurality of time points, that the one or more clones are sensitive clones or resistant clones, wherein the time-series evolution tree is based on sequence data for a tumor from a subject at a plurality of time points, wherein each time point in the time-series evolution tree represents an event in the subject's cancer treatment, and wherein a clone is a collection of gene alterations; based at least in part on determining that the one or more clones that are the sensitive or resistant clones, determining, by the processor, a geneset composition of the one or more clones that are the sensitive or resistant clones; and based at least in part on determining the geneset composition, determining by the processor, a further treatment for the subject.

Functional analysis of time-series phylogenetic tumor evolution tree

A computer-implemented method, computer program product, and computer processing system are provided for metagenomic pattern classification. The method includes pre-processing, by a processor, a taxonomy tree associated with a genome database to extract taxonomy related information therefrom. The genome database includes a plurality of genome sequences. The method further includes building, by the processor, a suffix tree on the genome database. The method also includes annotating, by the processor, nodes in the suffix tree, using a plurality of right maximal patterns derived from the extracted taxonomy related information as annotations, such that each of the plurality of right maximal patterns in the suffix tree points to a respective one of a plurality of nodes in the taxonomy tree and such that a leaf node in the taxonomy tree represents a respective sample organism. The annotations are configured to function as classifications for the plurality of genome sequences.

Taxonomic annotation of variable length metagenomic patterns

A computer-implemented method includes inputting, to a processor, an N&#xd7;K SSV frequency matrix M and an error tolerance &#x3b4;&#x2265;0, wherein N is a number of SSVs and K is a number of time points, wherein matrix M comprises a plurality of time-resolved mutation frequencies for each SSV; clustering, by the processor, matrix rows in M that satisfy the &#x3b4; to provide a plurality of SSV clusters; assigning, by the processor, a mean cluster frequency to each SSV within each SSV cluster; calculating errors for removing low frequency rows, for rounding rows to 1 or 0; assigning a root node for all SSV clusters of frequency 1; and calculating, by the processor, a &#x3b4;-compliant time-series evolution tree with error &#x2264;&#x3b4; comprising the root node and a plurality time-stratified nodes, wherein calculating includes assigning a clonal configuration, optionally re-configuring the clonal configuration, and calculating error for re-configuring.

Phylogenetic tumor evolution trees with distribution of variants in cell populations

Embodiments of the present invention include method, systems and computer program products for algebraic phasing of polyploids. Aspects of the invention include receiving a matrix including a set of two or more single-nucleotide poloymorphisms (SNPs) for two or more sample organisms. Each row of the matrix is set to a ploidy based on a number of ploidies present in the two or more sample organisms. Each allele in the set of two or more SNPs is represented as a binary number. A set of algebraic rules is received, wherein the set of algebraic rules include an algebraic phasing algorithm. And the set of algebraic rules are applied to the matrix to determine a haplotype of a parent of the two or more sample organisms.

Algebraic phasing of polyploids

Embodiments are directed to computer implemented method of assessing a relevancy of a pathway to a disease of interest, the pathway having a source and a target. The method includes developing an impact of the source on the pathway. The method further includes developing a value of targeting, based at least in part on an alteration of the pathway, the pathway with a drug of interest. The method further includes identifying a relationship between the source and the target within the pathway. The method further includes combining: the impact of the source on the pathway; the value of targeting, based at least in part on the alteration of the pathway, the pathway with a drug of interest; and the relationship between the source and the target within the pathway, wherein the combining results in an assessment that represents the relevancy of the pathway to the disease of interest.

Simplified visualization and relevancy assessment of biological pathways

Various embodiments reconstruct haplotypes from genotype data. In one embodiment, a set of progeny genotype data comprising n progenies encoded with m genetic markers is accessed. A first set of parent haplotypes associated with a first parent of the n progenies and a second set of parent haplotypes associated with a second parent of the n progenies are identified based on at least the set of progeny genotype data. A total minimum number of observable crossovers in the n progenies is determined. An agglomerate data structure comprising a collection of sets of haplotype sequences characterizing the n progenies is constructed based on the set of progeny genotype data and the first and second sets of parent haplotypes. Each set of haplotype sequences includes a number of crossovers equal to the total minimum number of observable crossovers in the n progenies.

Exact haplotype reconstruction of F2 populations

A system for reconstructing haplotypes from genotype data includes a memory, a processor, and a reconstruction module. The reconstruction module is configured to access a set of progeny genotype data including n progenies encoded with m genetic markers. A first set of parent haplotypes associated with a first parent of the n progenies and a second set of parent haplotypes associated with a second parent of the n progenies are identified based on at least the set of progeny genotype data. An agglomerate data structure including a collection of sets of haplotype sequences characterizing the n progenies is constructed based on the set of progeny genotype data and the first and second sets of parent haplotypes. Each set of haplotype sequences includes a number of crossovers equal to a total minimum number of observable crossovers in the n progenies.

Growing community of inventors

Pleasantville, NY, United States of America

Filippo Utro