Many biological functions are accomplished by altering the expression of various genes through transcriptional (e.g., through control of initiation, provision of RNA precursors, RNA processing, etc.) and/or translational control. For example, fundamental biological processes such as cell cycle, cell differentiation and cell death, are often characterized by the variations in the expression levels of groups of genes.
Changes in gene expression also are associated with pathogenesis. For example, the lack of sufficient expression of functional tumor suppressor genes and/or the over expression of oncogene/protooncogenes could lead to tumorgenesis or hyperplastic growth of cells (Marshall, Cell 64:313-326, 1991; Weinberg, Science, 254:1138-1146, 1991). Thus, changes in the expression levels of particular genes (e.g., oncogenes or tumor suppressors) serve as signposts for the presence and progression of various diseases.
Monitoring changes in gene expression may also provide certain advantages during drug screening and development. Often drugs are pre-screened for the ability to interact with a major target without regard to other effects the drugs have on cells. Often such other effects cause toxicity in the whole animal, which prevent the development and use of the potential drug.
Using pairs of samples from subjects, applicants have examined samples from diseased but non-cancerous liver tissue and from cancerous liver tissue to identify global changes in gene expression between tumor biopsies and surrounding non-cancerous tissue. Diseased but non-cancerous liver tissue was either inflamed tissue from chronic viral hepatitis patients or fibrotic tissue from liver cirrhosis patients. Non-cancerous tissue was removed from a point in the liver adjacent to a tumor biopsy site. These global changes in gene expression, also referred to as expression profiles, provide useful markers for diagnostic uses as well as markers that can be used to monitor disease states, disease progression, drug toxicity, drug efficacy and drug metabolism.
The gene expression profiles described herein were derived from diseased liver biopsy samples from Korean patients 34-65 years old. These patients had been diagnosed with chronic hepatitis or cirrhosis and, in each case, had subsequently developed liver cancer. The disease state associated with each sample is indicated in Table 2.
The present invention provides compositions and methods to detect the level of expression of genes that may be differentially expressed dependent upon the state of the cell, i.e., non-cancerous versus cancerous. These expression profiles of genes provide molecular tools for evaluating toxicity, drug efficacy, drug metabolism, development, and disease monitoring. Changes in the expression profile from a baseline profile can be used as an indication of such effects. Those skilled in the art can use any of a variety of known techniques to evaluate the expression of one or more of the genes and/or gene fragments identified in the instant application in order to observe changes in the expression profile in a tissue or sample of interest.
In the description that follows, numerous terms and phrases known to those skilled in the art are used. In the interest of clarity and consistency of interpretation, the definitions of certain terms and phrases are provided.
As used herein, the phrase “detecting the level of expression” includes methods that quantify expression levels as well as methods that determine whether a gene of interest is expressed at all. Thus, an assay which provides a yes or no result without necessarily providing quantification of an amount of expression is an assay that requires “detecting the level of expression” as that phrase is used herein.
As used herein, oligonucleotide sequences that are complementary to one or more of the genes described herein, refers to oligonucleotides that are capable of hybridizing under stringent conditions to at least part of the nucleotide sequence of said genes. Such hybridizable oligonucleotides will typically exhibit at least about 75% sequence identity at the nucleotide level to said genes, preferably about 80% or 85% sequence identity or more preferably about 90% or 95% or more nucleotide sequence identity to said genes.
“Bind(s) substantially” refers to complementary hybridization between a probe nucleic acid and a target nucleic acid and embraces minor mismatches that can be accommodated by reducing the stringency of the hybridization media to achieve the desired detection of the target polynucleotide sequence.
The terms “background” or “background signal intensity” refer to hybridization signals resulting from non-specific binding, or other interactions, between the labeled target nucleic acids and components of the oligonucleotide array (e.g., the oligonucleotide probes, control probes, the array substrate, etc.). Background signals may also be produced by intrinsic fluorescence of the array components themselves. A single background signal can be calculated for the entire array, or a different background signal may be calculated for each target nucleic acid. In a preferred embodiment, background is calculated as the average hybridization signal intensity for the lowest 5% to 10% of the probes in the array, or, where a different background signal is calculated for each target gene, for the lowest 5% to 10% of the probes for each gene. Of course, one of skill in the art will appreciate that where the probes to a particular gene hybridize well and thus appear to be specifically binding to a target sequence, they should not be used in a background signal calculation. Alternatively, background may be calculated as the average hybridization signal intensity produced by hybridization to probes that are not complementary to any sequence found in the sample (e.g., probes directed to nucleic acids of the opposite sense or to genes not found in the sample such as bacterial genes where the sample is mammalian nucleic acids). Background can also be calculated as the average signal intensity produced by regions of the array that lack any probes at all.
The phrase “hybridizing specifically to” refers to the binding, duplexing or hybridizing of a molecule substantially to or only to a particular nucleotide sequence or sequences under stringent conditions when that sequence is present in a complex mixture (e.g., total cellular) DNA or RNA.
Assays and methods of the invention may utilize available formats to simultaneously screen at least about 100, preferably about 1000, more preferably about 10,000 and most preferably about 1,000,000 or more different nucleic acid hybridizations.
The terms “mismatch control” or “mismatch probe” refer to a probe whose sequence is deliberately selected not to be perfectly complementary to a particular target sequence. For each mismatch (MM) control in a high-density array there typically exists a corresponding perfect match (PM) probe that is perfectly complementary to the same particular target sequence. The mismatch may comprise one or more bases that are not complementary to the corresponding bases of the target sequence.
While the mismatch(s) may be located anywhere in the mismatch probe, terminal mismatches are less desirable as a terminal mismatch is less likely to prevent hybridization of the target sequence. In a particularly preferred embodiment, the mismatch is located at or near the center of the probe such that the mismatch is most likely to destabilize the duplex with the target sequence under the test hybridization conditions.
The term “perfect match probe” refers to a probe that has a sequence that is perfectly complementary to a particular target sequence. The test probe is typically perfectly complementary to a portion (subsequence) of the target sequence. The perfect match (PM) probe can be a “test probe”, a “normalization control” probe, an expression level control probe and the like. A perfect match control or perfect match probe is, however, distinguished from a “mismatch control” or “mismatch probe.”
As used herein a “probe” is defined as a nucleic acid, preferably an oligonucleotide, capable of binding to a target nucleic acid of complementary sequence through one or more types of chemical bonds, usually through complementary base pairing, usually through hydrogen bond formation. As used herein, a probe may include natural (i.e., A, G, U, C or T) or modified bases (7-deazaguanosine, inosine, etc.). In addition, the bases in probes may be joined by a linkage other than a phosphodiester bond, so long as it does not interfere with hybridization. Thus, probes may be peptide nucleic acids in which the constituent bases are joined by peptide bonds rather than phosphodiester linkages.
The term “stringent conditions” refers to conditions under which a probe will hybridize to its target subsequence, but with only insubstantial hybridization to other sequences or to other sequences such that the difference may be identified. Stringent conditions are sequence-dependent and will be different in different circumstances. Longer sequences hybridize specifically at higher temperatures. Generally, stringent conditions are selected to be about 5° C. lower than the thermal melting point (Tm) for the specific sequence at a defined ionic strength and pH.
Typically, stringent conditions will be those in which the salt concentration is at least about 0.01 to 1.0 M sodium ion concentration (or other salts) at pH 7.0 to 8.3 and the temperature is at least about 30° C. for short probes (e.g., 10 to 50 nucleotide). Stringent conditions may also be achieved with the addition of destabilizing agents such as formamide.
The “percentage of sequence identity” or “sequence identity” is determined by comparing two optimally aligned sequences or subsequences over a comparison window or span, wherein the portion of the polynucleotide sequence in the comparison window may optionally comprise additions or deletions (i.e., gaps) as compared to the reference sequence (which does not comprise additions or deletions) for optimal alignment of the two sequences. The percentage is calculated by determining the number of positions at which the identical subunit (e.g., nucleic acid base or amino acid residue) occurs in both sequences to yield the number of matched positions, dividing the number of matched positions by the total number of positions in the window of comparison and multiplying the result by 100 to yield the percentage of sequence identity. Percentage sequence identity when calculated using the programs GAP or BESTFIT (see below) is calculated using default gap weights.
Homology or identity may be determined by BLAST (Basic Local Alignment Search Tool) analysis using the algorithm employed by the programs blastp, blastn, blastx, tblastn and tblastx (Karlin et al., Proc Natl Acad Sci USA 87:2264-2268, 1990 and Altschul, J Mol Evol 36:290-300, 1993, fully incorporated by reference) which are tailored for sequence similarity searching. The approach used by the BLAST program is to first consider similar segments between a query sequence and a database sequence, then to evaluate the statistical significance of all matches that are identified and finally to summarize only those matches which satisfy a preselected threshold of significance. For a discussion of basic issues in similarity searching of sequence databases, see Altschul et al., (Nature Genet 6:119-129, 1994) which is fully incorporated by reference. The search parameters for histogram, descriptions, alignments, expect (i.e., the statistical significance threshold for reporting matches against database sequences), cutoff, matrix and filter are at the default settings. The default scoring matrix used by blastp, blastx, tblastn, and tblastx is the BLOSUM62 matrix (Henikoff et al., Proc Natl Acad Sci USA 89:10915-10919, 1992, fully incorporated by reference). Four blastn parameters were adjusted as follows: Q=10 (gap creation penalty); R=10 (gap extension penalty); wink=1 (generates word hits at every winkth position along the query); and gapw=16 (sets the window width within which gapped alignments are generated). The equivalent Blastp parameter settings were Q=9; R=2; wink=1; and gapw=32. A Bestfit comparison between sequences, available in the GCG package version 10.0, uses DNA parameters GAP=50 (gap creation penalty) and LEN=3 (gap extension penalty) and the equivalent settings in protein comparisons are GAP=8 and LEN=2.
The present invention identifies those genes differentially expressed between cancerous and non-cancerous liver tissue. One of skill in the art can select one or more of the genes identified as being differentially expressed in Table 1 and use the information and methods provided herein to interrogate or test a particular sample. For a particular interrogation of two conditions or sources, it may be desirable to select those genes which display a great deal of difference in the expression pattern between the two conditions or sources. In other instances, it may be appropriate to select genes whose expression changes only slightly between the two conditions. At least a 1.5-fold difference may be desirable, but a three-fold, five-fold or ten-fold difference may be preferred in some instances. The data are subjected to statistical evaluation to ensure that the observed differences and-the disease association are statistically significant. Interrogations of the genes or proteins can be performed to yield different information.
As described herein, the genes and gene expression information provided in Table 1 may be used as diagnostic markers for the prediction or identification of a disease state of liver tissue. For instance, a liver tissue sample or other sample from a patient may be assayed by any of the methods known to those skilled in the art, and the expression levels from one or more genes from Table 1 may be compared to the expression levels found in non-cancerous liver tissue, cancerous liver tissue or both. Expression profiles generated from the tissue or other samples that substantially resemble an expression profile from non-cancerous or cancerous liver tissue may be used, for instance, to aid in disease diagnosis. Comparison of the expression data, as well as available sequence or other information, may be done by a researcher or diagnostician or may be done with the aid of a computer and databases as described herein.
Molecular expression markers for liver disease can be used to confirm the type and progression of disease made on the basis of morphological criteria. For example, non-cancerous liver tissue could be distinguished from cancerous tissue based on the level and type of genes expressed in a tissue sample. In some situations, identifications of cell type or source is ambiguous based on classical criteria. In these situations, the molecular expression markers of the present invention are useful for identifying the region of the liver from which a sample came, as well as whether or not normal levels of gene expression have been altered (signs of metabolic disturbances).
In addition, progression of hepatic carcinoma to new areas of the liver can be monitored by following the expression patterns of the involved genes using the molecular expression markers of the present invention. Monitoring of the efficacy of certain drug regimens can also be accomplished by following the expression patterns of the molecular expression markers.
As described above, the genes and gene expression information provided in Table 1 may also be used as markers for the direct monitoring of disease progression, for instance, the development of liver cancer. A liver tissue sample or other sample from a patient may be assayed by any of the methods known to those of skill in the art, and the expression levels in the sample from a gene or genes from Table 1 may be compared to the expression levels found in non-cancerous liver tissue, tissue from a hepatic carcinoma or both. Comparison of the expression data, as well as available sequence or other information may be done by a researcher or diagnostician or may be done with the aid of a computer and databases as described herein.
According to the present invention, potential drugs can be screened to determine if application of the drug alters the expression of one or more of the genes identified herein. This may be useful, for example, in determining whether a particular drug is effective in treating a particular patient with liver disease. In the case where a gene's expression is affected by the potential drug such that its level of expression returns to normal, the drug is indicated in the treatment of liver cancer. Similarly, a drug which causes expression of a gene which is not normally expressed by healthy liver cells may be contra-indicated in the treatment of liver cancer.
According to the present invention, the genes identified in Table 1 may also be used as markers to evaluate the effects of a candidate drug or agent on a cell, particularly a cell undergoing malignant transformation, for instance, a liver cancer cell or tissue sample. A candidate drug or agent can be screened for the ability to stimulate the transcription or expression of a given marker or markers (drug targets) or to down-regulate or inhibit the transcription or expression of a marker or markers. According to the present invention, one can also compare the specificity of a drug's effects by looking at the number of markers affected by the drug and comparing them to the number of markers affected by a different drug. A more specific drug will affect fewer transcriptional targets. Similar sets of markers identified for two drugs indicates a similarity of effects.
Assays to monitor the expression of a marker or markers as defined in Table 1 may utilize any available means of monitoring for changes in the expression level of the nucleic acids of the invention. As used herein, an agent is said to modulate the expression of a nucleic acid of the invention if it is capable of up- or down-regulating expression of the nucleic acid in a cell.
Agents that are assayed in the above methods can be randomly selected or rationally selected or designed. As used herein, an agent is said to be randomly selected when the agent is chosen randomly without considering the specific sequences involved in the association of a protein of the invention alone or with its associated substrates, binding partners, etc. An example of randomly selected agents is the use a chemical library or a peptide combinatorial library, or a growth broth of an organism.
As used herein, an agent is said to be rationally selected or designed when the agent is chosen on a nonrandom basis which takes into account the sequence of the target site and/or its conformation in connection with the agents action. Agents can be selected or designed by utilizing the peptide sequences that make up these sites. For example, a rationally selected peptide agent can be a peptide whose amino acid sequence is identical to or a derivative of any functional consensus site.
The agents of the present invention can be, as examples, peptides, small chemical molecules, vitamin derivatives, as well as carbohydrates, lipids, oligonucleotides and covalent and non-covalent combinations thereof. Dominant negative proteins, DNA encoding these proteins, antibodies to these proteins, peptide fragments of these proteins or mimics of these proteins may be introduced into cells to affect function. “Mimic” as used herein refers to the modification of a region or several regions of a peptide molecule to provide a structure chemically different from the parent peptide but topographically and functionally similar to the parent peptide (see Grant, in Molecular Biology and Biotechnology, Meyers (ed.), VCH Publishers, 1995). A skilled artisan can readily recognize that there is no limit as to the structural nature of the agents of the present invention.
Agents that up- or down-regulate or modulate the expression of the nucleic acid molecules of Table 1, or at least one activity of a protein encoded by the nucleic acid molecules of Table 1, such as agonists or antagonists, may be used to modulate biological and pathologic processes associated with the function and activity of the proteins encoded by these nucleic acid molecules. The agents can be the nucleic acid molecules of Table 1 themselves, the encoded proteins, or portions of these molecules, such as all or part of the open reading frames of these nucleic acid molecules.
Anti-sense oligonucleotide molecules derived from the nucleic acid sequences of Table 1 may also be used to down-regulate the expression of one or more of the genes in Table 1 that are expressed at elevated levels in liver cancer, the use of antisense gene therapy being an example. Down-regulation of expression of one or more of the genes of Table 1 is accomplished by administering an effective amount of antisense oligonucleotides. These antisense molecules can be fashioned from the DNA sequences of these genes or sequences containing various mutations, deletions, insertions or spliced variants. Isolated RNA or DNA sequences derived from these genes may also be used therapeutically in gene therapy. These agents may be used to induce gene expression in liver cancers associated with an absence of or considerably decreased expression of one or more of the proteins encoded by genes in Table 1.
As used herein, a subject can be any mammal, so long as the mammal is in need of modulation of a pathological or biological process mediated by a gene of the invention. The term “mammal” is defined as an individual belonging to the class Mammalia. The invention is particularly useful in the treatment of human subjects.
Pathological processes refer to a category of biological processes which produce a deleterious effect. For example, expression of a gene of the invention may be associated with hyperplasia in the liver, in particular malignant hyperplasia. As used herein, an agent is said to modulate a pathological process when the agent reduces the degree or severity of the process. For instance, liver cancer may be prevented or disease progression modulated by the administration of agents which up- or down-regulate or modulate in some way the expression or at least one activity of a gene of the invention.
The agents of the present invention can be provided alone, or in combination with other agents that modulate a particular pathological process. For example, an agent of the present invention can be administered in combination with other known drugs. As used herein, two agents are said to be administered in combination when the two agents are administered simultaneously or are administered independently in a fashion such that the agents will act at the same time.
The agents of the present invention can be administered via parenteral, subcutaneous, intravenous, intramuscular, intraperitoneal, transdermal, or buccal routes. Alternatively, or concurrently, administration may be by the oral route. The dosage administered will be dependent upon the age, health, and weight of the recipient, kind of concurrent treatment, if any, frequency of treatment, and the nature of the effect desired.
The present invention further provides compositions containing one or more agents which modulate expression or at least one activity of a protein of the invention. While individual needs vary, determination of optimal ranges of effective amounts of each component is within the skill of the art. Typical dosages comprise 0.1 to 100 μg/kg body wt. The preferred dosages comprise 0.1 to 10 μg/kg body wt. The most preferred dosages comprise 0.1 to 1 μg/kg body wt.
In addition to the pharmacologically active agent, the compositions of the present invention may contain suitable pharmaceutically acceptable carriers comprising excipients and auxiliaries which facilitate processing of the active compounds into preparations which can be used pharmaceutically for delivery to the site of action. Suitable formulations for parenteral administration include aqueous solutions of the active compounds in water-soluble form, for example, water-soluble salts. In addition, suspensions of the active compounds as appropriate oily injection suspensions may be administered. Suitable lipophilic solvents or vehicles include fatty oils, e.g., sesame oil, or synthetic fatty acid esters, e.g. ethyl oleate or triglycerides. Aqueous injection suspensions may contain substances which increase the viscosity of the suspension include, for example, sodium carboxymethyl cellulose, sorbitol, and/or dextran. Optionally, the suspension may also contain stabilizers. Liposomes can also be used to encapsulate the agent for delivery into the cell.
The pharmaceutical formulation for systemic administration according to the invention may be formulated for enteral, parenteral or topical administration. Indeed, all three types of formulations may be used simultaneously to achieve systemic administration of the active ingredient.
Suitable formulations for oral administration include hard or soft gelatin capsules, pills, tablets, including coated tablets, elixirs, suspensions, syrups or inhalations and controlled release forms thereof.
In practicing the methods of this invention, the compounds of this invention may be used alone or in combination, or in combination with other therapeutic or diagnostic agents. In certain preferred embodiments, the compounds of this invention may be coadministered along with other compounds typically prescribed for these conditions according to generally accepted medical practice. The compounds of this invention can be utilized in vivo, ordinarily in mammals, such as humans, rats, mice, dogs, cats, sheep, horses, cattle and pigs, or in vitro.
The genes identified as being differentially expressed in liver disease may be used in a variety of nucleic acid detection assays to detect or quantify the expression level of a gene or multiple genes in a given sample. For example, traditional Northern blotting, nuclease protection, RT-PCR and differential display methods may be used for detecting gene expression levels. In methods where small numbers of genes are assayed, such as 5-50 genes, high-throughput PCR may be used.
The protein products of the genes identified herein can also be assayed to determine the amount of expression. Methods for assaying for a protein include Western blot, immunoprecipitation and radioimmunoassay. In some methods, it is preferable to assay the mRNA as an indication of expression. Methods for assaying for mRNA include Northern blots, slot blots, dot blots, and hybridization to an ordered array of oligonucleotides. Any method for specifically and quantitatively measuring a specific protein or mRNA or DNA product can be used. However, methods and assays of the invention are most efficiently designed with array or chip hybridization-based methods for detecting the expression of a large number of genes.
Any hybridization assay format may be used, including solution-based and solid support-based assay formats. A preferred solid support is a high density array also known as a DNA chip or a gene chip. One variation of the DNA chip contains hundreds of thousands of discrete microscopic channels that pass completely through it. Probe molecules are attached to the inner surface of these channels, and molecules from the samples to be tested flow throughout the channels, coming into close proximity with the probes for hybridization. In one assay format, gene chips containing probes to at least two genes from Table 1 may be used to directly monitor or detect changes in gene expression in the treated or exposed cell as described herein.
The genes of the present invention may be assayed in any convenient sample form. For example, samples may be assayed in the form mRNA or reverse transcribed mRNA. Samples may be cloned or not, and the samples or individual genes may be amplified or not. The cloning itself does not appear to bias the representation of genes within a population. However, it may be preferable to use polyA+ RNA as a source, as it can be used with less processing steps. In some embodiments, it may be preferable to assay the protein or peptide expressed by the gene.
The sequences of the expression marker genes of Table 1 are available in the public databases. Table 1 provides the Accession number, Sequence Number ID and name for each of the sequences. The sequences of the genes in GenBank are herein expressly incorporated by reference in their entirety (see www.ncbi.nim.nih.gov).
Additional assay formats may be used to monitor the ability of the agent to modulate the expression of a gene identified in Table 1. For instance, as described above, mRNA expression may be monitored directly by hybridization of probes to the nucleic acids of the invention. Cell lines are exposed to an agent to be tested under appropriate conditions and time and total RNA or mRNA is isolated by standard procedures such those disclosed in Sambrook et al., Molecular Cloning—A Laboratory Manual, Third ed., Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y., 2001. In some embodiments, it may be desirable to amplify one or more of the RNA molecules isolated prior to application of the RNA to the gene chip. Using techniques well known in the art, the RNA may be reverse transcribed and amplified in the form of DNA or may be reverse transcribed into DNA and the DNA used as a template for transcription to generate recombinant RNA. Any method that results in the production of a sufficient quantity of nucleic acid to be hybridized effectively to the gene chip may be used.
In another format, cell lines that contain reporter gene fusions between the open reading frame and/or the 3′ or 5′ regulatory regions of a gene in Table 1 and any assayable fusion partner may be prepared. Numerous assayable fusion partners are known and readily available including the firefly luciferase gene and the gene encoding chloramphenicol acetyltransferase (Alam et al., Anal Biochem 188:245-254, 1990). Cell lines containing the reporter gene fusions are then exposed to the agent to be tested under appropriate conditions and time. Differential expression of the reporter gene between samples exposed to the agent and control samples identifies agents which modulate the expression of the nucleic acid.
In another assay format, cells or cell lines are first identified which express one or more of the gene products of the invention physiologically. Cells and/or cell lines so identified would preferably comprise the necessary cellular machinery to ensure that the transcriptional and/or translational apparatus of the cells would faithfully mimic the response of normal or cancerous liver tissue to an exogenous agent. Such machinery would likely include appropriate surface transduction mechanisms and/or cytosolic factors. Such cell lines may be, but are not required to be, derived from liver tissue. The cells and/or cell lines may then be contacted with an agent and the expression of one or more of the genes of interest may then be assayed. The genes may be assayed at the mRNA level and/or at the protein level.
In some embodiments, such cells or cell lines may be transduced or transfected with an expression vehicle (e.g., a plasmid or viral vector) containing an expression construct comprising an operable 5′-promoter containing end of a gene of interest identified in Table 1 fused to one or more nucleic acid sequences encoding one or more antigenic fragments. The construct may comprise all or a portion of the coding sequence of the gene of interest which may be positioned 5′- or 3′-to a sequence encoding an antigenic fragment. The coding sequence of the gene of interest may be translated or un-translated after transcription of the gene fusion. At least one antigenic fragment may be translated. The antigenic fragments are selected so that the fragments are under the transcriptional control of the promoter of the gene of interest and are expressed in a fashion substantially similar to the expression pattern of the gene of interest. The antigenic fragments may be expressed as polypeptides whose molecular weight can be distinguished from the naturally occurring polypeptides.
In some embodiments, gene products of the invention may further comprise an immunologically distinct tag. Such a process is well known in the art (see Sambrook et al., supra). Cells or cell lines transduced or transfected as outlined above are then contacted with agents under appropriate conditions; for example, the agent comprises a pharmaceutically acceptable excipient and is contacted with cells comprised in an aqueous physiological buffer such as phosphate buffered saline (PBS) at physiological pH, Eagles balanced salt solution (BSS) at physiological pH, PBS or BSS comprising serum or conditioned media comprising PBS or BSS and serum incubated at 37° C. Said conditions may be modulated as deemed necessary by one of skill in the art. Subsequent to contacting the cells with the agent, said cells will be disrupted and the polypeptides of the lysate are fractionated such that a polypeptide fraction is pooled and contacted with an antibody to be further processed by immunological assay (e.g., ELISA, immunoprecipitation or Western blot). The pool of proteins isolated from the “agent-contacted” sample will be compared with a control sample where only the excipient is contacted with the cells and an increase or decrease in the immunologically generated signal from the “agent-contacted” sample compared to the control will be used to distinguish the effectiveness of the agent.
Another embodiment of the present invention provides methods for identifying agents that modulate the levels, concentration or at least one activity of a protein(s) encoded by the genes in Table 1. Such methods or assays may utilize any means of monitoring or detecting the desired activity.
In one format, the relative amounts of a protein of the invention produced in a cell population that has been exposed to the agent to be tested may be compared to the amount produced in an unexposed control cell population. In this format, probes such as specific antibodies are used to monitor the differential expression of the protein in the different cell populations. Cell lines or populations are exposed to the agent to be tested under appropriate conditions and time. Cellular lysates may be prepared from the exposed cell line or population and a control, unexposed cell line or population. The cellular lysates are then analyzed with the probe, such as a specific antibody.
Probes based on the sequences of the genes described herein may be prepared by any commonly available method. Oligonucleotide probes for assaying the tissue or cell sample are preferably of sufficient length to specifically hybridize only to appropriate, complementary genes or transcripts. Typically the oligonucleotide probes will be at least 10, 12, 14, 16, 18, 20 or 25 nucleotides in length. In some cases longer probes of at least 30, 40, or 50 nucleotides will be desirable.
One of skill in the art will appreciate that an enormous number of array designs are suitable for the practice of this invention. The high density array will typically include a number of probes that specifically hybridize to the sequences of interest. See WO 99/32660 for methods of producing probes for a given gene or genes. In addition, in a preferred embodiment, the array will include one or more control probes.
High density array chips of the invention include “test probes.” Test probes may be oligonucleotides that range from about 5 to about 500 or about 5 to about 50 nucleotides, more preferably from about 10 to about 40 nucleotides and most preferably from about 15 to about 40 nucleotides in length. In other particularly preferred embodiments, the probes are about 20 or 25 nucleotides in length. In another preferred embodiment, test probes are double or single strand DNA sequences. DNA sequences may be isolated or cloned from natural sources or amplified from natural sources using natural nucleic acid as templates. These probes have sequences complementary to particular subsequences of the genes whose expression they are designed to detect. Thus, the test probes are capable of specifically hybridizing to the target nucleic acid they are to detect.
In addition to test probes that bind the target nucleic acid(s) of interest, the high density array can contain a number of control probes. The control probes fall into three categories referred to herein as (1) normalization controls; (2) expression level controls; and (3) mismatch controls.
Normalization controls are oligonucleotide or other nucleic acid probes that are complementary to labeled reference oligonucleotides or other nucleic acid sequences that are added to the nucleic acid sample. The signals obtained from the normalization controls after hybridization provide a control for variations in hybridization conditions, label intensity, “reading” efficiency and other factors that may cause the signal of a perfect hybridization to vary between arrays. In a preferred embodiment, signals (e.g., fluorescence intensity) read from all other probes in the array are divided by the signal (e.g., fluorescence intensity) from the control probes thereby normalizing the measurements.
Virtually any probe may serve as a normalization control. However, it is recognized that hybridization efficiency varies with base composition and probe length. Preferred normalization probes are selected to reflect the average length of the other probes present in the array, however, they can be selected to cover a range of lengths. The normalization control(s) can also be selected to reflect the (average) base composition of the other probes in the array, however in a preferred embodiment, only one or a few probes are used and they are selected such that they hybridize well (i.e., no secondary structure) and do not match any target-specific probes.
Expression level controls are probes that hybridize specifically with constitutively expressed genes in the biological sample. Virtually any constitutively expressed gene provides a suitable target for expression level controls. Typical expression level control probes have sequences complementary to subsequences of constitutively expressed “housekeeping genes” including, but not limited to the β-actin gene, the transferrin receptor gene, the GAPDH gene, and the like.
Mismatch controls may also be provided for the probes to the target genes, for expression level controls or for normalization controls. Mismatch controls are oligonucleotide probes or other nucleic acid probes identical to their corresponding test or control probes except for the presence of one or more mismatched bases. A mismatched base is a base selected so that it is not complementary to the corresponding base in the target sequence to which the probe would otherwise specifically hybridize. One or more mismatches are selected such that under appropriate hybridization conditions (e.g., stringent conditions) the test or control probe would be expected to hybridize with its target sequence, but the mismatch probe would not hybridize (or would hybridize to a significantly lesser extent). Preferred mismatch probes contain a central mismatch. Thus, for example, where a probe is a twenty-mer, a corresponding mismatch probe may have the identical sequence except for a single base mismatch (e.g., substituting a G, a C or a T for an A) at any of positions 6 through 14 (the central mismatch).
Mismatch probes thus provide a control for non-specific binding or cross hybridization to a nucleic acid in the sample other than the target to which the probe is directed. Mismatch probes also indicate whether a hybridization is specific or not. For example, if the target is present the perfect match probes should be consistently brighter than the mismatch probes. In addition, if all central mismatches are present, the mismatch probes can be used to detect a mutation. The difference in intensity between the perfect match and the mismatch probe (I(PM)-I(MM)) provides a good measure of the concentration of the hybridized material.
As is apparent to one of ordinary skill in the art, nucleic acid samples used in the methods and assays of the invention may be prepared by any available method or process. Methods of isolating total mRNA are also well known to those of skill in the art. For example, methods of isolation and purification of nucleic acids are described in detail in Chapter 3 of Laboratory Techniques in Biochemistry and Molecular Biology, Vol. 24, Hybridization With Nucleic Acid Probes: Theory and Nucleic Acid Probes, P. Tijssen (ed.) Elsevier Press, New York, 1993. Such samples include RNA samples, but also include cDNA synthesized from a mRNA sample isolated from a cell or tissue of interest. Such samples also include DNA amplified from the cDNA, and an RNA transcribed from the amplified DNA. One of skill in the art would appreciate that it may be desirable to inhibit or destroy RNase present in homogenates before homogenates can be used.
Biological samples may be of any biological tissue or fluid or cells from any organism as well as cells raised in vitro, such as cell lines and tissue culture cells. Frequently the sample will be a “clinical sample” which is a sample derived from a patient. Typical clinical samples include, but are not limited to, liver tissue biopsy, sputum, blood, blood-cells (e.g., white cells), tissue or fine needle biopsy samples, urine, peritoneal fluid, and pleural fluid, or cells therefrom. Biological samples may also include sections of tissues, such as frozen sections or formalin fixed sections taken for histological purposes.
Solid supports containing oligonucleotide probes for differentially expressed genes can be any solid or semisolid support material known to those skilled in the art. Suitable examples include, but are not limited to, membranes, filters, tissue culture dishes, polyvinyl chloride dishes, beads, test strips, silicon or glass based chips and the like. Suitable glass wafers and hybridization methods are widely available, for example, those disclosed by Beattie (WO 95/11755). Any solid surface to which oligonucleotides can be bound, either directly or indirectly, either covalently or non-covalently, can be used. In some embodiments, it may be desirable to attach some oligonucleotides covalently and others non-covalently to the same solid support.
A preferred solid support is a high density array or DNA chip. These contain a particular oligonucleotide probe in a predetermined location on the array. Each predetermined location may contain more than one molecule of the probe, but each molecule within the predetermined location has an identical sequence. Such predetermined locations are termed features. There may be, for example, from 2, 10, 100, 1000 to 10,000, 100,000 or 400,000 of such features on a single solid support. The solid support, or the area within which the probes are attached may be on the order of a square centimeter.
Oligonucleotide probe arrays for expression monitoring can be made and used according to any techniques known in the art (see for example, Lockhart et al., Nat Biotechnol 14:1675-1680, 1996; McGall et al., Proc Nat Acad Sci USA 93: 13555-13460, 1996). Such probe arrays may contain at least two or more oligonucleotides that are complementary to or hybridize to two or more of the genes described herein. Such arrays may also contain oligonucleotides that are complementary or hybridize to at least 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 50, 70 or more the genes described herein.
Methods of forming high density arrays of oligonucleotides with a minimal number of synthetic steps are known. The oligonucleotide analogue array can be synthesized on a solid substrate by a variety of methods, including, but not limited to, light-directed chemical coupling, and mechanically directed coupling (see Pirrung et al., (1992) U.S. Pat. No. 5,143,854; Fodor et al., (1998) U.S. Pat. No. 5,800,992; Chee et al., (1998) U.S. Pat. No. 5,837,832).
In brief, the light-directed combinatorial synthesis of oligonucleotide arrays on a glass surface proceeds using automated phosphoramidite chemistry and chip masking techniques. In one specific implementation, a glass surface is derivatized with a silane reagent containing a functional group, e.g., a hydroxyl or amine group blocked by a photolabile protecting group. Photolysis through a photolithographic mask is used selectively to expose functional groups which are then ready to react with incoming 5′ photoprotected nucleoside phosphoramidites. The phosphoramidites react only with those sites which are illuminated (and thus exposed by removal of the photolabile blocking group). Thus, the phosphoramidites only add to those areas selectively exposed from the preceding step. These steps are repeated until the desired array of sequences have been synthesized on the solid surface. Combinatorial synthesis of different oligonucleotide analogues at different locations on the array is determined by the pattern of illumination during synthesis and the order of addition of coupling reagents.
In addition to the foregoing, additional methods which can be used to generate an array of oligonucleotides on a single substrate are described in Fodor et al. WO 93/09668. High density nucleic acid arrays can also be fabricated by depositing pre-made or natural nucleic acids in predetermined positions. Synthesized or natural nucleic acids are deposited on specific locations of a substrate by light directed targeting and oligonucleotide directed targeting. Another embodiment uses a dispenser that moves from region to region to deposit nucleic acids in specific spots.
Nucleic acid hybridization simply involves contacting a probe and target nucleic acid under conditions where the probe and its complementary target can form stable hybrid duplexes through complementary base pairing (see Lockhart et al., (1999) WO 99/32660). The nucleic acids that do not form hybrid duplexes are then washed away leaving the hybridized nucleic acids to be detected, typically through detection of an attached detectable label. It is generally recognized that nucleic acids are denatured by increasing the temperature or decreasing the salt concentration of the buffer containing the nucleic acids. Under low stringency conditions (e.g., low temperature and/or high salt) hybrid duplexes (e.g., DNA-DNA, RNA-RNA or RNA-DNA) will form even where the annealed sequences are not perfectly complementary. Thus, specificity of hybridization is reduced at lower stringency. Conversely, at higher stringency (e.g., higher temperature or lower salt) successful hybridization requires fewer mismatches. One of skill in the art will appreciate that hybridization conditions may be selected to provide any degree of stringency. In a preferred embodiment, hybridization is performed at low, stringency, in this case in 6×SSPE-T at 37° C. (0.005% Triton x-100) to ensure hybridization and then subsequent washes are performed at higher stringency (e.g., 1×SSPE-T at 37° C.) to eliminate mismatched hybrid duplexes. Successive washes may be performed at increasingly higher stringency (e.g., down to as low as 0.25×SSPET at 37° C. to 50° C.) until a desired level of hybridization specificity is obtained. Stringency can also be increased by addition of agents such as formamide. Hybridization specificity may be evaluated by comparison of hybridization to the test probes with hybridization to the various controls that can be present (e.g., expression level control, normalization control, mismatch controls, etc.).
In general, there is a tradeoff between hybridization specificity (stringency) and signal intensity. Thus, in a preferred embodiment, the wash is performed at the highest stringency that produces consistent results and that provides a signal intensity greater than approximately 10% of the background intensity. Thus, in a preferred embodiment, the hybridized array may be washed at successively higher stringency solutions and read between each wash. Analysis of the data sets thus produced will reveal a wash stringency above which the hybridization pattern is not appreciably altered and which provides adequate signal for the particular oligonucleotide probes of interest.
The hybridized nucleic acids are typically detected by detecting one or more labels attached to the sample nucleic acids. The labels may be incorporated by any of a number of means well known to those of skill in the art (see Lockhart et al., (1999) WO 99/32660).
The present invention includes relational databases containing sequence information, for instance for one or more of the genes of Table 1, as well as gene expression information in various liver tissue samples. Databases may also contain information associated with a given sequence or tissue sample such as descriptive information about the gene associated with the sequence information, descriptive information concerning the clinical status of the tissue sample, or information concerning the patient from which the sample was derived. The database may be designed to include different parts, for instance a sequence database and a gene expression database. The databases of the invention may be stored on any available computer-readable medium. Methods for the configuration and construction of such databases are widely available, for instance, see Akerblom et al., (U.S. Pat. No. 5,953,727), which is specifically incorporated herein by reference in its entirety.
The databases of the invention may be linked to an outside or external database. In a preferred embodiment, as described in Table 1, the external database is GenBank and the associated databases maintained by the National Center for Biotechnology Information or NCBI (http://www.ncbi.nlm.nih.gov/Entrez/). Other external databases that may be used in the invention include those provided by Chemical Abstracts Service (http://stnweb.cas.org/) and Incyte Genomics (http://www.incyte.com/sequence/index.shtml).
Any appropriate computer platform may be used to perform the necessary comparisons between sequence information, gene expression information and any other information in the database or provided as an input. For example, a large number of computer workstations are available from a variety of manufacturers, such has those available from Silicon Graphics. Client-server environments, database servers and networks are also widely available and appropriate platforms for the databases of the invention.
The databases of the invention may be used to produce, among other things, electronic Northern blots (E-Northerns) to allow the user to determine the cell type or tissue in which a given gene is expressed and to allow determination of the abundance or expression level of a given gene in a particular tissue or cell. The E-northern analysis can be used as a tool to discover tissue specific candidate therapeutic targets that are not over-expressed in tissues such as the liver, kidney, or heart. These tissue types often lead to detrimental side effects once drugs are developed and a first-pass screen to eliminate these targets early in the target discovery and validation process would be beneficial.
The databases of the invention may also be used to present information identifying the expression level in a tissue or cell of a set of genes comprising at least one gene in Table 1, comprising the step of comparing the expression level of at least one gene in Table 1 in the tissue to the level of expression of the gene in the database. Such methods may be used to predict the physiological state of a given tissue by comparing the level of expression of a gene or genes in Table 1 from a sample to the expression levels found in normal liver tissue, tissue from liver carcinomas or both. Such methods may also be used in the drug or agent screening assays as described herein.
Without further description, it is believed that one of ordinary skill in the art can, using the preceding description and the following illustrative examples, make and utilize the compounds of the present invention and practice the claimed methods. The preceding working examples therefore, are illustrative only and should not be construed as limiting in any way the scope of the invention.
The patient tissue samples were derived from ten Korean patients, aged 34 to 65, and classified into two groups of five patients each. Each group contained samples from four men and one woman. One group of consisted of patients who had been diagnosed with chronic viral hepatitis B and who later developed hepatic carcinomas. The second group of patients had been diagnosed with cirrhosis of the liver. These people also later developed hepatic carcinomas. For each patient, tissue was obtained from two areas of the liver to produce a set of biopsy samples. In the first patient group (cancer/hepatitis), samples were removed from liver tumors and from the non-cancerous surrounding area composed of inflamed tissue (inflammation due to hepatitis). In the second group (cancer/cirrhosis), liver tissue was removed from tumors and from the non-cancerous surrounding area composed of fibrotic tissue (areas of fibrosis due to cirrhosis).
Histological analysis of each of the tissue samples was performed and samples were segregated into either non-cancerous or cancerous categories.
With minor modifications, the sample preparation protocol followed the Affymetrix GeneChip Expression Analysis Manual. Frozen tissue was first ground to powder using the Spex Certiprep 6800 Freezer Mill. Total RNA was then extracted using Trizol (Life Technologies). The total RNA yield for each sample (average tissue weight of 300 mg) was 200-500 μg. Next, mRNA was isolated using the Oligotex mRNA Midi kit (Qiagen). Since the mRNA was eluted in a final volume of 400 μl, an ethanol precipitation step was required to bring the concentration to 1 μg/μl. Using 1-5 μg of mRNA, double stranded cDNA was created using the SuperScript Choice system (Gibco-BRL). First strand cDNA synthesis was primed with a T7-(dT24) oligonucleotide. The cDNA was then phenol-chloroform extracted and ethanol precipitated to a final concentration of 1 μg/μl.
From 2 μg of cDNA, cRNA was synthesized according to standard procedures. To biotin label the cRNA, nucleotides Bio-11-CTP and Bio-16-UTP (Enzo Diagnostics) were added to the reaction. After a 37° C. incubation for six hours, the labeled cRNA was cleaned up according to the Rneasy Mini kit protocol (Qiagen). The cRNA was then fragmented (5× fragmentation buffer: 200 mM Tris-Acetate (pH 8.1), 500 mM KOAc, 150 mM MgOAc) for thirty-five minutes at 94° C.
55 μg of fragmented cRNA was hybridized on the human and the Human Genome U95 set of arrays for twenty-four hours at 60 rpm in a 45° C. hybridization oven. The chips were washed and stained with Streptavidin Phycoerythrin (SAPE) (Molecular Probes) in Affymetrix fluidics stations. To amplify staining, SAPE solution was added twice with an anti-streptavidin biotinylated antibody (Vector Laboratories) staining step in between. Hybridization to the probe arrays was detected by fluorometric scanning (Hewlett Packard Gene Array Scanner). Following hybridization and scanning, the microarray images were analyzed for quality control, looking for major chip defects or abnormalities in hybridization signal. After all chips passed QC, the data was analyzed using Affymetrix GeneChip software (v3.0), and Experimental Data Mining Tool (EDMT) software (v1.0).
All samples were prepared as described and hybridized onto the Affymetrix Human Genome U95 array. Each chip contains 16-20 oligonucleotide probe pairs per gene or cDNA clone. These probe pairs include perfectly matched sets and mismatched sets, both of which are necessary for the calculation of the average difference. The average difference is a measure of the intensity difference for each probe pair, calculated by subtracting the intensity of the mismatch from the intensity of the perfect match. This takes into consideration variability in hybridization among probe pairs and other hybridization artifacts that could affect the fluorescence intensities. Using the average difference value that has been calculated, an absolute call for each gene is made.
The absolute call of present, absent or marginal is used to generate a Gene Signature, a tool used to identify those genes that are commonly present or commonly absent in a given sample set, according to the absolute call.
The Gene Signature Curve is a graphic view of the number of genes consistently present in a given set of samples as the sample size increases, taking into account the genes commonly expressed among a particular set of samples, and discounting those genes whose expression is variable among those samples. The curve is also indicative of the number of samples necessary to generate an accurate Gene Signature. As the sample number increases, the number of genes common to the sample set decreases. The curve is generated using the positive Gene Signatures of the samples in question, determined by adding one sample at a time to the Gene Signature, beginning with the sample with the smallest number of present genes and adding samples in ascending order. The curve displays the sample size required for the most consistency and the least amount of expression variability from sample to sample. The point where this curve begins to level off represents the minimum number of samples required for the Gene Signature. Graphed on the x-axis is the number of samples in the set, and on the y-axis is the number of genes in the positive Gene Signature. As a general rule, the acceptable percent of variability in the number of positive genes between two sample sets should be less than 5%.
For the purposes of this study, the following statistical methods were used for the data analysis. A gene set consists of genes that have a certain percentage of present calls in at least one group of samples. These genes are analyzed, and others are excluded. For example, a gene having 40% present calls (2 out of 5 samples) in at least in one sample group, cancerous cells from either hepatitis or cirrhosis patients, or non-cancerous cells from either type of patient, is included in the analysis if 40% is above the lower limit for percent present calls. Also, the genes are divided into two groups depending on their expression values across samples. For the genes in the high expression group, the average difference value is transformed to log scale before the analysis. For the genes in the low expression group, the original values are used in the analysis. An Analysis of Variance (ANOVA) method is used for data analysis (Steel et al., Principles and Procedures of Statistics: A Biometrical Approach, Third Ed., McGraw-Hill, 1997). Prior to the final analysis, a leave-one-out approach is used for outlier detection. One sample is left out of the ANOVA analysis to see whether omitting a specific sample from the analysis has any significant effect on the final result. If so, that particular sample is excluded from the final analysis. After outlier detection, the final analysis produces a list of genes that are differentially expressed with a p-value ≦0.001 as determined by the contrast from the ANOVA.
Differentially expressed genes were discovered by comparing biopsy samples from cancerous and non-cancerous regions of the same liver in patients with chronic viral hepatitis (CH) or liver cirrhosis (LC) who went on to develop primary liver cancer (hepatocellular carcinoma or HCC). Genes which showed no difference in expression level between a the cancerous and non-cancerous samples were not included in Table 1. Group 1 of Table 1 (23 genes) lists the genes that were found to be differentially expressed when the level in liver tumor cells was compared to the level in non-cancerous cells from inflamed areas or from fibrotic areas. Group 2 (12 genes) lists the genes whose expression level differed in liver tumor cells compared to cells from areas of inflammation, and group 3 contains those genes whose expression level differed in liver tumor cells compared to cells from fibrotic regions of the liver (74 genes).
The data was first filtered to exclude all genes that showed no expression in any of the samples. The ratio (cancerous/non-cancerous, HCC/CH or HCC/LC) was calculated by comparing the mean expression value for each gene in the cancerous sample set against the mean expression value of that gene in the non-cancerous sample set. Genes were included in the analysis if they had a fold change ≧1.5 in either direction, and a p-value <0.0007 as determined by an Analysis of Variance Test (ANOVA). According to the criteria of the test, differences having p-values below 0.0007 were determined to be statistically significant. Out of the ˜60,000 genes surveyed by the Human Genome U95 set, 109 genes were present in the overall fold change analysis. In Table 1, numbers representing a comparison, or fold change, between the level of expression of a gene in two disease state liver biopsy samples can be positive or negative. Positive values indicate a higher expression level in the cancerous sample compared to the non-cancerous sample (up-regulation), while negative values indicate a lower expression level in the cancerous sample compared to the non-cancerous sample (down-regulation).
Using the above described methods, genes that were predominantly over-expressed in liver cancer, or predominantly under-expressed in liver cancer, were identified. Genes with consistent differential expression patterns provide potential targets for broad range diagnostics and therapeutics.
Table 1 lists the genes determined to be differentially expressed in cancerous liver tissue compared to non-cancerous liver tissue, with the fold change value for each gene. More specifically, the level of expression of the genes of Table 1 in liver cancer cells was compared to the level of expression in tissue from inflamed and/or fibrotic areas of the liver. The set of genes in each group, along with their relative expression levels, creates a profile for the diseases examined, chronic hepatitis with hepatic carcinoma and cirrhosis with hepatic carcinoma.
These genes or subsets of these genes confirm an overall liver disease gene expression profile. The genes in Table 1 may be used alone, or in combination with the methods, compositions, databases and computer systems of the invention.
Table 1 lists the members of diagnostic subsets of genes selected by p-value in groups 2 and 3 (12 and 74 genes, respectively). In addition to their diagnostic, monitoring, drug screening and therapeutic uses, these groups of genes can be used to differentiate between liver tumor samples from subjects with chronic hepatitis and liver tumor samples from subjects with cirrhosis. Assays measuring the expression level of these genes are capable of distinguishing between carcinomas arising in chronic hepatitis patients versus carcinomas arising in cirrhosis patients.
The gene subsets of Table 1 can, therefore, be used to identify the presence of a malignant tumor in liver tissue from chronic hepatitis or cirrhosis patients, to monitor the progression of the tumor (e.g., during cancer treatment or combined disease treatments), to evaluate the effects of therapeutic agents for treating the tumor or to distinguish the origin or predisposing condition of the tumor.
Although the present invention has been described in detail with reference to examples above, it is understood that various modifications can be made without departing from the spirit of the invention. Accordingly, the invention is limited only by the following claims. All cited patents and publications referred to in this application are herein incorporated by reference in their entirety.
This application claims the benefit of U.S. Provisional Applications 60/341,815 and 60/343,185, both of which are herein incorporated by reference in their entirety.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/US02/40718 | 12/20/2002 | WO | 00 | 2/24/2005 |
Number | Date | Country | |
---|---|---|---|
60343185 | Dec 2001 | US | |
60341815 | Dec 2001 | US |