US20170260590A1 - Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results - Google Patents

Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results Download PDF

Info

Publication number
US20170260590A1
US20170260590A1 US15/431,395 US201715431395A US2017260590A1 US 20170260590 A1 US20170260590 A1 US 20170260590A1 US 201715431395 A US201715431395 A US 201715431395A US 2017260590 A1 US2017260590 A1 US 2017260590A1
Authority
US
United States
Prior art keywords
cancer
genetic
tumor
polynucleotides
subject
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
US15/431,395
Other languages
English (en)
Inventor
Helmy Eltoukhy
AmirAli Talasaz
Bahram Ghaffarzadeh Kermani
Nnamdi IHUEGBU
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guardant Health Inc
Original Assignee
Guardant Health Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=56284976&utm_source=google_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=US20170260590(A1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by Guardant Health Inc filed Critical Guardant Health Inc
Priority to US15/431,395 priority Critical patent/US20170260590A1/en
Assigned to GUARDANT HEALTH, INC. reassignment GUARDANT HEALTH, INC. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: KERMANI, BAHRAM GHAFFARZADEH, ELTOUKHY, HELMY, IHUEGBU, Nnamdi, TALASAZ, AmirAli
Publication of US20170260590A1 publication Critical patent/US20170260590A1/en
Priority to US17/000,010 priority patent/US20210040565A1/en
Priority to US17/462,906 priority patent/US20210395837A1/en
Priority to US17/699,968 priority patent/US20240141432A9/en
Pending legal-status Critical Current

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • C12Q1/6886Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material for cancer
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12NMICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N15/00Mutation or genetic engineering; DNA or RNA concerning genetic engineering, vectors, e.g. plasmids, or their isolation, preparation or purification; Use of hosts therefor
    • C12N15/09Recombinant DNA-technology
    • C12N15/11DNA or RNA fragments; Modified forms thereof; Non-coding nucleic acids having a biological activity
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6827Hybridisation assays for detection of mutation or polymorphism
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • G06F19/22
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B30/00ICT specially adapted for sequence analysis involving nucleotides or amino acids
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/118Prognosis of disease development
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/156Polymorphic or mutational markers

Definitions

  • Cancer refers to any disorder of various malignant neoplasms characterized by the proliferation of anaplastic cells that tend to invade surrounding tissue and metastasize to new body sites and the pathological conditions characterized by such growths.
  • cancers are constantly changing and mutating. Cancers can accumulate genetic variants through, e.g., somatic cell mutation. Such variants include, for example, sequence variants and copy number variants. Analysis of tumors has indicated that different cells in a tumor can bear different genetic variants. Such differentiation between tumor cells has been referred to as tumor heterogeneity.
  • Cancers can evolve over time, becoming resistant to a therapeutic intervention. Certain variants are known to correlate with responsiveness or resistance to specific therapeutic interventions. More effective treatments for cancers exhibiting tumor heterogeneity would be beneficial. Such cancers may be treated with a second, different, therapeutic intervention to which the cancer responds.
  • DNA sequencing methods allow detection of genetic variants in DNA from tumor cells. Cancer tumors continually shed their unique genomic material into the bloodstream. Unfortunately, these telltale genomic “signals” are so weak that current genomic analysis technologies, including next-generation sequencing, may only detect such signals sporadically or in patients with terminally high tumor burden. The main reason for this is that such technologies are plagued by error rates and bias that can be orders of magnitude higher than what is required to reliably detect de novo genomic alterations associated with cancer.
  • a method comprising: (a) sequencing polynucleotides from cancer cells from a biological sample of a subject; (b) identifying and quantifying somatic mutations in the polynucleotides; (c) developing a profile of tumor heterogeneity in the subject indicating the presence and relative quantity of a plurality of the somatic mutations in the polynucleotides, wherein different relative quantities indicates tumor heterogeneity; and (d) determining a therapeutic intervention for a cancer exhibiting the tumor heterogeneity, wherein the therapeutic intervention is effective against a cancer having the profile of tumor heterogeneity determined.
  • the cancer cells are spatially distinct.
  • the therapeutic intervention is more effective against a cancer presenting with the plurality of somatic mutations than it is against a cancer presenting with any one, but not all, of the somatic mutations.
  • the method further comprises: (e) monitoring changes in tumor heterogeneity in the subject over time and determining different therapeutic interventions over time based on the changes.
  • the method further comprises: (e) displaying the therapeutic intervention.
  • the method further comprises: (e) implementing the therapeutic intervention.
  • the method further comprises: (e) generating a phylogeny of tumor evolution based on the tumor profile; wherein determining the therapeutic intervention takes into account the phylogeny.
  • determining is performed with the aid of computer-executed algorithm.
  • sequence reads generated by sequencing are subject to noise reduction before identifying and quantifying.
  • noise reduction comprises molecular tracking of sequences generated from a single polynucleotide in the sample.
  • determining a therapeutic intervention takes into account the relative frequencies of the tumor-related genetic alterations.
  • the therapeutic intervention comprises administering, in combination or in series, a plurality of drugs, wherein each drug is relatively more effective against a cancer presenting with a different one of somatic mutations that occur at different relative frequency.
  • a drug that is relatively more effective against a cancer presenting with a somatic mutation occurring at higher relative frequency is administered in higher amount.
  • the drugs are delivered at doses that are stratified to reflect the relative amounts of the variants in the DNA.
  • cancers presenting with at least one of the genetic variants is resistant to at least one of the drugs.
  • determining a therapeutic intervention takes into account the tissue of origin of the cancer. In some embodiments, the therapeutic intervention is determined based on a database of interventions shown to be therapeutic for cancers having tumor heterogeneity characterized by each of the somatic mutations.
  • the polynucleotides comprise cfDNA from a blood sample. In some embodiments, the polynucleotides comprise polynucleotides from spatially distinct cancer cells. In some embodiments, the polynucleotides comprise polynucleotides from different metastatic tumor sites. In some embodiments, the polynucleotides comprise polynucleotides from a solid tumor or a diffuse tumor. In some embodiments, the polynucleotides are comprised in a blood sample or in solid tumor biopsy.
  • identifying comprises generating a plurality of sequence reads for parent polynucleotides from the sample, and collapsing the sequence reads to generate consensus calls for bases in each parent polynucleotide.
  • quantifying comprises determining frequency at which the somatic mutations are detected in the population of polynucleotides from the biological sample.
  • the biological sample comprises biological molecules from non-disease cells.
  • the biological sample comprises biological molecules from a plurality of different tissues.
  • the biomolecules are comprised in one biological sample.
  • the biomolecules are comprised in a plurality of biological samples.
  • the plurality of biological samples are tumors from a plurality of metastases.
  • sequencing comprises sequencing all or part of a subset of genes in the subject's genome.
  • the somatic mutations are selected from single nucleotide variations (SNVs), insertions, deletions, inversions, transversions, translocations, copy number variations (CNVs) (e.g., aneuploidy, partial aneuploidy, polyploidy), chromosomal instability, chromosomal structure alterations, gene fusions, chromosome fusions, gene truncations, gene amplification, gene duplications, chromosomal lesions, DNA lesions, abnormal changes in nucleic acid chemical modifications, abnormal changes in epigenetic patterns and abnormal changes in nucleic acid methylation.
  • genetic loci are selected from single nucleotides, genes and chromosomes.
  • the cancer is selected from carcinomas, sarcomas, leukemias, lymphomas, myelomas and central nervous system cancers (e.g., breast cancer, prostate cancer, colorectal cancer, brain cancer, esophageal cancer, head and neck cancer, bladder cancer, gynecological cancer, liposarcoma, and multiple myeloma).
  • cancer cells of the tumor are derived from a common parent disease cell.
  • cancer cells of the tumor are derived from different parent cancer cells of the same or different cancer type.
  • the method further comprises determining a measure of the somatic mutations to one or more control references to determine the relative quantity.
  • the polynucleotides are sourced from both circulating cancer polynucleotides and from solid tumor biopsy. In some embodiments, profiles are separately developed for polynucleotides sourced from the circulating cancer polynucleotides and from the solid tumor biopsy.
  • a method comprising providing a therapeutic intervention for a subject having a cancer having a tumor profile from which tumor heterogeneity can be inferred, wherein the therapeutic intervention is effective against cancers with the tumor profile.
  • the tumor profile indicates relative frequency of a plurality of more somatic mutations.
  • the method further comprises monitoring changes in the relative frequencies in the subject over time and determining different therapeutic interventions over time based on the changes.
  • the therapeutic intervention is more effective against a cancer presenting with each of the somatic mutations than it is against a cancer presenting with any one, but not all, of the somatic mutations.
  • the therapeutic intervention comprises administering, in combination or in series, a plurality of drugs, wherein each drug is relatively more effective against a cancer presenting with a different one of somatic mutations that occur at different relative frequency.
  • a drug that is relatively more effective against a cancer presenting with a somatic mutation occurring at higher relative frequency is administered in higher amount.
  • the drugs are delivered at doses that are stratified to reflect the relative amounts of the variants in the DNA.
  • cancers presenting with at least one of the genetic variants is resistant to at least one of the drugs.
  • the cancer is selected from carcinomas, sarcomas, leukemias, lymphomas, myelomas and central nervous system cancers (e.g., breast cancer, prostate cancer, colorectal cancer, brain cancer, esophageal cancer, head and neck cancer, bladder cancer, gynecological cancer, liposarcoma, and multiple myeloma).
  • carcinomas e.g., breast cancer, prostate cancer, colorectal cancer, brain cancer, esophageal cancer, head and neck cancer, bladder cancer, gynecological cancer, liposarcoma, and multiple myeloma.
  • a method comprising administering to a subject a therapeutic intervention that is effective against a tumor exhibiting tumor heterogeneity, wherein the therapeutic intervention is based on a profile of tumor heterogeneity in the subject indicating the presence and relative quantity of a plurality of the somatic mutations in the polynucleotides, wherein different relative quantities indicates tumor heterogeneity.
  • a system comprising a computer readable medium comprising machine-executable code that, upon execution by a computer processor, implements a method comprising: (a) receiving into memory sequence reads of polynucleotides mapping to a genetic locus; (b) determining, among said sequence reads, identity of bases that are different than a base of a reference sequence at the locus of the total number of sequence reads mapping to a locus; (c) reporting the identity and relative quantity of the determined bases and their location in the genome; and (d) inferring heterogeneity of a given sample based on information in (c).
  • the method implemented further comprises receiving into memory sequence reads derived from samples at a plurality of different times and calculating a difference in relative amount and identity of a plurality of bases between the two samples.
  • kits comprising a first pharmaceutical drug and a second pharmaceutical drug, wherein a combination of the first drug and the second drug is more therapeutically effective against a cancer presenting with a first and a second somatic mutation than it is against a cancer presenting with any one, but not all, of the somatic mutations.
  • the combination is contained in a mixture or each drug is contained in a separate container.
  • a method comprising: (a) performing biomolecular analysis of biomolecular polymers from disease cells (e.g., spatially distinct disease cells) from a subject; (b) identifying and quantifying biomolecular variants in the biomolecular macromolecules; (c) developing a profile of disease cell heterogeneity in the subject indicating the presence and relative quantity of a plurality of the variants in the biomolecular macromolecules, wherein different relative quantities indicates disease cell heterogeneity; and (d) determining a therapeutic intervention for a disease exhibiting the disease cell heterogeneity, wherein the therapeutic intervention is effective against a disease having the profile of disease cell heterogeneity determined.
  • the disease cells are spatially distinct disease cells.
  • the therapeutic intervention is determined based on a database of interventions shown to be therapeutic for cancers having tumor heterogeneity characterized by each of the somatic mutations.
  • a method of detecting disease cell heterogeneity in a subject comprising: a) quantifying polynucleotides that bear a sequence variant at each of a plurality of genetic loci in polynucleotides from a sample from the subject, wherein the sample comprises polynucleotides from somatic cells and from disease cells; b) determining for each locus a measure of copy number variation (CNV) for polynucleotides bearing the sequence variant; c) determining for each locus a weighted measure of quantity of polynucleotides bearing a sequence variant at the locus as a function of CNV at the locus; and d) comparing the weighted measures at each of the plurality of loci, wherein different weighted measures indicate disease cell heterogeneity.
  • the disease cells are tumor cells.
  • polynucleotides comprise cfDNA.
  • a method comprising: a) subjecting a subject to one or more pulsed therapy cycles, each pulsed therapy cycle comprising: (i) a first period during which one or more drugs is administered at a first amount and (ii) a second period during which the one or more drugs is administered at a second, reduced (e.g., completely not administered) amount; wherein: (A) the first period is characterized by a tumor burden detected above a first clinical level; and (B) the second period is characterized by a tumor burden detected below a second clinical level.
  • tumor burden is measured as a function of a quantity of a selected somatic variant in tumor polynucleotides.
  • one or more drugs is a plurality of drugs and each amount of each drug in each cycle is determined as a function of tumor burden measured as a function of a quantity of each of a plurality of different selected somatic variants in tumor polynucleotides.
  • the method comprises subjecting the subject to a plurality of pulsed therapy cycles.
  • the method further comprises: b) when the subject exhibits resistance to the one or more drugs, subjecting the subject to one or more pulsed therapy cycles, each pulsed therapy cycle comprising: (i) a first period during which a different one or more drugs is administered at a first amount and (ii) a second period during which the different one or more drugs is administered at a second, reduced (e.g., completely not administered) amount; wherein: (A) the first period is characterized by a tumor burden detected above a first clinical level; and (B) the second period is characterized by a tumor burden detected below a second clinical level.
  • a method comprising: (a) sequencing polynucleotides from cancer cells from a subject; (b) identifying and quantifying somatic mutations in the polynucleotides; and (c) developing a profile of tumor heterogeneity in the subject for use in determining a therapeutic intervention effective for a cancer exhibiting tumor heterogeneity, wherein the profile indicates the presence and relative quantity of a plurality of the somatic mutations in the polynucleotides, wherein different relative quantities indicates tumor heterogeneity.
  • a method comprising providing a therapeutic intervention for a subject wherein the therapeutic intervention is determined from a profile of disease cell heterogeneity in the subject, wherein the profile indicates the presence and relative quantity of a plurality of the somatic mutations in the polynucleotides, wherein different relative quantities indicates disease cell heterogeneity; and wherein the therapeutic intervention is effective against a disease having the profile of disease cell heterogeneity determined, e.g., more effective against a disease presenting with the plurality of somatic mutations than it is against a disease presenting with any one, but not all, of the somatic mutations.
  • a method comprising: a) determining a measure of deviation from a value of central tendency (e.g., standard deviation, variance) of copy number in polynucleotides in a sample across a region of at least 1 kb, at least 10 kb, at least 100 kb, at least 1 mb, at least 10 mb or at least 100 mb of a genome; b) inferring a measure of burden of DNA from cells undergoing cell division in the sample based on the measure of deviation.
  • the value of central tendency is mean, median or mode.
  • determining comprises partitioning the region into a plurality of non-overlapping intervals, determining a measure of copy number at each interval and determining the measure of deviation based on measures of copy number at each interval.
  • the interval is no more than any of 1 base, 10 bases, 100 bases, 1 kb bases or 10 kb.
  • a method of inferring a measure of burden of DNA from cells undergoing cell division in a sample comprising measuring copy number variation induced by proximity of one or more genomic loci to cells' origins of replication, wherein increased CNV indicates cells undergoing cell division.
  • the burden is measured in cell-free DNA.
  • the measure of burden relates to the fraction of tumor cells or genome-equivalents of DNA from tumor cells in the sample.
  • CNV due to proximity to origins of replication is inferred from a set of control samples or cell-lines.
  • a hidden-markov model, regression model, principal component analysis-based model, or genotype-modified model is used to approximate variations due to origins of replications.
  • the measure of burden is presence or absence of cells undergoing cell division. In some embodiments, proximity is within 1 kb of an origin of replication.
  • the method comprises measuring CNV at a locus, determining amount of CNV due to proximity of the locus to an origin of replication, and correcting the measured CNV to reflect genomic CNV, e.g., by subtracting amount of CNV attributable to cell division.
  • the genomic data is obtained from cell-free DNA.
  • the measure of burden relates to the fraction of tumor cells or genome-equivalents of DNA in a sample.
  • variations due to origins of replication are inferred from a set of control samples or cell-lines.
  • a hidden-markov model, regression model, principal component analysis-based model, or genotype-modified model is used to approximate variations due to origins of replications.
  • a method comprising: a) determining a baseline measure of copies of DNA molecules at one or more loci from one or more control samples, wherein one or more of the loci includes an origin of replication, each containing DNA from cells undergoing a predetermined level of cell division; b) determining a test measure of DNA molecules in a test sample; wherein the measure in test sample is from one or more loci partitioned into one or more partitions and wherein one or more of the loci includes an origin of replication; c) comparing the test measure and the baseline measure, wherein a test measure above a baseline measure indicates DNA in the test sample from cells dividing at a rate faster than cells providing DNA to the control sample.
  • the measure is selected from molecule count, a measure of central tendency of molecule count across partitions or a measure of variation of molecule count across partitions.
  • a method comprising: (a) administering to a subject an intervention that increases an amount of tumor-derived DNA in the subject's circulation; and (b) when said amount is increased, collecting from the subject a sample containing tumor-derived DNA.
  • the intervention preferentially kills tumor cells.
  • the intervention comprises exposing the subject or suspected diseased areas of the subject to radiation.
  • the intervention comprises exposing the subject or suspected diseased areas of subject to ultrasound.
  • the intervention comprises exposing the subject or suspected diseased areas of subject to physical agitation.
  • the intervention comprises administering to the subject a low dose of chemotherapy.
  • the method comprises administering the intervention to the subject within 1 week before collecting the sample.
  • the sample is selected from blood, plasma, serum, urine, saliva, cerebral spinal fluid, vaginal secretion, mucous and semen.
  • a method comprising compiling a database, wherein the database includes, for each of a plurality of subjects having cancer, tumor genomic testing data, including somatic alterations, collected at two or more time intervals per subject, one or more therapeutic interventions administered to each of the subjects at one or more times and efficacy of the therapeutic interventions, wherein the database is useful to infer efficacy of the therapeutic interventions in subjects with a tumor genomic profile.
  • the plurality is at least 50, at least 500 or at least 5000.
  • the tumor genomic testing data is collected via serial biopsy, cell-free DNA, cell-free RNA or circulating tumor cells.
  • relative frequencies of detected genetic variants are used to classify treatment efficacy.
  • additional information is used to help classify treatment efficacy, including but not limited to, weight, adverse treatment effects, histological testing, blood testing, radiographic information, prior treatments, and cancer type.
  • treatment response per patient is collected and classified quantitatively through additional testing.
  • the additional testing is blood or urine based testing.
  • a method comprising use of a database to identify one or more effective therapeutic interventions for a subject having cancer, wherein the database includes, for each of a plurality of subjects having cancer, tumor genomic testing data, including somatic alterations, collected at two or more time intervals per subject, one or more therapeutic interventions administered to each of the subjects at one or more times and efficacy of the therapeutic interventions.
  • identified therapeutic interventions are stratified by efficacy.
  • quantitative bounds on predicted therapeutic interventions efficacy or lack thereof are reported.
  • the therapeutic interventions use information of predicted tumor genomic evolution or acquired resistance mechanisms in similar patients in response to treatment.
  • the method comprises classifying effectiveness of treatment using a classification algorithm, e.g., linear regression processes (e.g., multiple linear regression (MLR), partial least squares (PLS) regression and principal components regression (PCR)), binary decision trees (e.g., recursive partitioning processes such as CART—classification and regression trees), artificial neural networks such as back propagation networks, discriminant analyses (e.g., Bayesian classifier or Fischer analysis), logistic classifiers, and support vector classifiers (e.g., support vector machines).
  • linear regression processes e.g., multiple linear regression (MLR), partial least squares (PLS) regression and principal components regression (PCR)
  • binary decision trees e.g., recursive partitioning processes such as CART—classification and regression trees
  • artificial neural networks such as back propagation networks
  • discriminant analyses e.g., Bayesian classifier or Fischer analysis
  • logistic classifiers e.g., support vector machines
  • a method to report results of one or more genetic tests comprising: capturing genetic information including genetic variants and quantitative measures thereof over one or more test points using a genetic analyzer; normalizing the quantitative measures for rendering with the one or more test points and generating a scaling factor; applying the scaling factor to render a tumor response map; and generating a summary of genetic variants.
  • the method comprises analyzing non-CNV (copy number variation) mutant allele frequencies.
  • the method comprises transforming an absolute value into a relative metric for rendering the tumor response map.
  • the method comprises multiplying a mutant allele frequency by a predetermined value and taking a log thereof.
  • the method comprises: multiplying the scaling factor by a transformed value for each gene to determine a quantity indicator to be rendered on the tumor response map; and assigning a unique visual indicator for each alteration in a visual panel.
  • the method comprises Y-centering or vertically centering the quantity indicator in a contiguously placed panel that indicates continuity.
  • the assigning further comprises providing a unique color for each alteration.
  • the method comprises analyzing genetic information from another test point or test time. In some embodiments, wherein a new test result does not differ from a prior test result, the method comprises rendering the prior visual panel. In some embodiments, wherein if alterations remain the same, but quantities have changed, the method comprises: maintaining the order and unique visual indicator for each alteration; and determining a new quantity indicator and generating a new visual panel for all test points. In some embodiments, the method comprises determining a new alteration in the genetic information and adding the alteration to the top of existing alterations. In some embodiments, the method comprises determining a new alteration in the genetic information and determining new transform values and scaling factor and assigning a unique visual indicator for each new alteration.
  • the method comprises determining a new alteration in the genetic information and re-generating the tumor response map including alterations from a prior test point that are still detected in current test point and the new alteration. In some embodiments, the method comprises determining if a prior alteration is no longer present and if so, comprising using a height of zero when rendering the quantity of the alteration of the prior alteration for subsequent test points. In some embodiments, the method comprises determining if a prior alteration is no longer present and if so, reserving the unique visual indicator associated with the prior alternation from future use.
  • the method comprises analyzing CNV mutant allele frequencies and methylation mutant allele frequencies. In some embodiments, the method comprises grouping of maximum mutant allele frequencies for rendering first on the tumor response map. In some embodiments, the method comprises rendering alterations for the gene in decreasing mutant allele frequency order of alterations. In some embodiments, the method comprises rendering alterations for the gene in a decreasing order. In some embodiments, the method comprises selecting a next gene with next highest mutant allele frequency.
  • the method comprises generating a trend indicator for the alteration over the different test points. In some embodiments, the method comprises generating a summary of alterations. In some embodiments, the method comprises generating a summary of treatment options. In some embodiments, the method comprises generating a summary of mutant allele frequency, cell free amplification, clinical approval indication, and clinical trial. In some embodiments, the method comprises generating a panel based on a biological pathway. In some embodiments, the method comprises generating a panel based on an evidence level. In some embodiments, the genetic information includes one or more of single-nucleotide variations, copy number variations, insertions and deletions, and gene rearrangements. In some embodiments, the method comprises generating a clinical relevance report on detected alterations. In some embodiments, the method comprises generating a therapy result summary.
  • a method to generate a genetic report comprising: generating non-copy number variation (CNV) data using a genetic analyzer; determining a scaling factor for each non-CNV mutant allele frequency; for a first test, generating a visual panel each non-CNV alteration using the scaling factor; and for each subsequent test, generating changes in the non-CNV alteration for the visual panel using the scaling factor.
  • CNV non-copy number variation
  • the method comprises transforming an absolute value into a relative metric for rendering. In some embodiments, the method comprises multiplying a mutant allele frequency by a predetermined value and taking a log of the predetermined value. In some embodiments, the method comprises determining a scaling factor using a maximum observed value. In some embodiments, for each non-CNV alteration, the method comprises multiplying a scaling factor by a transformed value for each gene variant as a quantity indicator for visualizing the gene variant.
  • the method comprises assigning a unique visual indicator for each alteration.
  • the method comprises using the visual panel if the test result is unchanged.
  • the method comprises maintaining the order and unique visual indicator for each alteration; and recomputing a quantity indicator for visualizing that variant and re-rendering updated values in existing panel(s) and new panel for the latest test.
  • the method comprises adding the alterations to the top of all existing alterations; computing transform values and the scaling factor; and assigning a unique visual indicator for each new alterations.
  • the method comprises: re-rendering alterations in the prior test point and the new alteration; and vertically centering an image of the alterations in a contiguously placed panel that indicates continuity.
  • the method comprises using a height of zero as the quantity of the alteration for a subsequent rendering.
  • the method comprises rendering subject or intervention information associated with alteration changes.
  • the method comprises identifying an alteration with the maximum Mutant Allele Frequency.
  • the method comprises: reporting alterations for that gene in decreasing mutant allele frequency order of non-CNV alterations; and reporting CNV alterations for that gene in decreasing order of CNV value. In some embodiments, the method comprises selecting the next gene with next highest non-CNV mutant allele frequency and reporting alterations for that gene in decreasing mutant allele frequency order of non-CNV alterations; and reporting CNV alterations for that gene in decreasing order of CNV value.
  • the method comprises rendering a trend indicator for an alteration over different test dates.
  • the method comprises grouping of maximum mutant allele frequencies and generating annotations including biological pathways or evidence level.
  • the method comprises generating a panel based on an evidence level.
  • the method comprises generating a panel based on a biological pathway.
  • the genetic information includes one or more of single-nucleotide variations, copy number variations, insertions and deletions, and gene rearrangements.
  • a method comprising: a) providing a plurality of nucleic acid samples from a subject, the samples collected at serial time points; b) sequencing polynucleotides from the samples to generate sequences; c) determining a quantitative measure of each of a plurality of genetic variants among the polynucleotides in each sample; d) graphically representing by computer relative quantities of genetic variants at each serial time point for those somatic mutations present at a non-zero quantity at least one of the serial time points.
  • the quantitative measure is the frequency of the genetic variant among all sequences mapping to the same genetic locus.
  • the relative quantities are represented as a stacked area graph.
  • the relative quantities are stacked, at the earliest time point, highest to lowest from the bottom to the top of the graph, and wherein a genetic variant first appearing at a non-zero quantity at a later time point is stacked at the top of the graph.
  • the areas are represented by different colors.
  • the graphical representation further indicates, for each time point, the quantitative measure of the predominant genetic variant.
  • the graphical representation further includes a key identifying genetic variants represented on the graph.
  • graphically representing comprises normalizing and scaling the quantitative measures.
  • the polynucleotides comprise cfDNA.
  • the loci are located in oncogenes.
  • the plurality of the genetic variants maps to a different gene in the genome.
  • the plurality of the genetic variants maps to the same gene in the genome. In some embodiments, at least 10 different oncogenes are sequenced.
  • determining comprises receiving the sequences into computer memory and using a computer processor to execute software to determine the quantitative measurement.
  • graphically representing comprises using a computer processor to execute software that transforms the quantitative measures into a graphical format and representing the graphical format on an electronic graphical user interface, e.g., a display screen.
  • a method to generate a paper or electronic patient test report from data generated by a genetic analyzer comprising: a) summarizing data from two or more testing time points, whereby a union of all non-zero testing results are reported at each subsequent test point after the first test; and b) rendering the testing results on the paper or electronic patient test report.
  • summarizing and rendering are performed on a computer by executing code with a computer processor to (i) identify all non-zero testing results, (ii) generate the test report and (iii) display the test report on a graphical user interface.
  • a method of graphically representing evolution of genetic variants of a tumor in a subject from data generated by a genetic analyzer comprising: a) generating by computer a stacked representation of genetic variants detected at each of a plurality of time points in the subject, wherein a height or width of each layer in the stack that corresponds to a genetic variant represents a quantitative contribution of the genetic variant to the a total quantity of genetic variants at each time point; and b) displaying the stacked representation on a computer monitor or a paper report.
  • the method further comprises using a combination of a magnitude of detected genetic variants in a body-fluid based test to infer a disease burden.
  • the method further comprises using allele fractions of detected mutations, allelic imbalances, gene-specific coverage to infer the disease burden.
  • an overall stack height is representative of overall disease burden or a disease burden score in the subject.
  • a distinct color is used to represent each genetic variant.
  • only a subset of detected genetic variants is plotted. In some embodiments, the subset is chosen based on likelihood of being a driver alteration or association with increased or reduced response to treatment.
  • the method comprises producing a test report for a genomic test.
  • a non-linear scale is used for representing the heights or widths of each represented genetic variant.
  • a plot of previous test points is depicted on the report.
  • the method comprises estimating a disease progression or remission based on rate of change and/or quantitative precision of each testing result.
  • the method comprises displaying a therapeutic intervention between intervening testing points.
  • displaying comprises: a) receiving data representing the detected tumor genetic variants into computer memory; b) executing code with a computer processor to graphically represent the quantitative contribution of each genetic variant at a time point as a line or area proportional to the relative contribution; and c) displaying the graphical representation on a graphical user interface.
  • FIG. 1 shows a flow chart of an exemplary method of determination and use of a therapeutic intervention.
  • FIG. 2 shows a flow chart of an exemplary method of determining frequency of variants in a sample corrected based on CNV at a locus.
  • FIG. 3 shows a flow chart of an exemplary method of providing pulsed therapy cycles which can delay drug resistance.
  • FIG. 4 shows a flow chart of an exemplary method of detecting tumor burden using CNV at origins of replication to detect DNA from dividing cells.
  • FIG. 5 shows an exemplary computer system.
  • FIG. 6 shows an exemplary scan of CNV across a region of a genome from samples containing cells in a resting state and in a state of cell division.
  • No genomic CNV is seen in loci a and b, but locus c shows gene duplication.
  • copy number is relatively equal in all intervals in the region, except those intervals overlapping the locus of gene duplication.
  • copy number appears to increase immediately after origins of replication, providing variance in CNV over the region. Deviation is particularly dramatic at a locus exhibiting CNV at an origin of replication (c).
  • FIG. 7 shows an exemplary course of monitoring and treatment of disease in a subject.
  • FIG. 8 shows an exemplary panel of 70 genes that exhibit genetic variation in cancer.
  • FIG. 9A shows an exemplary system for communicating cancer test results.
  • FIG. 9B shows an exemplary process to reduce error rates and bias in DNA sequence readings and generate genetic reports for users.
  • FIG. 10A-10C show exemplary processes for reporting genetic test results to users.
  • FIG. 10D-10I-2 show pages from an exemplary genetic test report.
  • FIG. 10J-10P shows various exemplary modified streamgraph.
  • FIG. 11A-11B shows exemplary processes for detecting mutation and reporting test results to users.
  • Methods of the present disclosure can detect biomolecular mosaicism (e.g., genetic mosaicism) in a biological sample, such as a heterogeneous genomic population of cells or deoxyribonucleic acid (DNA).
  • Genetic mosaicism can exist at the organismal level. For example, genetic variants that arise early in development can result in different somatic cells having different genomes.
  • An individual can be a chimera, e.g., produced by the fusion of two zygotes.
  • Organ transplant from an allogeneic donor can result in genetic mosaics, which also can be detected by examining polynucleotides shed into the blood from the transplanted organ.
  • Disease cell heterogeneity in which diseased cells have different genetic variants, is another form of genetic mosaicism.
  • Methods provided herein can detect mosaicism and, in the case of disease, provide therapeutic intervention.
  • this disclosure provides methods for performing body-wide profiling of biomolecular mosaicism through the use of circulating polynucleotides, which may derive or otherwise originate from cells in diverse locations of the body of a subject.
  • Diseased cells such as tumors
  • Diseased cells may evolve over time, resulting in different clonal sub-populations having new genetic and phenotypic characteristics. This may result from natural mutations as the cells divide, or it may be driven by treatments that target certain clonal sub-populations, allowing clones more resistant to the treatment to proliferate by negative selection.
  • the existence of sub-populations of diseased cells that bear different genotypic or phenotypic characteristics is referred to herein as disease cell heterogeneity, or, in the case of cancer, tumor heterogeneity.
  • cancers are treated based on mutant forms found in a cancer biopsy.
  • the finding of Her2+ in even small amounts of breast cancer cells may be indicative of breast cancer, which may be followed through with a treatment using an anti-Her2 30 therapy.
  • a colorectal cancer in which a KRAS mutant is found in small amounts may be treated with a therapy for which KRAS is responsive.
  • This information can be used by a health care provider, e.g., a physician, to develop therapeutic interventions.
  • a subject that has a heterogeneous tumor can be treated as if they had two tumors, and a therapeutic intervention can treat each of the tumors.
  • the therapeutic intervention could include, for example, a combination therapy including a first drug effective against the first tumor type and a second drug effective against the second tumor type.
  • the drugs can be given in amounts that reflect the relative amounts of the mutant forms detected. For example, a drug to treat the mutant form that is found in higher relative amounts can be delivered at greater dose than a drug to treat the mutant form in lesser relative amount. Or, treatment for the mutant in the lesser relative amount can be delayed or staggered with respect the mutant in greater amount.
  • therapeutic intervention can be calibrated to an evolving tumor. For example, analysis may show increasing amounts of polynucleotides bearing drug resistance mutants.
  • the therapeutic intervention can be modified to decrease the amount of drug effective to treat a tumor that does not bear the resistance mutant and increase administration of a drug that does treat a tumor bearing the resistance marker.
  • Therapeutic interventions can be determined by a healthcare provider or by a computer algorithm, or a combination of the two.
  • a database can contain the results of therapeutic interventions against diseases having various profiles of disease cell heterogeneity. The database can be consulted in determining a therapeutic intervention for a disease with a particular profile.
  • This present disclosure provides, among other things, methods of determining a therapeutic intervention for a subject having a disease, such as cancer, that exhibits disease cell heterogeneity, e.g., tumor heterogeneity.
  • the method involves analyzing biological macromolecules (e.g., sequencing polynucleotides) of disease cells (e.g., spatially distinct disease cells) from a subject having the disease.
  • a profile of disease cell heterogeneity is developed that indicates the existence of genetic variants specific to the disease cells and the amount of these variants relative to each other. This information, in turn, is used to determine a therapeutic intervention that takes the profile into account.
  • a subject of the methods of this disclosure is any multicellular organism. More specifically, the subject can be a plant or an animal, a vertebrate, a mammal, a mouse, a primate, a simian or a human. Animals include, but are not limited to, farm animals, sport animals, and pets.
  • a subject can be a healthy individual, an individual that has or is suspected of having a disease or a pre-disposition to the disease, or an individual that is in need of therapy or suspected of needing therapy.
  • a subject can be a patient, e.g., a subject under the care of a professional heathcare provider.
  • the subject can have a pathological condition (disease).
  • Disease cells exhibiting pathology of disease are referred to herein as disease cells.
  • the disease can be a cancer.
  • Cancer is a condition characterized by abnormal cells that divide out of control. Cancers include, without limitation, carcinomas, sarcomas, leukemias, lymphomas, myelomas and central nervous system cancers. More specific examples of cancers are breast cancer, prostate cancer, colorectal cancer, brain cancer, esophageal cancer, head and neck cancer, bladder cancer, gynecological cancer, liposarcoma, and multiple myeloma.
  • cancers include, for example, acute lymphoblastic leukemia (ALL), acute myeloid leukemia (AML), adrenocortical carcinoma, Kaposi Sarcoma, anal cancer, basal cell carcinoma, bile duct cancer, bladder cancer, bone cancer, osteosarcoma, malignant fibrous histiocytoma, brain stem glioma, brain cancer, craniopharyngioma, ependymoblastoma, ependymoma, medulloblastoma, medulloeptithelioma, pineal parenchymal tumor, breast cancer, bronchial tumor, Burkitt lymphoma, Non-Hodgkin lymphoma, carcinoid tumor, cervical cancer, chordoma, chronic lymphocytic leukemia (CLL), chronic myelogenous leukemia (CML), colon cancer, colorectal cancer, cutaneous T-cell lymphoma, ductal carcinoma in situ, endometrial cancer
  • a tumor is a collection of cancer cells (cancer disease cells). This includes, for example, a collection of cells in a single mass of cells (e.g., a solid tumor), a collection of cells from different metastatic tumor sites (metastatic tumors), and diffuse tumors (e.g., circulating tumor cells).
  • a tumor can include cells of a single cancer (e.g., colorectal cancer), or multiple cancers (e.g., colorectal cancer and pancreatic cancer).
  • a tumor can include cells originating from a single original somatic cell or from different somatic cells.
  • disease cells in the subject are spatially distinct.
  • Disease cells are spatially distinct if the cells are located at least 1 cm, at least 2 cm, at least 5 cm or at least 10 cm apart in a body, e.g, in different tissues or organs, or the same tissue or organ.
  • examples of spatially distinct cancer cells include cancer cells from diffuse cancers (such as leukemias), cancer cells at different metastatic sites, and cancer cells from the same mass of tumor cells that are separated by at least 1 cm.
  • Disease cell burden is a quantitative measure of the amount of disease cells in a subject.
  • One measure of disease cell burden is the fraction of total biological macromolecules in a sample that are disease biological macromolecules, e.g., the relative amount of tumor polynucleotides in a sample of cell free polynucleotides. For example, if cfDNA from a first subject has 10% cancer polynucleotides, the subject may be said to have a cell-free tumor burden of 10%, If cfDNA from a second subject has 5% cancer polynucleotides, the a second subject may be said to have half the cell-free tumor burden of the first subject.
  • Polynucleotides to be sequenced can be sourced from spatially distinct sites. This includes polynucleotides sourced from biopsies of different locations in a single tumor mass. It also includes polynucleotides sourced from cells at different metastatic tumor sites. Cells shed polynucleotides into the blood where it is detectable as cell free polynucleotides (e.g., circulating tumor DNA). Cell free polynucleotides also can be found in other bodily fluids such as urine. Therefore, cfDNA provides a more accurate profile of tumor heterogeneity across the entire disease cell population than DNA sourced from a single tumor location. DNA sampled from cells across the disease cell population in a body is referred to as “disease burden DNA” or, in the case of cancer, “tumor burden DNA”.
  • Disease cells such as tumors
  • tumors may share one, two, three or more genetic variants.
  • Such variants may share the same stratification, for example highest frequency, second highest frequency, etc.
  • Profiles can also share similar disease cell burdens, e.g., cfDNA burdens, e.g., within 15%, within 10%, within 5% or within 2%.
  • a macromolecule is a molecule formed from monomeric subunits.
  • Monomeric subunits forming biological macromolecules include, for example, nucleotides, amino acids, monosaccharides and fatty acids.
  • Biological macromolecules include, for example, biopolymers and non-polymeric macromolecules.
  • a polynucleotide is a macromolecule comprising a polymer of nucleotides.
  • Polynucleotides include, for example, polydeoxyribonucleotides (DNA) and polyribonucleotides (RNA).
  • a polypeptide is a macromolecule comprising a polymer of amino acids.
  • a polysaccharide is a macromolecule comprising a polymer of monosaccharides.
  • Lipids are a diverse group of organic compounds including, for example, fats, oils and hormones that share the functional characteristic of not interacting appreciably with water.
  • a triglyceride is a fat formed from three fatty acid chains.
  • a cancer polynucleotide is a polynucleotide (e.g., DNA) derived from a cancer cell.
  • Cancer DNA and/or RNA can be extracted from tumors, from isolated cancer cells or from biological fluids (e.g., saliva, serum, blood or urine) in the form of cell free DNA (cfDNA) or cell free RNA.
  • Cell free DNA is DNA located outside of a cell in a bodily fluid, e.g., in blood or urine.
  • Circulating nucleic acids are nucleic acids found in the blood stream.
  • Cell free DNA in the blood is a form of circulating nucleic acid.
  • Cell free DNA is believed to arise from dying cells that shed their DNA into the blood. Because spatially distinct cancer cells will shed DNA into bodily fluids, such as blood, cfDNA of cancer subjects typically comprises cancer DNA from spatially distinct cancer cells.
  • Analytes for analysis in the methods of this disclosure can derive from a biological sample, e.g., a sample comprising a biological macromolecule.
  • a biological sample can be derived from any organ, tissue or biological fluid.
  • a biological sample can comprise, for example, a bodily fluid or a solid tissue sample.
  • An example of a solid tissue sample is a tumor sample, e.g., from a solid tumor biopsy.
  • Bodily fluids include, for example, blood, serum, tumor cells, saliva, urine, lymphatic fluid, prostatic fluid, seminal fluid, milk, sputum, stool and tears.
  • Bodily fluids are particularly good sources of biological macromolecules from spatially distinct disease cells, as such cells from many locations in a body can shed these molecules into the bodily fluid.
  • blood and urine are good sources of cell free polynucleotides. Macromolecules from such sources can provide a more accurate profile of the diseased cells than macromolecules derived from a localized disease cell mass.
  • Amounts of disease polynucleotides in a bodily fluid sample can be increased. Such increases can increase sensitivity of detection of disease polynucleotides.
  • an intervention such as a therapeutic intervention, is administered to a subject that causes disease cells to lyse, emptying their DNA into the surrounding fluid.
  • Such interventions can include administration of chemotherapy. It also can include administering radiation or ultrasound to the whole body of a subject, or to a portion of the body of a subject, such as being directed to a tumor or a diseased organ.
  • a fluid sample is collected for analysis. The interval between administration of the intervention and collection can be long enough for the disease polynucleotides to increase, but not so long that they are cleared from the body. For example, a low dose of chemotherapy can be administered about a week before collection of the sample.
  • Genomic analysis can be performed by, for example, a genetic analyzer, e.g., using DNA sequencing. Methylation analysis can be performed by, for example, conversion of methylated bases followed by DNA sequencing. RNA expression analysis can be performed by, for example, polynucleotide array hybridization. Proteomic analysis can be performed by, for example, mass spectrometry.
  • the term “genetic analyzer” refers to a system including a DNA sequencer for generating DNA sequence information and a computer comprising software that performs bioinformatic analysis on the DNA sequence information.
  • Bioinformatic analysis can include, without limitation, assembling sequence data, detecting and quantifying genetic variants in a sample, including either of germline variants (e.g., heterozygosity) and somatic cell variants (e.g., cancer cell variants).
  • Analytic methods can include generating and capturing genetic information.
  • Genetic information can include genetic sequence information, ploidy states, the identity of one or more genetic variants, as well as a quantitative measure of the variants.
  • quantitative measure refers to any measure of quantity including absolute and relative measures.
  • a quantitative measure can be, for example, a number (e.g., a count), a percentage, a frequency, a degree or a threshold amount.
  • Polynucleotides can be analyzed by any method known in the art.
  • the DNA sequencer will employ next generation sequencing (e.g., Illumina, 454, Ion torrent, SOLiD).
  • Sequence analysis can be performed by massively parallel sequencing, that is, simultaneously (or in rapid succession) sequencing any of at least 100,000, 1 million, 10 million, 100 million, or 1 billion polynucleotide molecules.
  • Sequencing methods may include, but are not limited to: high-throughput sequencing, pyrosequencing, sequencing-by-synthesis, single-molecule sequencing, nanopore sequencing, semiconductor sequencing, sequencing-by-ligation, sequencing-by-hybridization, RNA-Seq (Illumina), Digital Gene Expression (Helicos), Next generation sequencing, Single Molecule Sequencing by Synthesis (SMSS) (Helicos), massively-parallel sequencing, Clonal Single Molecule Array (Solexa), shotgun sequencing, Maxam-Gilbert or Sanger sequencing, primer walking, sequencing using PacBio, SOLiD, Ion Torrent, Genius (GenapSys) or Nanopore (e.g., Oxford Nanopore) platforms and any other sequencing methods known in the art.
  • SMSS Single Molecule Sequencing by Synthesis
  • Solexa Solexa
  • shotgun sequencing Maxam-Gilbert or Sanger sequencing
  • primer walking sequencing using PacBio, SOLiD, Ion
  • the DNA sequencer can apply Gilbert's sequencing method based on chemical modification of DNA followed by cleavage at specific bases, or it can apply Sanger's technique which is based on dideoxynucleotide chain termination.
  • the Sanger method became popular due to its increased efficiency and low radioactivity.
  • the DNA sequencer can use techniques that do not require DNA amplification (polymerase chain reaction—PCR), which speeds up the sample preparation before sequencing and reduces errors.
  • PCR polymerase chain reaction
  • sequencing data is collected from the reactions caused by the addition of nucleotides in the complementary strand in real time.
  • the DNA sequencers can utilize a method called Single-molecule real-time (SMRT), where sequencing data is produced by light (captured by a camera) emitted when a nucleotide is added to the complementary strand by enzymes containing fluorescent dyes.
  • SMRT Single-molecule real-time
  • Sequencing of the genome can be selective, e.g., directed to portions of the genome of interest. For example, many genes (and mutant forms of these genes) are known to be associated with various cancers. Sequencing of select genes, or portions of genes may suffice for the analysis desired. Polynucleotides mapping to specific loci in the genome that are the subject of interest can be isolated for sequencing by, for example, sequence capture or site-specific amplification.
  • a nucleotide sequence (e.g., DNA sequence) can refer to raw sequence reads or processed sequence reads, such as unique molecular counts inferred from raw sequence reads.
  • Sequence reads generated from sequencing are subject to analysis including, for example, identifying genetic variants. This can include identifying sequence variants and quantifying numbers of base calls at each locus. Quantifying can involve, for example, counting the number of reads mapping to a particular genetic locus. Different numbers of reads at different loci can indicate copy number variation (CNV).
  • CNV copy number variation
  • Sequencing and bioinformatics methods that reduce noise and distortion are particularly useful when the number of target polynucleotides in a sample is small compared with non-target polynucleotides.
  • the signal from the target may be weak. This can be the case, for example, in the case of cell free DNA, where a small number of tumor polynucleotides may be mixed with a much larger number of polynucleotides from healthy cells.
  • Molecular tracking methods can be useful in such situations. Molecular tracking involves tracking sequence reads from a sequencing protocol back to molecules in an original sample (e.g., before amplification and/or sequencing) from which the reads are derived.
  • Certain methods involve tagging molecules in such a way that multiple sequence reads produced from original molecules can be grouped into families of sequences derived from original molecules. In this way, base calls representing noise can be filtered out.
  • Such methods are described in more detail in, for example, WO 2013/142389 (Schmitt et al.), US 2014/0227705 (Vogelstein et al.) and WO 2014/149134 (Talasaz et al.). Up-sampling methods also are useful to more accurately determine counts of molecules in a sample.
  • up-sampling methods involve determining a quantitative measure of individual DNA molecules for which both strands (Watson and Crick strands) are detected; determining a quantitative measure of individual DNA molecules for which only one of the DNA strands is detected; inferring from these measures a quantitative measure of individual DNA molecules for which neither strand was detected; and using these measures to determine the quantitative measure indicative of a number of individual double-stranded DNA molecules in the sample.
  • This method is described in more detail in PCT/US2014/072383, filed Dec. 24, 2014.
  • Methods of the present disclosure can be used in the detection of genetic variants (also referred to a “gene alterations”).
  • Genetic variants are alternative forms at a genetic locus. In the human genome, approximately 0.1% of nucleotide positions are polymorphic, that is, exist in a second genetic form occurring in at least 1% of the population. Mutations can introduce genetic variants into the germ line, and also into disease cells, such as cancer.
  • Reference sequences such as hg19 or NCBI Build 37 or Build 38, intend to represent a “wild type” or “normal” genome. However, to the extent they have a single sequence, they do not identify common polymorphisms which may also be considered normal.
  • Genetic variants include sequence variants, copy number variants and nucleotide modification variants.
  • a sequence variant is a variation in a genetic nucleotide sequence.
  • a copy number variant is a deviation from wild type in the number of copies of a portion of a genome.
  • Genetic variants include, for example, single nucleotide variations (SNPs), insertions, deletions, inversions, transversions, translocations, gene fusions, chromosome fusions, gene truncations, copy number variations (e.g., aneuploidy, partial aneuploidy, polyploidy, gene amplification), abnormal changes in nucleic acid chemical modifications, abnormal changes in epigenetic patterns and abnormal changes in nucleic acid methylation.
  • SNPs single nucleotide variations
  • insertions insertions, deletions, inversions, transversions, translocations
  • gene fusions chromosome fusions
  • gene truncations copy number variations (e.g
  • Genetic variants can be detected by comparing sequences from polynucleotides in a sample to a reference, e.g., to a reference genome sequence, to an index or to a database of known mutations.
  • the reference sequence is a publicly available reference sequence, such as the human genome sequence HG-19 or NCBI Build 37.
  • the reference sequence is a sequence in a non-public database.
  • the reference sequence is a germ line sequence of an organism inferred or determined from sequencing polynucleotides from the organism.
  • a somatic mutation or somatic alteration is a genetic variant that arises in a somatic cell.
  • Somatic mutations are distinguished from mutations that arise in the genome of a germ line cell (i.e., sperm or egg) or a zygote, of an individual.
  • Somatic mutations e.g., those found in cancer cells, are distinguishable from the germ line genome of a subject in which the cancer arose. They also can be detected by comparing the cancer genome with the germ line genome or with a reference genome.
  • a database of SNVs in human cancer can be found at the website: cancer.sanger.ac.uk/cancergenome/projects/cosmic/.
  • FIG. 8 shows genes known, in cancer, to exhibit point mutations, amplifications, fusions and indels.
  • a diploid cell having 2N chromosomes with replicated DNA may correspond to about 4 ⁇ DNA content, whereas a diploid cell having 2N chromosomes without replicated DNA may correspond to about 2 ⁇ DNA content.
  • Replication proceeds from origins of replication. In mammals, origins of replication are spaced at intervals of about 15 kb to 300 kb. During this period, portions of the genome exist in polyploid form. Those areas between origins of replication and the position of the polymerase are duplicated, while those areas beyond the position of the polymerase (or just before the origin of replication) are still in single copy number in the strand undergoing replication.
  • One method to detect tumor burden involves determining copy-number variation due to proximity of examined locus or loci to various origins of replication. Regions that include a replication origin will have very close to 4 copies of DNA in that locus (in a diploid cell), while regions that are far removed from a replication origin will have closer to 2 copies (in a diploid cell).
  • the examined locus or loci include, at least 1kb, at least 10 kb, at least 100 kb, at least 1 mb, at least 10 mb, at least 100 mb, across an entire chromosome or across an entire genome.
  • a measure of replication origin CNV (ROCNV) across the region is determined. This can be, for example, a measure of deviation in copy number from a value of central tendency.
  • the value of central tendency can be, for example, mean, median or mode.
  • the measure of deviation can be for example, variance or standard deviation. This measure can be compared with a measure of ROCNVs across the same region in a control sample, e.g., from a healthy individual or cells in resting state.
  • ROCNVs can be determined by partitioning the region or regions analyzed into non-overlapping partitions of various lengths and taking a measure of CNV in this partition. This measure of CNV can be derived from the number of reads or fragments determined to map to those regions after sequencing.
  • the partitions can have various sizes, to produce various levels of resolution, e.g., a single base level (base-per-base), 10 bases, 100 bases, 1 kb, 10 kb or 100 kb.
  • Deviations that are greater than a control indicate the presence of DNA undergoing replication, which, in turn, indicates malignancy. The greater the degree of deviation, the greater the amount of DNA from cells undergoing cell division in the sample.
  • Various methods can be used to calculate true genetic copy number variations that differ from replication origin based distortion.
  • heterozygous SNP positions at affected CNV loci can be used to infer copy number variation by calculating the deviation from 50% or the allelic imbalance at those loci. Distortion due to replication origin proximity should not affect this imbalance since both copies would generally be copied at similar time intervals and thus self-normalizing (although allelic changes could conceivably change the replication of origin between the two allelic variants).
  • duplication of a chromosome segment containing a SNP could be detected in around 67% of reads, while duplication resulting from ROCNV would be detected in about 50% of reads.
  • counting-based techniques that use the density of detected fragments or reads at a certain locus are used to calculate relative copy number. These techniques are generally limited by poisson noise and systematic bias due to DNA sample preparation and sequencing bias. A combination of these methods may also be to obtain even greater accuracy.
  • ROCNV can be calculated for a given sample and be used to give a value on cell-free tumor burden despite lack of detection of traditional somatic variants, such as, SNVs, gene-specific CNVs, genomic rearrangements, epigenetic variants, loss of heterozygosity, etc.
  • ROCNVs can also be used to subtract distortion for a given sample to increase sensitivity and/or specificity of a given CNV detection/estimation method by removing variation that is related to replication origin proximity rather than due to true copy number changes in a cell.
  • Cell-lines with known or no copy number changes over a reference can also be used as a reference of ROCNVs for use in estimating its contribution to a given sample.
  • the method involves determining a baseline level of copies of DNA molecules at one or more loci from one or more control samples, each containing DNA from cells undergoing a predetermined level of cell division, e.g., cells in resting state or rapidly dividing tumor cells.
  • a measure of copies of DNA molecules in a test sample is also determined.
  • the measure in test samples can be from one or more loci partitioned into one or more partitions. In each case, a plurality of loci each include an origin or replication.
  • the measure of copies from the test sample can be an average across all partitions, or a level of variance across loci.
  • a measure of central tendency or of variation e.g., variance or standard deviation
  • a measure that is greater in a test sample than in a control of cells in resting state, or slowly dividing indicates that cells generating the DNA in the test sample are dividing more rapidly than cells providing DNA to the control sample, e.g., are cancerous.
  • measures that are similar between a test sample and a control of cells in actively dividing state indicates that cells generating the DNA in the test sample are dividing at a rate similar to the rapidly dividing cells, e.g., are cancerous.
  • Disease cell heterogeneity e.g., tumor heterogeneity
  • Disease cell heterogeneity is the occurrence of diseased cells having different genetic variants.
  • Disease cell heterogeneity can be determined by examination of polynucleotides isolated from diseased cells and detection of differences in their genomes.
  • Disease cell heterogeneity also can be inferred from examination of polynucleotides from a sample containing polynucleotides from both diseased and healthy cells based on differences in relative frequency of somatic mutations.
  • cancer is characterized by changes at the genetic level, e.g., through the accumulation of somatic mutations in different clonal groups of cells. These changes can contribute to unregulated growth of the cancer cells, or function as markers of responsiveness or non-responsiveness to various therapeutic interventions.
  • Tumor heterogeneity is a condition in which a tumor characterized by cancer cells containing different combinations of genetic variants, e.g., different combinations of somatic mutations. That is, the tumor can have different cells containing alterations in different genes, or containing different alterations in the same gene.
  • a first cell could include a mutant form of BRAF
  • a second cell could include mutant forms of both BRAF and ERBB2.
  • a first cancer cell could include the single nucleotide polymorphism EGRF 55249063 G>A
  • a second cell could include the single nucleotide polymorphism EGRF 55238874 T>A. (Numbers refer to nucleotide position in genomic reference sequence.)
  • an original tumor cell can include a genetic variant in a gene, e.g., an oncogene.
  • a genetic variant in a gene e.g., an oncogene.
  • some progeny cells, which carry the original mutation may independently develop genetic variants in other genes or in different parts of the same gene. In subsequent divisions, tumor cells can accumulate still more genetic variants.
  • the profile includes information from polynucleotides from spatially distinct disease cells.
  • the profile is a whole body profile containing information from cells distributed throughout the body. Analysis of polynucleotides in cfDNA allows sampling of DNA across the entire geographic extent of a tumor, in contrast with sampling of a localized area of a tumor. In particular, it allows sampling of diffuse and metastatic tumors. This contrasts with methods that detect the mere existence of tumor heterogeneity through the localized sampling of a tumor.
  • the profile can indicate the exact nucleotide sequence of the variant, or may simply indicate a gene bearing the somatic mutation.
  • a profile of disease cell heterogeneity such as tumor cell heterogeneity
  • the profile identifies genetic variations and the relative amounts of each variant. From this information, one can infer possible distributions of the variants in different cell sub-population. For example, a cancer may begin with a cell bearing somatic mutation X. As a result of clonal evolution, some progeny of this cell may develop variant Y. Other progeny may develop variant Z.
  • the tumor may be characterized as 50% X, 35% XY and 15% XZ.
  • the profile may indicate 100% X, 35% Y and 15% Z.
  • Tumor heterogeneity can be detected from analysis of sequences of cancer polynucleotides, based on the existence of genomic variations at different loci occurring at different frequencies. For example, in a sample of cell free DNA (which is likely to contain germ line DNA as well as cancer DNA), it may be found that a sequence variant of BRAF occurs at a frequency of 17%, a sequence variant of CDKN2A occurs at a frequency of 6%, a sequence variant of ERBB2 occurs at a frequency of 3% and a sequence variant of ATM occurs at a frequency of 1%. These different frequencies of sequence variants indicate tumor heterogeneity. Similarly, genetic sequences exhibiting different amounts of copy number variation also indicate tumor heterogeneity. For example, analysis of a sample may show different levels of amplification for the EGFR and CCNE1 genes. This also indicates tumor heterogeneity.
  • somatic mutations can be made by comparing base calls in the sample to a reference sequence or, internally, as less frequent base calls to more common base calls, presumed to be in the germ line sequence. In either case, the existence of sub-dominant forms (e.g., less than 40% of total base calls) at different loci and at different frequency indicates disease cell heterogeneity.
  • Cell free DNA typically comprises a preponderance of DNA from normal cells having the germ line genome sequence and, in the case of a disease, such as cancer, a small percentage of DNA from cancer cells and having a cancer genome sequence.
  • Sequences generated from polynucleotides in a sample of cfDNA can be compared with a reference sequence to detect differences between the reference sequence and the polynucleotides in the cfDNA.
  • all or nearly all of the polynucleotides from a test sample may be identical to a nucleotide in the reference sequence.
  • a nucleotide detected at nearly 100% frequency in a sample may be different than a nucleotide in the reference sequence.
  • a first nucleotide that matches a reference nucleotide is detected at about 50% and a second nucleotide that is different than a reference nucleotide is detected at about 50%, this most likely indicates normal heterozygosity.
  • Heterozygosity may present at allele ratios divergent from 50:50, e.g., 60:40 or even 70:30.
  • the sample comprises a nucleotide detectable above noise at a frequency below (of above) an unambiguously heterozygote range (for example, less than about 45%, less than 40%, less than 30%, less than 20%, less than 10% or less than 5%), this can be attributed to the existence of somatic mutations in a percentage of the cells contributing DNA to the cfDNA population. These may come from disease cells, e.g., cancer cells. (The exact percentage is a function of tumor load.) If the frequency of somatic mutations at two different genetic loci are different, e.g., 16% at one locus and 5% at another locus, this indicates that the disease cells, e.g., the cancer cells, are heterogeneous.
  • somatic mutations In the case of DNA from solid tumors, which is expected to predominantly comprise tumor DNA, somatic mutations also can be detected by comparison to a reference sequence. Detection of somatic mutations that exist in 100% of the tumor cells may require reference to a standard sequence or information about known mutants to. However, the existence of sub-dominant sequences among the polynucleotide pool at different loci and at different relative frequencies, indicates tumor heterogeneity.
  • the profile may include genetic variants in genes that are known to be actionable. Knowledge of such variants can contribute to selecting therapeutic interventions, as therapies can be targeted to such variants. In the case of cancer, many actionable genetic variants are already known.
  • the copy number state of a gene should be reflected in the frequency of a genetic form of the gene in the sample.
  • a sequence variant may be detected at a frequency consistent with homozygosity or heterozygosity (e.g., about 100% or about 50%, respectively) with no copy number variation. This is consistent with a germ line polymorphism or mutation.
  • a sequence variant is detected at a level consistent with homozygosity (e.g., about 100%) but at amounts consistent with copy number variation, this is more likely to reflect the presence of disease cell polynucleotides having undergone gene amplification.
  • a sequence variant is detected at a level not inconsistent with heterozygosity (e.g., deviating somewhat from 50%) but at amounts consistent with copy number variation, this also is more likely to reflect the presence of disease cell polynucleotides; the diseased polynucleotides create some level of imbalance in allele frequency away from 50:50.
  • sequence variants in a gene detected at levels ideally consistent with heterozygosity in the germ line are more probably the product of a somatic mutation in disease cells if copy number variation also is detected in that gene.
  • tumor heterogeneity can be inferred when two genes are detected at different frequency but their copy number is relatively equal.
  • tumor homogeneity can be inferred when the difference in frequency between two sequence variants is consistent with difference in copy number for the two genes.
  • an EGFR variant is detected at 11% and a KRAS variant is detected at 5%, and no CNV is detected at these genes, the difference in frequency likely reflects tumor heterogeneity (e.g., all tumor cells carry an EGFR mutant and half the tumor cells also carry a KRAS mutant).
  • the EGFR gene carrying the mutant is detected at increased copy number, one consistent interpretation is a homogenous population of tumor cells, each cell carrying a mutant in the EGFR and KRAS genes, but in which the KRAS gene is duplicated. Accordingly, both the frequency of a sequence variant and a measure of CNV at the locus of the sequence variant in a sample can be determined. The frequency can then be corrected to reflect the relative number of cells bearing the variant by weighing the frequency based on dose per cell determined from the measure of CNV. This result is now more comparable in terms of number of cells carrying the variant to a sequence variant that does not vary in copy number.
  • a report of results from genetic variant analysis may be provided by a report generator, for example to a healthcare practitioner, e.g., a physician, to aid the interpretation of the test results (e.g., data) and selection of treatment options.
  • a report generated by a report generator may provide additional information, such as clinical lab results, that may be useful for diagnosing disease and selecting treatment options.
  • the report generator system can be a central data processing system configured to establish communications directly with: a remote data site or lab 2 , a medical practice/healthcare provider (treating professional) 4 , and/or a patient/subject 6 through communication links.
  • the lab 2 can be medical laboratory, diagnostic laboratory, medical facility, medical practice, point-of-care testing device, or any other remote data site capable of generating subject clinical information.
  • Subject clinical information includes but it is not limited to laboratory test data, e.g., analysis of genetic variants; imaging and X-ray data; examination results; and diagnosis.
  • the healthcare provider or practice 6 may include medical services providers, such as doctors, nurses, home health aides, technicians and physician's assistants, and the practice may be any medical care facility staffed with healthcare providers. In certain instances the healthcare provider/practice is also a remote data site. Where cancer is a disease to be treated, the subject may be afflicted with cancer, among other possible diseases or disorders.
  • Other clinical information for a cancer subject 6 can include the results of laboratory tests, e.g., analysis of genetic variants, metabolic panel, complete blood count, etc.; medical imaging data; and/or medical procedures directed to diagnosing the condition, providing a prognosis, monitoring the progression of the disease, determining relapse or remission, or combinations thereof.
  • the list of appropriate sources of clinical information for cancer includes, but it is not limited to, CT scans, MRI scans, ultrasound scans, bone scans, PET Scans, bone marrow test, barium X-ray, endoscopies, lymphangiograms, IVU (Intravenous urogram) or IVP (IV pyelogram), lumbar punctures, cystoscopy, immunological tests (anti-malignin antibody screen), and cancer marker tests.
  • the subject 6 's clinical information may be obtained from the lab 2 manually or automatically. Where simplicity of the system is desired, the information may be obtained automatically at predetermined or regular time intervals.
  • a regular time interval can refer to a time interval at which the collection of the laboratory data is carried out automatically by the methods and systems described herein based on a measurement of time such as hours, days, weeks, months, years etc.
  • the collection of data and processing is carried out at least once a day.
  • the transfer and collection of data is carried out about any of monthly, biweekly, weekly, several times a week or daily.
  • the retrieval of information may be carried out at predetermined time intervals, which may not be regular time intervals. For instance, a first retrieval step may occur after one week and a second retrieval step may occur after one month.
  • the transfer and collection of data can be customized according to the nature of the disorder that is being managed and the frequency of required testing and medical examinations of the subjects.
  • FIG. 9B shows an exemplary process to generate genetic reports, including a tumor response map and associated summary of alterations.
  • a tumor response map is a graphical representation of genetic information indicating changes over time in genetic information from a tumor, e.g., qualitative and quantitative changes. Such changes can reflect response of a subject to a therapeutic intervention.
  • This process can reduce error rates and bias that may be orders of magnitude higher than what is required to reliably detect de novo genetic variants associated with cancer.
  • the process can comprise first capturing genetic information by collecting body fluid samples as sources of genetic material (e.g., blood, saliva, sweat, urine, etc). Then, the process can comprise sequencing the materials ( 11 ). For example, polynucleotides in a sample can be sequenced, producing a plurality of sequence reads.
  • the tumor burden in a sample that comprises polynucleotides can be estimated as the relative number of sequence reads bearing a variant to the total number of sequence reads generated from the sample. Where copy number variants are analyzed, the tumor burden can be estimated as the relative excess (e.g., in the case of gene duplication) or relative deficit (e.g., in the case of gene elimination) of the total number of sequence reads at test and control loci. For example, a run may produce 1000 reads mapping to an oncogene locus of which 900 correspond to wild type and 100 correspond to a cancer mutant, indicating a copy number variant at this gene. More details on exemplary specimen collection and sequencing of the genetic materials are discussed below in FIGS. 10-11 .
  • genetic information can be processed ( 12 ). Genetic variants can then be identified.
  • the process can comprise determining the frequency of genetic variants in the sample containing the genetic material.
  • the process can comprise separating information from noise ( 13 ) if this process is noisy.
  • the sequencing methods for genetic analysis may have error rates.
  • the mySeq system of Illumina can produce percent error rates in the low single digits.
  • about 50 reads about 5% may be expected to include errors.
  • Certain methodologies, such as those described in WO 2014/149134 can significantly reduce the error rate. Errors create noise that can obscure signals from cancer present at low levels in a sample. For example, if a sample has a tumor burden at a level around the sequencing system error rate, e.g., around 0.1%-5%, it may be difficult to distinguish a signal corresponding to a genetic variant due to cancer from one due to noise.
  • Analysis of genetic variants may be used for diagnosing in the presence of noise.
  • the analysis can be based on the frequency of Sequence Variants or Level of CNV ( 14 ) and a diagnosis confidence indication or level for detecting genetic variants in the noise range can be established ( 15 ).
  • the process can comprise increasing the diagnosis confidence. This can be done using a plurality of measurements to increase confidence of diagnosis ( 16 ), or alternatively using measurements at a plurality of time points to determine whether cancer is advancing, in remission or stabilized ( 17 ).
  • the diagnostic confidence can be used to identify disease states.
  • cell free polynucleotides taken from a subject can include polynucleotides derived from normal cells, as well as polynucleotides derived from diseased cells, such as cancer cells. Polynucleotides from cancer cells may bear genetic variants, such as somatic cell mutations and copy number variants. When cell free polynucleotides from a sample from a subject are sequenced, these cancer polynucleotides are detected as sequence variants or as copy number variants.
  • Measurements of a parameter may be provided with a confidence interval. Tested over time, one can determine whether a cancer is advancing, stabilized or in remission by comparing confidence intervals over time. When confidence intervals overlap, one may not be able to tell whether disease is increasing or decreasing, because there is no statistically significant difference between the measures. However, where the confidence intervals do not overlap, this indicates the direction of disease. For example, comparing the lowest point on a confidence interval at one time point and the highest point on a confidence interval at a second time point indicates the direction.
  • the process can comprise generating genetic Report/Diagnosis.
  • the process can comprise generating genetic graph for a plurality of measurements showing mutation trend ( 18 ) and generating report showing treatment results and options ( 19 ).
  • FIGS. 10A-10C show in more details one embodiment for generating genetic reports and diagnosis (e.g., Report/Diagnosis).
  • FIG. 10C shows an exemplary pseudo-code executed by the system of FIG. 9A to process non-CNV reported mutant allele frequencies.
  • the system can process CNV reported mutant allele frequencies as well.
  • Samples comprising genetic material can be collected from a subject at a plurality of time points, that is, serially.
  • the genetic material can be sequenced, e.g., using a high-throughput sequencing system. Sequencing can target loci of interest to detect genetic variants, such genes bearing somatic mutations, genes that undergo copy number variation, or genes involved in gene fusions, for example, in cancer.
  • a quantitative measure of the genetic variants found can be determined.
  • the quantitative measure can be the frequency or percentage of a genetic variant among polynucleotides mapping to a locus, or the absolute number of sequence reads or polynucleotides mapping to a locus.
  • variants having a non-zero quantity at at least one time point can then be represented graphically through all time points. For example, in a collection of 1000 sequences, variant 1 may be found at time points 1, 2 and 3 in amounts of 50, 30 and 0, respectively. Variant 2 may be found in amounts 0, 10 and 20 at these time points. These amounts can be normalized, for variant 1, to 5%, 3% and 0%, and, for variant 2, 0%, 1% and 2%. A graphical representation showing the union of all non-zero results can indicate these amounts for both variants at all of the time points. The normalized amounts can be scaled so that each percentage is represented by a layer, for example, having height 1 mm.
  • the heights would be at time point 1: heights 5 mm (variant 1) and 0 mm (variant 2); at time point 2: heights 3 mm (variant 1) and 1 mm (variant 2), at time point 3: heights 0 mm (variant 1) and 2 mm (variant 2).
  • the graphical representation can be in the form of a stacked area graph, such as a streamgraph. A “zero” time point (before the first time point) can be represented by a point, with all values at 0.
  • the height of the quantity of the variants in the graphical representation can be, for example, relative or proportional to each other.
  • a variant frequency 5% at one time point could be represented with a height of twice that of a variant with frequency of 2.5% at the same time point.
  • the order of stacking can be chosen for ease of understanding.
  • variants can be stacked in order of quantity high to low from bottom to top. Or, they can be stacked in a streamgraph with the variant of largest initial amount in the middle, and other variants of decreasing quantity on either side.
  • the areas can be color coded based on variant.
  • Variants in the same gene can be shown in different hues of the same color.
  • KRAS mutants can be shown in different shades of blue, EGFR mutants in different shades of red.
  • the process can comprise receiving genetic information from a DNA sequencer ( 30 ).
  • the process can then comprise determining specific gene alterations and quantities thereof ( 32 ).
  • a tumor response map is generated.
  • the process can comprise normalizing the quantities for each gene alteration for rendering across all test points and then generates a scaling factor ( 34 ).
  • the term “normalize” generally refers to means adjusting values measured on different scales to a notionally common scale. For example, data measured at different points are converted/adjusted so that all values can be resized to a common scale.
  • the process can comprise rendering information on a tumor response map ( 36 ).
  • the process can comprise rendering alterations and relative heights using the determined scaling factor ( 38 ) and assigns a unique visual indicator for each alteration ( 40 ).
  • the process can comprise generating a summary of alterations and treatment options ( 42 ). Also, information from clinical trials that may help the particular genetic alterations and other helpful treatment suggestions is presented, along with explanations of terminology, test methodology, and other information is added to the report and rendered for the user.
  • the copy number variation may be reported as graph, indicating various positions in the genome and a corresponding increase or decrease or maintenance of copy number variation at each respective position. Additionally, copy number variation may be used to report a percentage score indicating how much disease material (or nucleic acids having a copy number variation) exists in the cell free polynucleotide sample.
  • the report includes annotations to help physicians interpret the results and recommend treatment options.
  • the annotating can include annotating a report for a condition in the NCCN Clinical Practice Guidelines in OncologyTM or the American Society of Clinical Oncology (ASCO) clinical practice guidelines.
  • the annotating can include listing one or more FDA-approved drugs for off-label use, one or more drugs listed in a Centers for Medicare and Medicaid Services (CMS) anti-cancer treatment compendia, and/or one or more experimental drugs found in scientific literature, in the report.
  • CMS Centers for Medicare and Medicaid Services
  • the annotating can include connecting a listed drug treatment option to a reference containing scientific information regarding the drug treatment option.
  • the scientific information can be from a peer-reviewed article from a medical journal.
  • the annotating can include providing a link to information on a clinical trial for a drug treatment option in the report.
  • the annotating can include presenting information in a pop-up box or fly-over box near provided drug treatment options in an electronic based report.
  • the annotating can include adding information to a report selected from the group consisting of one or more drug treatment options, scientific information concerning one or more drug treatment options, one or more links to scientific information regarding one or more drug treatment options, one or more links to citations for scientific information regarding one or more drug treatment options, and clinical trial information regarding one or more drug treatment options.
  • FIG. 10B shows an exemplary process to generate a tumor response map pathway which may be used by a healthcare practitioner, e.g., physician, for example to make patient care decisions.
  • the process can comprise first determining a global scaling factor ( 43 ).
  • the process can comprise transforming the absolute value into a relative metric/scale that may be more amenable for plotting (e.g. Multiply mutant allele frequency by 100 and take log of that value) and determines a global scaling factor using maximum observed value.
  • the process then involves visualizing information from the earliest test dataset ( 44 ).
  • Visualizing can comprise graphically representing the information on a user interface (e.g., a computer screen) or in tangible form (e.g., on a piece of paper).
  • a user interface e.g., a computer screen
  • tangible form e.g., on a piece of paper.
  • the process can comprise multiplying the scaling factor by a transformed value for each gene and use as a quantity indicator for plotting that variant, and then assigns a color/unique visual indicator for each alteration. Then the process can comprise visualizing information for subsequent test points ( 45 ) using the following pseudo-code:
  • Each subsequent panel denoting a test date may also include additional patient or intervention information that may correlate with the alteration changes seen in the remainder of the map. Similar scaling, plotting, and transformation may be also implemented on CNV and other types of DNA alterations (e.g. methylation) to display these quantities in separate or combined charts. These additional annotations may themselves also be quantifiable and similarly plotted on the map.
  • the process can then comprise determining a summary of alterations and treatment options ( 46 ).
  • the following actions are done:
  • the process can comprise including a trend indicator for that alteration over the different test date points.
  • Grouping of maximum mutant allele frequencies may also extend beyond just the genes they are harbored in to greater encapsulating annotations such as biological pathways, evidence level, etc.
  • FIGS. 10D-10I show one exemplary report generated by the system of FIG. 9A .
  • a patient identification section 52 provides patient information, reporting date, and physician contact information.
  • a tumor response map 54 includes a modified streamgraph 56 that shows tumor activities with unique colors for each mutant gene.
  • the graph 56 has accompanying summary explanation textbox 58 . More details are provided in a summary of alterations and treatment option section 60 .
  • the alterations 62 and 64 are presented in section 60 , along with mutation trend, mutant allele frequency, cell-free amplification, FDA Approved Drug Indication, FDA Approved Drugs with other Indications, and Clinical Drug Trial information.
  • FIGS. 10D-1, 10D-2, and 10D-3 provide enlarged views of FIG. 10D .
  • FIG. 10E shows an exemplary report section providing definitions, comments, and interpretation of the tests.
  • FIGS. 10E-1 and 10E-2 provide enlarged views of FIG. 10E .
  • FIG. 10F shows an exemplary detailed therapy result portion of the report.
  • FIGS. 10E-1 and 10E-2 provide enlarged views of FIG. 10F .
  • FIG. 10G shows an exemplary discussion of the clinical relevance of detected alterations.
  • FIGS. 10G-1 and 10G-2 provide enlarged views of FIG. 10G .
  • FIG. 10H shows potentially available medications that are going through clinical trials.
  • FIG. 10I shows the test methods and limitations thereof.
  • FIGS. 10I-1 and 10I-2 provide enlarged views of FIG. 10I .
  • FIG. 10J-10P shows various exemplary modified streamgraph 56 .
  • a streamgraph, or stream graph is a type of stacked area graph which is displaced around a central axis, resulting in a flowing, organic shape.
  • Streamgraphs are a generalization of stacked area graphs where the baseline is free. By shifting the baseline, it is possible to minimize the change in slope (or “wiggle”) in individual series, thereby making it easier to perceive the thickness of any given layer across the data.
  • FIG. 10J shows seven layers representing at least 8 mutants over three time periods, and a “0” time point (all values “0”).
  • FIG. 10K shows a single mutant over 4 time periods. No mutants are detected at the second, third and fourth time points.
  • FIG. 10L indicates frequency of dominant allele at each time point.
  • FIG. 10M shows a single time point with a total of four mutants in two genes. Mutants are identified by amino acid at a position changed (i.e., EGFR T790M).
  • One embodiment renders a streamgraph so that it is not x-axis reflective.
  • the modified graph applies a unique scaling to denote proportional attributes.
  • the graph can indicate the addition of new attributes over time.
  • the presence or absence of a mutation may be reflected in graphical form, indicating various positions in the genome and a corresponding increase or decrease or maintenance of a frequency of mutation at each respective position.
  • mutations may be used to report a percentage score indicating how much disease material exists in the cell free polynucleotide sample.
  • a confidence score may accompany each detected mutation, given known statistics of typical variances at reported positions in non-disease reference sequences. Mutations may also be ranked in order of abundance in the subject or ranked by clinically actionable importance.
  • the mapping of genome positions and copy number variation for the subject with cancer can indicate that a particular cancer is aggressive and resistant to treatment.
  • the subject may be monitored for a period and retested. If at the end of the period, the copy number variation profile, e.g., as depicted in a tumor response map, begins to increase dramatically, this may indicate that the current treatment is not working.
  • a comparison can also done with genetic profiles of other subjects. For example, if it is determined that this increase in copy number variation indicates that the cancer is advancing, then the original treatment regimen as prescribed is no longer treating the cancer and a new treatment is prescribed.
  • These reports can be submitted and accessed electronically via the internet. Analysis of sequence data may occur at a site other than the location of the subject. The report can be generated and transmitted to the subject's location. Via an internet enabled computer, the subject may access the reports reflecting his tumor burden.
  • an exemplary process receives genetic materials from blood sample or other body samples ( 1102 ).
  • the process can comprise converting the polynucleotides from the genetic materials into tagged parent nucleotides ( 1104 ).
  • the tagged parent nucleotides are amplified to produce amplified progeny polynucleotides ( 1106 ).
  • a subset of the amplified polynucleotides is sequenced to produce sequence reads ( 1108 ), which are grouped into families, each generated from a unique tagged parent nucleotide ( 1110 ).
  • the process can comprise assigning each family a confidence score for each family ( 1112 ).
  • a consensus is determined using prior readings. This is done by reviewing prior confidence score for each family, and if consistent prior confidence scores exists, then the current confidence score is increased ( 1114 ). If there are prior confidence scores, but they are inconsistent, the current confidence score is not modified in one embodiment ( 1116 ). In other embodiments, the confidence score is adjusted in a predetermined manner for inconsistent prior confidence scores. If this is a first time the family is detected, the current confidence score can be reduced as it may be a false reading ( 1118 ). The process can comprise inferring the frequency of the family at the locus in the set of tagged parent polynucleotides based on the confidence score. Then genetic test reports are generated as discussed above ( 1120 ).
  • the historical comparison can be used in conjunction with other consensus sequences mapping to a particular reference sequence to detect instances of genetic variation.
  • Consensus sequences mapping to particular reference sequences can be measured and normalized against control samples. Measures of molecules mapping to reference sequences can be compared across a genome to identify areas in the genome in which copy number varies, or heterozygosity is lost.
  • Consensus methods include, for example, linear or non-linear methods of building consensus sequences (e.g., voting, averaging, statistical, maximum a posteriori or maximum likelihood detection, dynamic programming, Bayesian, hidden Markov or support vector machine methods, etc.) derived from digital communication theory, information theory, or bioinformatics.
  • a stochastic modeling algorithm is applied to convert the normalized nucleic acid sequence read coverage for each window region to the discrete copy number states.
  • this algorithm may comprise one or more of the following: Hidden Markov Model, dynamic programming, support vector machine, Bayesian network, trellis decoding, Viterbi decoding, expectation maximization, Kalman filtering methodologies and neural networks.
  • cell free DNAs are extracted and isolated from a readily accessible bodily fluid such as blood, sweat, saliva, urine, etc.
  • a readily accessible bodily fluid such as blood, sweat, saliva, urine, etc.
  • cell free DNAs can be extracted using a variety of methods known in the art, including but not limited to isopropanol precipitation and/or silica based purification.
  • Cell free DNAs may be extracted from any number of subjects, such as subjects without cancer, subjects at risk for cancer, or subjects known to have cancer (e.g. through other means).
  • any of a number of different sequencing operations may be performed on the cell free polynucleotide sample.
  • Samples may be processed before sequencing with one or more reagents (e.g., enzymes, unique identifiers (e.g., barcodes), probes, etc.).
  • reagents e.g., enzymes, unique identifiers (e.g., barcodes), probes, etc.
  • the samples or fragments of samples may be tagged individually or in subgroups with the unique identifier. The tagged sample may then be used in a downstream application such as a sequencing reaction and individual molecules may be tracked to parent molecules.
  • the cell free polynucleotides can be tagged or tracked in order to permit subsequent identification and origin of the particular polynucleotide.
  • the assignment of an identifier to individual or subgroups of polynucleotides may allow for a unique identity to be assigned to individual sequences or fragments of sequences. This may allow acquisition of data from individual samples and is not limited to averages of samples.
  • nucleic acids or other molecules derived from a single strand may share a common tag or identifier and therefore may be later identified as being derived from that strand.
  • all of the fragments from a single strand of nucleic acid may be tagged with the same identifier or tag, thereby permitting subsequent identification of fragments from the parent strand.
  • gene expression products may be tagged in order to quantify expression.
  • a barcode or barcode in combination with sequence to which it is attached can be counted.
  • the systems and methods can be used as a PCR amplification control. In such cases, multiple amplification products from a PCR reaction can be tagged with the same tag or identifier. If the products are later sequenced and demonstrate sequence differences, differences among products with the same identifier can then be attributed to PCR error. Additionally, individual sequences may be identified based upon characteristics of sequence data for the read themselves.
  • the detection of unique sequence data at the beginning (start) and end (stop) portions of individual sequencing reads may be used, alone or in combination, with the length, or number of base pairs of each sequence read to assign unique identities to individual molecules. Fragments from a single strand of nucleic acid, having been assigned a unique identity, may thereby permit subsequent identification of fragments from the parent strand. This can be used in conjunction with bottlenecking the initial starting genetic material to limit diversity.
  • unique sequence data at the beginning (start) and end (stop) portions of individual sequencing reads and sequencing read length may be used, alone or combination, with the use of barcodes.
  • the barcodes may be unique as described herein. In other cases, the barcodes themselves may not be unique. In this case, the use of non-unique barcodes, in combination with sequence data at the beginning (start) and end (stop) portions of individual sequencing reads and sequencing read length may allow for the assignment of a unique identity to individual sequences. Similarly, fragments from a single strand of nucleic acid having been assigned a unique identity may thereby permit subsequent identification of fragments from the parent strand.
  • Sequencing methods may include, but are not limited to: high-throughput sequencing, pyrosequencing, sequencing-by-synthesis, single-molecule sequencing, nanopore sequencing, semiconductor sequencing, sequencing-by-ligation, sequencing-by-hybridization, RNA-Seq (Illumina), Digital Gene Expression (Helicos), Next generation sequencing, Single Molecule Sequencing by Synthesis (SMSS) (Helicos), massively-parallel sequencing, Clonal Single Molecule Array (Solexa), shotgun sequencing, Maxim-Gilbert sequencing, primer walking, and any other sequencing methods known in the art.
  • SMSS Single Molecule Sequencing by Synthesis
  • Solexa Single Molecule Array
  • Sequencing methods typically involve sample preparation, sequencing of polynucleotides in the prepared sample to produce sequence reads and bioinformatic manipulation of the sequence reads to produce quantitative and/or qualitative genetic information about the sample.
  • Sample preparation typically involves converting polynucleotides in a sample into a form compatible with the sequencing platform used. This conversion can involve tagging polynucleotides.
  • the tags comprise polynucleotide sequence tags. Conversion methodologies used in sequencing may not be 100% efficient. For example, it is not uncommon to convert polynucleotides in a sample with a conversion efficiency of about 1-5%, that is, about 1-5% of the polynucleotides in a sample are converted into tagged polynucleotides.
  • Polynucleotides that are not converted into tagged molecules are not represented in a tagged library for sequencing. Accordingly, polynucleotides having genetic variants represented at low frequency in the initial genetic material may not be represented in the tagged library and, therefore may not be sequenced or detected. By increasing conversion efficiency, the probability that a polynucleotide in the initial genetic material will be represented in the tagged library and, consequently, detected by sequencing is increased. Furthermore, rather than directly address the low conversion efficiency issue of library preparation, most protocols to date call for greater than 1 microgram of DNA as input material. However, when input sample material is limited or detection of polynucleotides with low representation is desired, high conversion efficiency can efficiently sequence the sample and/or to adequately detect such polynucleotides.
  • mutation detection may be performed on selectively enriched regions of the genome or transcriptome purified and isolated ( 1302 ).
  • specific regions which may include but are not limited to genes, oncogenes, tumor suppressor genes, promoters, regulatory sequence elements, non-coding regions, miRNAs, snRNAs and the like may be selectively amplified from a total population of cell free polynucleotides. This may be performed as herein described.
  • multiplex sequencing may be used, with or without barcode labels for individual polynucleotide sequences.
  • sequencing may be performed using any nucleic acid sequencing platforms known in the art. This step generates a plurality of genomic fragment sequence reads ( 1304 ).
  • a reference sequence is obtained from a control sample, taken from another subject.
  • the control subject may be a subject known to not have known genetic aberrations or disease.
  • these sequence reads may contain barcode information. In other examples, barcodes are not utilized.
  • a quality score may be a representation of reads that indicates whether those reads may be useful in subsequent analysis based on a threshold. In some cases, some reads are not of sufficient quality or length to perform the subsequent mapping step. Sequencing reads with a quality score at least 90%, 95%, 99%, 99.9%, 99.99% or 99.999% may be filtered out of the data set. In other cases, sequencing reads assigned a quality scored at least 90%, 95%, 99%, 99.9%, 99.99% or 99.999% may be filtered out of the data set. In step 1306 , the genomic fragment reads that meet a specified quality score threshold are mapped to a reference genome, or a reference sequence that is known not to contain mutations.
  • mapping score may be a representation or reads mapped back to the reference sequence indicating whether each position is or is not uniquely mappable.
  • reads may be sequences unrelated to mutation analysis. For example, some sequence reads may originate from contaminant polynucleotides. Sequencing reads with a mapping score at least 90%, 95%, 99%, 99.9%, 99.99% or 99.999% may be filtered out of the data set. In other cases, sequencing reads assigned a mapping scored less than 90%, 95%, 99%, 99.9%, 99.99% or 99.999% may be filtered out of the data set.
  • bases that do not meet the minimum threshold for mappability, or low quality bases may be replaced by the corresponding bases as found in the reference sequence.
  • the frequency of variant bases may be calculated as the number of reads containing the variant divided by the total number of reads 1308 after ascertaining read coverage and identifying variant bases relative to the control sequence in each read. This may be expressed as a ratio for each mappable position in the genome.
  • the frequencies of all four nucleotides, cytosine, guanine, thymine, adenine can be analyzed in comparison to the reference sequence.
  • a stochastic or statistical modeling algorithm can be applied to convert the normalized ratios for each mappable position to reflect frequency states for each base variant.
  • this algorithm may comprise one or more of the following: Hidden Markov Model, dynamic programming, support vector machine, Bayesian or probabilistic modeling, trellis decoding, Viterbi decoding, expectation maximization, Kalman filtering methodologies, and neural networks.
  • the discrete mutation states of each base position can be utilized to identify a base variant with high frequency of variance as compared to the baseline of the reference sequence.
  • the baseline might represent a frequency of at least 0.0001%, 0.001%, 0.01%, 0.1%, 1.0%, 2.0%, 3.0%, 4.0% 5.0%, 10%, or 25%.
  • the baseline might represent a frequency of at least 0.0001%, 0.001%, 0.01%, 0.1%, 1.0%, 2.0%, 3.0%, 4.0% 5.0%. 10%, or 25%.
  • all adjacent base positions with the base variant or mutation can be merged into a segment to report the presence or absence of a mutation.
  • various positions can be filtered before they are merged with other segments.
  • the variant with largest deviation for a specific position in the sequence derived from the subject as compared to the reference sequence can be identified as a mutation.
  • a mutation may be a cancer mutation.
  • a mutation might be correlated with a disease state.
  • a mutation or variant may comprise a genetic aberration that includes, but is not limited to a single base substitution, or small indels, transversions, translocations, inversion, deletions, truncations or gene truncations.
  • a mutation may be at most 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15 or 20 nucleotides in length. On other cases a mutation may be at least 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15 or 20 nucleotides in length.
  • a consensus is determined using prior readings. This is done by reviewing prior confidence score for the corresponding bases, and if consistent prior confidence scores exists, then the current confidence score is increased ( 1314 ). If there are prior confidence scores, but they are inconsistent, the current confidence score is not modified in one embodiment ( 1316 ). In other embodiments, the confidence score is adjusted in a predetermined manner for inconsistent prior confidence scores. If this is a first time the family is detected, the current confidence score can be reduced as it may be a false reading ( 1318 ). The process can comprise then converting the frequency of variance per each base into discrete variant states for each base position ( 1320 ).
  • Cancers cells as most cells, can be characterized by a rate of turnover, in which old cells die and are replaced by newer cells. Generally dead cells, in contact with vasculature in a given subject, may release DNA or fragments of DNA into the blood stream. This is also true of cancer cells during various stages of the disease. Cancer cells may also be characterized, dependent on the stage of the disease, by various genetic aberrations such as copy number variation as well as mutations. This phenomenon may be used to detect the presence or absence of cancers individuals using the methods and systems described herein.
  • blood from subjects at risk for cancer may be drawn and prepared as described herein to generate a population of cell free polynucleotides.
  • this might be cell free DNA.
  • the systems and methods of the disclosure may be employed to detect mutations or copy number variations that may exist in certain cancers present. The method may help detect the presence of cancerous cells in the body, despite the absence of symptoms or other hallmarks of disease.
  • the types and number of cancers that may be detected may include but are not limited to blood cancers, brain cancers, lung cancers, skin cancers, nose cancers, throat cancers, liver cancers, bone cancers, lymphomas, pancreatic cancers, skin cancers, bowel cancers, rectal cancers, thyroid cancers, bladder cancers, kidney cancers, mouth cancers, stomach cancers, solid state tumors, heterogeneous tumors, homogenous tumors and the like.
  • the system and methods may be used to detect any number of genetic aberrations that may cause or result from cancers. These may include but are not limited to mutations, mutations, indels, copy number variations, transversions, translocations, inversion, deletions, aneuploidy, partial aneuploidy, polyploidy, chromosomal instability, chromosomal structure alterations, gene fusions, chromosome fusions, gene truncations, gene amplification, gene duplications, chromosomal lesions, DNA lesions, abnormal changes in nucleic acid chemical modifications, abnormal changes in epigenetic patterns, abnormal changes in nucleic acid methylation infection and cancer.
  • the systems and methods described herein may also be used to help characterize certain cancers.
  • Genetic data produced from the system and methods of this disclosure may allow practitioners to help better characterize a specific form of cancer. Often times, cancers are heterogeneous in both composition and staging. Genetic profile data may allow characterization of specific sub-types of cancer that may be important in the diagnosis or treatment of that specific sub-type. This information may also provide a subject or practitioner clues regarding the prognosis of a specific type of cancer.
  • the systems and methods provided herein may be used to monitor already known cancers, or other diseases in a particular subject. This may allow either a subject or practitioner to adapt treatment options in accord with the progress of the disease.
  • the systems and methods described herein may be used to construct genetic profiles of a particular subject of the course of the disease. In some instances, cancers can progress, becoming more aggressive and genetically unstable. In other examples, cancers may remain benign, inactive or dormant. The system and methods of this disclosure may be useful in determining disease progression.
  • the systems and methods described herein may be useful in determining the efficacy of a particular treatment option.
  • successful treatment options may actually increase the amount of copy number variation or mutations detected in subject's blood if the treatment is successful as more cancers may die and shed DNA. In other examples, this may not occur.
  • certain treatment options may be correlated with genetic profiles of cancers over time. This correlation may be useful in selecting a therapy.
  • the systems and methods described herein may be useful in monitoring residual disease or recurrence of disease.
  • the methods and systems described herein may not be limited to detection of mutations and copy number variations associated with only cancers.
  • Various other diseases and infections may result in other types of conditions that may be suitable for early detection and monitoring.
  • genetic disorders or infectious diseases may cause a certain genetic mosaicism within a subject. This genetic mosaicism may cause copy number variation and mutations that could be observed.
  • the system and methods of the disclosure may also be used to monitor the genomes of immune cells within the body. Immune cells, such as B cells, may undergo rapid clonal expansion upon the presence certain diseases. Clonal expansions may be monitored using copy number variation detection and certain immune states may be monitored. In this example, copy number variation analysis may be performed over time to produce a profile of how a particular disease may be progressing.
  • systems and methods of this disclosure may also be used to monitor systemic infections themselves, as may be caused by a pathogen such as a bacteria or virus.
  • Copy number variation or even mutation detection may be used to determine how a population of pathogens are changing during the course of infection. This may be particularly important during chronic infections, such as HIV/AIDs or Hepatitis infections, whereby viruses may change life cycle state and/or mutate into more virulent forms during the course of infection.
  • transplanted tissue undergoes a certain degree of rejection by the body upon transplantation.
  • the methods of this disclosure may be used to determine or profile rejection activities of the host body, as immune cells attempt to destroy transplanted tissue. This may be useful in monitoring the status of transplanted tissue as well as altering the course of treatment or prevention of rejection.
  • a disease may be heterogeneous. Disease cells may not be identical.
  • some tumors are known to comprise different types of tumor cells, some cells in different stages of the cancer.
  • heterogeneity may comprise multiple foci of disease. Again, in the example of cancer, there may be multiple tumor foci, perhaps where one or more foci are the result of metastases that have spread from a primary site.
  • the methods of this disclosure may be used to generate or profile, fingerprint or set of data that is a summation of genetic information derived from different cells in a heterogeneous disease.
  • This set of data may comprise copy number variation and mutation analyses alone or in combination.
  • systems and methods of the disclosure may be used to diagnose, prognose, monitor or observe cancers or other diseases of fetal origin. That is, these methodologies may be employed in a pregnant subject to diagnose, prognose, monitor or observe cancers or other diseases in a unborn subject whose DNA and other polynucleotides may co-circulate with maternal molecules.
  • these reports are submitted and accessed electronically via the internet. Analysis of sequence data occurs at a site other than the location of the subject. The report is generated and transmitted to the subject's location. Via an internet enabled computer, the subject accesses the reports reflecting his tumor burden.
  • the annotated information can be used by a health care provider to select other drug treatment options and/or provide information about drug treatment options to an insurance company.
  • the method can include annotating the drug treatment options for a condition in, for example, the NCCN Clinical Practice Guidelines in OncologyTM or the American Society of Clinical Oncology (ASCO) clinical practice guidelines.
  • the drug treatment options that are stratified in a report can be annotated in the report by listing additional drug treatment options.
  • An additional drug treatment can be an FDA-approved drug for an off-label use.
  • OBRA Omnibus Budget Reconciliation Act
  • the drugs used for annotating lists can be found in CMS approved compendia, including the National Comprehensive Cancer Network (NCCN) Drugs and Biologics CompendiumTM, Thomson Micromedex DrugDex®, Elsevier Gold Standard's Clinical Pharmacology compendium, and American Hospital Formulary Service—Drug Information Compendium®.
  • the drug treatment options can be annotated by listing an experimental drug that may be useful in treating a cancer with one or more molecular markers of a particular status.
  • the experimental drug can be a drug for which in vitro data, in vivo data, animal model data, pre-clinical trial data, or clinical-trial data are available.
  • the drug treatment options can be annotated by providing a link on an electronic based report connecting a listed drug to scientific information regarding the drug.
  • a link can be provided to information regarding a clinical trial for a drug (clinicaltrials.gov). If the report is provided via a computer or computer website, the link can be a footnote, a hyperlink to a website, a pop-up box, or a fly-over box with information, etc.
  • the report and the annotated information can be provided on a printed form, and the annotations can be, for example, a footnote to a reference.
  • the information for annotating one or more drug treatment options in a report can be provided by a commercial entity that stores scientific information.
  • a health care provider can treat a subject, such as a cancer patient, with an experimental drug listed in the annotated information, and the health care provider can access the annotated drug treatment option, retrieve the scientific information (e.g., print a medical journal article) and submit it (e.g., a printed journal article) to an insurance company along with a request for reimbursement for providing the drug treatment.
  • Physicians can use any of a variety of Diagnosis-related group (DRG) codes to enable reimbursement.
  • DSG Diagnosis-related group
  • a drug treatment option in a report can also be annotated with information regarding other molecular components in a pathway that a drug affects (e.g., information on a drug that targets a kinase downstream of a cell-surface receptor that is a drug target).
  • the drug treatment option can be annotated with information on drugs that target one or more other molecular pathway components.
  • the identification and/or annotation of information related to pathways can be outsourced or subcontracted to another company.
  • the annotated information can be, for example, a drug name (e.g., an FDA approved drug for off-label use; a drug found in a CMS approved compendium, and/or a drug described in a scientific (medical) journal article), scientific information concerning one or more drug treatment options, one or more links to scientific information regarding one or more drugs, clinical trial information regarding one or more drugs (e.g., information from clinicaltrials.gov/), one or more links to citations for scientific information regarding drugs, etc.
  • a drug name e.g., an FDA approved drug for off-label use; a drug found in a CMS approved compendium, and/or a drug described in a scientific (medical) journal article
  • scientific information concerning one or more drug treatment options e.g., an FDA approved drug for off-label use; a drug found in a CMS approved compendium, and/or a drug described in a scientific (medical) journal article
  • scientific information concerning one or more drug treatment options e.g.,
  • the annotated information can be inserted into any location in a report.
  • Annotated information can be inserted in multiple locations on a report.
  • Annotated information can be inserted in a report near a section on stratified drug treatment options.
  • Annotated information can be inserted into a report on a separate page from stratified drug treatment options.
  • a report that does not contain stratified drug treatment options can be annotated with information.
  • the system can also include reports on the effects of drugs on sample (e.g. tumor cells) isolated from a subject (e.g. cancer patient).
  • sample e.g. tumor cells
  • An in vitro culture using a tumor from a cancer patient can be established using techniques known to those skilled in the art.
  • the system can also include high-throughput screening of FDA approved off-label drugs or experimental drugs using said in vitro culture and/or xenograft model.
  • the system can also include monitoring tumor antigen for recurrence detection.
  • the system can provide internet enabled access of reports of a subject with cancer.
  • the system can use a handheld DNA sequencer or a desktop DNA sequencer.
  • the DNA sequencer is a scientific instrument used to automate the DNA sequencing process. Given a sample of DNA, a DNA sequencer is used to determine the order of the four bases: adenine, guanine, cytosine, and thymine. The order of the DNA bases is reported as a text string, called a read.
  • Some DNA sequencers can be also considered optical instruments as they analyze light signals originating from fluorochromes attached to nucleotides.
  • the data is sent by the DNA sequencers over a direct connection or over the internet to a computer for processing.
  • the data processing aspects of the system can be implemented in digital electronic circuitry, or in computer hardware, firmware, software, or in combinations of them.
  • Data processing apparatus of the invention can be implemented in a computer program product tangibly embodied in a machine-readable storage device for execution by a programmable processor; and data processing method steps of the invention can be performed by a programmable processor executing a program of instructions to perform functions of the invention by operating on input data and generating output.
  • the data processing aspects of the invention can be implemented advantageously in one or more computer programs that are executable on a programmable system including at least one programmable processor coupled to receive data and instructions from and to transmit data and instructions to a data storage system, at least one input device, and at least one output device.
  • Each computer program can be implemented in a high-level procedural or object-oriented programming language, or in assembly or machine language, if desired; and, in any case, the language can be a compiled or interpreted language.
  • Suitable processors include, by way of example, both general and special purpose microprocessors. Generally, a processor will receive instructions and data from a read-only memory and/or a random access memory.
  • Storage devices suitable for tangibly embodying computer program instructions and data include all forms of nonvolatile memory, including by way of example semiconductor memory devices, such as EPROM, EEPROM, and flash memory devices; magnetic disks such as internal hard disks and removable disks; magneto-optical disks; and CD-ROM disks. Any of the foregoing can be supplemented by, or incorporated in, ASICs (application-specific integrated circuits).
  • the invention can be implemented using a computer system having a display device such as a monitor or LCD (liquid crystal display) screen for displaying information to the user and input devices by which the user can provide input to the computer system such as a keyboard, a two-dimensional pointing device such as a mouse or a trackball, or a three-dimensional pointing device such as a data glove or a gyroscopic mouse.
  • the computer system can be programmed to provide a graphical user interface through which computer programs interact with users.
  • the computer system can be programmed to provide a virtual reality, three-dimensional display interface.
  • the methods of this disclosure allow one to provide therapeutic interventions more precisely directed to the form of a disease in a subject, and to calibrate these therapeutic interventions over time.
  • This precision reflects, in part, the precision by which one is able to profile the whole body tumor status of a subject as reflected in tumor heterogeneity.
  • the therapeutic intervention is more effective against cancers with this profile than against cancers with any single one of these variants.
  • a therapeutic intervention is an intervention that produces a therapeutic effect, (e.g., is therapeutically effective).
  • Therapeutically effective interventions prevent, slow the progression of, improve the condition of (e.g., causes remission of), or cure a disease, such as a cancer.
  • a therapeutic intervention can include, for example, administration of a treatment, such as chemotherapy, radiation therapy, surgery, immunotherapy, administration of a pharmaceutical or a nutraceutical, or, a change in behavior, such as diet.
  • One measure of therapeutic effectiveness is effectiveness for at least 90% of subjects undergoing the intervention over at least 100 subjects.
  • a therapeutic intervention is determined that takes into account both the type of genetic variants found in the disease cells and their relative amounts (e.g., proportion).
  • the therapeutic intervention can treat the subject as if each clonal variant were a different cancer to be treated independently.
  • these variants may be left out of the therapeutic intervention until they rise to a clinical threshold or significant relative frequency (e.g., greater than the threshold stated above).
  • a therapeutic intervention can include treatments effective against diseases with each of the genetic variants.
  • genetic variants such as mutant forms of a gene or gene amplification, may be detected in several genes (e.g., a major clone and a minor clone).
  • Each of these forms may be actionable, that is, a treatment may be known for which cancers with the particular variant are responsive.
  • the profile of tumor heterogeneity may indicate that one of the variants is present in the polynucleotides at, for example, five times the level of each of the other two variants.
  • a therapeutic intervention can be determined that involves delivering three different drugs to the subject, each drug relatively more effective against cancers bearing each of the variants.
  • the drugs can be delivered as a cocktail, or sequentially.
  • the drugs can be administered in doses stratified to reflect the relative amounts of the variants in the DNA.
  • a drug effective against the most common variant can be administered in greater amount than drugs effective against the two less common variants.
  • the profile of tumor heterogeneity can show the presence of a sub-population of cancer cells bearing a genetic variant that is resistant to a drug to which the disease typically responds.
  • the therapeutic intervention can involve including both a first drug effective against tumor cells without the resistance variant and a second drug effective against tumor cells with the resistant variant.
  • doses can be stratified to reflect relative amounts of each variant detected in the profile.
  • changes in the profile of tumor heterogeneity are examined over time, and therapeutic interventions are developed to treat the changing tumor.
  • disease heterogeneity can be determined at a plurality of different times.
  • profiling methods of this disclosure more precise inferences can be made about tumor evolution. This allows the practitioner to monitor the evolution of the disease, in particular as new clonal sub-populations emerge after remission effected by a first wave of therapy.
  • therapeutic interventions can be calibrated over time to treat the changing tumor.
  • a profile may show that a cancer has a form that is responsive to a certain treatment. The treatment is delivered and the tumor burden is seen to decrease over time.
  • a genetic variant is found in the tumor indicating the presence of a population of cancer cells that is not responsive to the treatment.
  • a new therapeutic intervention is determined that targets the cells bearing the marker of non-responsiveness.
  • a dominant tumor form can eventually give way through Darwinian selection to cancer cells carrying mutants that render the cancer unresponsive to the therapy regimen. Appearance of these resistance mutants can be delayed through methods of this disclosure.
  • a subject is subjected to one or more pulsed therapy cycles, each pulsed therapy cycle comprising a first period during which a drug is administered at a first amount and a second cycle during which the drug is administered at a second, reduced amount.
  • the first period is characterized by a tumor burden detected above a first clinical level.
  • the second period is characterized by a tumor burden detected below a second clinical level.
  • First and second clinical levels can be different in different pulsed therapy cycles. So, for example, the first clinical level can be lower in succeeding cycles.
  • a plurality of cycles can include at least 2, 3, 4, 5, 6, 7, 8 or more cycles.
  • the BRAF mutant V600E may be detected in disease cell polynucleotides at an amount indicating a tumor burden of 5% in cfDNA.
  • Chemotherapy can commence with dabrafenib. Subsequent testing can show that the amount of the BRAF mutant in the cfDNA falls below 0.5% or to undetectable levels. At this point, dabrafenib therapy can stop or be significantly curtailed. Further subsequent testing may find that DNA bearing the BRAF mutation has risen to 2.5% of polynucleotides in cfDNA. At this point, dabrafenib therapy is re-started, e.g., at the same level as the initial treatment. Subsequent testing may find that DNA bearing the BRAF mutation has decreased to 0.5% of polynucleotides in cfDNA. Again, dabrafenib therapy is stopped or reduced. The cycle can be repeated a number of times.
  • FIG. 7 shows an exemplary course of monitoring and treatment of disease in a subject.
  • a subject tested at the time of blood draw 1 has a tumor burden of 1.4% and presents with genetic alterations in genes 1 , 2 and 3 .
  • the subject is treated with Drug A. After a time, treatment is discontinued.
  • a second blood draw shows the cancer in remission.
  • a third blood draw indicates that the cancer has recurred, in this instance, presenting with a genetic variant in Gene 4 .
  • the subject is now put on a course of Drug B, to which cancers having this variant are responsive.
  • a therapeutic intervention can be changed upon detection of the rise of a mutant form resistant to an original drug.
  • cancers with the EGFR mutation L858R respond to therapy with erlotinib.
  • cancers with the EGFR mutation T790M are resistant to erlotinib.
  • they are responsive to ruxolitinib.
  • a method of this disclosure involves monitoring changes in tumor profile and changing a therapeutic intervention when a genetic variant associated with drug resistance rises to a predetermined clinical level.
  • a database is built in which genetic information from serial samples collected from cancer patients is recorded.
  • This database may also contain intervening treatment and other clinically relevant information, such as, weight, adverse effects, histological testing, blood testing, radiographic information, prior treatments, cancer type, etc.
  • Serial test results can be used to infer efficacy of treatment, especially when used with blood samples, which can give a more unbiased estimate of tumor burden than self-reporting or radiographic reporting by a medical practitioner.
  • Treatment efficacy can be clustered by those with similar genomic profiles and vice versa. Genomic profiles can be organized around, for example, primary genetic alteration, secondary genetic alteration(s), relative amounts of these genetic alterations, and tumor load. This database can be used for decision support for subsequent patients.
  • Both germline and somatic alterations can be used for determining treatment efficacy as well.
  • Acquired resistance alterations that can also be inferred from the database when treatments that were effective initially begin to fail. This failure can be detected through radiographic, blood or other means.
  • the primary data used for inference of acquired resistance mechanisms are genomic tumor profiles collected after treatment per patient. This data can also be used to place quantitative bounds on likely treatment response as well as predict time to treatment failure. Based on likely acquired resistance alterations for a given treatment and tumor genomic profile, a treatment regimen can be modified to suppress acquisition of most likely resistance alterations.
  • FIG. 5 shows a computer system 1501 that is programmed or otherwise configured to implement the methods of the present disclosure.
  • the computer system 1501 includes a central processing unit (CPU, also “processor” and “computer processor” herein) 1505 .
  • the computer system 1501 also includes memory or memory location 1510 (e.g., random-access memory, read-only memory, flash memory), electronic storage unit 1515 (e.g., hard disk), communication interface 1520 (e.g., network adapter) for communicating with one or more other systems, and peripheral devices 1525 , such as cache, other memory, data storage and/or electronic display adapters.
  • memory or memory location 1510 e.g., random-access memory, read-only memory, flash memory
  • electronic storage unit 1515 e.g., hard disk
  • communication interface 1520 e.g., network adapter
  • peripheral devices 1525 such as cache, other memory, data storage and/or electronic display adapters.
  • the memory 1510 , storage unit 1515 , interface 1520 and peripheral devices 1525 are in communication with the CPU 1505 through a communication bus (solid lines).
  • the storage unit 1515 can be a data storage unit (or data repository) for storing data.
  • the computer system 1501 can be operatively coupled to a computer network (“network”) 1530 with the aid of the communication interface 1520 .
  • the network 1530 can be the Internet, an internet and/or extranet, or an intranet and/or extranet that is in communication with the Internet.
  • the network 1530 in some cases is a telecommunication and/or data network.
  • the network 1530 can include one or more computer servers, which can enable distributed computing, such as cloud computing.
  • the CPU 1505 can execute a sequence of machine-readable instructions, which can be embodied in a program or software.
  • the instructions may be stored in a memory location, such as the memory 1510 .
  • the storage unit 1515 can store files, such as drivers, libraries and saved programs.
  • the computer system 1501 can communicate with one or more remote computer systems through the network 1530 .
  • Methods as described herein can be implemented by way of machine (e.g., computer processor) executable code stored on an electronic storage location of the computer system 1501 , such as, for example, on the memory 1510 or electronic storage unit 1515 .
  • the machine executable or machine readable code can be provided in the form of software. Aspects of the systems and methods provided herein, such as the computer system 1501 , can be embodied in programming.
  • the computer system 1501 can include or be in communication with an electronic display that comprises a user interface (UI) for providing, for example, one or more results of sample analysis.
  • UI user interface
  • Nucleotide positions (e.g., loci) in the genome can be designated by number, as depicted in FIG. 2 . Positions at which about 100% of the base calls are identical to the reference sequence or at which about 100% of the base calls are different than the reference sequence are inferred to represent homozygosity of the cfDNA (presumed normal). Positions at which about 50% of the base calls are identical to the reference sequence are inferred to represent heterozygosity of the cfDNA (also presumed normal). Positions at which the percentage of base calls at a locus are substantially below 50% and above the detection limit of the base calling system are inferred to represent tumor-associated genetic variants.
  • the sample is subjected to proteinase K digestion.
  • DNA is precipitated with isopropanol.
  • DNA is captured on a DNA purification column (e.g., a QIAamp DNA Blood Mini Kit) and eluted in 100 ⁇ l solution.
  • DNAs below 500 bp are selected with Ampure SPRI magnetic bead capture (PEG/salt).
  • HGE haploid genome equivalents
  • High-efficiency DNA tagging (>80%) is performed by end repair, A-tailing and sticky-end ligation with 2 different octomers (i.e., 4 combinations) with overloaded hairpin adaptors.
  • 2.5 ng DNA i.e. approximately 800 HGE
  • Each hairpin adaptor comprises a random sequence on its non-complementary portion. Both ends of each DNA fragment are attached with hairpin adaptors.
  • Each tagged fragment can be identified by a combination of the octomer sequence on the hairpin adaptors and endogenous portions of the insert sequence.
  • Tagged DNA is amplified by 12 cycles of PCR to produce about 1-7 ⁇ g DNA that contain approximately 500 copies of each of the 800 HGE in the starting material.
  • Buffer optimization, polymerase optimization and cycle reduction may be performed to optimize the PCR reactions.
  • Amplification bias e.g., non-specific bias, GC bias, and/or size bias are also reduced by optimization.
  • Noise(s) e.g., polymerase-introduced errors are reduced by using high-fidelity polymerases.
  • Sequences may be enriched as follow: DNAs with regions of interest (ROI) are captured using biotin-labeled bead with probe to ROIs. The ROIs are amplified with 12 cycles of PCR to generate a 2000 times amplification.
  • ROI regions of interest
  • 0.1 to 1% of the sample are used for sequencing.
  • the resulting DNA is then denatured and diluted to 8 pM and loaded into an Illumina sequencer.
  • Sequence reads are grouped into families, with about 10 sequence reads in each family. Families are collapsed into consensus sequences by voting (e.g., biased voting) each position in a family. A base is called for consensus sequence if 8 or 9 members agree. A base is not called for consensus sequence if no more than 60% of the members agree.
  • the resulting consensus sequences are mapped to a reference genome, such as hg19.
  • a reference genome such as hg19.
  • Each base in a consensus sequence is covered by about 3000 different families.
  • a quality score for each sequence is calculated and sequences are filtered based on their quality scores.
  • Base calls at each position in a consensus sequence are compared with the HG-19 reference sequence. At each position at which a base call differs from the reference sequence, the identity of the different base or bases, and their percentage as a function of total base calls at the locus is determined and reported.
  • Sequence variation is detected by counting distribution of bases at each locus. If 98% of the reads have the same base (homozygous) and 2% have a different base, the locus is likely to have a sequence variant, presumably from cancer DNA.
  • CNV is detected by counting the total number of sequences (bases) mapping to a locus and comparing with a control locus.
  • CNV analysis is performed specific regions, including regions on ALK, APC, BRAF, CDKN2A, EGFR, ERBB2, FBXW7, KRAS, MYC, NOTCH1, NRAS, PIK3CA, PTEN, RB1, TP53, MET, AR, ABL1, AKT1, ATM, CDH1, CSF1R, CTNNB1, ERBB4, EZH2, FGFR1, FGFR2, FGFR3, FLT3, GNA11, GNAQ, GNAS, HNF1A, HRAS, IDH1, IDH2, JAK2, JAK3, KDR, KIT, MLH1, MPL, NPM1, PDGFRA, PROC, PTPN11, RET,SMAD4, SMARCB1, SMO, SRC, STK11, VHL, TERT, CCND1, CDK4, CDKN2
  • fragments are amplified and the sequences of amplified fragments are read and aligned, the fragments are subjected to base calling. Variations in the number of amplified fragments and unseen amplified fragments can introduce errors in base calling. These variations are corrected by calculating the number of unseen amplified fragments.
  • sequence readouts can come from two types of fragments: double-strand fragments and single-strand fragments. The following is a theoretical example of calculating the total number of unseen molecules in a sample.
  • N is the total number of molecules in the sample.
  • P is the probability of seeing a strand.
  • Q is the probability of not detecting a strand.
  • An assay is used to analyze a panel of genes to identify genetic variants in cancer-associated somatic variants with high sensitivity.
  • Cell-free DNA is extracted from plasma of a patient and amplified by PCR. Genetic variants are analyzed by massively parallel sequencing of the amplified target genes. For one set of genes, all exons are sequenced as such sequencing coverage had shown to have clinically utility (Table 3). For another set of genes, sequencing coverage included those exons with a previously reported somatic mutation (Table 4). The minimum detectable mutant allele (limit of detection) is dependent on the patient's sample cell-free DNA concentration, which varied from less than 10 to over 1,000 genomic equivalents per mL of peripheral blood. Amplification may not be detected in samples with lower amounts of cell-free DNA and/or low-level gene copy amplification. Certain sample or variant characteristics resulted in reduced analytic sensitivity, such as low sample quality or improper collection.
  • the percentage of genetic variants found in cell-free DNA circulating in blood is related to the unique tumor biology of this patient. Factors that affected the amount/percentages of detected genetic variants in circulating cell-free DNA in blood include tumor growth, turn-over, size, heterogeneity, vascularization, disease progression or treatment. Table 5 annotates the percentage, or allele frequency, of altered circulating cell-free DNA (% cfDNA) detected in this patient. Some of the detected genetic variants are listed in descending order by % cfDNA.
  • Genetic variants are detected in the circulating cell-free DNA isolated from this patient's blood specimen. These genetic variants are cancer-associated somatic variants, some of which have been associated with either increased or reduced clinical response to specific treatment. “Minor Alterations” are defined as those alterations detected at less than 10% the allele frequency of “Major Alterations”. A Major Alteration is the predominant alteration at a locus. The detected allele frequencies of these alterations (Table 5) and associated treatments for this patient are annotated.
  • a nucleotide detected at at least 98.8% frequency in the sample is different than a nucleotide in the reference sequence, indicating homozygosity at these loci.
  • T was detected rather than reference nucleotide C in 100% of cases.
  • a nucleotide detected at between 41.4% and 55% frequency in the sample is different than a nucleotide in the reference sequence, indicating heterozygosity at these loci.
  • G was detected rather than reference nucleotide A in 50% of cases.
  • nucleotide detected at less than 9% frequency is different than a nucleotide in the reference sequence.
  • BRAF (140453136 A>T, 8.9%)
  • NRAS 115256530 G>T 2.6%)
  • JAK2 5073770 G>T 1.5%). They are presumed to be somatic mutations from cancer DNA.
  • the relative amounts of tumor-associated genetic variants are calculated.
  • the ratio of amounts of BRAF:NRAS:JAK2 is 8.9:2.6:1.5, or 1:0.29: 0.17. From this result one can infer the presence of tumor heterogeneity. For example, one possible interpretation is that 100% of tumor cells contain a variant in BRAF, 83% contain variants in BRAF and NRAS, and 17% contain variants in BRAF, NRAS and JAK2. However, analysis of CNV may show amplification of BRAF, in which case 100% of tumor cells may have variants in both BRAF and NRAS.
  • Double-stranded cell-free DNA is isolated from the plasma of a patient.
  • the cell-free DNA fragments are tagged using 16 different bubble-containing adaptors, each of which comprises a distinctive barcode.
  • the bubble-containing adaptors are attached to both ends of each cell-free DNA fragment by ligation. After ligation, each of the cell-free DNA fragment can be distinctly identified by the sequence of the distinct barcodes and two 20 bp endogenous sequences at each end of the cell-free DNA fragment.
  • the tagged cell-free DNA fragments are amplified by PCR.
  • the amplified fragments are enriched using beads comprising oligonucleotide probes that specifically bind to a group of cancer-associated genes. Therefore, cell-free DNA fragments from the group of cancer-associated genes are selectively enriched.
  • Sequencing adaptors each of which comprises a sequencing primer binding site, a sample barcode, and a cell-flow sequence, are attached to the enriched DNA molecules.
  • the resulting molecules are amplified by PCR.
  • each bubble-containing adaptor comprises a non-complementary portion (e.g., the bubble)
  • the sequence of the one strand of the bubble-containing adaptor is different from the sequence of the other strand (complement). Therefore, the sequence reads of amplicons derived from the Watson strand of an original cell-free DNA can be distinguished from amplicons from the Crick strand of the original cell-free DNA by the attached bubble-containing adaptor sequences.
  • sequence reads from a strand of an original cell-free DNA fragment are compared to the sequence reads from the other strand of the original cell-free DNA fragment. If a variant occurs in only the sequence reads from one strand, but not other strand, of the original cell-free DNA fragment, this variant will be identified as an error (e.g., resulted from PCR and/or amplification), rather than a true genetic variant.
  • sequence reads are grouped into families. Errors in the sequence reads are corrected.
  • the consensus sequence of each family is generated by collapsing.
  • a therapeutic intervention is determined to treat the cancer.
  • Cancers with BRAF mutants respond to treatment with vemurafenib, regorafenib, tranetinib and dabrafenib.
  • Cancers with NRAS mutants respond to treatment with trametinib.
  • Cancers with JAK2 mutants respond to treatment with ruxolitinib.
  • a therapeutic intervention including administration of trametinib and ruxolitinib is determined to be more effective against this cancer than treatment with any one of the aforementioned drugs alone.
  • the subject is treated with a combination of trametinib and ruxolitinib at a dose ratio of 5:1.
  • the cfDNA from the subject is tested again for the presence of tumor heterogeneity.
  • Results show that the ratio of the BRAF:NRAS:JAK2 is now about 4:2:1.5. This indicates that the therapeutic intervention has reduced the number of cells with the BRAF and NRAS mutants, and has halted growth of cells with JAK2 mutants.
  • a second therapeutic intervention is determined in which trametinib and ruxolitinib are determined to be effective in a dose ratio of 1:1.
  • the subject is given a course of chemotherapy at amounts at this ratio.
  • Subsequent testing shows that BRAF, NRAS and JAK2 mutants are present in cfDNA at amounts below 1%.
  • a blood sample is collected from an individual with melanoma pre-treatment and the patient is determined to have a BRAF V600E mutation at a concentration of 2.8% and no detectable NRAS mutations using cell-free DNA analysis.
  • the patient is put on an anti-BRAF therapy (dabrafenib). After 3 weeks, another blood sample is collected and tested.
  • the BRAF V600E level is determined to have dropped to 0.1%.
  • the therapy is stopped and the test repeated every 2 weeks.
  • the BRAF V600E level rises again and therapy is reinitiated when the BRAF V600E level rises to 1.5%. Therapy is again stopped when the level drops down to 0.1% again. This cycle is repeated.
  • Copy number variations in a patient sample are determined. Methods for determining can include molecular tracking and upsampling, as described above.
  • a hidden-markov model based on expected locations of origins of replication is used to remove the effect of replication origin proximity from the estimated copy number variations in the patient sample. The standard deviation of copy-number variations for each gene is subsequently reduced by 40%.
  • the replication origin proximity model is also used to infer cell-free tumor burden in the patient.
  • the level of cell-free tumor derived may be low or below the detection limit of a particular technology. This can be the case when the number of human genome equivalents of tumor derived DNA in plasma is below 1 copy per 5 mL.
  • Radiation and chemotherapies have been shown to affect rapidly dividing cells more than stable, healthy cells, hence their efficacy in treating advanced cancer patients.
  • a procedure with minimal adverse effects is administered to a patient pre-blood collection to preferentially increase the fraction of tumor-derived DNA collected.
  • a low dose of chemotherapy could be administered to the patient and a blood sample could be collected within 24 hours, 48 hours, 72 hours or less than 1 week.
  • this blood sample contains higher concentrations of cell-free tumor-derived DNA due to potentially higher rates of cell-death of cancer cells.
  • low-dose radiation therapy is applied via a whole-body radiographic instrument or locally to the affected regions instead of low-dose chemotherapy.
  • Other procedures are envisioned, including subjecting a patient to ultrasound, sound waves, exercise, stress, etc.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Engineering & Computer Science (AREA)
  • Zoology (AREA)
  • Wood Science & Technology (AREA)
  • Physics & Mathematics (AREA)
  • Genetics & Genomics (AREA)
  • Analytical Chemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Immunology (AREA)
  • General Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Microbiology (AREA)
  • Biochemistry (AREA)
  • Pathology (AREA)
  • Hospice & Palliative Care (AREA)
  • Oncology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Plant Pathology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
  • Medical Treatment And Welfare Office Work (AREA)
  • Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
US15/431,395 2014-12-31 2017-02-13 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results Pending US20170260590A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US15/431,395 US20170260590A1 (en) 2014-12-31 2017-02-13 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/000,010 US20210040565A1 (en) 2014-12-31 2020-08-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/462,906 US20210395837A1 (en) 2014-12-31 2021-08-31 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/699,968 US20240141432A9 (en) 2014-12-31 2022-03-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US201462098426P 2014-12-31 2014-12-31
US201562155763P 2015-05-01 2015-05-01
PCT/US2015/067717 WO2016109452A1 (en) 2014-12-31 2015-12-28 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US15/431,395 US20170260590A1 (en) 2014-12-31 2017-02-13 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2015/067717 Continuation WO2016109452A1 (en) 2014-12-31 2015-12-28 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US17/000,010 Continuation US20210040565A1 (en) 2014-12-31 2020-08-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results

Publications (1)

Publication Number Publication Date
US20170260590A1 true US20170260590A1 (en) 2017-09-14

Family

ID=56284976

Family Applications (4)

Application Number Title Priority Date Filing Date
US15/431,395 Pending US20170260590A1 (en) 2014-12-31 2017-02-13 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/000,010 Pending US20210040565A1 (en) 2014-12-31 2020-08-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/462,906 Pending US20210395837A1 (en) 2014-12-31 2021-08-31 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/699,968 Pending US20240141432A9 (en) 2014-12-31 2022-03-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results

Family Applications After (3)

Application Number Title Priority Date Filing Date
US17/000,010 Pending US20210040565A1 (en) 2014-12-31 2020-08-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/462,906 Pending US20210395837A1 (en) 2014-12-31 2021-08-31 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
US17/699,968 Pending US20240141432A9 (en) 2014-12-31 2022-03-21 Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results

Country Status (9)

Country Link
US (4) US20170260590A1 (de)
EP (3) EP3766986B1 (de)
JP (4) JP6783768B2 (de)
CN (2) CN107406876B (de)
AU (3) AU2015374259B2 (de)
CA (1) CA2972433A1 (de)
ES (2) ES2828279T3 (de)
GB (1) GB2552267B (de)
WO (1) WO2016109452A1 (de)

Cited By (27)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10240209B2 (en) 2015-02-10 2019-03-26 The Chinese University Of Hong Kong Detecting mutations for cancer screening
US10364467B2 (en) 2015-01-13 2019-07-30 The Chinese University Of Hong Kong Using size and number aberrations in plasma DNA for detecting cancer
WO2019204360A1 (en) * 2018-04-16 2019-10-24 Grail, Inc. Systems and methods for determining tumor fraction in cell-free nucleic acid
US10535434B2 (en) * 2017-04-28 2020-01-14 4D Path Inc. Apparatus, systems, and methods for rapid cancer detection
WO2020064390A1 (en) * 2018-09-28 2020-04-02 Siemens Healthcare Gmbh A noise measure for copy number analysis on targeted panel sequencing data
US10633713B2 (en) 2017-01-25 2020-04-28 The Chinese University Of Hong Kong Diagnostic applications using nucleic acid fragments
US10636512B2 (en) 2017-07-14 2020-04-28 Cofactor Genomics, Inc. Immuno-oncology applications using next generation sequencing
US20200227136A1 (en) * 2019-01-10 2020-07-16 Travera LLC Identifying cancer therapies
WO2020150656A1 (en) 2017-08-07 2020-07-23 The Johns Hopkins University Methods for assessing and treating cancer
US10741270B2 (en) 2012-03-08 2020-08-11 The Chinese University Of Hong Kong Size-based analysis of cell-free tumor DNA for classifying level of cancer
US11062791B2 (en) 2016-09-30 2021-07-13 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
US11062789B2 (en) 2014-07-18 2021-07-13 The Chinese University Of Hong Kong Methylation pattern analysis of tissues in a DNA mixture
CN113227468A (zh) * 2018-11-21 2021-08-06 卡里乌斯公司 感染性疾病的检测和预测
US11180803B2 (en) 2011-04-15 2021-11-23 The Johns Hopkins University Safe sequencing system
WO2021243303A1 (en) * 2020-05-29 2021-12-02 Children's National Medical Center Genetic diagnostic tool for facioscapulohumeral muscular dystrophy (fshd)
US11211147B2 (en) 2020-02-18 2021-12-28 Tempus Labs, Inc. Estimation of circulating tumor fraction using off-target reads of targeted-panel sequencing
US11211144B2 (en) 2020-02-18 2021-12-28 Tempus Labs, Inc. Methods and systems for refining copy number variation in a liquid biopsy assay
US11286531B2 (en) 2015-08-11 2022-03-29 The Johns Hopkins University Assaying ovarian cyst fluid
CN114267445A (zh) * 2021-12-23 2022-04-01 山东众阳健康科技集团有限公司 一种诊断一致性检验方法、系统、设备及介质
US11410750B2 (en) 2018-09-27 2022-08-09 Grail, Llc Methylation markers and targeted methylation probe panel
US11435339B2 (en) 2016-11-30 2022-09-06 The Chinese University Of Hong Kong Analysis of cell-free DNA in urine
US11475981B2 (en) 2020-02-18 2022-10-18 Tempus Labs, Inc. Methods and systems for dynamic variant thresholding in a liquid biopsy assay
US11525163B2 (en) 2012-10-29 2022-12-13 The Johns Hopkins University Papanicolaou test for ovarian and endometrial cancers
US11581063B2 (en) 2015-07-23 2023-02-14 The Chinese University Of Hong Kong Analysis of fragmentation patterns of cell-free DNA
US11643693B2 (en) 2019-01-31 2023-05-09 Guardant Health, Inc. Compositions and methods for isolating cell-free DNA
WO2023183750A1 (en) * 2022-03-23 2023-09-28 Foundation Medicine, Inc. Methods and systems for determining tumor heterogeneity
US12024750B2 (en) 2018-04-02 2024-07-02 Grail, Llc Methylation markers and targeted methylation probe panel

Families Citing this family (30)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3561072A1 (de) 2012-12-10 2019-10-30 Resolution Bioscience, Inc. Verfahren zur gezielten genomanalyse
CN108474040B (zh) 2015-10-09 2023-05-16 夸登特健康公司 使用无细胞dna的基于群体的治疗推荐
ES2856598T3 (es) 2015-11-11 2021-09-27 Resolution Bioscience Inc Construcción de alta eficiencia de bibliotecas de ADN
CN108603228B (zh) 2015-12-17 2023-09-01 夸登特健康公司 通过分析无细胞dna确定肿瘤基因拷贝数的方法
RU2019108294A (ru) 2016-08-25 2020-09-25 Резолюшн Байосайенс, Инк. Способы обнаружения изменений количества геномных копий в образцах днк
US9850523B1 (en) 2016-09-30 2017-12-26 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
CA3049926A1 (en) 2017-01-17 2018-07-26 Heparegenix Gmbh Protein kinase inhibitors for promoting liver regeneration or reducing or preventing hepatocyte death
KR102083501B1 (ko) * 2017-02-09 2020-03-02 사회복지법인 삼성생명공익재단 종양 치료를 위한 표적 유전자 판별 방법
CN118711654A (zh) * 2017-05-16 2024-09-27 夸登特健康公司 无细胞dna的体细胞来源或种系来源的鉴定
JP7072825B2 (ja) * 2017-09-13 2022-05-23 三菱電機ソフトウエア株式会社 コピー数計測装置、コピー数計測プログラムおよびコピー数計測方法
EP3685386A1 (de) * 2017-09-20 2020-07-29 Guardant Health, Inc. Verfahren und systeme zur differenzierung von somatischen und keimbahnvarianten
WO2019090156A1 (en) * 2017-11-03 2019-05-09 Guardant Health, Inc. Normalizing tumor mutation burden
US11961589B2 (en) 2017-11-28 2024-04-16 Grail, Llc Models for targeted sequencing
CA3067229A1 (en) * 2017-12-01 2019-06-06 Illumina, Inc. Methods and systems for determining somatic mutation clonality
US11597967B2 (en) 2017-12-01 2023-03-07 Personal Genome Diagnostics Inc. Process for microsatellite instability detection
CN108491689B (zh) * 2018-02-01 2019-07-09 杭州纽安津生物科技有限公司 基于转录组的肿瘤新抗原鉴定方法
CA3090951C (en) * 2018-02-12 2023-10-17 F.Hoffmann-La Roche Ag Method of predicting response to therapy by assessing tumor genetic heterogeneity
KR20210009299A (ko) * 2018-02-27 2021-01-26 코넬 유니버시티 게놈-와이드 통합을 통한 순환 종양 dna의 초민감 검출
CA3092352A1 (en) * 2018-02-27 2019-09-06 Cornell University Systems and methods for detection of residual disease
WO2019195769A1 (en) * 2018-04-06 2019-10-10 The Brigham And Women's Hospital, Inc. Methods of diagnosing and treating aggressive cutaneous t-cell lymphomas
US20240254547A1 (en) * 2018-04-27 2024-08-01 Kao Corporation Highly accurate sequencing method
JP7274504B2 (ja) * 2018-05-08 2023-05-16 エフ. ホフマン-ラ ロシュ アーゲー 多様度指数を確立することで腫瘍バリアント多様度を評価することによるがん予後診断の方法
KR20210038577A (ko) * 2018-07-23 2021-04-07 가던트 헬쓰, 인크. 종양 분율 및 커버리지에 의해 종양 돌연변이 부담을 조정하기 위한 방법 및 시스템
KR20210089240A (ko) * 2018-11-13 2021-07-15 미리어드 제네틱스, 인크. 체세포성 돌연변이를 위한 방법 및 시스템 그리고 그 용도
CN109712671B (zh) * 2018-12-20 2020-06-26 北京优迅医学检验实验室有限公司 基于ctDNA的基因检测装置、存储介质及计算机系统
CN110428905B (zh) * 2019-07-02 2022-03-29 江南大学附属医院 一种肿瘤生长趋势预测方法
CN110895963A (zh) * 2019-10-31 2020-03-20 深圳兰丁医学检验实验室 一种基于人工智能的细胞dna定量检测系统
JP7365656B1 (ja) 2022-02-18 2023-10-20 Dic株式会社 硫化オレフィンの製造方法
CN116403644B (zh) * 2023-03-03 2023-12-05 深圳吉因加信息科技有限公司 一种用于癌症风险预测的方法及装置
CN117219162B (zh) * 2023-09-12 2024-07-02 四川大学 针对肿瘤组织str图谱进行身源鉴定的证据强度评估方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7385605B2 (en) * 2003-12-04 2008-06-10 International Business Machines Corporation Computer display system for dynamically modifying stacked area line graphs to change the order or presence of a set of stacked areas in the graph respectively representative of the proportions contributed to a total by each of a set of time dependent variables
US20060199189A1 (en) * 2005-03-07 2006-09-07 Bradford Sherry A Low-dose, sequenced, individualized chemotherapy dosing method
NZ597655A (en) * 2009-07-08 2013-05-31 Worldwide Innovative Network Method for predicting efficacy of drugs in a patient
AU2011291599B2 (en) * 2010-08-18 2015-09-10 Caris Life Sciences Switzerland Holdings Gmbh Circulating biomarkers for disease
US9476095B2 (en) 2011-04-15 2016-10-25 The Johns Hopkins University Safe sequencing system
CN103003447B (zh) * 2011-07-26 2020-08-25 维里纳塔健康公司 用于确定样品中存在或不存在不同非整倍性的方法
US20130110407A1 (en) * 2011-09-16 2013-05-02 Complete Genomics, Inc. Determining variants in genome of a heterogeneous sample
SG190466A1 (en) * 2011-11-18 2013-06-28 Agency Science Tech & Res Methods for diagnosis and/or prognosis of ovarian cancer
DK2828218T3 (da) 2012-03-20 2020-11-02 Univ Washington Through Its Center For Commercialization Methods of lowering the error rate of massively parallel dna sequencing using duplex consensus sequencing
US11261494B2 (en) * 2012-06-21 2022-03-01 The Chinese University Of Hong Kong Method of measuring a fractional concentration of tumor DNA
CA2883901C (en) * 2012-09-04 2023-04-11 Guardant Health, Inc. Systems and methods to detect rare mutations and copy number variation
EP3421613B1 (de) * 2013-03-15 2020-08-19 The Board of Trustees of the Leland Stanford Junior University Identifikation und verwendung von zirkulierenden nukleinsäure-tumormarkern
CN105408496A (zh) 2013-03-15 2016-03-16 夸登特健康公司 检测稀有突变和拷贝数变异的系统和方法
US10435740B2 (en) * 2013-04-01 2019-10-08 University Of Florida Research Foundation, Incorporated Determination of methylation state and chromatin structure of target genetic loci
WO2014191938A1 (en) * 2013-05-31 2014-12-04 Novartis Ag Combination therapy containing a pi3k-alpha inhibitor and fgfr kinase inhibitor for treating cancer

Cited By (45)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11773440B2 (en) 2011-04-15 2023-10-03 The Johns Hopkins University Safe sequencing system
US11453913B2 (en) 2011-04-15 2022-09-27 The Johns Hopkins University Safe sequencing system
US11459611B2 (en) 2011-04-15 2022-10-04 The Johns Hopkins University Safe sequencing system
US11180803B2 (en) 2011-04-15 2021-11-23 The Johns Hopkins University Safe sequencing system
US12006544B2 (en) 2011-04-15 2024-06-11 The Johns Hopkins University Safe sequencing system
US11031100B2 (en) 2012-03-08 2021-06-08 The Chinese University Of Hong Kong Size-based sequencing analysis of cell-free tumor DNA for classifying level of cancer
US10741270B2 (en) 2012-03-08 2020-08-11 The Chinese University Of Hong Kong Size-based analysis of cell-free tumor DNA for classifying level of cancer
US11525163B2 (en) 2012-10-29 2022-12-13 The Johns Hopkins University Papanicolaou test for ovarian and endometrial cancers
US11062789B2 (en) 2014-07-18 2021-07-13 The Chinese University Of Hong Kong Methylation pattern analysis of tissues in a DNA mixture
US11984195B2 (en) 2014-07-18 2024-05-14 The Chinese University Of Hong Kong Methylation pattern analysis of tissues in a DNA mixture
US10364467B2 (en) 2015-01-13 2019-07-30 The Chinese University Of Hong Kong Using size and number aberrations in plasma DNA for detecting cancer
US10240209B2 (en) 2015-02-10 2019-03-26 The Chinese University Of Hong Kong Detecting mutations for cancer screening
US11168370B2 (en) 2015-02-10 2021-11-09 The Chinese University Of Hong Kong Detecting mutations for cancer screening
US11615865B2 (en) 2015-07-23 2023-03-28 The Chinese University Of Hong Kong Analysis of fragmentation patterns of cell-free DNA
US11605445B2 (en) 2015-07-23 2023-03-14 The Chinese University Of Hong Kong Analysis of fragmentation patterns of cell-free DNA
US11581063B2 (en) 2015-07-23 2023-02-14 The Chinese University Of Hong Kong Analysis of fragmentation patterns of cell-free DNA
US11286531B2 (en) 2015-08-11 2022-03-29 The Johns Hopkins University Assaying ovarian cyst fluid
US11062791B2 (en) 2016-09-30 2021-07-13 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
US12094573B2 (en) 2016-09-30 2024-09-17 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
US11817177B2 (en) 2016-09-30 2023-11-14 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
US12100482B2 (en) 2016-09-30 2024-09-24 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
US11817179B2 (en) 2016-09-30 2023-11-14 Guardant Health, Inc. Methods for multi-resolution analysis of cell-free nucleic acids
US11435339B2 (en) 2016-11-30 2022-09-06 The Chinese University Of Hong Kong Analysis of cell-free DNA in urine
US11479825B2 (en) 2017-01-25 2022-10-25 The Chinese University Of Hong Kong Diagnostic applications using nucleic acid fragments
US10633713B2 (en) 2017-01-25 2020-04-28 The Chinese University Of Hong Kong Diagnostic applications using nucleic acid fragments
US10535434B2 (en) * 2017-04-28 2020-01-14 4D Path Inc. Apparatus, systems, and methods for rapid cancer detection
US11901078B2 (en) 2017-04-28 2024-02-13 4D Path Inc. Apparatus, systems, and methods for rapid cancer detection
US10636512B2 (en) 2017-07-14 2020-04-28 Cofactor Genomics, Inc. Immuno-oncology applications using next generation sequencing
WO2020150656A1 (en) 2017-08-07 2020-07-23 The Johns Hopkins University Methods for assessing and treating cancer
US12024750B2 (en) 2018-04-02 2024-07-02 Grail, Llc Methylation markers and targeted methylation probe panel
WO2019204360A1 (en) * 2018-04-16 2019-10-24 Grail, Inc. Systems and methods for determining tumor fraction in cell-free nucleic acid
US11725251B2 (en) 2018-09-27 2023-08-15 Grail, Llc Methylation markers and targeted methylation probe panel
US11795513B2 (en) 2018-09-27 2023-10-24 Grail, Llc Methylation markers and targeted methylation probe panel
US11410750B2 (en) 2018-09-27 2022-08-09 Grail, Llc Methylation markers and targeted methylation probe panel
US11685958B2 (en) 2018-09-27 2023-06-27 Grail, Llc Methylation markers and targeted methylation probe panel
WO2020064390A1 (en) * 2018-09-28 2020-04-02 Siemens Healthcare Gmbh A noise measure for copy number analysis on targeted panel sequencing data
CN113227468A (zh) * 2018-11-21 2021-08-06 卡里乌斯公司 感染性疾病的检测和预测
US20200227136A1 (en) * 2019-01-10 2020-07-16 Travera LLC Identifying cancer therapies
US11643693B2 (en) 2019-01-31 2023-05-09 Guardant Health, Inc. Compositions and methods for isolating cell-free DNA
US11211147B2 (en) 2020-02-18 2021-12-28 Tempus Labs, Inc. Estimation of circulating tumor fraction using off-target reads of targeted-panel sequencing
US11211144B2 (en) 2020-02-18 2021-12-28 Tempus Labs, Inc. Methods and systems for refining copy number variation in a liquid biopsy assay
US11475981B2 (en) 2020-02-18 2022-10-18 Tempus Labs, Inc. Methods and systems for dynamic variant thresholding in a liquid biopsy assay
WO2021243303A1 (en) * 2020-05-29 2021-12-02 Children's National Medical Center Genetic diagnostic tool for facioscapulohumeral muscular dystrophy (fshd)
CN114267445A (zh) * 2021-12-23 2022-04-01 山东众阳健康科技集团有限公司 一种诊断一致性检验方法、系统、设备及介质
WO2023183750A1 (en) * 2022-03-23 2023-09-28 Foundation Medicine, Inc. Methods and systems for determining tumor heterogeneity

Also Published As

Publication number Publication date
US20210395837A1 (en) 2021-12-23
AU2023206196A1 (en) 2023-08-10
EP3240911A4 (de) 2018-11-21
EP3766986B1 (de) 2022-06-01
AU2020264326A1 (en) 2020-11-26
CN107406876A (zh) 2017-11-28
EP3240911A1 (de) 2017-11-08
EP3766986A1 (de) 2021-01-20
ES2923602T3 (es) 2022-09-28
US20220213562A1 (en) 2022-07-07
JP2024029174A (ja) 2024-03-05
AU2015374259A1 (en) 2017-08-10
WO2016109452A1 (en) 2016-07-07
GB2552267A (en) 2018-01-17
JP7458360B2 (ja) 2024-03-29
JP2018507682A (ja) 2018-03-22
US20240141432A9 (en) 2024-05-02
CN113930507A (zh) 2022-01-14
AU2020264326B2 (en) 2023-05-18
JP2020127408A (ja) 2020-08-27
EP4123032A1 (de) 2023-01-25
US20210040565A1 (en) 2021-02-11
GB2552267B (en) 2020-06-10
JP7145907B2 (ja) 2022-10-03
ES2828279T3 (es) 2021-05-25
JP2022024040A (ja) 2022-02-08
AU2015374259B2 (en) 2020-08-13
CA2972433A1 (en) 2016-07-07
EP3240911B1 (de) 2020-08-26
GB201712299D0 (en) 2017-09-13
JP6783768B2 (ja) 2020-11-11
CN107406876B (zh) 2021-09-07

Similar Documents

Publication Publication Date Title
AU2020264326B2 (en) Detection and treatment of disease exhibiting disease cell heterogeneity and systems and methods for communicating test results
Pleasance et al. Pan-cancer analysis of advanced patient tumors reveals interactions between therapy and genomic landscapes
JP2022532897A (ja) マルチラベルがん分類のためのシステムおよび方法
JP5955557B2 (ja) 膵臓腫瘍形成の根底にある経路および遺伝性の膵癌遺伝子
AU2016293025A1 (en) System and methodology for the analysis of genomic data obtained from a subject
CN113151474A (zh) 用于癌症检测的血浆dna突变分析
CN114026646A (zh) 用于评估肿瘤分数的系统和方法
US20210257055A1 (en) Estimation of circulating tumor fraction using off-target reads of targeted-panel sequencing
US20230162815A1 (en) Methods and systems for accurate genotyping of repeat polymorphisms
US20240071628A1 (en) Database for therapeutic interventions
US20240052419A1 (en) Methods and systems for detecting genetic variants
Williams et al. Tracking clonal evolution of drug resistance in ovarian cancer patients by exploiting structural variants in cfDNA
WO2023183750A1 (en) Methods and systems for determining tumor heterogeneity
WO2023183751A1 (en) Characterization of tumor heterogeneity as a prognostic biomarker

Legal Events

Date Code Title Description
AS Assignment

Owner name: GUARDANT HEALTH, INC., CALIFORNIA

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:ELTOUKHY, HELMY;TALASAZ, AMIRALI;KERMANI, BAHRAM GHAFFARZADEH;AND OTHERS;SIGNING DATES FROM 20161116 TO 20170428;REEL/FRAME:043120/0192

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: FINAL REJECTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

STPP Information on status: patent application and granting procedure in general

Free format text: SPECIAL NEW

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER

STPP Information on status: patent application and granting procedure in general

Free format text: FINAL REJECTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: FINAL REJECTION MAILED

STPP Information on status: patent application and granting procedure in general

Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION

STPP Information on status: patent application and granting procedure in general

Free format text: NON FINAL ACTION MAILED