The invention generally relates to a molecular classification of disease and particularly to molecular markers for cancer and methods of use thereof.
Cancer is a major health challenge. Nearly 560,000 people die from cancer annually in the United States alone, representing almost 23% of all deaths. Despite recent advances in molecular and imaging diagnostics, one of the most vexing aspects of cancer remains early detection. In fact, for certain types of cancer—e.g., pancreatic adenocarcinoma—detection often occurs so late as to practically preclude any good prognosis. Thus there is an urgent need for sensitive methods of detecting cancer.
Mutations in certain genes are associated with cancer in general and with specific cancer types. For example, inactivating mutations in the TP53 gene are found in approximately 50% of all solid tumors and activating mutations in the KRAS or BRAF genes are often found in colorectal cancer. It has been discovered that screening patients for mutations in certain genes can detect and classify cancer. More specifically, it has been determined that (a) screening certain genes (e.g., APC, EGFR, KRAS, PTEN, and TP53) for mutations will detect nearly 95% of all cancers, while (b) screening certain genes (e.g., AIM1, APC, CDKN2A, EGFR, FBN2, FBXW7, FLJ13479, IDH1, KRAS, PIK3CA, PIK3R1, PTEN, RB1, SMAD4, TGFBR2, TNN, and TP53) for mutations can accurately classify the cancer (e.g., as breast cancer, colon cancer, glioblastoma, pancreatic cancer, etc.).
Thus the invention generally provides methods comprising analyzing panels of genes from a sample obtained from a patient (e.g., mRNA or cDNA synthesized therefrom) and determining the mutational status of the genes in the panel, wherein the presence of a particular mutational status in particular genes in the panel indicates (a) the patient has cancer and/or (b) the patient has a particular cancer.
One aspect of the invention provides a method of detecting mutations comprising: (1) analyzing in a bodily fluid sample from a human subject a panel of genes consisting of between 5 and 5,000 genes, wherein said panel comprises at least five genes chosen from the group consisting of the genes listed in Table 1; and (2) determining whether any of the genes in Table 1 harbors a mutation.
In some embodiments the panel comprises the APC, EGFR, KRAS, PTEN, and TP53 genes. In some embodiments the panel comprises the genes listed in Table 3. In some embodiments the panel comprises the genes listed in Table 2. In some embodiments the panel comprises the genes listed in Table 1.
One aspect of the invention provides a method of detecting cancer comprising: (1) analyzing a panel of genes comprising the APC, EGFR, KRAS, PTEN, and TP53 genes in a bodily fluid sample; and (2) determining whether any of the APC, EGFR, KRAS, PTEN, or TP53 genes harbors a mutation; wherein said mutation indicates the presence of cancer.
In some embodiments the panel comprises the genes listed in Table 3. In some embodiments the panel comprises the genes listed in Table 2. In some embodiments the panel comprises the genes listed in Table 1. In some embodiments the mutation is selected from the group consisting of those listed in Table 7 and/or Table 8.
One aspect of the invention provides a method of determining the likelihood a patient has cancer c1 comprising: (1) analyzing in a fluid sample a panel of genes comprising the genes listed in Table 3; (2) detecting a mutation in at least one of said genes listed in Table 3; and (3) calculating a likelihood said patient has cancer c1 using the formula: P(c1|g1, g2, . . . , gn)=P0(c1) Πi M(gi|c1)/Σt P0(t) Πi M(gi|t); wherein the product is taken over all genes in said panel mutated in the sample (i=1, 2, . . . , n), the sum is taken over all cancer types t, M(g|c1) is the frequency of somatic mutations in gene g in cancer type c1, and P0(c1) is the a priori probability of cancer c1 given that the patient has a cancer.
In some embodiments such method further comprises calculating a likelihood said patient has a second cancer c2 using the formula: P(c2|g1, g2, . . . , gn)=P0(c2) Πi M(gi|c2)/Σt P0(t) Πi M(gi|t); wherein the product is taken over all genes in said panel mutated in the sample (i=1, 2, . . . , n), the sum is taken over all cancer types t, M(g|c2) is the frequency of somatic mutations in gene g in cancer type c2, and P0(c2) is the a priori probability of cancer type c2 given that the patient has a cancer.
Some embodiments further comprise recommending, prescribing, ordering, or performing a test for the presence of cancer c1 in said patient. In some embodiments the test for the presence of cancer c1 is recommended, prescribed, ordered, or performed if the calculated likelihood said patient has said cancer c1 is above a threshold value (e.g., 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 96%, 97%, 98%, 99%, or 100%).
In some embodiments the test for the presence of cancer c1 is recommended, prescribed, ordered, or performed if the calculated likelihood said patient has said cancer c1 is higher than the calculated likelihood said patient has cancer c2. In some embodiments the method further comprises recommending, prescribing, ordering, or performing a test for the presence of cancer c2 in said patient if said test for the presence of cancer c1 does not indicate the presence of cancer c1.
In various embodiments of the invention the said bodily fluid sample is a blood sample. In some embodiments the blood sample is a plasma sample. In some embodiments the blood sample is a serum sample.
In some embodiments detecting a mutation or determining whether a gene harbors a mutation comprises analyzing an mRNA molecule from a sample or analyzing a DNA molecule synthesized using the mRNA molecule as a template. In some embodiments detecting a mutation or determining whether a gene harbors a mutation comprises analyzing a nucleic acid from a sample by a technique chosen from resequencing, TaqMan™, microarray analysis, and FISH.
In some embodiments nucleic acids to be analyzed are derived from an extracellular vesicle. In some embodiments such extracellular vesicle is an exosome.
One aspect of the invention provides a kit comprising reagents for analyzing a panel of genes consisting of between 5 and 5,000 genes, said kit comprising reagents for detecting mutations in at least five genes selected from the group consisting of the genes listed in Table 1. In some embodiments the kit comprises reagents for detecting mutations in the APC, EGFR, KRAS, PTEN, and TP53 genes.E3. In some embodiments the kit comprises reagents for detecting mutations in the genes listed in Table 3. In some embodiments the kit comprises reagents for detecting mutations in the genes listed in Table 2. In some embodiments the kit comprises reagents for detecting mutations in the genes listed in Table 1.
Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described below. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.
Other features and advantages of the invention will be apparent from the following detailed description, and from the claims.
Mutations in certain genes are associated with cancer in general and with specific cancer types. For example, inactivating mutations in the TP53 gene are found in approximately 50% of all solid tumors and activating mutations in the KRAS or BRAF genes are often found in colorectal cancer.
The invention is based in part on the discovery that analyzing patient samples for mutations in a relatively small number of genes can (a) detect the vast majority of cancers and (b) specify in which tissue the cancer is located. More specifically, it has been determined that (a) screening certain genes (e.g., the genes listed in Table 4 below) for mutations will detect cancer (e.g., nearly 95% of all cancers), while (b) screening certain genes (e.g., the genes listed in Tables 2 & 3 below) for mutations can accurately classify the cancer (e.g., as breast cancer, colon cancer, glioblastoma, pancreatic cancer, etc.).
Thus the invention provides a method of detecting mutations comprising (1) analyzing a panel of genes consisting of between 5 and 5,000 genes in a bodily fluid sample from a human subject, wherein said panel comprises at least five genes chosen from the group consisting of the genes listed in Table 1; and (2) determining whether at least one of said five genes harbors a mutation.
In some embodiments the panel consists of between 5 and 4,500, between 5 and 4,000, between 5 and 3,500, between 5 and 3,000, between 5 and 2,500, between 5 and 2,000, between 5 and 1,500, between 5 and 1,000, between 5 and 500, between 5 and 400, between 5 and 300, between 5 and 200, between 5 and 150, between 5 and 100, between 5 and 75, or between 5 and 50 genes. In some embodiments the genes chosen from Table 1 comprise at least 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 96%, 97%, 98%, 99% or 100% of the panel.
It has been discovered that screening patient samples for mutations in the genes listed in Table 1 below will detect the vast majority of cancers. In Example 2, for instance, screening for mutations in the APC, EGFR, KRAS, PTEN and TP53 genes is shown to detect nearly 95% of cancers (
In some embodiments of this aspect the panel consists of between 5 and 4,500, between 5 and 4,000, between 5 and 3,500, between 5 and 3,000, between 5 and 2,500, between 5 and 2,000, between 5 and 1,500, between 5 and 1,000, between 5 and 500, between 5 and 400, between 5 and 300, between 5 and 200, between 5 and 150, between 5 and 100, between 5 and 75, or between 5 and 50 genes. In some embodiments the genes chosen from Table 1 comprise at least 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 96%, 97%, 98%, 99% or 100% of the panel.
It has further been discovered that one can detect a particular cancer c in a patient by screening for somatic mutations in n genes g1, g2, . . . , gn in the sample and applying the following equation:
P(c|g1,g2, . . . ,gn)=P0(c)ΠiM(gi|c)/ΣtP0(t)ΠiM(gi|t) (1)
where the product is taken over all genes mutated in the sample (i=1, 2, . . . , n) and the sum is taken over all cancer types t. See Example 1, infra. M(g|c) is the frequency of somatic mutations in gene g in cancer type c. See, e.g.,
Note that the reference values discussed herein (e.g., frequency of mutations in any particular gene in any particular cancer type and probability of a particular cancer type given the patient has cancer) may be tailored to suit the needs of the skilled artisan. For example, mutation frequencies and the relative prevalence of particular cancer types may vary between, e.g., ethnic populations, countries, regions, etc.
Thus one aspect of the invention provides a method of determining the likelihood a patient has a particular cancer c1 comprising:
As used herein, the “a priori probability of cancer c given that the patient has a cancer” refers to the general incidence of the particular cancer c in the relevant cancer patient population (e.g., males or females). In other words, this is the relative proportion of all cancers in the relevant population represented by the particular cancer c. Such incidences may be gathered from various sources—e.g., yearly American Cancer Society reports on cancer incidence (as in Example 1, infra), which often give detailed breakdowns of specific cancer incidence in relevant patient subpopulations such as male vs. female, race or ethnicity, etc.
In some embodiments it is concluded that the patient has a particular cancer c1 only if the calculated likelihood said patient has said cancer c1 is above a threshold value. This threshold value may be arbitrarily chosen (e.g., 95% probability is good enough) or determined empirically (e.g., patients with a calculated probability above 80% have ended up with the particular cancer with enough frequency to validate this as a good threshold). In some embodiments said threshold value is chosen from the group consisting of 35%, 40%, 45%, 50%, 55%, 60%, 65%, 70%, 75%, 80%, 85%, 90%, 95%, 96%, 97%, 98%, and 99%.
Since some organs can develop cancers of different types (such as adenocarcinoma and squamous cell carcinoma in lung), one may calculate the probability P(o) that the cancer has developed in organ o:
P(o|g1,g2, . . . ,gn)=ΣcP(c|g1,g2, . . . ,gn) (2)
where the sum is over all cancer types c of the organ o. Using Equation (2), the probabilities are calculated for each organ o, and the organ with the highest probability is the most likely cancer site in the patient. The patient may then optionally be examined by additional diagnostic techniques to confirm cancer site. If the most likely cancer site is not confirmed, the organ with the second highest probability may then be examined and so on.
Thus one aspect of the invention provides a method of diagnosing cancer in a particular organ o1 comprising:
P(o|g1,g2, . . . ,gn)=ΣcP(c|g1,g2, . . . ,gn)
P(c1|g1,g2, . . . ,gn)=P0(c1)ΠiM(gi|c1)/ΣtP0(t)ΠiM(gi|t)
When screening a patient for cancer (e.g., early detection), it will often be desirable to calculate the probabilities of several different cancers (e.g., the most prevalent cancers in the relevant patient population or the cancers listed in Tables 3 & 4) so as to allow comparison to determine which of a plurality of cancers is the most likely. Thus another aspect of the invention provides a method of determining the likelihood a patient has a particular cancer c1 comprising:
P(c1|g1,g2, . . . ,gn)=P0(c1)ΠiM(gi|c1)/ΣtP0(t)ΠiM(gi|t)
P(c2|g1,g2, . . . ,gn)=P0(c2)ΠiM(gi|c2)/ΣtP0(t)ΠiM(gi|t)
This may be repeated and the various probabilities compared to give the desired confidence that the patient has any particular cancer. In some embodiments the method further comprises concluding the patient has c1 if P(c1) is higher than P(c2), P(c3), P(c4), . . . , P(cx), where P(c2) through P(cx) represent the calculated probabilities of each cancer (e.g., major cancers such as those listed in Tables 3 & 4) other than c1.
It will often be useful to know what particular cancer is present. Thus one aspect of the invention provides a method of diagnosing cancer comprising:
P(c1|g1,g2, . . . ,gn)=P0(c1)ΠiM(gi|c1)/ΣtP0(t)ΠiM(gi|t)
As mentioned above, screening the five genes in Table 4 can detect nearly 95% of solid tumor types and the genes in Tables 2 & 3 can classify the cancer. Thus in some embodiments the presence of a mutation in any one of the genes listed in Table 4 is used as a pan-cancer screen to determine for which patients additional analysis should be done on a panel comprising at least one of the genes listed in Table 2 or 3. In some embodiments a mutation in any one of the genes in the first panel indicates the patient has cancer and application of the second panel classifies which type.
In some circumstances somatic mutations are the most informative mutations (e.g., as in Example 1). In such cases one may determine the mutational status of the panel genes in both germline and somatic tissue to confirm that the mutation detected in the mutation screen is in fact somatic. In some embodiments this may be done with a single patient blood sample since germline mutational status may be determined from circulating blood cells while the somatic mutational analysis can be done with, e.g., circulating tumor cells, exosomes derived from tumor cells, or circulating nucleic acids derived from tumor cells.
Calculating a patient's likelihood of having a particular cancer can be useful in various clinical settings. For example, if the calculated probability of the patient having a particular cancer is high enough one may diagnose the particular cancer, prescribe a treatment for the specific cancer, etc. If the patient is at particularly high risk of a specific cancer (e.g., BRCA mutation carrier), then even a lower calculated likelihood of breast or ovarian cancer might be sufficient to make a diagnosis. A high likelihood of a particular cancer may alternatively prompt the doctor to recommend, prescribe, order, or perform an additional test (e.g., biopsy, MRI, CT scan, digital rectal exam, mammography, etc.) to confirm the cancer.
Thus in aspects comprising calculating the likelihood of cancer c1, some embodiments further comprise recommending, prescribing, ordering, or performing a test to confirm the presence of cancer c1. In some embodiments the test is prescribed, ordered, recommended, or performed if the calculated likelihood exceeds some threshold value. In aspects comprising calculating the likelihood of cancer c1 and the likelihood of cancer c2, some embodiments further comprise recommending, prescribing, ordering, or performing a test to confirm the presence of cancer c1 in said patient if the calculated likelihood said patient has said cancer c1 is higher than the calculated likelihood said patient has cancer c2. In some embodiments the test is prescribed, ordered, recommended, or performed if the calculated likelihood of c1 exceeds that of c2 and also exceeds some threshold value.
As used herein, a “panel of genes” is a plurality of genes. In some embodiments the panel consists of between 2 and 500, between 3 and 500, between 4 and 500, between 5 and 500, between 6 and 500, between 7 and 500, between 8 and 500, between 9 and 500, between 10 and 500, between 11 and 500, between 12 and 500, between 13 and 500, between 14 and 500, between 15 and 500, between 16 and 500, between 17 and 500, between 18 and 500, between 19 and 500, between 20 and 500, between 25 and 500, between 30 and 500, between 35 and 500, between 40 and 500, between 45 and 500, between 50 and 500, between 55 and 500, between 60 and 500, between 65 and 500, between 70 and 500, between 75 and 500, between 80 and 500, between 85 and 500, between 90 and 500, between 95 and 500, between 100 and 500, between 2 and 400, between 2 and 350, between 2 and 300, between 2 and 250, between 2 and 200, between 2 and 150, between 2 and 100, between 2 and 90, between 2 and 80, between 2 and 70, between 2 and 60, between 2 and 50, between 2 and 45, between 2 and 40, between 2 and 35, between 2 and 30, between 2 and 25, between 2 and 20, between 2 and 19, between 2 and 18, between 2 and 17, between 2 and 16, between 2 and 15, between 2 and 14, between 2 and 13, between 2 and 12, between 2 and 11, between 2 and 10, between 2 and 9, between 2 and 8, between 2 and 7, between 2 and 6, between 2 and 5, between 2 and 4, or between 2 and 3 genes and comprises at least one of the gene listed in Table 1 or a subset of the genes in Table 1. As used in the context of ranges, “between” includes the end of the range (i.e., “between 2 and 500” includes both 2 and 500).
In some embodiments of the invention the panel comprises genes listed in Table 1 below:
In some embodiments the panel comprises subsets (e.g., at least 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 40 or more) of the genes in Table 1. In some embodiments the panel comprises APC, EGFR, KRAS, PTEN, and TP53. In some embodiments the panel comprises AIM1, APC, CDKN2A, EGFR, FBN2, FBXW7, FLJ13479, IDH1, KRAS, PIK3CA, PIK3R1, PTEN, RB1, SMAD4, TGFBR2, TNN, and TP53. In some embodiments the panel comprises APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CTNNB1, EGFR, FBXW7, FGFR3, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MSH2, MSH6, NE1, NF2, PIK3CA, PRKDC, PTEN, RB1, RET, SMAD4, SMO, STK11, TAF1L, TP53, TRRAP, and VHL. In some embodiments the panel comprises the genes listed in Table 4. In some embodiments the panel comprises the genes listed in Table 3. In some embodiments the panel comprises the genes listed in Table 1.
Mutations useful in the methods of the invention include missense mutations, deletions, insertions, frameshifts, copy number variations, and loss of heterozygosity. Deleterious mutations (i.e., mutations that reduce or abolish gene and/or protein function) are particularly relevant in the context of tumor suppressors (e.g., APC, TP53, PTEN). Activating mutations (i.e., mutations that increase gene and/or protein function) are particularly relevant in the context of oncogenes (e.g., KRAS, EGFR). Those skilled in the art are familiar with various deleterious and activating mutations for the genes listed in Tables 1, 3, and 4 (e.g., codons 12 and 13 in KRAS). Skilled artisans are also familiar with various techniques for determining whether a particular mutation is in fact deleterious or activating. For example, frameshift mutations resulting in early truncation of a tumor suppressor gene are generally expected to be deleterious. Table 7 includes mutations found in some of the genes listed in Table 1. Those skilled in the art are familiar with various resources and databases cataloguing mutations in the genes listed in Table 1. For example, the COSMIC [Catalogue of Somatic Mutations in Cancer] database currently contains over 26,000 entries for these genes. Those skilled in the art will be able to use these entries in the methods of the invention for detecting and classifying cancer.
As used herein, determining the “mutational status” of a gene means determining at least one of the following: (a) whether the gene (or any of its products) harbors a sequence mutation (including point mutations, deletions, insertions, copy number variants, etc.), (b) the prevalence of such mutations in a sample, or (c) whether such a sequence mutation is activating or inactivating. Thus a particular mutational status includes, but is not limited to, the presence or absence of a mutation, a relatively high or relatively low prevalence of a mutation, an inactivating mutation, an activating mutation, etc. In some embodiments determining the mutational status of a gene comprises assaying some marker whose status itself is correlated with the mutational status of the gene of interest. Determining the mutational status of a panel of genes means determining the mutational status of each gene in the panel.
Mutational status of a gene may be determined by any of several techniques familiar to those skilled in the art. Exemplary techniques include resequencing (either of selected regions of the gene or of the entire gene), allele-specific amplification (e.g., TaqMan™ using mutant allele-specific probes), microarray analysis (e.g., arrays for CNV or arrays containing mutant allele-specific probes), etc. In some embodiments of the invention the method comprises physically amplifying and/or isolating nucleic acid of a panel of genes from a sample obtained from a patient. As used herein, “amplifying a nucleic acid” and “isolating nucleic acid” have their conventional meanings in the art. Thus in some embodiments the method further comprises amplifying nucleic acid of a panel of genes (e.g., comprising the genes listed in Table 3) from a sample obtained from a patient, determining the mutational status of each gene in the panel, and calculating the likelihood of a particular cancer as discussed above and below.
“Sample” as used herein refers to any biological specimen, including any tissue or fluid, that can be obtained from, or derived from a specimen obtained from, a human subject. Such samples include but are not limited to healthy or tumor tissue, bodily fluids (e.g., blood), waste matter (e.g., urine, stool), etc. “Bodily fluid sample” as used herein means any fluid that can be extracted or collected from a human body. In some embodiments of each aspect of the invention the bodily fluid sample is blood or a blood derivative. Examples of blood derivatives include, but are not limited to, plasma and serum. In some embodiments the bodily fluid sample is urine, stool, pleural effusion, lacrimal effusion, saliva, sputum, etc. As used herein, “analyzing genes in a sample” refers to analyzing nucleic acids corresponding to those genes in a sample or any substance derived from that sample. For example, analyzing the APC, EGFR, KRAS, PTEN and TP53 genes in blood includes analyzing PCR™ amplified portions of these genes in a patient blood sample (including plasma or serum), or in DNA or RNA isolated (i.e., derived) from such a sample. In some embodiments such a nucleic acid is chosen from the group consisting of genomic DNA (including PCR™ amplified copies of genomic DNA), mRNA, cDNA, and a portion of any of these.
The cancer screening and classification methods of the inventions will often involve analyzing nucleic acids from bodily fluids since these are often the least invasive samples to obtain from patients. For example, the method of the invention may involve isolating nucleic acids from circulating tumor cells from the blood. This may involve capturing circulating tumor cells (e.g., using tumor-specific capture antibodies) and subsequent analysis of the DNA or RNA contained in the cell. Alternatively, the methods of the invention may isolate and analyze nucleic acids that float freely in the bodily fluid. As discussed in more detail below, the methods of the invention may also isolate nucleic acids from extracellular vesicles found in the bodily fluid sample.
Mutations in some of genes are associated with particular cancer types. As used herein, “cancer type” and “type of cancer” mean a cancer in or originating from a particular tissue or organ and/or a cancer with a particular molecular or clinical feature. Often, the specificity of the “cancer type” varies with the application, including tissue type (e.g., squamous versus cuboidal), organ type (e.g., breast versus lung), and clinical subtype (e.g., triple-negative breast cancer). Thus another aspect of the invention provides a method of classifying cancer comprising isolating nucleic acids corresponding to a panel of genes from a sample obtained from a patient and determining the mutational status of each such nucleic acid, wherein a particular mutational status in particular genes in the panel indicates the patient has a particular cancer. Those skilled in the art will appreciate that methods according to this aspect may simultaneously detect and classify cancer. In some embodiments the panel comprises the AIM1, APC, CDKN2A, EGFR, FBN2, FBXW7, FLJ13479, IDH1, KRAS, PIK3CA, PIK3R1, PTEN, RB1, SMAD4, TGFBR2, TNN, and TP53 genes or a subset (e.g., at least 3, 4, 5, 6, 7, 8, 9, 10 or 15 or more) thereof. In other embodiments the panel comprises the APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CTNNB1, EGFR, FBXW7, FGFR3, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MSH2, MSH6, NE1, NF2, PIK3CA, PRKDC, PTEN, RB1, RET, SMAD4, SMO, STK11, TAF1L, TP53, TRRAP, and VHL genes or a subset (e.g., at least 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, or 30 or more) thereof. In still other embodiments the panel comprises the AIM1, APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CD95, CTNNB1, EGFR, FBN2, FBXW7, FLJ13479, FGFR3, IDH1, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MLH2, MSH1, MSH2, NE1, NF2, PIK3CA, PIK3R1, PRKDC, PTEN, PMS1, PMS2, RB1, RET, SMAD4, SMO, STK11, TAF1L, TGFBR2, TNN, TP53, TRRAP, and VHL genes or a subset (e.g., at least 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 40 or more) thereof.
As used herein, “classifying a cancer” and “cancer classification” refer to determining one or more clinically-relevant features of a cancer. Thus “classifying a cancer” includes, but is not limited to: (i) determining the tissue type or organ of origin of the cancer (e.g., cancer type); (ii) determining clinical subtype of cancer (e.g., EGFR amplified); (iii) evaluating metastatic potential, potential to metastasize to specific organs, risk of recurrence, and/or course of the tumor; (iv) evaluating tumor stage; (v) determining patient prognosis in the absence of treatment of the cancer; (vi) determining prognosis of patient response (e.g., tumor shrinkage or progression-free survival) to treatment (e.g., chemotherapy, radiation therapy, surgery to excise tumor, etc.); (vii) diagnosis of actual patient response to current and/or past treatment; (viii) determining a preferred course of treatment for the patient; (ix) prognosis for patient relapse after treatment (either treatment in general or some particular treatment); (x) prognosis of patient life expectancy (e.g., prognosis for overall survival), etc. The methods of the invention are particularly suited to determining tumor origin.
The cancer screening and cancer classification aspects of the invention may also be combined to provide a method for diagnosing specific cancer types. This will often involve screening a patient for the presence of cancer generally and, if it is present, classifying the cancer. Thus this aspect of the invention provides a method of diagnosing cancer comprising (1) isolating nucleic acids corresponding to a first panel of genes from a sample obtained from a patient; and (2) determining the mutational status of each nucleic acid corresponding to a gene in the first panel, wherein a particular mutational status in particular genes in the first panel indicates the patient has cancer; (3) isolating nucleic acids corresponding to a second panel of genes from the sample; (4) determining the mutational status of each nucleic acid corresponding to a gene in the second panel, wherein a particular mutational status in particular genes in the second panel indicates the patient has a particular cancer type. As used herein, “cancer type” refers to tissue, tissue type or organ of origin for a cancer.
In some embodiments the isolating steps (1) and (3) are performed sequentially. This allows for a relatively less expensive, quicker initial assessment of the general presence of cancer which can, if necessary, be followed with further analysis of more genes to determine cancer type. Alternatively, in other embodiments the isolating steps (1) and (3) are done at the same time—i.e., they are in essence collapsed into a single step that isolates and/or analyzes nucleic acids from both panels simultaneously. Isolation and analysis may be performed on the same patient sample or on different samples.
In some embodiments the first panel comprises the APC, EGFR, KRAS, PTEN, and TP53 genes and the second panel comprises the AIM1, APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CD95, CTNNB1, EGFR, FBN2, FBXW7, FLJ13479, FGFR3, IDH1, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MLH2, MSH1, MSH2, NE1, NF2, PIK3CA, PIK3R1, PRKDC, PTEN, PMS1, PMS2, RB1, RET, SMAD4, SMO, STK11, TAF1L, TGFBR2, TNN, TP53, TRRAP, and VHL genes or subsets thereof (e.g., AIM1, APC, CDKN2A, EGFR, FBN2, FBXW7, FLJ13479, IDH1, KRAS, PIK3CA, PIK3R1, PTEN, RB1, SMAD4, TGFBR2, TNN, and TP53 or APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CTNNB1, EGFR, FBXW7, FGFR3, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MSH2, MSH6, NE1, NF2, PIK3CA, PRKDC, PTEN, RB1, RET, SMAD4, SMO, STK11, TAF1L, TP53, TRRAP, and VHL).
Knowing that a patient has cancer can be valuable in various clinical settings beyond diagnosis. Thus other aspects of the invention provide methods of detecting cancer in a patient identified as being at heightened risk of having or developing cancer, methods of monitoring cancer therapy (e.g., for recurrence or progression), methods of determining whether a patient is a candidate for biopsy or other further testing, methods of determining drug response, etc. These methods will generally comprise isolating nucleic acids corresponding to a panel of genes from a patient sample and determining the mutational status of each such nucleic acid, wherein a particular mutational status in particular genes in the panel will indicate some particular clinical feature (e.g., desirability of biopsy, desirability of a particular treatment, etc.). For example, a panel of genes comprising KRAS may be assayed to determine that a patient has colon cancer, with knowledge of an activating mutation in KRAS further indicating a decreased likelihood of response to anti-EGFR therapy.
Thus one aspect of the invention provides a method of screening for cancer in a patient comprising identifying a patient at risk of having, or in need of screening for, cancer and determining the mutational status of a panel of genes in a sample obtained from the patient, wherein a particular mutational status in the sample indicates the presence of cancer. Patients may be identified as at risk of having, or in need of screening for, cancer in a variety of ways and based on numerous clinical and/or molecular characteristics. One class of patients at risk of having cancer and in need of screening is those patients known to carry a germline deleterious mutation in a tumor suppressor gene. Examples include, but are not limited to, BRCA1 (breast or ovarian), BRCA2 (breast or ovarian), PTEN (glioma), p16 (melanoma), MLH1 (colorectal), MSH6 (colorectal), APC (colorectal), MYH (colorectal), etc. In such patients, cancer-type specificity is often less crucial since, for example, a BRCA1-mutant patient whose mutational status in a panel of predictive genes (e.g., APC, EGFR, KRAS, PTEN, and TP53) indicates cancer would be expected have breast or ovarian cancer rather than some other type of cancer. The relatively non-invasive nature of serum detection (i.e., simple blood draw) makes such widespread screening attractive and practical.
Thus in some embodiments the invention provides a method of detecting cancer comprising identifying a patient having a mutation in a gene selected from the group consisting of BRCA1, BRCA2, PTEN, p16, MLH1, MSH6, APC, and MYH; and determining the mutational status of a panel of genes in a sample obtained from the patient; wherein a particular mutational status indicates the presence of cancer. In some such embodiments the method further comprises additional tests to determine/confirm which type of cancer is present.
Another aspect of the invention provides a method of detecting recurrence in a cancer patient comprising determining the mutational status of a panel of genes in a sample obtained from the patient, wherein a particular mutational status indicates recurrence. Because it is difficult to remove or kill all cancerous cells, one of the main challenges in cancer treatment is making sure a cancer removed by surgery and/or treated with drugs has not returned. Thus this aspect of the invention is particularly useful in monitoring cancer patients following treatment. Much like the at-risk patients discussed above, cancer-type specificity is often not crucial: If a lung cancer patient is found to have a particular mutational status in his serum several months or years after treatment, then the new cancer is likely to be a return of the former lung cancer. As above, in some embodiments further testing (e.g., imaging) to confirm the type of cancer or to characterize the cancer (e.g., stage) is encompassed by the invention. In some embodiments mutational status is measured soon after treatment (e.g., to determine a post-treatment baseline) and then monitored at regular intervals there after in order to catch any significant change (e.g., from this baseline).
Yet another promising way in which the invention may be used clinically is to identify patients who need further testing to confirm the existence, location, and/or character of a cancer. Biopsies, for example, are generally quite invasive, involving substantial discomfort and risk (e.g., infection). Imaging tests (e.g., MRI, CT scan, etc.) are generally less invasive, but are very expensive and some a priori idea of the location of a tumor is generally needed. By indicating which patients are likely to have cancer in a particular organ or tissue, the methods of the present invention may be used to identify patients who are good candidates for biopsy or imaging. For example, the invention provides a method of diagnosing cancer comprising isolating nucleic acids corresponding to a panel of genes from a sample obtained from a patient; determining the mutational status of each such nucleic acid, wherein a particular mutational status in particular genes in the first panel indicates the patient has a particular cancer; and recommending, prescribing or performing further testing to confirm the presence, location or character of the cancer. In some embodiments the further testing comprises a biopsy or an imaging test. In some embodiments, especially if the genetic screen indicates cancer in a large organ like the lung, further testing may involve an imaging test to better pinpoint the location of any mass and then biopsy to further analyze the mass (e.g., to confirm malignancy). In the case of patients already identified as at-risk for particular cancers (e.g., BRCA mutation carriers), a simple pan-cancer screen according to the present invention may give the information necessary to prompt further testing of the at-risk area (e.g., breasts or ovaries).
Nucleic acids (e.g., mRNA) for analysis according to the present invention may come from any suitable source, especially those likely to be enriched for tumor nucleic acids. One example may be tumor tissue itself (e.g., unknown metastasis for which origin is to be determined). In another example, the blood (or serum or plasma) of a patient may be treated to isolate mRNA or DNA for mutation analysis since such body fluids carry circulating mRNA and DNA. This nucleic acid may come from circulating tumor cells or it may be free circulating nucleic acid. Techniques for isolating and analyzing nucleic acids from blood and blood derivatives are known to those skilled in the art. See, e.g., U.S. Pat. No. 7,442,507. Thus in some embodiments of the invention the sample is a bodily fluid (e.g., blood, pleural fluid, urine, etc.). In some embodiments the bodily fluid is blood. In some embodiments the sample is a blood derivative such as serum or plasma.
An additional source of nucleic acids is small extracellular vesicles, including exosomes, which are abundant in the blood (and serum and plasma) of cancer patients due to increased production by tumor cells. This is especially true of epithelial cancers (e.g., those of the lung, colon, breast, prostate, ovaries, endometrium, etc). Exosomes carry important biomolecules on their surface (e.g., protein) and within their interior (e.g., mRNA). Because exosomes are often derived from tumor cells, the biomolecules they carry can provide valuable information regarding the tumor cells from which they are derived. Thus, circulating exosomes, by generally yielding a relatively high concentration of tumor-derived mRNA, can provide an enriched snapshot or non-invasive “virtual biopsy” of tumor cells. This is especially helpful in general cancer screening, where minimal invasiveness is particularly advantageous. mRNA from exosomes may be isolated and analyzed to determine the mutational status of genes.
Thus in some embodiments of the invention nucleic acids are isolated from exosomes obtained from a patient blood (or blood derivative) sample. Several techniques for isolating nucleic acids from exosomes and for isolating exosomes themselves are known in the art. See, e.g., U.S. Pat. No. 7,198,923. Examples include differential centrifugation, immunoseparation, bead-assisted centrifugation, fluorescence-assisted cell sorting (FACS), affinity chromatography, etc. At times it will be desirable to differentiate tumor-derived exosomes from exosomes derived from some other cell, especially since normal immune cells in the blood release exosomes. This can be done, e.g., by FACS or immunocentrifugation using a surface marker specific for cancer or a marker specific for non-immune cells (e.g., epithelial membrane antigen [EMA] or EpCAM).
Other information may be combined with mutational status in some aspects of the invention. For example, expression levels of certain genes often differ between cancer and non-cancer and among different cancer types and subtypes. Thus some embodiments provide methods as described below further comprising determining the expression level of a gene, wherein a particular mutational status and a high expression level indicate cancer, a particular cancer type, etc. Examples of such genes whose expression level is often informative include, but are not limited to, EGFR, HER2, PSA, CA125, CEA, etc. Determining the expression level of a gene can include determining the amount of mRNA and/or protein products of the gene. In some embodiments the level (including presence, absence, or qualitative amount) of a marker is used not so much to indicate cancer or cancer type, but instead simply to indicate tissue or organ type from which the nucleic acid (e.g., by way of an exosome) is derived. Examples include EpCAM, 34βE12, Ae1/3, AFP, B72.3, CA-125, Calictonin, Calretinin, CAM5.2, CD10, CD15, CD56, CEA, Chromogranin, CK19, CK5/6, cytokeratin 20, cytokeratin 7, EMA, GCDFP-15, HBME-1, HepPar1, HER2, Leu, Leu7, M1, Mesothelin, Mucicarmine, NCAM, PSA, PSAP, PSMA, RCC, Synaptophysin, Thyroglobulin, UroplakinIII, Villin, Vimentin, etc.
In some embodiments the panel of tissue markers comprises two or more markers shown in
In further embodiments the status of individual markers in the panel is tested in a certain order in order to narrow down which specific cancer type is present. One example is illustrated in
If CK7 is absent and CK20 is present as in
If CK7 is present and CK20 is absent as in
If CK7 and CK20 are both present as in
As mentioned above, some embodiments of the invention involve mutational analysis combined with more traditional diagnostic techniques. For example, physical examination (e.g., digital rectal exam for prostate cancer), imaging (e.g., mammography), and/or biopsy may be used to confirm a diagnosis indicated by mutational analysis according to the invention. Alternatively, such techniques may be combined with mutational analysis (and optionally exosome surface marker analysis) to yield a more comprehensive diagnosis. As an illustrative example, a mutational screen may indicate the presence of cancer in a patient and exosomes may be found to be CK7+/CK20− and have the marker CK5/6 associated with them. One may not be able to conclusively determine based solely on this information whether the cancer is squamous cell carcinoma of the lung, cervical cancer, or mesothelioma at some unknown organ (see
As used herein in the context of biomarkers and their expression, the “level” of something in a sample has its conventional meaning in the art. Determining a “level” herein includes quantitative determinations—e.g., mg/mL, fold change, etc. Determining a “level” herein also includes qualitative determinations—e.g., determining the presence or absence of a marker or determining whether the level of the marker is “high,” “low” or even “present” relative to some index value.
In one embodiment, in determining the level of expression in accordance with the present invention the amount of expression is measured within one or more samples and compared to some index value. The index value may represent the average expression level of a marker in a plurality of training patients (e.g., both diseased and healthy patients). For example, a “cancer index value” can be generated from a plurality of training patients characterized as having cancer. A “cancer-free index value” can be generated from a plurality of training patients defined as not having cancer. Thus, a cancer index value of expression may represent the average level of expression in patients having cancer, whereas a cancer-free index value of expression may represent the average level of expression in patients not having cancer. Thus, when the level of expression is more similar to the cancer index value than to the cancer-free index value, then it can be concluded that the patient has or is likely to have cancer. On the other hand, if the level of expression is more similar to the cancer-free index value than to the cancer index value, then it can be concluded that the patient does not have or has no increased likelihood of having cancer.
The results of these and any other analyses according to the invention will often be communicated to physicians, genetic counselors and/or patients (or other interested parties such as researchers) in a transmittable form that can be communicated or transmitted to any of the above parties. Such a form can vary and can be tangible or intangible. The results can be embodied in descriptive statements, diagrams, photographs, charts, images or any other visual forms. For example, graphs showing mutational status information for various genes can be used in explaining the results. Diagrams showing such information for additional target gene(s) are also useful in indicating some testing results. The statements and visual forms can be recorded on a tangible medium such as papers, computer readable media such as floppy disks, compact disks, etc., or on an intangible medium, e.g., an electronic medium in the form of email or website on internet or intranet. In addition, results can also be recorded in a sound form and transmitted through any suitable medium, e.g., analog or digital cable lines, fiber optic cables, etc., via telephone, facsimile, wireless mobile phone, internet phone and the like.
Thus, the information and data on a test result can be produced anywhere in the world and transmitted to a different location. As an illustrative example, when an assay is conducted outside the United States, the information and data on a test result may be generated, cast in a transmittable form as described above, and then imported into the United States. Accordingly, the present invention also encompasses a method for producing a transmittable form of information on at least mutational status for a panel of genes for at least one patient sample. The method comprises the steps of (1) determining mutational status as described above according to methods of the present invention; and (2) embodying the result of the determining step in a transmittable form. The transmittable form is the product of such a method. Thus the processing of physical samples may be temporally and physically separated from their analysis in the methods of the invention. Indeed, mutational status may be determined in a blood sample for some other purpose and the stored mutational data from an earlier assay may be applied to the methods of the invention in diagnosing cancer.
Techniques for analyzing mutational status or expression (indeed any data obtained according to the invention) will often be implemented using hardware, software or a combination thereof in one or more computer systems or other processing systems capable of effectuating such analysis. The computer-based analysis function can be implemented in any suitable language and/or browsers. For example, it may be implemented with C language and preferably using object-oriented high-level programming languages such as Visual Basic, SmallTalk, C++, and the like. The application can be written to suit environments such as the Microsoft Windows™ environment including Windows™ 98, Windows™ 2000, Windows™ NT, and the like. In addition, the application can also be written for the MacIntosh™, SUN™, UNIX or LINUX environment. In addition, the functional steps can also be implemented using a universal or platform-independent programming language. Examples of such multi-platform programming languages include, but are not limited to, hypertext markup language (HTML), JAVA™, JavaScript™, Flash programming language, common gateway interface/structured query language (CGI/SQL), practical extraction report language (PERL), AppleScript™ and other system script languages, programming language/structured query language (PL/SQL), and the like. Java™- or JavaScript™-enabled browsers such as HotJava™, Microsoft™ Explorer™, or Netscape™ can be used. When active content web pages are used, they may include Java™ applets or ActiveX™ controls or other active content technologies.
The analysis function can also be embodied in computer program products and used in the systems described above or other computer- or internet-based systems. Accordingly, another aspect of the present invention relates to a computer program product comprising a computer-usable medium having computer-readable program codes or instructions embodied thereon for enabling a processor to carry out gene status analysis. These computer program instructions may be loaded onto a computer or other programmable apparatus to produce a machine, such that the instructions which execute on the computer or other programmable apparatus create means for implementing the functions or steps described above. These computer program instructions may also be stored in a computer-readable memory or medium that can direct a computer or other programmable apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory or medium produce an article of manufacture including instruction means which implement the analysis. The computer program instructions may also be loaded onto a computer or other programmable apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions or steps described above.
Thus in some embodiments the invention provides a method comprising: accessing mutational status information derived from a patient sample and stored in a computer-readable medium; querying this information to determine whether the patient has a particular mutational status for a panel of genes; calculating the likelihood of the patient having a particular cancer type based on the mutational status of the panel; outputting [or displaying] the likelihood of the patient having a particular cancer type based on the mutational status of the panel. A similar computer-implemented diagnostic method may use a panel of genes to indicate likelihood of the presence of cancer generally. Yet another method may combine the pan-cancer screen and the cancer type-specific screen described above. For example, one embodiment provides a method comprising: accessing mutational status information on a first panel of genes derived from a patient sample and stored in a computer-readable medium; querying this information to determine whether the patient has a particular mutational status for the first panel; calculating the likelihood of the patient having cancer based on the mutational status of the first panel; accessing mutational status information on a second panel of genes derived from a patient sample and stored in a computer-readable medium; querying this information to determine whether the patient has a particular mutational status for the second panel; calculating the likelihood of the patient having a particular cancer based on the mutational status of the second panel; outputting [or displaying] the likelihood of the patient having a particular cancer type based on the mutational status of the second panel. One may optionally also output [or display] the likelihood of the patient having cancer generally, either before analyzing the mutational status information for the second panel or together with the output of the likelihood of the patient having a particular cancer type. Some embodiments further comprise displaying the mutational status information.
As used herein in the context of computer-implemented embodiments of the invention, “displaying” means communicating any information by any sensory means. Examples include, but are not limited to, visual displays, e.g., on a computer screen or on a sheet of paper printed at the command of the computer, and auditory displays, e.g., computer generated or recorded auditory expression of a patient's genotype.
The practice of the present invention may also employ conventional biology methods, software and systems. Computer software products of the invention typically include computer readable media having computer-executable instructions for performing the logic steps of the method of the invention. Suitable computer readable medium include floppy disk, CD-ROM/DVD/DVD-ROM, hard-disk drive, flash memory, ROM/RAM, magnetic tapes and etc. Basic computational biology methods are described in, for example, Setubal et al., I
The present invention may also make use of various computer program products and software for a variety of purposes, such as probe design, management of data, analysis, and instrument operation. See U.S. Pat. Nos. 5,593,839; 5,795,716; 5,733,729; 5,974,164; 6,066,454; 6,090,555; 6,185,561; 6,188,783; 6,223,127; 6,229,911 and 6,308,170. Additionally, the present invention may have embodiments that include methods for providing genetic information over networks such as the Internet as shown in U.S. Ser. Nos. 10/197,621 (U.S. Pub. No. 20030097222); 10/063,559 (U.S. Pub. No. 20020183936), 10/065,856 (U.S. Pub. No. 20030100995); 10/065,868 (U.S. Pub. No. 20030120432); 10/423,403 (U.S. Pub. No. 20040049354).
Another aspect of the invention provides microarrays and kits (including a microarray kit) for practicing the methods of the invention. The kit may include a carrier for its various components. The carrier can be a container or support, in the form of, e.g., bag, box, tube, rack, and is optionally compartmentalized. The carrier may define an enclosed confinement for safety purposes during shipment and storage.
Microarrays and kits (including microarray kits) of the invention may comprise reagents for determining the mutational status of a panel of genes consisting of between 5 and 5,000 genes and comprising at least one gene chosen from the group consisting of: AIM1, APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CD95, CTNNB1, EGFR, FBN2, FBXW7, FLJ13479, FGFR3, IDH1, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MLH2, MSH1, MSH2, NE1, NF2, PIK3CA, PIK3R1, PRKDC, PTEN, PMS1, PMS2, RB1, RET, SMAD4, SMO, STK11, TAF1L, TGFBR2, TNN, TP53, TRRAP, and VHL. In some embodiments the panel comprises subsets of these genes, e.g., APC, EGFR, KRAS, PTEN, and TP53; or AIM1, APC, CDKN2A, EGFR, FBN2, FBXW7, FLJ13479, IDH1, KRAS, PIK3CA, PIK3R1, PTEN, RB1, SMAD4, TGFBR2, TNN, and TP53; or APC, ATM, BRAF, BRCA1, BRCA2, CDKN2A, CTNNB1, EGFR, FBXW7, FGFR3, KIT, KRAS, HRAS, NRAS, MAP2K4, MET, MLH1, MSH2, MSH6, NE1, NF2, PIK3CA, PRKDC, PTEN, RB1, RET, SMAD4, SMO, STK11, TAF1L, TP53, TRRAP, and VHL.
Those skilled in the art are familiar with various reagents that may be used for determining whether a particular gene harbors a mutation. For example, one may use oligonucleotide probes (e.g., probes specific for a mutant allele) and/or primers (e.g., PCR primers in RT-PCR reactions) to determine mutational status. In some embodiments the invention provides the use of such reagents for the manufacture of an invitro diagnostic kit.
Kits of the invention may further comprise reagents (e.g., antibodies) for assessing the status (e.g., presence, absence, level) of various additional markers, e.g., those given in
Consider a sample from a patient with some type of cancer. The mutation screening of this sample identifies somatic mutations in n genes g1, g2, . . . , gn. Assuming that somatic mutations occur independently, the probability that this patient has cancer of type c is given by the following equation:
P(c|g1,g2, . . . ,gn)=P0(c)ΠiM(gi|c)/ΣtP0(t)ΠiM(gi|t) (1)
where the product is taken over all genes mutated in the sample (i=1, 2, . . . , n) and the sum is taken over all cancer types t. M(g|c) is the frequency of somatic mutations in gene g in cancer type c. See
Using Equation (1), the probabilities were calculated for each cancer type c, and the cancer with the highest probability was designated the most likely cancer type in the patient. Such a patient may be examined by available diagnostic techniques for this cancer type. If the most likely cancer type is not confirmed, the cancer type with the second highest probability should be examined and so on.
Since some organs can develop cancers of different types (such as adenocarcinoma and squamous cell carcinoma in lung), one may calculate the probability P(o) that the cancer has developed in organ o:
P(o|g1,g2, . . . ,gn)=ΣcP(c|g1,g2, . . . ,gn) (2)
where the sum is over all cancer types c of the organ o. Using Equation (2), the probabilities are calculated for each organ o, and the organ with the highest probability is the most likely cancer site in the patient. The patient may optionally be examined by additional diagnostic techniques to confirm this cancer site. If the most likely cancer site is not confirmed, the organ with the second highest probability may then be examined and so on.
In order to evaluate the power of using mutations to determine the tumor site, we used three published studies (PMID: 17932254, PMID: 18772397, PMID: 18772396) in which over 20,000 genes were sequenced in samples representing four cancers: 11 breast ductal carcinoma samples, 11 colon adenocarcinoma samples, 22 glioblastoma samples, and 24 pancreatic ductal carcinoma samples. We used these datasets as a validation dataset for our approach. In order to calculate the probabilities given by Equation (1) for these samples, we used two sets of genes.
The first set of genes represents all the genes with mutation frequency above 5% in one of 29 common cancer types. Using COSMIC database we identified 33 such genes:
Using this set of genes the following results were obtained:
“Percent Correct 1” is the percent of samples for which the cancer type with highest predicted probability coincided with the true cancer type of the sample, “Percent Correct 2” is the percent of samples for which cancer type with the second highest predicted probability coincided with the true cancer type of the sample, and “Percent Wrong” is the percent of samples for which cancers types with neither highest nor second highest predicted probabilities coincided with the true cancer type of the sample.
The second set of genes was based on the validation dataset. The set was composed of genes which satisfied the following conditions:
17 genes satisfied these conditions:
Using this list of genes thus gave better prediction accuracy, as shown in the following table:
Some modifications to the above approach may be applied individually or in combination to improve results or under certain circumstances.
P(c|g1,g2, . . . ,gn)=P0(c)ΠiM(gi|c)Πj(1−M(gi|c))/ΠtP0(t)ΠiM(gi|t)Πj(1−M(gi|c))
Since somatic mutations are very specific to cancer or pre-cancerous conditions, the main performance measure of using mutation screening of a set of genes is its sensitivity. The sensitivity of screening for any cancer depends on sensitivities within individual cancers as well as on the incidences of the cancers. The sensitivity was defined by the following equation:
S=Σ
t
P
0(t)S(t) (4)
where S(t) was the sensitivity within cancer type t. S(t) was defined as the percentage of patients with somatic mutations in one or more genes within a predefined set of one or more genes.
The following algorithm was used to define a small set of genes with high sensitivity:
The same validation dataset described above was used. The list of genes in the order they were define by the above algorithm is shown below, with the cumulative sensitivity as a function of the number of the genes in the list presented in
Some modifications to the above approach may be applied individually or in combination to improve results or under certain circumstances.
To confirm our ability to detect cancer-related mutations in serum exosomes, cell culture supernatants (1-10 ml from ovarian and colon cancer cell lines) or ovarian and colon cancer patient serum samples (1-3 ml) were used to prepare exosomes by high-speed centrifugation. Total RNA was extracted from exosomal pellets and converted to cDNA by standard methods. PCR amplicons for a set of mutation hot spots in TP53, KRAS, EGFR and APC were designed and optimized for multiplexing. Exosomal cDNA was pre-amplified with a multiplex of all amplicons. The pre-amplification product was split into separate reactions and re-amplified with the individual target amplicons. Re-amplification primers were synthesized with tails for dye-primer sequencing. Individual PCR products were sequenced by dye-primer chemistry to identify particular mutations.
Mutations were found in exosomes harvested from cell lines as follows:
Mutations were found in cancer serum samples as follows (gels showing mutations in ovarian cancer serum shown in
Examples of important mutations in genes listed in Table 1 are shown below in Table 7:
Examples of important mutations found in cancer serum exosomes in Table 1 genes are shown below in Table 8:
Although the foregoing invention has been described in some detail by way of illustration and example for purposes of clarity of understanding, it will be clear to those skilled in the art that certain changes and modifications may be practiced within the scope of the appended claims.
This application claims the priority benefit of PCT/US10/037,659, filed Jun. 7, 2010, which claims the priority benefit of U.S. Provisional Application Ser. No. 61/184,685 (filed on Jun. 5, 2009), which is hereby incorporated by reference in its entirety.