Patent Publication Number: US-2003224382-A1

Title: Reduced complexity nucleic acid targets and methods of using same

Description:
[0001] This application claims the benefit of priority of provisional application serial No. 60/083,331, filed Apr. 27, 1998, No. 60/098,070, filed Aug. 27, 1998, and No. 60/118,624, filed Feb. 4, 1999, each of which is incorporated herein by reference. 
    
    
     [0002] This invention was made with government support under grant number CA68822, NS33377, AI34829 awarded by the National Institutes of Health and under grant number BC961294 awarded by the Department of Defense. The government has certain rights in the invention. 
    
    
     
       BACKGROUND OF THE INVENTION  
       [0003] The present invention relates generally to methods of measuring nucleic acid molecules in a target and more specifically, to methods of detecting differential gene expression.  
       [0004] Every living organism requires genetic material, deoxyribonucleic acid (DNA), which contains genes that impart a unique collection of characteristics to the organism. DNA is composed of two strands of complementary sequences of nucleotide building blocks. The two strands bind, or hybridize, with the complementary sequence to form a double helix. Genes are discreet segments of the DNA and provide the information required to generate a new organism and to give that organism its unique characteristics. Even simple organisms, such as bacteria, contain thousands of genes, and the number is many fold greater in complex organisms such as humans. Understanding the complexities of the development and functioning of living organisms requires knowledge of these genes.  
       [0005] For many years, scientists have searched for and identified a number of genes important in the development and function of living organisms. The search for new genes has greatly accelerated in recent years due to directed projects aimed at identifying genetic information with the ultimate goal being the determination of the entire genome of an organism and its encoded genes, termed genomic studies. One of the most ambitious of these genomic projects has been the Human Genome Project, with the goal of sequencing the entire human genome. Recent advances in sequencing technology have led to a rapid accumulation of genetic information, which is available in both public and private databases. These newly discovered genes as well as those genes soon to be discovered provide a rich resource of potential targets for the development of new drugs.  
       [0006] Despite the rapid pace of gene discovery, there remains a formidable task of characterizing these genes and determining the biological function of these genes. The characterization of newly discovered genes is often a time consuming and laborious undertaking, sometimes taking years to determine the function of a gene or its gene product, particularly in complex higher organisms.  
       [0007] Another level of complexity arises when complex interactions between genes and their gene products are contemplated. To understand how an organism works, it is important not only to understand what role a gene, its transcript and its gene product plays in the workings of an organism, it is also important to understand potentially complex interactions between the gene, its transcript, or its gene product and other genes and their gene products.  
       [0008] A number of approaches have been used to assess gene expression in a particular cell or tissue of an organism. These approaches have been used to characterize gene expression under various conditions, including looking at differences in expression under differing conditions. However, most of these methods are useful for detecting transcripts that are abundant transcripts but have proven less useful for detecting transcripts that are of low abundance, particularly when looking at the expression of a number of genes rather than a selected few genes. Since genes expressed at low levels often regulate the physiological pathways in a cell, it is desirable to detect transcripts having at low abundance.  
       [0009] Thus, a need exists for a method to characterize the expression pattern of genes under a given set of conditions and to detect low abundance transcripts. The present invention satisfies this need and provides related advantages as well.  
       SUMMARY OF THE INVENTION  
       [0010] The invention provides a method of measuring the level of two or more nucleic acid molecules in a target by contacting a probe with a target comprising two or more nucleic acid molecules, wherein the nucleic acid molecules are arbitrarily sampled and wherein the arbitrarily sampled nucleic acid molecules comprise a subset of the nucleic acid molecules in a population of nucleic acid molecules; and detecting the amount of specific binding of the target to the probe. The invention also provides a method of measuring the level of two or more nucleic acid molecules in a target by contacting a probe with a target comprising two or more nucleic acid molecules, wherein the nucleic acid molecules are statistically sampled and wherein the statistically sampled nucleic acid molecules comprise a subset of the nucleic acid molecules in a population of nucleic acid molecules; and detecting the amount of specific binding of the target to the probe. 
     
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
     [0011]FIG. 1 shows differential hybridization to clone arrays. Each image is an autoradiogram that spans about 4000 double spotted  E. coli  colonies, each carrying a different EST clone. Panel A shows the binding of a total target made from 1 μg of polyA +  RNA from confluent human keratinocytes that was radiolabeled during reverse transcription. Panels B and C show RAP-PCR fingerprint with a pair of arbitrary primers that was performed on cDNA from oligo(dT) primed cDNA of confluent human keratinocytes that were untreated (Panel B) and treated with epidermal growth factor (EGF) (Panel C). The two radiolabeled colonies from one differentially expressed cDNA are indicated with an arrow. Panel D shows a RAP-PCR fingerprint with a different pair of arbitrary primers that was performed on RNA from confluent human keratinocytes.  
     [0012]FIG. 2 shows RAP-PCR fingerprints resolved on a polyacrylamide-urea gel. Reverse transcription was performed with an oligo-dT primer on 250, 125, 62.5 and 31.25 ng RNA in lanes 1, 2, 3, and 4 respectively. RNA was from untreated, TGF-β and EGF treated HaCaT cells, as indicated. RAP-PCR was performed with two sets of primers, primers GP14 and GP16 (Panel A) or Nucl+ and OPN24 (Panel B). Molecular weight markers are indicated on the left of each panel, and the sizes of the two differentially amplified RAP-PCR-products are indicated with arrows (317 and 291).  
     [0013]FIG. 3 shows hybridization of targets generated by RAP-PCR to arrays. Shown are autoradiograms of the bottom half of duplicates of the same filter (Genome Systems) hybridized with radiolabeled DNA. Panels A and B show hybridization of two RAP-PCR reactions generated using the same primers and derived from untreated (Panel A) and EGF treated (Panel B) HaCaT cells. Three double-spotted clones that show differential hybridization signals are marked on each array. The GenBank accession numbers of the clone and the corresponding genes are H10045 and H10098, corresponding to vav-3 and AF067817 (square); H28735, gene unknown, similar to $heparan sulfate 3-O-sulfotransferase-1, AF019386 (circle); R48633, gene unknown (diamond). Panel C shows an array hybridized with a RAP-PCR target generated using the same RNA as in panel A but with a different pair of primers. Panel D shows an array hybridized with cDNA target generated by reverse transcription of 1 μg poly(A) + -selected mRNA. Panel E shows an array hybridized with human genomic DNA labeled using random priming.  
     [0014]FIG. 4 shows resolution of RT-PCR products on polyacrylamide-urea gels and confirmation of differential regulation in response to EGF using low stringency RT-PCR. Reverse transcription was performed at two RNA concentrations (500 ng, left column; 250 ng, right column) at different cycle numbers. Shown are bands for the control (22 cycles); for GenBank accession number H11520 (22 cycles); for TSC-22, corresponding to GenBank accession numbers H11073 and H11161 (19 cycles); and for R48633 (19 cycles).  
     [0015]FIG. 5 shows differential display of untreated and EGF treated HaCaT cells. Panel A shows differential display reactions performed at four different starting concentrations of total RNA (designated 1, 2, 3 and 4 and corresponding to 800, 400, 200 and 100 ng, respectively), which was then used for PCR. An anchored oligo(dT) primer, H-T 11 C or H-T 11 A, was used in combination with one of two different arbitrary primers, H-AP3 or H-AP4, which are indicated above the lanes. Panel B shows differential display using the arbitrary primer KA2 with three different anchored oligo(dT) primers, T 13 V, AT 15 A and GT 15 G, used at four different starting concentrations of RNA (designated 1, 2, 3 and 4 and corresponding to 1000, 500, 250 and 125 ng, respectively), which was then used for PCR.  
     [0016]FIG. 6 shows hybridization of differential display reactions to cDNA arrays. Differential display products generated with the primers GT 15 G and KA2 from untreated (Panel A) and EGF treated (Panel B) HaCaT cells were labeled by random priming and hybridized to cDNA arrays. A section representing less than 5% of a membrane is shown with a differentially regulated gene indicated by an arrow. Panel C shows hybridization of differential display products generated with the primers AT 15 A and KA2 from untreated HaCaT cells.  
     [0017]FIG. 7 shows confirmation of differential regulation of genes by EGF using low stringency RT-PCR. Reverse transcription was performed at twofold different RNA concentrations, and low stringency PCR was performed at different cycle numbers. The amount of input RNA used for initial first strand cDNA synthesis and used in each RAP-PCR reaction was 125 ng, left column and 250 ng, right column. The RT-PCR products from 19 cycle reactions were resolved on polyacrylamide-urea gels. Shown are the products for the control (unregulated) and genes exhibiting≧1.6-fold regulation in response to EGF, corresponding to GenBank accession numbers R72714, H14529, H27389, H05545, H27969, R73247, and H21777.  
     [0018]FIG. 8 shows the nucleotide sequence for GenBank accession number H11520 (SEQ ID NO:1).  
     [0019]FIG. 9 shows the nucleotide sequence for GenBank accession number H11161 (SEQ ID NO:2).  
     [0020]FIG. 10 shows the nucleotide sequence for GenBank accession number H11073 (SEQ ID NO:3).  
     [0021]FIG. 11 shows the nucleotide sequence for GenBank accession number U35048 (SEQ ID NO:4).  
     [0022]FIG. 12 shows the nucleotide sequence for GenBank accession number R48633 (SEQ ID NO:5).  
     [0023]FIG. 13 shows the nucleotide sequence for GenBank accession number H28735 (SEQ ID NO:6).  
     [0024]FIG. 14 shows the nucleotide sequence for GenBank accession number AF019386 (SEQ ID NO:7).  
     [0025]FIG. 15 shows the nucleotide sequence for GenBank accession number H25513 (SEQ ID NO:8).  
     [0026]FIG. 16 shows the nucleotide sequence for GenBank accession number H25514 (SEQ ID NO:9).  
     [0027]FIG. 17 shows the nucleotide sequence for GenBank accession number M13918 (SEQ ID NO:10).  
     [0028]FIG. 18 shows the nucleotide sequence for GenBank accession number H12999 (SEQ ID NO:11).  
     [0029]FIG. 19 shows the nucleotide sequence for GenBank accession number H05639 (SEQ ID NO:12).  
     [0030]FIG. 20 shows the nucleotide sequence for GenBank accession number L49207 (SEQ ID NO:13).  
     [0031]FIG. 21 shows the nucleotide sequence for GenBank accession number H15184 (SEQ ID NO:14).  
     [0032]FIG. 22 shows the nucleotide sequence for GenBank accession number H15124 (SEQ ID NO:15) FIG. 23 shows the nucleotide sequence for GenBank accession number X79781 (SEQ ID NO:16).  
     [0033]FIG. 24 shows the nucleotide sequence for GenBank accession number H25195 (SEQ ID NO:17).  
     [0034]FIG. 25 shows the nucleotide sequence for GenBank accession number H24377 (SEQ ID NO:18).  
     [0035]FIG. 26 shows the nucleotide sequence for GenBank accession number M31627 (SEQ ID NO:19).  
     [0036]FIG. 27 shows the nucleotide sequence for GenBank accession number H23972 (SEQ ID NO:20).  
     [0037]FIG. 28 shows the nucleotide sequence for GenBank accession number H27350 (SEQ ID NO:21).  
     [0038]FIG. 29 shows the nucleotide sequence for GenBank accession number AB000712 (SEQ ID NO:22).  
     [0039]FIG. 30 shows the nucleotide sequence for GenBank accession number R75916 (SEQ ID NO:23).  
     [0040]FIG. 31 shows the nucleotide sequence for GenBank accession number X85992 (SEQ ID NO:24).  
     [0041]FIG. 32 shows the nucleotide sequence for GenBank accession number R73021 (SEQ ID NO:25).  
     [0042]FIG. 33 shows the nucleotide sequence for GenBank accession number R73022 (SEQ ID NO:26).  
     [0043]FIG. 34 shows the nucleotide sequence for GenBank accession number U66894 (SEQ ID NO:27).  
     [0044]FIG. 35 shows the nucleotide sequence for GenBank accession number H10098 (SEQ ID NO:28).  
     [0045]FIG. 36 shows the nucleotide sequence for GenBank accession number H10045 (SEQ ID NO:29).  
     [0046]FIG. 37 shows the nucleotide sequence for GenBank accession number AF067817 (SEQ ID NO:30).  
     [0047]FIG. 38 shows the nucleotide sequence for GenBank accession number R72714 (SEQ ID NO:31).  
     [0048]FIG. 39 shows the nucleotide sequence for GenBank accession number X52541 (SEQ ID NO:32).  
     [0049]FIG. 40 shows the nucleotide sequence for GenBank accession number H14529 (SEQ ID NO:33).  
     [0050]FIG. 41 shows the nucleotide sequence for GenBank accession number M10277 (SEQ ID NO:34).  
     [0051]FIG. 42 shows the nucleotide sequence for GenBank accession number H27389 (SEQ ID NO:35).  
     [0052]FIG. 43 shows the nucleotide sequence for GenBank accession number D89092 (SEQ ID NO:36).  
     [0053]FIG. 44 shows the nucleotide sequence for GenBank accession number D89678 (SEQ ID NO:37).  
     [0054]FIG. 45 shows the nucleotide sequence for GenBank accession number H05545 (SEQ ID NO:38).  
     [0055]FIG. 46 shows the nucleotide sequence for GenBank accession number J03804 (SEQ ID NO:39).  
     [0056]FIG. 47 shows the nucleotide sequence for GenBank accession number H27969 (SEQ ID NO:40).  
     [0057]FIG. 48 shows the nucleotide sequence for GenBank accession number R73247 (SEQ ID NO:41).  
     [0058]FIG. 49 shows the nucleotide sequence for GenBank accession number U51336 (SEQ ID NO:42).  
     [0059]FIG. 50 shows the nucleotide sequence for GenBank accession number H21777 (SEQ ID NO:43).  
     [0060]FIG. 51 shows the nucleotide sequence for GenBank accession number K00558 (SEQ ID NO:44).  
     [0061]FIG. 52 shows the nucleotide sequence for GenBank accession number D31765 (SEQ ID NO:45). 
    
    
     DETAILED DESCRIPTION OF THE INVENTION  
     [0062] The invention provides methods for measuring the level of two or more nucleic acid molecules in a target by contacting a probe with an arbitrarily sampled target or a statistically sampled target and detecting the amount of specific binding to the probe. The invention also provides methods of identifying two or more differentially expressed nucleic acid molecules associated with a condition by measuring the level of two or more nucleic acid molecules in a target and comparing the expression levels to expression levels of the nucleic acid molecules in a second target. The methods of the invention are useful for obtaining a profile of nucleic acid molecules expressed in a target under a given set of conditions. The methods of the invention are particularly useful for comparing the relative abundance of low abundance nucleic acid molecules between two or more targets. The methods of the invention are advantageous in that a profile of nucleic acid molecule abundance can be determined and correlated with a given set of conditions or compared to another target to determine if the original target was exposed to a particular set of conditions, thereby providing information useful for assessing the diagnosis or treatment of a disease.  
     [0063] The invention provides a method of measuring the abundance of two or more nucleic acid molecules in a target. The method of the invention includes the steps of contacting a probe with a target comprising two or more nucleic acid molecules, wherein the nucleic acid molecules are arbitrarily sampled and wherein the arbitrarily sampled nucleic acid molecules comprise a subset of the nucleic acid molecules in a population of nucleic acid molecules; and detecting the amount of specific binding of the target to the probe.  
     [0064] As used herein, the term “nucleic acid molecule” refers to a nucleic acid of two or more nucleotides. A nucleic acid molecule can be RNA or DNA. For example, a nucleic acid molecule can include messenger RNA (mRNA), transfer RNA (tRNA) or ribosomal RNA (rRNA). A nucleic acid molecule can also include, for example, genomic DNA or cDNA. A nucleic acid molecule can be synthesized enzymatically, either in vivo or in vitro, or the nucleic acid molecule can be chemically synthesized by methods well known in the art. A nucleic acid molecule can also contain modified bases, for example, the modified bases found in tRNA such as inosine, methylinosine, dihyrouridine, ribothymidine, pseudouridine, methylguanosine and dimethylguanosine. Furthermore, a chemically synthesized nucleic acid molecule can incorporate derivatives of nucleotide bases.  
     [0065] As used herein, the term “population of nucleic acid molecules” refers to a group of two or more different nucleic acid molecules. A population of nucleic acid molecules can also be 3 or more, 5 or more, 10 or more, 20 or more, 50 or more, 100 or more, 1000 or more or even 10,000 or more different nucleic acid molecules. The nucleic acid molecules can differ, for example, by a single nucleotide or by modification of a single base. Generally, a population of nucleic acid molecules is obtained from a target sample, for example, a cell, tissue or organism. In such a case, the population of nucleic acid molecules contains the nucleic acid molecules of the target sample.  
     [0066] A population of nucleic acid molecules has characteristics that can differentiate one population of nucleic acid molecules from another. These characteristics are based on the number and nature of individual nucleic acid molecules comprising the population. Such characteristics include, for example, the abundance of nucleic acid molecules in the population. The abundance of an individual nucleic acid molecule can be an absolute amount in a given target sample or can be the amount relative to other nucleic acid molecules in the target sample. In a population of nucleic acid molecules obtained from a target, individual nucleic acid molecules can be more abundant or less abundant relative to other nucleic acid molecules in the sample target. A less abundant sequence can also be relative abundance between two samples.  
     [0067] As used herein, a less abundant nucleic acid molecule can be, for example, less than about 10% as abundant as the most abundant nucleic acid molecule in a population. A less abundant nucleic acid molecule can also be less than about 1% as abundant, less than about 0.1% as abundant or less than about 0.01% as abundant as the most abundant nucleic acid molecule in a population. For example, a low abundance nucleic acid molecule can be less than about 10 copies per cell, or even as low as 1 copy per cell.  
     [0068] Another characteristic of a population of nucleic acid molecules is the complexity of the population. As used herein, “complexity” refers to the number of nucleic acid molecules having different sequences in the population. For example, a population of nucleic acid molecules representative of the mRNA in a bacterial cell has lower complexity than a population of nucleic acid molecules representative of the mRNA in a eukaryotic cell, a tissue or an organism because a smaller number of genes are expressed in a bacterial cell relative to a eukaryotic cell, tissue or organism.  
     [0069] A population of nucleic acid molecules can also be characterized by the properties of individual nucleic acid molecules in the population. For example, the length of individual nucleic acid molecules contributes to the characteristics of a population of nucleic acid molecules. Similarly, the sequence of individual nucleic acid molecules in the population contributes to the characteristics of the population of nucleic acid molecules, for example, the G+C content of the nucleic acid sequences and any secondary structure that can form due to complementary stretches of nucleotide sequence that can undergo intrastrand hybridization.  
     [0070] As used herein, the term “subset of nucleic acids” means less than all of a set of nucleic acid molecules. For example, a subset of nucleic acid molecules of a target sample population would be less than all of the nucleic acid molecules in the target sample population. Specifically excluded from a subset of nucleic acid molecules is a group of nucleic acid molecules representative of all the nucleic acid molecules in a sample target, for example, a target generated using total cDNA or total mRNA.  
     [0071] As used herein, the term “target” refers to one or more nucleic acid molecules to which binding of a probe is desired. A target is detectable when bound to a probe. A target of the invention generally comprises two or more different nucleic acid molecules. A target can be derived from a population of nucleic acid molecules from a cell, tissue or organism. A target can also contain 3 or more, 5 or more, 10 or more, 20 or more, 30 or more, 50 or more, 100 or more, 200 or more, 500 or more, 1000 or more, 2000 or more, 5000 or more, or even 10,000 or more different nucleic acid molecules. A target can have a detectable moiety associated with it such as a radioactive label, a fluorescent label or any label that is detectable. When a target is labeled, for example, with a radioactive label, the target can be used “to probe” or hybridize with other nucleic acid molecules. Methods of making a target are disclosed herein.  
     [0072] A method of detection that directly measures binding of the target to a probe, without the need for a detectable moiety attached to the target, can also be used. In such a case, the nucleic acid molecules are directly detectable without modification of a nucleic acid molecule of the target, for example, by attaching a detectable moiety. An example of such a detection method using a target without a detectable moiety is detection of binding of a target using mass spectrometry. Another example of a method using a target containing nucleic acid molecules without an attached detectable moiety is binding the target to a probe that contains molecules having a detectable moiety. In such a case, the binding of a target to the probe containing molecules having a detectable moiety is detected and, as such, the target is detectable when bound to the probe. An example is the “molecular beacon,” where probe binding causes separation of a fluorescent tag from a fluorescence quencher.  
     [0073] As used herein, the term “specific binding” means binding that is measurably different from a non-specific interaction. Specific binding can be measured, for example, by determining binding of a molecule compared to binding of a control molecule, which generally is a molecule of similar structure that does not have binding activity. For example, specific binding of a target to a probe can be determined by comparing binding of the target with binding control nucleic acids not included in the target. Specific binding can also be determined by competition with a control molecule that is similar to the target, for example, an excess of non-labeled target. In this case, specific binding is indicated if the binding of a labeled target to a probe is competitively inhibited by excess unlabeled target.  
     [0074] The term “specific binding,” as used herein, includes both low and high affinity specific binding. Specific binding can be exhibited, for example, by a low affinity molecule having a Kd of at least about 10 −4  M. Specific binding also can be exhibited by a high affinity molecule, for example, a molecule having a Kd of at least about of 10 −7  M, at least about 10 −8  M, at least about 10 −9  M, at least about 10 −10  M, or can have a Kd of at least about 10 −11  M or 10 −12  M or greater.  
     [0075] In the case of a probe comprising an array of nucleic acid molecules, binding of a specific nucleic acid molecule of the probe to another nucleic acid molecule is also known as hybridizing or hybridization. As used herein, the term “hybridizing” or “hybridization” refers to the ability of two strands of nucleic acid molecules to hydrogen bond in a sequence dependent manner. Under appropriate conditions, complementary nucleotide sequences can hybridize to form double stranded DNA or RNA, or a double stranded hybrid of RNA and DNA. Nucleic acid molecules with similar but non-identical sequences can also hybridize under appropriate conditions.  
     [0076] As used herein, the term “probe” refers to a population of two or more molecules to which binding of a target is desired. The molecules of a probe include nucleic acid molecules, oligonucleotides and polypeptide-nucleic acid molecules. A probe can additionally be an array of molecules.  
     [0077] In general, a probe is comprised of molecules immobilized on a solid support and the target is in solution. However, it is understood that a target can be bound to a solid support and a probe can be in solution. Furthermore, both the probe and the target can be in solution. It is understood that the configuration of the probe and target can be in solution or bound to a solid support, so long as the probe and target can bind to each other. When bound to a solid support, the binding of the probe or target to the support can be covalent or non-covalent, so long as the bound probe or target remains bound under conditions of contacting the solid support with a probe or target in solution and washing of the solid support. If the probe and target hybridize or otherwise specifically interact, the probe or target bound to a solid support remains bound during the hybridization and washing steps.  
     [0078] As used herein, the term “sampled” or “samples,” when used in reference to a nucleic acid molecule, refers to a nucleic acid molecule to which specific binding can be detected. A nucleic acid molecule that samples another molecule is capable of specifically binding to that molecule and being detected. For example, a probe can sample molecules in a target by detectably binding to molecules in the target. Those molecules in the target to which nucleic acid molecules in the probe specifically bind are therefore sampled.  
     [0079] As used herein, the term “arbitrarily sampled” or “arbitrarily sampled nucleic acid molecule” means that a nucleic acid molecule is sampled by binding based on its sequence without sampling based on a particular site where a molecule will bind. When generating a target comprising arbitrarily sampled nucleic acid molecules from a population of nucleic acid molecules, the target is generated without prior reference to the sequences of nucleic acid molecules in the population. Thus, it is not necessary to have previous knowledge of the nucleotide sequence of nucleic acid molecules in the population to arbitrarily sample the population. It is understood that knowledge of a nucleotide sequence of a nucleic acid molecule in the population does not preclude the ability to arbitrarily sample the population so long as the nucleotide sequence is not referenced before sampling the population. Methods for generating a probe containing arbitrarily sampled nucleic acid molecules are disclosed herein (see below and Examples I to III).  
     [0080] An arbitrarily sampled probe containing arbitrarily sampled nucleic acid molecules can be generated using one or more arbitrary oligonucleotides. As used herein, the term “arbitrary oligonucleotide” means that the oligonucleotide is a sequence that is selected randomly and is not selected based on its complementarity to any known sequence. As such, an arbitrary oligonucleotide can be used to arbitrarily sample a population of nucleic acid molecules.  
     [0081] An arbitrarily sampled nucleic acid molecule is sampled based on its sequence and is not based on binding to a predetermined sequence. For example, arbitrary oligonucleotides are oligonucleotides having an arbitrary sequence and, as such, will bind to a given nucleic acid molecule because the complementary sequence of the arbitrary oligonucleotide occurs by chance in the nucleic acid molecule. Because the oligonucleotides can bind to a nucleic acid molecule based on the presence of a complementary sequence, the sampling of the nucleic acid molecule is based on that sequence. However, the binding of the arbitrary oligonucleotide to any particular nucleic acid molecule in a population is not determined prior to the binding of the oligonucleotide, for example, by comparing the sequence of the arbitrary oligonucleotides to known nucleic acid sequences and selecting the oligonucleotides based on previously known nucleic acid sequences. The use of arbitrary oligonucleotides as primers for amplification is well known in the art (Liang and Pardee,  Science  257:967-971 (1992)).  
     [0082] As used herein, the term “oligonucleotide” refers to a nucleic acid molecule of at least 2 and less than about 1000 nucleotides. An oligonucleotide can be, for example, at least about 5 nucleotides and less than about 100 nucleotides, for example less than about 50 nucleotides.  
     [0083] The invention also provides a method of measuring the level of two or more nucleic acid molecules in a target by contacting a probe with a target comprising two or more nucleic acid molecules, wherein the nucleic acid molecules are statistically sampled and wherein the statistically sampled nucleic acid molecules comprise a subset of the nucleic acid molecules in a population of nucleic acid molecules; and detecting the amount of specific binding of the target to the probe.  
     [0084] As used herein, the term “statistically sampled nucleic acid molecule” means that a nucleic acid sequence is sampled based on its sequence with prior reference to its nucleotide sequence by predetermining the statistical occurrence of a nucleotide sequence in two or more nucleic acid molecules. Thus, to obtain a statistically sampled nucleic acid molecule, it is necessary to have previous knowledge of the nucleotide sequence of at least two nucleic acid molecules in the population.  
     [0085] A statistically sampled nucleic acid molecule is sampled based on the sequence of a nucleic acid molecule with prior reference to its nucleotide sequence but without prior reference to a preselected portion of its nucleotide sequence. A group of oligonucleotides can be identified without prior reference to a preselected portion of a nucleotide sequence, for example, by determining a group of arbitrary oligonucleotides. The arbitrary oligonucleotides can then be referenced to known nucleotide sequences by determining which of the arbitrary primers match the known nucleotide sequences. Such arbitrary oligonucleotides referenced to known nucleotide sequences are selected based on the known sequences and thus become statistical primers. This method is in contrast to a method where a preselected site in a known nucleotide sequence is identified and an oligonucleotide is specifically designed to match that preselected site.  
     [0086] Statistical sampling is advantageous because a set of oligonucleotides can be determined based on the presence in a group of known sequences of a sequence complementary to the oligonucleotides. The oligonucleotides can further be ranked based on complexity binding. Complexity binding means that a given oligonucleotide binds to more than one nucleic acid molecule. The larger the number of molecules to which an oligonucleotide can bind, the higher the “complexity binding.” Statistical selection can be used to enhance for complexity binding by ranking oligonucleotides based on the number of sequences to which the oligonucleotides will bind and selecting those that bind to the highest number (see, for example, WO 99/11823). Statistical sampling can be based, for example, on the binding of an oligonucleotide to 5 or more nucleic acid molecules, and can be based on the binding to 10 or more, 50 or more, 100 or more, 200 or more, 500 or more, 1000 or more, or even 10,000 or more nucleic acid molecules.  
     [0087] In addition, statistical sampling can enhance for the highest complexity binding for a given oligonucleotide, for example, by selecting the above average ranked oligonucleotides that are complementary to above the average number of nucleic acid molecules. The oligonucleotides can be selected for the any range of complexity binding, for example, the top 10% of highest ranked complexity binding, the top 20% of highest ranked complexity binding, or the top 50% of highest ranked complexity binding.  
     [0088] Furthermore, statistical selection can be used to exclude undesirable nucleotide sequences, including conserved sequences in a family of related nucleic acid molecules (WO 99/11823). A statistical oligonucleotide can be about 5 nucleotides in length to about 1000 nucleotides in length, for example, about 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 18, 20, 25, 30 or 50 nucleotides in length. A set of statistical primers can contain degenerate bases, for example, more than one nucleotide at any given position.  
     [0089] A sampled nucleic acid molecule obtained using a preselected portion of a nucleotide sequence is specifically excluded from the meaning of the term “statistically sampled nucleic acid molecule.” For example, if a portion of a known nucleotide sequence is identified and an oligonucleotide that matches the identified portion is generated to sample a nucleic acid molecule, such a sampled nucleic acid molecule would not be a statistically sampled nucleic acid molecule. However, if a group of oligonucleotides is first identified and then compared to two or more known nucleotide sequences in a population of nucleic acid molecules to determine oligonucleotides statistically present in or similar to the known nucleotide sequences, such statistically identified oligonucleotides can be used to obtain a statistically sampled nucleic acid molecule. Methods for generating a target containing statistically sampled nucleic acid molecules are disclosed herein.  
     [0090] A statistically sampled target containing statistically sampled nucleic acid molecules can be generated using one or more statistical oligonucleotides. As used herein, the term “statistical oligonucleotide” means that an oligonucleotide is a sequence that is selected based on its statistical occurrence of complementarity in more than one known nucleic acid molecule. As such, a statistical oligonucleotide can be used to statistically sample a population of nucleic acid molecules.  
     [0091] The methods of the invention detect specific binding of a target to a probe. A target can be generated, for example, by amplifying nucleic acid molecules. As used herein, the term “amplified target” refers to a target generated by enzymatically copying a nucleic acid molecule to generate more than one copy of the nucleic acid molecules in a population of nucleic acid molecules. An amplified nucleic acid target can be generated, for example, using an amplification method such as polymerase chain reaction (PCR). A target having a single copy of each nucleic acid molecule in a target sample from which the target sample is derived, which would have identical abundance and complexity as the original population, would not be considered an amplified target. An amplified target can be useful, for example, if nucleic acid molecules sampled by the probe are in limited quantities in the target. A nucleic acid molecule that is to be sampled and which is present in very low quantities would be difficult to detect without amplification and increasing the mass of the nucleic acid molecules in the probe. However, a limited complexity target, in which the complexity or number of different molecules is limited, need not be amplified.  
     [0092] Other methods for generating an amplified target include, for example, the ligase chain reaction (LCR); self-sustained sequence replication (3SR); beta replicase reaction, for example, using Q-beta replicase; phage terminal binding protein reaction; strand displacement amplification (SDA); nucleic acid sequence based amplification (NASBA); cooperative amplification by cross hybridization (CATCH); rolling circle amplification (RCA) and AFLP (Trippler et al.,  J. Viral. Hepat.  3:267 (1996); Hofler et al.,  Lab. Invest.  73:577 (1995); Tyagi et al.,  Proc. Natl. Acad. Sci. USA  93:5395 (1996); Blanco et al.,  Proc. Natl. Acad. Sci. USA  91:12198 (1994); Spears et al.,  Anal. Biochem.  247:130 (1997); Spargo et al.,  Mol. Cell. Probes  10:247 (1996); Gobbers et al.,  J. Virol. Methods  66:293 (1997); Uyttendaele et al.,  Int. J. Food Microbiol.  37:13 (1997); and Leone et al.,  J. Virol. Methods  66:19 (1997); Ellinger et al.,  Chem. Biol.  5:729-741 (1998); Ehricht et al.,  Nucleic Acids Res.  25:4697-4699 (1997); Ehricht et al.,  Eur. J. Biochem.  243:358-364 (1997); Lizardi et al.,  Nat. Genet.  19:225-232 (1998)).  
     [0093] The methods of the invention are useful for measuring the level of two or more nucleic acid molecules in a target. The methods of the invention can also be used to compare expression levels between two targets. In particular, the methods of the invention are useful for measuring differential expression of nucleic acid molecules (see below).  
     [0094] A total target, using the full complexity of the mRNA population for target preparation, can easily examine the top few hundred or a few thousand of the mRNAs in the cell (Pietu et al.,  Genome Res.  6:492-503 (1996)). However, a total labeled cDNA target from a mammalian cell typically has a complexity of over 100 million bases which complicates attempts to detect differential expression among the rarer mRNAs using differential hybridization. Recent advances in the use of fluorescence and confocal microscopy have led to improvements in the sensitivity and dynamic range of differential hybridization methods, with a dynamic range of detection of 10,000-fold and the detection of transcripts at a sensitivity approaching 1/500,000 (Marshall and Hodgson,  Nat. Biotechnol.  16:27-31 (1998); Ramsay,  Nat. Biotechnol.  16:40-44 (1998)). Despite the improvements in sensitivity, methods using total target remain biased toward more abundant mRNAs in a sample.  
     [0095] The standard method for differential screening, which typically uses targets derived from reverse transcription of total message and autoradiography or phosphoimaging, can be used to detect differential expression (Pietu, supra, 1996). However, the method is limited to the most abundant messages. Only abundant transcripts are represented highly enough to yield effective targets with a sensitivity of perhaps 1/15,000 (Boll,  Gene  50:41-53 (1986)). As disclosed herein, differential screening can be improved greatly by reducing the complexity of the target and by systematically increasing the amount of rarer nucleic acid molecules in the target. By enhancing the amount of less abundant nucleic acids in a target, differential screening is not confined to only the most abundant nucleic acid molecules, as observed using total target.  
     [0096] By reducing the complexity of the target, the ability to identify all mRNA species in a source simultaneously is sacrificed for improved kinetics and an improved signal to noise ratio. Complexity reduction methods generate a target having a subset of nucleic acid molecules in a population that allow a few rare mRNAs to contribute significantly to the final mass of the target, thereby enhancing the ability to observe differential gene expression among rare mRNAs in a source. Any method that generates a mixture of products that reliably enriches for only part of each mRNA or only a subset of the mRNA population is useful for generating a reduced complexity target.  
     [0097] There are two fundamentally different types of complexity reduction methods, methods that maintain the relative stoichiometry among the mRNAs they sample and methods that do not maintain stoichiometry. One class of methods yields nucleic acids representing a subset of the mRNA population and maintains the approximate stoichiometry of the input RNA. Such methods are exemplified by most amplified restriction fragment length polymorphism (AFLP) and restriction strategies that sample the 3′ end or internal fragments of mRNAs (Habu et al.,  Biochem. Biophys. Res. Commun.  234:516-521 (1997); Money et al.,  Nucleic Acids Res.  24:2616-2617 (1996); Bachem et al.,  Plant J.  9:745-753 (1996)). Another example is the use of size fractionated mRNAs to generate cDNA targets. All the mRNAs, for example, the 2.0 to 2.1 kb range can be used as a reduced complexity target. Stoichiometry among these mRNAs would be mostly preserved in the target (Dittmar et al.,  Cell Biol. Int.  21:383-391 (1997)).  
     [0098] A second class of methods for generating reduced complexity targets does not preserve the stoichiometry of the starting mRNAs, though it does preserve differences among individual RNAs between target samples from which targets are made. One method to generate a reduced complexity target that does not maintain stoichiometry is to use subtracted targets, which have shown sensitivity for rare messages comparable to chips, in particular methods based on representational difference analysis or suppression subtractive hybridization (Rhyner et al.,  J. Neurosci. Res.  16:167-181 (1986); Lisitsyn et al.,  Science  259:946-951 (1993); Lisitsyn &amp; Wigler,  Methods Enzymol.  254:291-304 (1995); Jin et al.,  Biotechniques  23:1084-1086 (1997)).  
     [0099] Particularly useful methods for generating a reduced complexity target that does not maintain stoichiometry are exemplified by using arbitrarily sampled targets or statistically sampled targets. Methods using arbitrarily sampled targets and statistically sampled targets are disclosed herein. The methods using arbitrarily sampled or statistically sampled targets allow detection of low abundance nucleic acid molecules in a target. The methods of the invention are advantageous because they enhance the ability to detect low abundance nucleic acid molecules in a target and also allow detection of nucleic acid molecules in a target derived from limited quantities of nucleic acid molecules, such as a few cells or even a single cell.  
     [0100] An arbitrarily sampled target or statistically sampled target can be generated, for example, by amplification. If an amplified target is generated using arbitrary oligonucleotides or statistical oligonucleotides, the amplified products reflect a function of both the starting abundance of each target nucleic acid molecule and the quality of the match of the oligonucleotide to the target nucleic acid molecule to be sampled. Thus, the final mixture of amplified products can include quite abundant amplified products that derive from low abundance nucleic acid molecules that have a good match with the oligonucleotide primers used and have favorable “amplifiability” after the initial priming events. Amplifiability includes effects such as secondary structure and product size.  
     [0101] A consequence of generating an amplified target using arbitrary oligonucleotides or statistical oligonucleotides is that the same nucleic acid molecules in two different targets experience an identical combination of primability and amplifiability so that changes in abundance for particular mRNAs are maintained, even as the relative abundances between different nucleic acid molecules within one target are profoundly changed. This is in contrast to methods that maintain stoichiometry, where less abundant nucleic acid molecules would be present as less abundant nucleic acid molecules in the target.  
     [0102] When generating an amplified target, there are generally no particular constraints on the oligonucleotide primers. The oligonucleotide primers preferably contain at least a few C or G bases. The oligonucleotide primers also preferably do not contain 3′ ends complementary with themselves or the other primer in the reaction, to avoid primer dimers. The oligonucleotide primers are also preferably chosen to have different sequences so that the same parts of mRNA are not amplified in different fingerprints.  
     [0103] As disclosed herein, methods of generating arbitrarily sampled targets or statistically sampled targets can be based on methods that have been traditionally used to “fingerprint” a target sample containing nucleic acid molecules. The fingerprints are characteristic of the expression of nucleic acid molecules in a target sample. To generate an arbitrarily sampled target, one method that can be used is based on RNA arbitrarily primed PCR (RAP-PCR) (see Examples I and II; Welsh et al.,  Nucleic Acids Res.  18:7213-7218 (1990); Welsh et al.,  Nucleic Acids Res.  20:4965-4970 (1992); Liang and Pardee,  Science  257:967-971 (1992)).  
     [0104] In RAP-PCR, both the abundance and the extent of match with the primers contribute to the prevalence of any particular product. Thus, rare mRNAs that happen to have excellent matches with the primers and are efficiently amplified are found among the more abundant RAP-PCR products, which makes a target generated by RAP-PCR non-stoichiometric. This is a very useful feature of RAP-PCR because it allows the sampling of mRNAs that are difficult to sample using other methods.  
     [0105] In a typical RAP-PCR fingerprint, about 50-100 cDNA fragments per lane are visible on a polyacrylamide gel, including products from relatively rare mRNAs that happen to have among the best matches with the arbitrary primers. If only 100 cDNA clones could be detected in an array by each target, then hybridization to arrays would be inefficient. However, RAP-PCR fingerprints contain many products that are too rare to visualize by autoradiography of a polyacrylamide gel. Nonetheless, these rarer products are reproducible and of sufficient abundance to serve as target for arrays when labeled at high specific activity.  
     [0106] As disclosed herein, a single target derived from RAP-PCR can detect about a thousand cDNAs on an array containing about 18,000 EST clones, a 10-20 fold improvement over the performance of fingerprints displayed on denaturing polyacrylamide gels. In addition, when a differentially regulated gene is detected on a cDNA array, a clone representing the transcript is immediately available, and often sequence information for the clone is also available. Furthermore, the clones are usually much longer than the usual RAP-PCR product. In contrast, the standard approaches to RNA fingerprinting require that the product be gel purified and sequenced before verification of differential expression can be performed. As disclosed herein, differentially amplified RAP-PCR products that are below the detection capabilities of the standard denaturing polyacrylamide gel and autoradiography methods can be detected using hybridization to cDNA arrays.  
     [0107] An arbitrarily sampled target generated by RAP-PCR can sample the top few thousand highest expressed nucleic acid molecules in a target sample and can sample different subsets of the nucleic acid molecules in a population, depending on the oligonucleotide primers used for amplification. Some of the rare nucleic acid molecules in a target are sufficiently represented to be easily detected on arrays of colonies (see Examples I and II).  
     [0108] To generate an arbitrarily sampled target using RAP-PCR, the RAP-PCR fingerprint is made by arbitrarily primed reverse transcription and PCR of nucleic acid molecules in a target sample, for example, messenger RNA (McClelland et al., in  Differential Display Methods and Protocols,  Liang and Pardee, eds., Humana Press (1997)). Alternatively, first strand cDNA can be primed with oligo dT or with random short oligomers, followed by arbitrary priming. Analysis of such a RAP-PCR “fingerprint” by gel electrophoresis reveals a complex fingerprint showing relative abundances of an arbitrary sample of about 100 transcripts (see Example II).  
     [0109] As disclosed herein, RAP-PCR fingerprints were converted to targets- to probe or hybridize human cDNA clones arrayed as  E. coli  colonies on nylon membranes (Example II). Each array contained 18,432 cDNA clones from the Integrated Molecular Analysis of Genomes and their Expression (I.M.A.G.E.) consortium. Hybridization to about 1000 cDNA clones was detected using each arbitrarily sampled target generated by RAP-PCR. Different RAP-PCR fingerprints gave hybridization patterns having very little overlap (&lt;3%) with each other, or with hybridization patterns from total cDNA targets. Consequently, repeated application of RAP-PCR targets allows a greater fraction of the message population to be screened on this type of array than can be achieved with a radiolabeled total cDNA target.  
     [0110] The arbitrarily sampled targets were generated from HaCaT keratinocytes treated with EGF. Two RAP-PCR targets hybridized to 2000 clones, from which 22 candidate differentially expressed genes were observed (Example II). Differential expression was tested for 15 of these clones using RT-PCR and 13 were confirmed. The use of this cDNA array to analyze RAP-PCR fingerprints allowed for an increase in detection of 10- to 20-fold over the conventional denaturing polyacrylamide gel approach to RAP-PCR or differential display. Throughput is vastly improved by the reduction in cloning and sequencing afforded by the use of arrays. Also, repeated cloning and sequencing of the same gene, or of genes already known to be regulated in the system of interest, is minimized.  
     [0111] The use of RAP-PCR to generate an arbitrarily sampled target is particularly useful because it allows very high throughput discovery of differentially regulated genes (see Examples II and III). The throughput using this method is about 20 times faster. Essentially, once a RAP-PCR fingerprint has been generated, instead of analyzing the product by gel electrophoresis, the RAP-PCR fingerprint is used as a target to probe or hybridize to nucleic acid molecules. Such an arbitrarily sampled target generated by RAP-PCR is particularly useful as a target for an array.  
     [0112] Parameters of the RAP-PCR reaction can be varied, for example, to optimize complexity of the target and enhance complexity binding. For example, to increase the complexity, Taq polymerase Stoffel fragment, which is more promiscuous than AMPLITAQ, can be used for amplification. The oligonucleotide primers used herein (Example II) were 10 or 11 bases in length and were not degenerate, having a single base at each position. Longer oligonucleotide primers used at the same temperature can give a more complex product, as would primers with some degeneracy. However, the greater the complexity of the target, the more closely it will resemble a total mRNA target, which loses the advantage of non-stoichiometric sampling. To further vary RAP-PCR parameters, the oligonucleotide primer length, degeneracy, and 3′ anchoring can be varied in the reverse transcription and PCR reactions. Various different polymerases can also be used.  
     [0113] The RAP-PCR fingerprint can be radiolabeled or labeled with fluorescent dyes, as described below, and used as a target to probe against dense arrays such as arrays of cDNA clones. Differences in the level of nucleic acid molecules between two targets can indicate, for example, differences in mRNA transcript levels, which usually reflects differences in gene expression levels. Differences in expression can also reflect degradation or post-translational processsing. Using an arbitrarily sampled target, each target is estimated to allow the detection of roughly 10% of the total complexity of the message population, and most importantly, this 10% very effectively includes the rare message class. The rare message class is included in the target because, while RAP-PCR reflects message abundance between target samples, the cDNAs selected for amplification in any particular RAP-PCR reaction is determined by sequence rather than abundance. When the sequence match between oligonucleotide primers and nucleic acid molecules is very good, even if the nucleic acid molecule is in low abundance, the low abundance nucleic acid molecules have a good chance of having a larger amount of the less abundant nucleic acid molecule relative to more abundant nucleic acid molecules in the final target.  
     [0114] To be suitable for either gel- or array-based analysis, RAP-PCR fingerprints should remain almost identical over an eight-fold dilution of the input RNA. Low quality RAP-PCR fingerprints are usually the consequence of poor control over RNA quality and concentration. Before proceeding with the array hybridization steps, the quality of the RAP-PCR products can be verified. Because the array method has such high throughput, this extra step is neither costly, nor time-consuming, and can greatly improve efficiency by reducing the number of false positives due to poor fingerprint reproducibility. The reproducibility of RAP-PCR fingerprints as targets is exemplified herein (see Example II).  
     [0115] The enhanced ability of the methods of the invention to detect low abundance nucleic acid molecules in a target sample provides a major improvement over previously used methods that have limited ability to detect rare messages. It is likely that the entire complexity of the message population of a cell could be examined in a short period of time, for example, in a few weeks.  
     [0116] For example, as disclosed in Example II, targets generated by RAP-PCR sample a population of mRNAs largely independent of message abundance. This is because the low abundance class of messages has much higher complexity than the abundant class, making it more likely that the arbitrary primers will find good matches. Unlike differential display, RAP-PCR demands two arbitrary priming events, possibly biasing RAP-PCR toward the complex class. It is likely that the majority of the mRNA population in a cell (&lt;20,000 mRNAs) can be found in as few as ten RAP-PCR fingerprints.  
     [0117] In addition to using RAP-PCR, differential display can also be used to generate an arbitrarily sampled target (see Example III). For differential display, first, reverse transcription uses a 3′ anchored primer such as an oligo(dT) primer. Next, second strand cDNA is primed with an arbitrary primer. Then PCR takes place between the arbitrary primer and the 3′ anchor.  
     [0118] As disclosed in Example III, a combination of one arbitrary and one oligo(dT) anchor primer was used to generate an arbitrarily sampled target for cDNA arrays. Both the RAP-PCR and differential display approaches to target preparation can use less than {fraction (1/200)}th of the amount of RNA used in some other array hybridization methods. Each fingerprint detected about 5-10% of the transcribed mRNAs, sampled almost independent of abundance, using inexpensive  E. coli  colony arrays of EST clones. The differential display protocol was modified to generate a sufficient mass of PCR products for use as a target to probe nucleic acid molecules. The use of different oligo(dT) anchor primers with the same arbitrary primer resulted in considerable overlap among the genes sampled by each target. Overlap of sampled genes can be avoided by using different arbitrary primers with each oligo(dT) anchor primer. Four genes not previously known to be regulated by EGF and three genes known to be regulated by EGF in other cell types were characterized using the arbitrarily sampled targets generated by differential display. The use of arbitrarily sampled targets generated by differential display is particularly useful for identification of differentially regulated genes.  
     [0119] A very large number of fingerprints that have been previously generated can be converted to effective targets to be probed by nucleic acid molecule arrays if the mass is increased by performing PCR on an aliquot of each fingerprint in the presence of sufficient dNTPs (100 μM) and primers (about 1 μM). Fingerprints can be reamplified, as previously shown (Ralph et al.  Proc. Natl. Acad. Sci. USA  90:10710-10714 (1993)). Thus, previously determined differential display samples can be used to generate targets to probe arrays, allowing additional information to be obtained.  
     [0120] As disclosed herein, differential display was used to generate targets based on the method of Liang and Pardee (supra, 1992). The use of targets derived from oligo(dT) anchoring has some potential advantages for certain types of arrays. For example, some arrays are generated by oligo(dT) primed reverse transcription, and these clones are 3′ biased. A target generated by an oligo(dT) anchored primer and an arbitrary primer should also be 3′ biased so that each PCR product can hybridize to the corresponding 3′ biased clone. In contrast, a target generated using arbitrary priming can sample regions internal to mRNAs. If the arbitrary product is located further 5′ in the mRNA than the 3′ truncated clone, the target cannot bind to the corresponding mRNA.  
     [0121] Arbitrarily sampled targets generated using differential display with 3′ anchored oligonucleotide primers are particularly useful for probing 3′ biased libraries and, in particular, 3′ biased ESTs. 3′ anchoring is not useful for sampling RNAs that do not have poly(A) tails, such as most bacterial RNAs. Targets generated using 3′ anchor primers would also not be suitable for PCR arrays based on internal products. 3′ biased targets are also less useful for random primed libraries.  
     [0122] Other methods for generating an arbitrarily sampled target can also be used. One such method is a variant of RAP-PCR, called complexity limited arbitrary sample sequencing (CLASS). CLASS was conceived as a solution to a well known and frustrating limitation of Serial Analysis of Gene Expression (SAGE) (Velculescu et al.,  Science  270:484-487 (1995)). SAGE is a method for generating small pieces of cDNA from two sources, linking them together, and sequencing them in large numbers. The average cell contains 200,000 mRNA transcripts, representing about 20,000 different sequences, and SAGE allows sequencing of about 40 at one time. Therefore, to compare two targets using a standard sequencing apparatus, a very large number of sequencing gels, about 100, would be required to obtain information on 400,000 mRNAs, representing 200,000 mRNAs from two populations being compared. Although the method is useful for obtaining information on expression of nucleic acid molecules, each additional RNA sample increases the number of gels needed by 50, which is very expensive and time consuming. The main problem is that all 100 gels have to be run to have confidence in the statistics on rare messages that have changed in expression from 1 to 10 copies per cell.  
     [0123] To solve this problem, CLASS was devised. CLASS is similar to RAP-PCR except that the oligonucleotide primers used have degenerate 3′ ends. The degeneracy causes the primers to prime often, generating short sequence tags. By choosing a short PCR extension time, the predominant products come only from a fraction of the total complexity of the mRNA, and the size of this fraction can be adjusted at will by varying the number of 3+ degenerate bases. These short tags can then be concatenated and sequenced, rapidly yielding reliable statistics on a subsample of the message complexity, similar to the ligation and sequencing strategy used in SAGE (Valculescu et al., supra, 1995). The CLASS products can also be used as a target to probe, for example, against arrays.  
     [0124] The CLASS method is advantageous because additional sets of primers having degenerate 3′ ends can be generated and used to obtain a different sampling of nucleic acid molecules. This iterative approach to determining nucleic acid molecule expression provides more information about a pattern of expression in a source of nucleic acid molecules than the holistic approach of SAGE (Velculescu et al., supra, 1995).  
     [0125] In contrast to SAGE, which requires nearly complete sequencing of the 100 gels to be certain of any of the rare messages, CLASS allows nucleic acid molecule populations to be partitioned into small groups so that, with 10% of the work, confidence is generated for the results of 10% of all of the genes in the cell. With one round of CLASS, no information is obtained on 90% of the rare messages in the first pass (10 gels), but there is high confidence in the results for 10% of the nucleic acid molecules in a target sample. The high confidence in 10% of the genes is preferable because, when hunting for differentially regulated genes, it is expected that a pattern or “type of behavior” occurs during differential gene regulation. It is seldom, if ever, that a single gene is activated without the coordinate regulation of others controlled by the same pathway. Thus, if one is seeking any one of 10 low abundance transcripts regulated, for example, by a topoisomerase inhibitor, SAGE would require running 100 sequencing gels that would yield all 10 low abundance genes. In contrast, CLASS allows running 10 gels, in one-tenth the time, to identify at least one gene, which can be sufficient to identify a pattern of gene expression. Furthermore, CLASS can be used iteratively using different primers to run additional gels, for example, 50 gels, to get information on five times as many genes, whereas running 50 gels with SAGE would reveal no statistically relevant information. Therefore, CLASS is a much more economic approach to identifying a gene expression pattern.  
     [0126] CLASS can be applied to any species, even those for which arrays are unavailable, and to mRNAs that have not yet been deposited on arrays. Thus, whereas use of targets generated by RAP-PCR on known arrays gives expression information on known genes, CLASS gives expression information on any gene, even if not previously encountered in libraries that have been arrayed. CLASS thus provides a low cost, relatively high throughput method for obtaining information on gene expression.  
     [0127] The invention also provides methods of measuring the level of nucleic acid molecules in a target using a statistically sampled target. Methods useful for generating a statistically sampled target have been previously described (WO 99/11823; McClelland et al., supra, 1997; Pesole et al.,  Biotechniques  25:112-123 (1998); Lopez-Nieto and Nigam,  Nature Biotechnology  14:857-861 (1996)). An exemplary method for generating a statistically sampled target is statistically primed PCR (SP-PCR). The main difference between a statistical priming method and RAP-PCR is that the primers are selected by a computer program to determine the statistical occurrence of a nucleotide sequence in a group of nucleic acid molecules, rather than selecting primers arbitrarily.  
     [0128] A method for generating a statistically sampled target can be a directed statistical selection. For example, a program called GeneUP has been devised that uses an algorithm to select primer pairs to sample sequences in a list of interest, for example, a list of human mRNA associated with apoptosis, while excluding sequences in another list, for example, a list of abundantly expressed mRNA in human cells and structural RNAs such as rRNAs, Alu repeats and mtDNA (Pesole et al., supra, 1998). A directed statistical method provides a systematic determination of whether any given oligonucleotide matches any given nucleotide sequence and the number of different nucleic acid molecules to which a given oligonucleotide can bind. Such a directed statistical method can be used to generate a statistically sampled target useful in the invention.  
     [0129] Another method for generating a statistically sampled target is a Monte-Carlo statistical selection method (Lopez-Nieto and Nigam, supra, 1996). A Monte-Carlo statistical selection method randomly pairs a set of primers using a Monte-Carlo method. A Monte-Carlo method approximates the solution of determining primers that can be used for amplification by simulating a random process of primer matching. A Monte-Carlo statistical method differs from a directed statistical method in that a directed statistical method provides a systematic determination of whether any given oligonucleotide matches any given nucleotide sequence and the number of different nucleic acid molecules to which a given oligonucleotide can bind.  
     [0130] In general, two arbitrarily sampled targets, generated using different pairs of arbitrary oligonucleotides, will hybridize to largely non-overlapping sets of nucleic acid molecules in a target sample. Similarly, two statistically sampled targets, generated using different pairs of statistical oligonucleotides, will hybridize to largely non-overlapping sets of nucleic acid molecules in a target. Generally, fewer than 100 products overlap among the most intensely hybridizing 2000 colonies in two differently primed reduced complexity target (see Example I). The pattern of expression is also almost entirely different from the pattern generated by directly labeling the whole mRNA population. However, as more nucleic acid molecules are sampled by additional arbitrary sampling of the RNA population or additional statistic sampling of the RNA population, the number of non-overlapping nucleic acid molecules sampled will decrease. To some extent, the efficiency of coverage of nucleic acid molecules can be improved by the use of statistically selected primers (Pesole et al., supra, 1998). Multiple arbitrarily sampled targets generated by RAP-PCR could supply sufficient targets to cover all genes.  
     [0131] The methods described above for generating arbitrarily sampled targets and statistically sampled targets can be modified. For example, a subtraction strategy can be used to generate arbitrarily sampled targets or statistically sampled targets enriched for differentially regulated nucleic acids. A target from one source of nucleic acid molecules (A) is labeled, then mixed with a few-fold excess of unlabeled target from the other source (B). The whole mixture is denatured and added to the hybridization solution for binding to the probe. The amplified nucleic acid products present in both targets form double stranded nucleic acid molecules, and the remaining available labeled target is primarily from the differences between the two targets. The same experiment can be-done with labeled target from source (B) and excess unlabeled target from source (A). The probe bound to both sets of subtracted targets are compared to detect differential gene expression. This procedure also partly quenches repeats present in the target cDNA mixtures. The use of such a subtraction method to generate an arbitrarily sampled target or statistically sampled target can thus be used to compare two conditions by using an unlabeled target from one condition to quench the labeled target from another condition.  
     [0132] A limitation of subtraction is that it can eliminate small differences in expression that can appear to be total absence of a mRNA. Furthermore, while subtraction is useful in a binary question, it is of limited utility in cases where a large number of conditions are to be compared, combinatorially.  
     [0133] Detection of specific binding is limited by background hybridization and incomplete blockage of repeats. Therefore, in addition to using the methods described above for generating reduced complexity targets, Cot 1  DNA can be used to quench nucleic acid repetitive elements. A Cot 1  DNA genomic fraction is enriched in repeats. A target that contains Cot 1  DNA is useful for looking at low abundance nucleic acid molecules that can be difficult to detect. Although low abundance sequences can be partly quenched by the use of total genomic DNA, Cot 1  DNA is useful for the more sophisticated arrays such as PCR-based arrays, where the signal to noise ratio is sufficiently high to be concerned about relatively poorly amplified products.  
     [0134] When generating an arbitrarily sampled target or a statistically sampled target, various promoters such as T7 polymerase, T3 polymerase, SP6 polymerase or others can be incorporated into a primer so that transcription with the corresponding polymerase is used to generate the target. Using transcription to generate the target has the advantage of generating a single stranded target. A primer comprising an RNA polymerase promoter can be used in combination with any other statistical or arbitrary primer.  
     [0135] An arbitrarily sampled target or a statistically sampled target can also be generated using digestion ligation. In this case, a population of nucleic acid molecules used to generate the target is digested with a restriction enzyme and an oligonucleotide primer is ligated to generate an amplified target. Ligation-mediated PCR is where a primer binding site or part of the primer binding site is placed on a template by ligation, for example, after site-specific cleavage.  
     [0136] Nested PCR can also be used to generate an arbitrarily sampled target or statistically sampled target. Nested PCR involves two PCR steps, with a first round of PCR performed using a first primer followed by PCR with a second primer that differs from the first primer in that it includes a sequence that extends one or more nucleotides beyond the first primer sequence.  
     [0137] Targets can be enriched for those that hybridize to a particular probe. Once a target generated by a particular arbitrary or statistically primed method has been used on a particular nucleic acid population and the resulting target used against a set of probes, then the set of targets that are detectably hybridized will be known. At that point it is possible to devise a new set of targets that includes only those that were detected or mostly those that were detected by that probe. For example, if a particular primer “A” is used for RAP-PCR using RNA from the human brain and the resulting target is hybridized to an array of cDNA clones, some of the clones will be detectably hybridized. It is then possible to make an array of only those probes that were hybridized by that particular target. Most of the cDNAs on the array can be expected to hybridize with a target developed from human brain RNA made with the same primer “A”.  
     [0138] In some cases, the sequences of the nucleic acids that are the basis of targets are known. Some targets hybridize detectably with a particular probe and others do not. The sequence information associated with the targets can be used to deduce the rules of arbitrary or statistical priming events that resulted in the target that hybridized to those probes. Such information will help to predict what sequences are likely to be sampled by a particular primer if that sequence occurs in the target. Such information can improve the estimates of which sequences are sampled efficiently and which sequences are sampled efficiently by a particular primer.  
     [0139] The methods of the invention are particularly useful for measuring the level of a molecule in a target using an array. As used herein, the term “array” or “array of molecules” refers to a plurality of molecules stably bound to a solid support. An array can comprise, for example, nucleic acid, oligonucleotide or polypeptide-nucleic acid molecules. It is understood that, as used herein, an array of molecules specifically excludes molecules that have been resolved electrophoretically prior to binding to a solid support and, as such, excludes Southern blots, Northern blots and Western blots of DNA, RNA and proteins, respectively.  
     [0140] As used herein, the term “non-dot blot” array refers to an array in which-the molecules of the array are attached to the solid support by a means other than vacuum filtration or spotting onto a nitrocellulose or nylon membrane in a configuration of at least about 2 spots per cm 2 .  
     [0141] As used herein, the term “peptide-nucleic acid” or “PNA” refers to a peptide and nucleic acid molecule covalently bound (Nielson, Current Opin. Biotechnol. 10:71-75 (1999)).  
     [0142] As used herein, the term “polypeptide,” when used in reference to PNA, means a peptide, polypeptide or protein of two or more amino acids. The term is similarly intended to refer to derivatives, analogues and functional mimetics thereof. For example, derivatives can include chemical modifications of the polypeptide such as alkylation, acylation, carbamylation, iodination, or any modification which derivatizes the polypeptide. Analogues can include modified amino acids, for example, hydroxyproline or carboxyglutamate, and can include amino acids that are not linked by peptide bonds. Mimetics encompass chemicals containing chemical moieties that mimic the function of the polypeptide regardless of the predicted three-dimensional structure of the compound. For example, if a polypeptide contains two charged chemical moieties in a functional domain, a mimetic places two charged chemical moieties in a spatial orientation and constrained structure so that the charged chemical function is maintained in three-dimensional space. Thus, all of these modifications are included within the term “polypeptide.” 
     [0143] The solid support for the arrays can be nylon membranes, glass, derivatized glass, silicon or other substrates. The arrays can be flat surfaces such as membranes or can be spheres or beads, if desired. The molecules can be attached as “spots” on the solid support and generally can be spotted at a density of at least about 5/cm 2  or 10/cm 2 , but generally does not exceed about 1000/cm 2 .  
     [0144] Various methods to manufacture arrays of DNA molecules have been described (reviewed in Ramsay, supra, 1998; Marshall and Hodgson, supra, 1998). Arrays are available containing nucleic acid molecules from various species, including yeast, mouse and human. The use of arrays is advantageous because differential expression of many genes can be determined in parallel.  
     [0145] One type of array contains thousands of PCR products per square centimeter. Arrays of PCR products from segments of mRNAs have been attached to glass, for example, and probed using cDNA populations from two sources. Each cDNA or cRNA population is labeled with a different fluorescent dye and hybridization is assessed using fluorescence (DeRisi et al.,  Nature Genet.  14:457-460 (1996); Schena et al.,  Science  270:467-470 (1995)). Arrays are also available containing over 5000 PCR products from selected I.M.A.G.E. clones. An array of PCR products also is available for every yeast ORF and for a subset of human ESTs.  
     [0146] Another type of array contains colonies of 18,432  E. coli  clones, each carrying a different I.M.A.G.E. EST plasmid, and each spotted twice on a 22×22 cm membrane (Genome Systems). One advantage of using the arrays from the I.M.A.G.E. consortium is that more than 80% of the clones have single pass sequence reads from the 5′ or 3′ end, or both, deposited in the GenBank database. Thus, it is usually not necessary to clone or sequence any DNA to determine if there is a known gene or other ESTs that share the same sequence. UniGene clustering of human and mouse ESTs that appear to be from the same gene greatly aids in this process (http://www.ncbi.nlm.nih.gov/UniGene/index.html). Mapping onto chromosomes at a resolution of a few centiMorgans is also available for most of these clusters at the same web site. The clones on these arrays are all available to be used to probe nucleic acid molecules or to complete the sequencing (www-bio.llnl.gov). It is often possible to identify a close homolog in other species. In contrast to PCR product arrays and oligonucleotide arrays, which are free of other DNAs, each spotted EST is associated with  E. coli  genomic DNA from the host. Thus, the clone arrays can have higher background than PCR arrays or oligonucleotide arrays.  
     [0147] If EST arrays are used, 5′ RACE can be used to extend beyond the ESTs currently available (Zhang and Frohman,  Methods Mol. Biol.  69:61-87 (1997)). When cDNA libraries that contain near full length clones are available and end sequenced, it will be possible to go from a differentially hybridized spot to a full length cDNA, directly.  
     [0148] Another class of arrays uses oligonucleotides that are either attached to a glass or silicon surface or manufactured by sequential photochemistry on the DNA chip (Chee et al.,  Science  274:610-614 (1996)). Such chips can contain tens of thousands of different oligonucleotide sequences per square centimeter. Arrays of oligonucleotide nucleic acid analogs such as peptide-nucleic acids, for example, can be prepared (Weiler et al.,  Nucleic Acids Res.  25:2792-2799 (1997)).  
     [0149] Hybridization of fingerprints to arrays has the huge advantage that there is generally no need to isolate, clone, and sequence the genes detected. In principle, all known human mRNAs will fit on three membranes (about 50,000 genes), or in a smaller area on glass arrays or other solid supports. At present, each fingerprint has a sufficient complexity to hybridize to over 2000 of the 50,000 known genes.  
     [0150] The use of arrays, which can have thousands of genes that can bind to a target, particular genes for further characterization can be selected based on desired criteria. For example, identified genes can be chosen that are already known and for which a new role in the condition of interest can be deduced. Alternatively, some of the genes can be family members of known genes with known functions for which a plausible role can be determined.  
     [0151] In addition to arrays, a number of cDNA libraries are available, for example, from the I.M.A.G.E. consortium (www-bio.llnl.gov/bbrp/image/image.html), including libraries available on nylon membranes, for example, from Research Genetics (Huntsville Ala.; www.resgen.com), Genome Systems (St. Louis Mo.; www.genomesystems.com), and the German Human Genome Project (www.rzpd.de). These libraries include clones from various human tissues, stages of development, disease states and other sources.  
     [0152] The methods of the invention include the step of detecting the amount of specific binding of the probe to the target. As disclosed herein, a variety of detection methods can be used. For example, if a detectable moiety is a radioactive moiety, the method of detection can be autoradiography or phosphoimaging. Phosphoimaging is advantageous for quantitation and shortened data collection time. If a detectable moiety is a fluorescent moiety, the method of detection can be fluorescence spectroscopy or confocal microscopy.  
     [0153] The methods of the invention use nucleic acid probes to measure the level of expression of a nucleic acid molecule in a target. If a radioactive moiety is attached to a target, for example, incorporation of the radioactive moiety can be by any enzymatic or chemical method that allows attachment of the radioactive moiety. For example, end-labeling can be used to attach a radioactive moiety to the end of a nucleic acid molecule. Alternatively, a radioactive nucleotide, in particular a  32 P-,  33 P-, or  35 S-labeled nucleotide, can be incorporated into the nucleic acid molecule during synthesis. The use of random primed synthesis is particularly useful for generating a high specific activity target. Generally, random primed synthesis generates approximately equal amounts of randomly primed nucleic acid molecules from both strands of double stranded PCR products, which will re-anneal to some degree during hybridization to the target (see Example I). If desired, the amount of re-annealing can be limited, for example, using exoII digestion.  
     [0154] When generating a labeled target or probe, it is generally preferable to incorporate a labeled nucleotide that is not ATP or dATP. The use of labeled dATP can cause an increase in the background because any poly-A sequences in the target or probe will become heavily labeled and will hybridize to the strands containing poly-T stretches complementary to the poly-A tails present in all of the clones. Similarly, the use of dTTP would heavily label poly-T stretches complementary to the polyA tails in mRNA.  
     [0155] A fluorescent dye can also be attached to or incorporated in the probe or target. If desired, a different fluor detectable at different wavelengths can be incorporated into different targets and used simultaneously on the same probe. The use of different fluors is advantageous since multiple targets can be bound to the same probe and detected. A fluorescently labeled target can be detected using, for example, a fluorescent scanner or confocal microscope. Measuring the relative abundance of two targets simultaneously on the same array rather than on two different arrays eliminates problems that arise due to differences in the hybridization conditions or the quantity of target PCR product on replicates of the same array. Nylon membranes are typically unsuitable for most commercially available fluorescent tags due to background fluorescence from the membrane itself.  
     [0156] Infrared dyes are also useful as detectable moieties for attachment to a probe or target. Infrared dyes are particularly useful with targets or probes such as arrays attached to nylon membranes, provided the membrane is free of protein.  
     [0157] When determining the level of a nucleic acid molecule in a target, some variation can occur, in particular for certain amplification products that are very sensitive to the amplification conditions. To control for variation in amplification products between nucleic acid targets, the target can be generated at two concentrations of nucleic acid molecules, differing by a factor of two or more. The use of various nucleic acid concentrations to generate a target to confirm differential expression is described herein (see Examples II and III).  
     [0158] The methods of the invention are directed to detecting specific binding of a target to a probe. When hybridizing a target to a probe, the specificity of binding is determined-by the stringency of the hybridization conditions. The length of oligonucleotide primers and the temperature of the amplification reaction contributes to the final product. The products are a function of both the starting abundance of each target nucleic acid molecule and the quality of the match between the oligonucleotide primer and the amplified nucleic acid target. For example, oligonucleotide primers of about 8 bases in length at reaction temperatures of about 60° C. can be used to generate a target. Hybridization conditions can range, for example, from about 32° C. in about 2×SSC to about 68° in about 0.1×SSC. The hybridization temperature can be, for example, about 40° C., about 45° C., about 50° C., about 55° C., about 60° C. or about 65° C. Furthermore, the SSC concentration (see below) can be, for example, about 0.2×,0.3×, 0.5×, 1× or 1.5×.  
     [0159] The invention additionally provides a method for determining the relative amounts of nucleic acid molecules in two targets by comparing the amount of specific binding of a probe to the target, wherein the amount of specific binding corresponds to an expression level of the nucleic acid molecules in the target, to an expression level of the nucleic acid molecules in a second target. For example, if desired, the expression level in a first target, which can be a target for which the level of expression is unknown, can be compared to the expression level in a second target. The expression level in the second target can be determined, for example, by binding the same probe to the second target and determining the level of expression in the second target. The expression level in the first and second target can then be compared.  
     [0160] The relative expression level in a first target can also be compared to the expression level in a second target, where the abundance in the second target is already known. As used herein, the term “known” when used in reference to expression level of a nucleic acid molecule means that an abundance of a nucleic acid molecule has been previously determined. It is understood that such a known abundance would apply to a particular set of conditions. It is also understood that, for the purpose of comparing the abundance of a nucleic acid molecule in an unknown target to a known abundance, the same method of measuring the abundance between the targets is used.  
     [0161] The invention also provides a method of identifying two or more differentially expressed nucleic acid molecules associated with a condition. The method includes the step of measuring the level of two or more nucleic acid molecules in a target, for example using an arbitrarily sampled target or a statistically sampled target, wherein the amount of specific binding of the target to the probe corresponds to an abundance of the nucleic acid molecules in the target. The method further includes the step of comparing the relative expression level of the nucleic acid molecules in the target to an expression level of the nucleic acid molecules in a second target, whereby a difference in expression level between the targets indicates a condition.  
     [0162] As used herein, the term “differentially expressed” means that the abundance of a molecule is expressed at different levels between two targets. Two targets can be from different cells or tissues, or the target can be from the same cell or tissue under different conditions. The condition can be, for example, associated with a disease state such as cancer, autoimmune disease, infection with a pathogen, including bacteria, virus, fungal, yeast, or single-celled and multi-celled parasites; associated with a treatment such as efficacy, resistance or toxicity associated with a treatment; or associated with a stimulus such as a chemical, for example, a drug or a natural product, for example, a growth factor.  
     [0163] The methods of the invention are useful for determining differential gene expression between two targets. The methods of the invention can be applied to any system where differential gene expression is thought to be of significance, including drug and hormone responses, normal development, abnormal development, inheritance of a genotype, disease states such as cancer or autoimmunge disease, aging, infectious disease, pathology, drug treatment, hormone activity, aging, cell cycle, homeostatic mechanisms, and others, including combinations of the above conditions.  
     [0164] As disclosed herein, the abundance of nucleic acid molecules in two targets can be compared to identify two or more differentially expressed nucleic acid molecules (see Examples I to III). Using arbitrarily sampled targets, targets treated with and without EGF were hybridized with probes and a number of genes regulated by EGF were identified. EGF-regulated genes were found that increased in response to EGF and decreased in response to EGF (see Tables 1 and 2 in Examples II and III, respectively). The methods of the invention can therefore be used to determine nucleic acid molecules that increase in response to a stimulus or decrease in response to a stimulus (see Example II).  
     [0165] The arbitrarily sampled targets and statistically sampled targets used in the invention can readily detect less abundant nucleic acid molecules in a population. Therefore, the methods of the invention are particularly useful for identifying differentially expressed nucleic acid molecules since differentially expressed nucleic acid molecules are often less abundant.  
     [0166] The methods of the invention can be applied to any two targets to determine differential gene expression. The methods of the invention can be used, for example, to diagnose a disease state. In such a case, a “normal” target is compared to a potential disease target to determine differential gene expression associated with the disease. A normal target can be a target sample of the same tissue nearby the diseased tissue from the patient. A normal target can also be a sample of the same tissue from a different individual. Using methods of the invention, a profile of normal expression can be established by determining a gene expression pattern in one to many normal target samples, which can then be used to compare to a potentially diseased target sample. Differential gene expression between the normal and diseased tissue can be used to diagnose or confirm a particular disease state. Furthermore, a collection of target samples obtained from known diseased tissue can similarly be determined to identify an abundance profile of the target reflecting gene expression associated with that disease. In such a case, comparison of a potential disease target sample to a known disease target sample with no differential gene expression would indicate that the potential disease target sample was associated with the disease.  
     [0167] The methods of the invention can also be used to assess treatment of an individual with a drug. The analysis of gene expression patterns associated with a particular drug treatment is also known as pharmacogenomics. The methods of the invention can be used to determine efficacy of a treatment, resistance to a treatment or toxicity associated with a treatment. For example, a gene expression profile can be determined on an individual prior to treatment and after treatment for a particular disease or condition. A difference in gene expression can then be correlated with the effectiveness of the treatment. For example, if an individual is found to be responsive to treatment and if that treatment is associated with differential gene expression, the identification of differential gene expression can be used to correlate with efficacy of that treatment. As described above, a gene expression pattern associated with an untreated individual can be determined in the individual prior to treatment or can be determined in a number of individuals who have not been given the treatment. Similarly, a change in expression pattern associated with efficacy of the treatment can be determined in a number of individuals for which the treatment was efficacious. In such a case, comparison of a treated target sample to a known target sample associated with efficacious treatment with no differential gene expression would indicate that the treatment was likely to be efficacious. A similar approach can be used to determine the association of a treatment with toxicity of the treatment or resistance to a treatment. Resistance to a treatment could be associated with a change in expression pattern from an untreated target sample or could be associated with no change in the expression pattern compared to an untreated target sample.  
     [0168] The methods of the invention can also be used to determine co-regulated genes that can be potential targets for drug discovery. For example, a cell or organism can be treated with a stimulus and differential gene expression between the untreated target sample and the target sample treated with a stimulus can be determined. The stimulus can be, for example, a drug or growth factor. A difference in the abundance of nucleic acid molecules between an untreated target sample and a target sample treated with a stimulus can be used to identify differential gene expression associated with the stimulus. Such a differential expression pattern can be used to determine if a target sample has been exposed to a stimulus. Additionally, the gene expression profile can be used to identify other chemicals that mimic the stimulus by screening for compounds that elicit the same gene expression profile as the original stimulus. Thus, the methods of the invention can be used to identify new drugs that have a similar effect as a known drug.  
     [0169] The methods of the invention are useful for identifying a marker for a pathway that correlates with a drug response by determining an abundance profile for a given target sample that reflects the expression profile of the source population of nucleic acids such as the source RNA. For example, the methods of the invention can be used to define the “neighborhood” of potential therapeutic targets by identifying several genes regulated in response to a drug, thereby providing “neighbors” in a pathway that are potential drug targets. The invention can also be used to define bad neighborhoods, for example, pathways that “failed” therapeutics, which can indicate that a particular pathway should not be perturbed. Additional insights into the function of a pathway can be obtained by sequencing any differentially expressed genes for which complete sequence information is unavailable. The methods are particularly useful for drug comparison. Correlation of gene expression patterns with a drug response can be used to determine why two similar drugs have a somewhat different spectrum of effects.  
     [0170] With knowledge of the correlation between gene expression and response to a drug, drugs can be tested in cell types that are of more relevance to a particular disease or condition. By knowing the pathways that are present in a cell type associated with a pathology, predictions can be made regarding the drug responses of the cell type and thereby allow choice of drugs from a tested panels of drugs that are most likely to affect the pathology. The correlation of information on drug response and gene expression also can aid in choosing drugs that would be synergistic, for example, drugs that hit non-overlapping pathways, or, for example, drugs that affect overlapping pathways when genes in the overlap are targeted.  
     [0171] The methods of the invention can be applied to determining the response to a stimulus, in particular to determining a response to a stimulus for drug discovery. One potential application is to use the methods of the invention on the 60 cell lines in the National Cancer Institute (NCI) drug screening panel. These 60 cell lines are maintained by the NCI and used to assess drug activity.  
     [0172] For example, each of the 60 cell lines of the NCI panel can be used as a complex measuring device that reports the single variable of cell growth and, secondarily, apoptosis. Changes in each cell type&#39;s growth upon treatment with a chemical such as a drug is determined. Studies of tens of thousands of drugs, when compared over all 60 cell lines, have shown that similar effects on growth have proven to share mechanisms of action. Comparing the response of the 60 cell lines to various drugs allows grouping of drugs according to their detailed chemical functionality. Consequently, the panel of cell lines has become one of the most important analytical tools for drug discovery.  
     [0173] The methods of the invention can be applied to analyzing drug response in the 60 cell lines of the NCI panel. As disclosed herein, the methods are applicable to determining differential gene expression, which can be correlated with the response of the cells to a particular drug. The methods can be used to identify many differentially expressed genes associated with a drug response. Therefore, an analysis of gene expression in untreated cells in the 60 cell line NCI drug screening panel can be used to determine a profile of gene expression, based on the presence or absence of mRNAs, that correlate with some of the many 10,000&#39;s of drugs that have been used on the panel.  
     [0174] Differential gene expression patterns are expected to correlate with drug response. Following identification of such a correlation in 30 of the cell lines, prediction of drug responses in the remaining 30 cell lines can be tested. This strategy circumvents the need to determine extensive expression profiles for all 60 cell lines for every new drug to find genes that correlate with the ability to respond to the drug. This strategy differs from previous methods in that differential expression of the gene after treatment does not need to occur. All that is necessary is that the gene be differentially regulated between cell types prior to treatment.  
     [0175] Each of the 60 cell lines has its characteristic response to drugs, and these responses depend on the cell&#39;s phenotype. The response of any cell to any drug depends on which genetic systems are operative in that cell. Once treated, the cell&#39;s genetic mechanisms are perturbed, leading to differential gene expression, differential protein modification, and a wide variety of other changes that can be subtle. Nonetheless, it is the ground state genetic pattern or profile of gene expression, before any exposure to drug, that determines how the cell responds to drugs.  
     [0176] The ground state of genetic profile is an important state to characterize for cells, for example, cells of the NCI panel. The ground state of the cell has predictive power for how a given cell will respond to a given drug. Furthermore, the ground state is the only unifying point of reference for the behavior of almost 100,000 different drugs and can be used to determine response to additional drugs.  
     [0177] For example, if two steroids and two alkylating agents are applied to the panel of 60 cell lines, and their growth spectra are compared, the average responses of the cell lines to the steroids tends to be similar, the average responses to the alkylating agents tend to be similar, but a comparison of responses to steroids versus alkylating agents show fewer similarities. This reflects the fact that steroids elicit their effects through naturally existing receptors, whereas alkylating agents elicit their effects by causing widespread damage. The signal transduction pathways for handling steroidal signals versus handling damage are largely different.  
     [0178] When a panel of steroids are used to challenge the 60 cell lines, some of the cells are growth accelerated, some growth inhibited, and some are indifferent to steroids. Much of this data is available on the NCI web site (http://www.nci.nih.gov/). An obvious next step is to examine gene responses to the steroids to see which genes are activated, which are inactivated, and which are indifferent. Each cell type&#39;s genes will respond differently, depending on which of about 30 steroid receptor genes are expressed in the cell type before steroid treatment.  
     [0179] The various responses of genes to steroids are cell type-dependent, in large part due to which receptors are present. By comparing the ground state gene expression of the NCI panel of cells, the spectrum of steroid receptor genes expressed in each cell type can be described, thereby explaining what is needed, in genetic terms, for a cell to be responsive to any particular steroid.  
     [0180] The drug-receptor, or hormone-receptor, relationship described above is one example of a correlation that can be drawn between the NCI panel baseline gene expression database and the NCI panel drug response database. Other drug responses can be readily determined. For example, drugs that induce apoptosis also induce gene expression, and different apoptotic responses correlating with cell type can be used to determine gene products that control apoptosis.  
     [0181] It is understood that methods of the invention can be applied to any cell type, in addition to the NCI panel of cells, for characterization of a response to a drug or other stimulus. The functional overlap between drugs is an important concern in drug discovery. A study of the responses of genes to drugs in different cell types is useful because gene expression determines the response of the cell to the drug. The methods of the invention can therefore be applied to determine the response of one or more cell lines to a particular drug.  
     [0182] The methods can also be applied to characterize the ground state of the NCI panel of cells. The methods described herein can be used to correlate the response of tens of thousands of drugs with genes in the pathways regulated by the drug. The methods of the invention can be applied to determine an expression profile for the &gt;80,000 drugs previously tested with the NCI panel of cells. The methods are applicable to determining coordinate mechanisms of drug action, likely pathways controlling drug activity, pathways that correlate with toxicity, apoptosis and other effects of drugs.  
     [0183] The invention also provides methods for the use of the patterns of gene expression by a panel of different untreated cells or tissues to correlate basal gene expression with susceptibility to a treatment, such as differences in the growth of cells, for example, the NCI panel of cells, in the presence of a drug, pathogen or other stimulus. The methods can be applied to determine genes and pathways that are present prior to treatment and also to correlate treatment with the phenotype induced by the treatment.  
     [0184] To obtain additional information on gene expression, the expression pattern of two different RNA populations from different conditions can be determined (McClelland et al.,  Nucleic Acids Res.  22:4419-4431 (1994); McClelland et al.,  Trends Genet.  11:242-246 (1995)). For example, if interested in apoptosis, using a target from a cell that has been stressed but which has not undergone apoptosis can be used to determine genes responsive to apoptosis, genes responsive to stress, and genes that respond to both. The identification of differentially regulated genes can be used to further characterize transcriptional activity of genes under various conditions. The genes can be further characterized to correlate promoters of regulated genes with signal transduction pathways that respond to a given condition.  
     [0185] When determining differential expression of a nucleic acid molecule, the determination that an RNA sampled in a target is differentially regulated is initially made by comparing differential abundance at two different concentrations of nucleic acid in the target sample. Abundance is determined for the nucleic acid molecules of the target sample for which no difference in abundance is observed at two different concentrations of RNA source. Only those hybridization events that indicate differential expression at both RNA concentrations in both RNA sources are used (see Examples II and III).  
     [0186] For hybridization to an array to determine differential expression, four membranes were used for radioactively labeled target, one for each of two concentrations of RNA for each of the two RNA samples compared (see Examples I to III). If two color fluorescence is used for detecting the target, then two membranes are used, one for each of the two concentrations of starting target sample nucleic acids, because the two targets with different detectable fluorescent markers can be mixed and applied to the same probe. If a subsequent verification step is employed, for example, RT-PCR, one marker can be used for each target sample.  
     [0187] Confirmation of differential expression does not need a full length sequence and can be confirmed using RT-PCR of the known region. In particular, low stringency PCR can be used to generate products a few hundred bases in length (Mathieu-Daude et al.,  Mol. Biochem. Parasitol.  92:15-28 (1998)). This method generates internal “control” PCR products that can be used to confirm the quality of the PCR reaction and the quality and quantity of the RNA used.  
     [0188] The invention additionally provides a profile of five or more stimulus-regulated nucleic acid molecules. As used herein, the term “profile” refers to a group of two or more nucleic acid molecules that are characteristic of a target under a given set of conditions. The invention provides a profile comprising a portion of a nucleotide sequence selected from the group consisting of the nucleotide sequences referenced as SEQ ID NOS:1-45. The profile includes a portion of a nucleotide sequence of the GenBank accession numbers H11520, H11161, H11073, U35048, R48633, H28735, AF019386, H25513, H25514, M13918, H12999, H05639, L49207, H15184, H15124, X79781, H25195, H24377, M31627, H23972, H27350, AB000712, R75916, X85992, R73021, R73022, U66894, H10098, H10045, AF067817, R72714, X52541, H14529, M10277, H27389, D89092, D89678, H05545, J03804, H27969, R73247, U51336, H21777, K00558, and D31765. The profile of the invention includes a portion of the nucleotide sequences encoding TSC-22, fibronectin receptor α-subunit, ray gene, X-box binding protein-1, CPE receptor, epithelium-restricted ets protein ESX and Vav-3.  
     [0189] The invention also provides a target comprising a portion of each of the nucleotide sequences referenced as SEQ ID NOS:1-45. The target includes a portion of a nucleotide sequence of the GenBank accession numbers H11520, H11161H11073, U35048, R48633, H28735, AF019386, H25513, H25514, M13918, H12999, H05639, L49207, H15184, H15124, X79781, H25195, H24377, M31627, H23972, H27350, AB000712, R75916, X85992, R73021, R73022, U66894, H10098, H10045, AF067817, R72714, X52541, H14529, M10277, H27389, D89092, D89678, H05545, J03804, H27969, R73247, U51336, H21777, K00558, and D31765. The invention also provides a probe comprising a portion of a nucleic acid sequence selected from the group consisting of SEQ ID NOS:1-45.  
     [0190] The invention further provides a substantially pure nucleic acid molecule comprising a nucleic acid sequence selected from the group consisting of SEQ ID NOS:1-45, or a functional fragment thereof, so long as the nucleic acid molecule does not include the exact SEQ ID NOS:1-45.  
     [0191] The invention additionally provides a method of measuring the amount of two or more nucleic acid molecules in a first target relative to a second target. The method includes the step of hybridizing a first amplified nucleic acid target comprising two or more nucleic acid molecules to a probe, wherein the target is amplified from a population of nucleic acid molecules using one or more oligonucleotides, wherein the oligonucleotide hybridizes by chance to a nucleic acid molecule in the population of nucleic acid molecules, wherein the amplification is not based on abundance of nucleic acids in the population of nucleic acid molecules, and wherein the amplified nucleic acids in the target are enhanced for less abundant nucleic acids in the population of nucleic acid molecules. Further included in the method are the steps of detecting the amount of hybridization of the first amplified nucleic acid target to the probe, wherein the amount of hybridization corresponds to an abundance of the nucleic acid molecules in the first target; and comparing the abundance of the nucleic acid molecules in the first target to the abundance of the nucleic acid molecules in a second target, wherein the amplified nucleic acid target comprises a subset of nucleic acids in the initial nucleic acid populations.  
     [0192] The invention further provides a method of measuring the amount of two or more nucleic acid molecules in a first target relative to a second target. The method includes the step of hybridizing a first amplified nucleic acid target comprising 50 or more nucleic acid molecules to a probe, wherein the target is amplified from a population of nucleic acid molecules, wherein the amplification is not based on abundance of nucleic acids in the population of nucleic acid molecules, and wherein the amplified nucleic acids in the target are enhanced for less abundant nucleic acids in the population of nucleic acid molecules. The method further includes the steps of detecting the amount of hybridization of the amplified nucleic acid target to the probe, wherein the amount of hybridization corresponds to an expression level of the nucleic acid molecules in the first target; and comparing the abundance of the nucleic acid molecules in the first target to an abundance of the nucleic acid molecules in a second target, wherein the amplified nucleic acid target comprises a subset of nucleic acids in each nucleic acid population such as an RNA population.  
     [0193] As used herein, the term “hybridizes by chance,” when referring to an oligonucleotide, means that hybridization of the oligonucleotide to a complementary sequence is based on the statistical frequency of the complementary sequence occurring in a given nucleic acid molecule. An oligonucleotide that hybridizes by chance is generated by determining the sequence of the oligonucleotide and subsequently determining if the oligonucleotide will hybridize to one or more nucleic acid molecules. The hybridization of such an oligonucleotide is not predetermined by the sequence of a known nucleic acid molecule and therefore occurs by chance. As such, an arbitrary oligonucleotide is considered to hybridize by chance since the oligonucleotides are determined without reference to the exact sequence to be amplified. In contrast, an oligonucleotide that does not hybridize by chance is one that is generated by first analyzing a known sequence and then identifying an exact sequence in the nucleic acid molecule that can be used as an oligonucleotide that will amplify an exact sequence between the oligonucleotides. The hybridization of such an oligonucleotide has been predetermined by the sequence of a known nucleic acid molecule and, therefore, does not occur by chance.  
     [0194] As used herein, the phrase “amplification is not based on abundance” means a target comprises nucleic acid molecules which are representative of the nucleic acid molecules in a population of nucleic acid molecules without regard to the relative amount of individual nucleic acid molecules in the population.  
     [0195] As used herein, the phrase “enhanced for less abundant nucleic acids” means that individual nucleic acid molecules that are less abundant in the population of nucleic acid molecules are amplified so that the amount of these less abundant nucleic acid molecules would be increased relative to the amount of these nucleic acid molecules in the original population of nucleic acid molecules. Thus, the relative proportion of nucleic acid molecules in the population of nucleic acid molecules would not be maintained in the target.  
     [0196] As used herein, the term “single sample” when used in reference to a target means that the target is generated using nucleic acid molecules from a single cell, tissue or organism sample that has not been previously exposed to another sample. For example, if a target was generated from a population of nucleic acid molecules that was determined by the exposure of one sample to another, for example, the subtraction of the nucleic acid molecules of one sample from another, such a target would not be considered as coming from a single sample.  
     [0197] The following examples are intended to illustrate but not limit the present invention.  
     EXAMPLE I  
     Generation and Use of Arbitrarily Sampled Targets to Probe a DNA Array  
     [0198] This example describes the generation of an arbitrarily sampled target having reduced complexity to probe a DNA array to determine mRNA expression.  
     [0199] A DNA fingerprint was generated using RAP-PCR and was converted to high specific activity probe using random hexamer oligonucleotides (Genosys Biotechnologies; The Woodlands Tex.). Up to 10 μg of PCR product from RAP-PCR was purified using a QIAQUICK PCR Purification Kit (Qiagen, Inc.; Chatsworth Calif.), which removes unincorporated bases, primers, and primer dimers smaller than 40 base pairs. The DNA was recovered in 100 μl of 10 mM Tris, pH 8.3. Random primed synthesis with incorporation of radioactive phosphorus from (α- 32 P)dCTP was used under standard conditions. 10% of the recovered fingerprint DNA (10 μl) was combined with 6 μg random hexamer oligonucleotide primer, and 1 μg of one of the fingerprint primers (Genosys) in a total volume of 28 μl, boiled for 3 min, then placed on ice. The hexamer/primer/DNA mix was mixed with 22 μl reaction mix to yield a 50 μl reaction containing a 0.05 mM concentration of three dNTP (dATP, dTTP and dGTP; minus dCTP), 100 μCi of 3000 Ci/mmol (α- 32 P) dCTP (10 μl), 1× Klenow fragment buffer (50 mM Tris-HCl, pH 8.0, 10 mM MgCl 2 , 50 mM NaCl) and 8 U Klenow fragment (3.82 U/μl; Gibco-BRL Life Technologies; Gaithersburg Md.). The reaction was performed at room temperature for 4 hr. For maximum target length, the reaction was chased by adding 1 μl of 2.5 mM dCTP and incubated for 15 min at room temperature followed by an additional 15 min incubation at 37° C. The unincorporated nucleotides and hexamers were removed with the Qiagen Nucleotide Removal Kit (Qiagen) and the purified products were eluted twice in 140 μl 10 mM Tris, pH 8.3.  
     [0200] For hybridization to the array, four membranes were used for radioactively labeled target, one for each of two concentrations of RNA for each of the two RNA samples to be compared. To prepare the cDNA filters (Genome Systems), the filters were prewashed in three changes of 2×SSC and 0.1% sodium dodecyl sulfate (SDS) in a horizontally shaking flat bottom container to reduce the residual bacterial debris. 20×SSC contains 3 M NaCl, 0.3 M Na 3 citrate-2H 2 O, pH 7.0. The first wash was carried out in 500 ml for 10 min at room temperature. The second and third washes were carried out in 1 liter of prewarmed (50° C.) prewash solution for 10 min each.  
     [0201] For prehybridization, the filters were transferred to roller bottles and prehybridized in 60 ml prewarmed (42° C.) prehybridization solution containing 6×SSC, 5× Denhardt&#39;s reagent, 0.5% SDS, 100 pg/ml fragmented, denatured salmon sperm DNA (Pharmacia; Piscataway N.J.) and 50% formamide (Aldrich; Milwaukee Wis.) for 1-2 hr at 42° C. 50× Denhardt&#39;s solution contains 1% Ficoll, 1% polyvinylpyrrolidone and 1% bovine serum albumin, sterile filtered.  
     [0202] For hybridization, the prehybridization solution was removed and 7 ml prewarmed (42° C.) hybridization solution, containing 6×SSC, 0.5% SDS, 100 μg/ml fragmented, denatured salmon sperm DNA and 50% formamide, was added. To decrease the background hybridization due to repeated sequences such as Alu repeats, long interspersed repetitive elements (LINE) or centromeric DNA repeats, sheared human genomic DNA (1 μg/ml stock concentration) was denatured in a boiling water bath for 10 min and immediately added to the hybridization solution to a final concentration of 10 μg/ml. Simultaneously, the labeled target (280 μl) was denatured in a boiling water bath for 4 min and immediately added to the hybridization solution. Hybridization was carried out at 42° C. for 2 to 48 hrs, typically 18 hr, in a hybridization oven using roller bottles or sealed in a plastic bag and incubated in a water bath.  
     [0203] For the washes, the temperature was set to 55° C. in the incubator oven (Techne HB-1D; VWR Scientific; San Francisco Calif.). The hybridization solution was poured off and the membrane was washed twice with 50 ml 2×SSC and 0.1% SDS for 5 min at room temperature. The membrane was then washed with 100 ml 0.1×SSC and 0.1% SDS and incubated for 10 min at room temperature. For the further washes, the wash solution, containing 0.1×SSC and 0.1% SDS, was prewarmed to 50° C. and the filter was washed for 40 min in a roller bottle with 100 ml wash solution. The filter was then transferred to a horizontally shaking flat bottom container and washed in 1 liter of the wash solution for 20 min under gentle agitation. The filter was transferred back to a roller bottle containing 100 ml prewarmed 0.1×SSC and 0.1% SDS and incubated for 1 hr. The final wash solution was removed and the filter briefly rinsed in 2×SSC at room temperature.  
     [0204] After washing, the membranes were lightly dried with 3MM paper and the slightly moist membranes were wrapped in SARAN wrap. The membranes were exposed to X-ray film.  
     [0205]FIG. 1 shows differential hybridization to clone arrays. All four images show a closeup of an autoradiogram for the same part of a larger membrane. Each image spans about 4000 double spotted  E. coli  colonies, each carrying a different EST clone. Panel A shows hybridization of 1 μg of polyA +  RNA from confluent human keratinocytes that was radiolabeled during reverse transcription. About 500 clearly hybridizing clones can be seen. Panels B and C show RAP-PCR fingerprints with a pair of arbitrary primers that was performed on cDNA from oligo(dT) primed cDNA of confluent human keratinocytes that were untreated (Panel B) or treated with EGF (Panel C). The pattern of hybridizing genes was almost identical in Panels B and C, but entirely different from that seen with total polyA+ RNA (compare to Panel A). The two radiolabeled colonies from one differentially expressed cDNA are indicated with an arrow. Differential expression of this gene was subsequently confirmed by specific RT-PCR (Trenkle et al.,  Nucl. Acids Res.  26:3883-3891 (1998)).  
     [0206]FIG. 1D shows a RAP-PCR fingerprint with a different pair of arbitrary primers that was performed on RNA from confluent human keratinocytes. This pattern of hybridization is almost entirely different from that found with the previous primer pair (Panel B) and with mRNA (Panel A), with very few overlapping spots between Panel D and Panels A and B.  
     [0207] These results demonstrate that arbitrarily sampled targets, which have reduced complexity, allow detection of mRNAs that are not detectable using total message as a target. Thus, unlike a total message target, which detects mRNAs based on their abundance, an arbitrarily sampled target can be used to detect less abundant mRNAs.  
     EXAMPLE II  
     An Arbitrarily Sampled Target Generated by RT-PCR Detects Genes Differentially Expressed in Response to EGF  
     [0208] This example describes the use of RT-PCR with arbitrary primers to generate an arbitrarily sampled target for detecting differential gene expression upon treatment of cells with EGF.  
     [0209] An arbitrarily sampled target generated by RT-PCR was used to probe arrays for differential gene expression (Trenkle et al.,  Nucleic Acids Res.  26:3883-3891 (1998)). For RNA preparation, the immortal human keratinocyte cell line HaCaT (Boukamp et al.,  Genes Chromosomes Cancer  19:201-214 (1997)) was grown to confluence and maintained at confluence for two days. The media, DMEM containing 10% fetal bovine serum (FBS) and penicillin/streptomycin was changed one day prior to experiments. EGF (Gibco-BRL) was added at 20 ng/ml, or TGF-β (R&amp;D Systems; Minneapolis Minn.) was added at 5 ng/ml. Treated and untreated cells were harvested after four hours by scraping the petri dishes in the presence of lysis buffer (RLT buffer; Qiagen) and homogenized through Qiashredder columns (Qiagen). On average, 7×10 6  cells, grown to confluency in a 100 mm diameter petri dish, yielded 40 μg of total RNA from the RNEASY total RNA purification kit (Qiagen). RNA, in 20 mM Tris, 10 mM MgCl 2  buffer, pH 8 was incubated with 0.08 U/μl of RNase free DNase and 0.32 U/μl of RNase inhibitor (both from Boehringer Mannheim Biochemicals; Indianapolis Ind.) for 40 min at 37° C. and cleaned again using the RNEASY kit, which is important for removing small amounts of genomic DNA that can contribute to the fingerprints. RNA quantity was measured by spectrophotometry, and RNA samples were adjusted to 400 ng/μl in water. RNA samples were checked for quality and concentration by agarose gel electrophoresis and stored at −20° C.  
     [0210] For RNA fingerprinting, RAP-PCR was performed using standard protocols (McClelland et al., supra, 1994; Reverse transcription was performed on total RNA using four concentrations per sample (1000, 500, 250 and 125 ng per reaction) and a oligo d(T) primer (15-mer) (Genosys). RNA (5 μl) was mixed with 5 μl of buffer for a 10 μl final reaction volume containing 50 mM Tris, pH 8.3, 75 mM KCl, 3 mM MgCl 2 , 20 mM dithiothreitol (DTT), 0.2 mM of each dNTP, 0.5 μM of primer, and 20 U of MuLV-reverse transcriptase (Promega; Madison Wis.). RNA samples are checked for DNA contaminants by including a reverse transcriptase-free control in initial RAP-PCR experiments. The reaction was performed at 37° C. for 1 hr, after a 5 min ramp from 25° C. to 37° C. The enzyme was inactivated by heating the samples at 94° C. for 5 min, and the newly synthesized cDNA was diluted 4-fold in water.  
     [0211] PCR was performed after the addition of a pair of two different 10- or 11-mer oligonucleotide primers of arbitrary sequence; pair A: GP14 (GTAGCCCAGC; SEQ ID NO:) plus GP16 (GCCACCCAGA; SEQ ID NO:), pair B: Nucl+ (ACGAAGAAGAAGAG; SEQ ID NO:) plus OPN24 (AGGGGCACCA; SEQ ID NO:). In general, there are no particular constraints on the primers except that they contain at least-a few C or G bases, that the 3′ ends are not complementary with themselves or the other primer in the reaction, to avoid primer dimers, and that primer sets are chosen that are different in sequence so that the same parts of mRNA are not amplified in different fingerprints.  
     [0212] Diluted cDNAs (10 μl) were mixed with the same volume of 2×PCR mixture containing 20 mM Tris, pH 8.3, 20 mM KCl, 6.25 mM MgCl 2 , 0.35 mM of each dNTP, 2 μM of each oligonucleotide primer, 2 μCi α-( 32 P)-dCTP (ICN; Irvine Calif.) and 5 U AMPLITAQ DNA polymerase Stoffel fragment, (Perkin-Elmer-Cetus; Norwalk Conn.) for a 20 μl final reaction volume. Thermocycling was performed using 35 cycles of 94° C. for 1 min, 35° C. for 1 min and 72° C. for 2 min.  
     [0213] A 3.5 μl aliquot of the amplification products was mixed with 9 μl of formamide dye solution, denatured at 85° C. for 4 min, and chilled on ice. 2.4 pl was loaded onto a 5% polyacrylamide, 43% urea gel prepared with 1×TBE buffer containing 0.09 M Tris-borate, 0.002 M ethylene diamine tetraacetic acid (EDTA). The PCR products resulting from the four different concentrations of the same RNA template were loaded side by side on the gel.  
     [0214] Electrophoresis was performed at 1,700 V or at a constant power of 50-70 Watts until the xylene cyanol tracking dye reached the bottom of the gel (approximately 4 h). The gel was dried under vacuum and placed on Kodak BioMax X-Ray film for 16 to 48 hours.  
     [0215] For labeling of RAP-PCR products for use as targets to probe arrays, up to 10 pg of PCR product from RAP-PCR was purified using a QIAQUICK PCR Purification Kit (QIAGEN) which removes unincorporated bases, primers, and primer dimers under 40 base pairs. The DNA was recovered in 50 μl of 10 mM Tris, pH 8.3.  
     [0216] Random primed synthesis with incorporation of α-( 32 P)-dCTP was performed essentially as described in Example I. Briefly, 10% of the recovered fingerprint DNA, typically about 100 ng in 5 μl, was combined with 3 μg random hexamer oligonucleotide primer and 0.3 μg of each of the fingerprint primers in a total volume of 14 μl, which was boiled for 3 min and then placed on ice.  
     [0217] The hexamer/primer/DNA mix was mixed with 11 μl reaction mix to yield a 25 μl reaction containing 0.05 mM of three dNTP (minus dCTP), 50 μCi of 3000 Ci/mmol α-( 32 P)-dCTP (5 μl), 1× Klenow fragment buffer, containing 50 mM Tris-HCl, 10 mM MgCl 2 , 50 mM NaCl, pH 8.0, and 4 U Klenow fragment (Gibco-BRL). The reaction was performed at room temperature for 4 hrs. For maximum target length, the reaction was chased by adding 1 pl of 1.25 mM dCTP and incubated for 15 min at 25° C., followed by an additional 15 min incubation at 37° C. The unincorporated nucleotides, hexamers and primers were removed with the Qiagen Nucleotide Removal Kit (Qiagen) and the purified products were eluted using two aliquots of 140 μl of 10 mM Tris, pH 8.3.  
     [0218] For labeling of poly(A) +  mRNA and genomic DNA for use as a target, random hexamers were used to label poly(A) + -selected mRNA and genomic DNA. Genomic DNA (150 ng) was labeled using the same protocol used for labeling the RAP-PCR products described above. Poly(A) +  mRNA (1 μg) and 9 μg random hexamer in a volume of 27 μl were incubated at 70° C. for 2 min and chilled on ice. The RNA/hexamer mix was mixed with 23 μl master mix, which contained 10 μl 5×AMV reaction buffer, containing 250 mM Tris-HCl, pH 8.5, 40 mM MgCl 2 , 150 mM KCl, 5 mM DTT, 1 μl three dNTP, each 33 mM (dATP, dTTP, dGTP; minus dCTP), 2 μl AMV reverse transcriptase (20 units; Boehringer Mannheim) and 10 μl 3000 Ci/mmol α-( 32 P)-dCTP in a final volume of 50 μl. The reaction was incubated at room temperature for 15 min, ramped for 1 hour to 47° C., held at 47° C. for 1 hr, and chased with 1 μl of 33 mM dCTP for another 30 min at 47° C. The labeled products were purified as described above.  
     [0219] For hybridization to the array, four membranes were used, one membrane for each of two concentrations of RNA for each of the two RNA samples to be compared. The cDNA filters (Genome Systems) were washed in three changes of 2×SSC and 0.1% SDS in a horizontally shaking flat bottom container to reduce the residual bacterial debris. The first wash was carried out in 500 ml for 10 min at room temperature. The second and third washes were carried out in 1 liter of prewash solution, prewarmed to 55° C., for 10 min each wash.  
     [0220] For prehybridization, the filters were transferred to roller bottles and prehybridized in 60 ml prehybridization solution, prewarmed to 42° C., containing 6×SSC, 5× Denhardt&#39;s reagent, 0.5% SDS, 100 pg/ml fragmented, denatured salmon sperm DNA, and 50% formamide for 1-2 hrs at 42° C. in a hybridization oven.  
     [0221] For hybridization, the prehybridization solution was removed and 7 ml hybridization solution, prewarmed to 42° C., containing 6×SSC, 0.5% SDS, 100 μg/ml fragmented, denatured salmon sperm DNA, and 50% formamide, was added. To decrease the background hybridization due to repeats such as Alu and Line elements, sheared human genomic DNA was denatured in a boiling water bath for 10 min and immediately added to the hybridization solution to a final concentration of 10 μg/ml. 10 ng/ml poly(dA) was added to block oligo d(T) stretches in the radiolabeled target. Simultaneously, the labeled target, in a total volume of 280 μl, was denatured in a boiling water bath for 4 min and immediately added to the hybridization solution. The hybridization was carried out at 42° C. for 2-48 hrs, typically 18 hrs, in large roller bottles.  
     [0222] For the washes, the incubator oven temperature was set to 68° C. The hybridization solution was poured off and the membrane was washed twice with 50 ml 2×SSC and 0.1% SDS at room temperature for 5 min. The wash solution was then replaced with 100 ml 0.1×SSC and 0.1% SDS and incubated for 10 min at room temperature. For the further washes, the wash solution, containing 0.1×SSC and 0.1% SDS, was prewarmed to 68° C. The membranes were incubated 40 min in 100 ml of wash solution in the roller bottles, then the filters were transferred to horizontally shaking flat bottom containers and washed in 1 liter for 20 min under gentle agitation. The filters were transferred back to the roller bottles containing 100 ml 0.1×SSC and 0.1% SDS, prewarmed to 68° C., and incubated for 1 hr. The final wash solution was removed and the filters are briefly rinsed in 2×SSC at room temperature.  
     [0223] After washing, the membranes were blotted with 3MM paper, wrapped in SARAN wrap while moist, and exposed to X-ray film. The membranes were usually sufficiently radioactive that a one-day exposure with a screen revealed the top 1000 products on an array of 18,432 bacterial colonies carrying EST clones. Weaker targets or fainter hybridization events were visualized using an intensifying screen at −70° C. for a few days.  
     [0224] For confirmation of differential expression, low stringency RT-PCR was used. The initial confirmation of differential expression was the use of two RNA concentrations per sample. Only those hybridization events that indicated differential expression at both RNA concentrations in both RNA samples were relied upon.  
     [0225] More than 70% of the I.M.A.G.E. consortium clones have single pass sequence reads from the 5′ or 3′ end, or both, deposited in the GenBank database. In cases where there is no prior sequence information available, the clones can be ordered from Genome Systems and sequenced. Sequences were used to derive PCR primers of 18 to 25 bases in length using MacVector 6.0 (Oxford Molecular Group; Oxford UK). Generally, primers were chosen to generate PCR products of 50 to 250 base pairs and have melting temperatures of at least 60° C.  
     [0226] Reverse transcription was performed under the same conditions as in the RAP-PCR protocol described above, using an oligo-d(T) primer or a mixture of random 9-mer primers (Genosys). The PCR reaction was performed using the two pairs of specific primers described below (18 to 25-mers). The PCR conditions were the same as in the RAP-PCR fingerprint protocol except that 1.5 μM of each primer was used. A low stringency thermal profile was used: 94° C. for 40 sec, 47° C. for 40 sec, and 72° C. for 1 min, for 19, 22 and 25 cycles in three separate reaction tubes. The reactions were carried out in three sets of tubes at different cycle numbers because the abundance of the transcripts, the performance of the primer pairs, and the amplifiability of the PCR products can vary. PCR products were run under the same conditions as above on a 5% polyacrylamide and 43% urea gel. The gel was dried and exposed to X-ray film for 18 to 72 hours. Invariance among the other arbitrary products in the fingerprint was used as an internal control to indicate the reliability of the relative quantitation.  
     [0227] Primer pairs (Genosys) were used for confirmation of differential expression. For GenBank accession number H11520 (90 nucleotide product); primer A, AATGAGGGGGACAAATGGGAAGC (SEQ ID NO:); primer B, GGAGAGCCCTTCCTCAGACATGAAG (SEQ ID NO:). For TSC-22 gene (GenBank accession numbers U35048, H11073, H11161; 179 nucleotide product); primer A, TGACAAAATGGTGACAGGTAGCTGG (SEQ ID NO:); primer B, AAGTCCACACCTCCTCAGACAGCC (SEQ ID NO).  
     [0228] For GenBank accession number R48633 (178 nucleotide product); primer A, CCCAGACACCCAAACAGCCGTG (SEQ ID NO); primer B, TGGAGCAGCCGTGTGTGCTG (SEQ ID NO:).  
     [0229] The array analyzed contains 18,432  E. coli  colonies, each carrying a different I.M.A.G.E. consortium EST plasmid (www-bio.llnl.gov/bbrp/image/image.html), spotted twice on a 22×22 cm membrane (Genome Systems). The Genome Systems arrays are advantageous in that they contain by far the largest number of ESTs per unit cost. RNA fingerprinting for target preparation.  
     [0230] RAP-PCR amplifications were performed to look for differential gene expression in keratinocytes (HaCaT) when treated with EGF or TGF-β for four hours (Boukamp et al., supra, 1997). These experiments were designed to detect genes differentially regulated by EGF and TGF-β treatment in confluent keratinocytes. Using RAP-PCR, about 1% of the genes in normal or immortal keratinocytes responded to EGF, and fewer responded to TGF-β in this time frame.  
     [0231] Shown in FIG. 2 are RAP-PCR fingerprints of RNA from confluent keratinocytes treated with TGF-β or EGF using multiple RNA concentrations and two sets of arbitrarily chosen primers. Reverse transcription was performed with an oligo-dT primer on 250, 125, 62.5 and 31.25 ng RNA in lanes 1, 2, 3, and 4, respectively. RNA was from untreated, TGF-β treated or EGF treated HaCaT cells, as indicated. RAP-PCR was performed with two sets of primers, GP14 and GP16 (Panel A) or Nucl+ and OPN24 (Panel B). The sizes of the two differentially amplified RAP-PCR products are indicated with arrows (317 and 291 nucleotides).  
     [0232] In the first fingerprint shown in FIG. 2A, two differentially regulated products were detected, which were cloned and sequenced. The sizes of these two products, 291 and 317 nucleotides, are indicated with arrows (see FIG. 2A). The Genome Systems arrays used were chosen based on the presence of these two clones. This fingerprint was used to demonstrate that differentially regulated genes in an array can be identified without isolating, cloning and sequencing the RAP-PCR products. The fingerprint shown in FIG. 2A and the second fingerprint shown in FIG. 2B, which displayed no differential regulation in response to the treatments, were also used to demonstrate that fainter differentially regulated products not visible on the fingerprint gel could, nevertheless, be observed by the array approach.  
     [0233] The results obtained were highly reproducible. Using gel electrophoresis, there were no differences among the ˜100 bands visible in any of the fingerprints from a single treatment condition performed at different RNA concentrations (see FIG. 2). Similarly, more than 99% of the top 1000 clones hybridized by the targets derived from the fingerprint in FIG. 2A were visible at both input RNA concentrations. Furthermore, more than 98% of the products were the same between the two treatment conditions, plus and minus EGF, at a single RNA concentration. These results indicated high reproducibility among the top 1000 PCR products in the RAP-PCR amplification.  
     [0234] The untreated control and EGF-treated samples were further characterized. RAP-PCR fingerprints shown in FIG. 2 were converted into high specific activity radioactive targets by random primed synthesis using α-( 32 P)-dCTP as described above. For each of the two conditions, EGF treated and untreated, fingerprints generated from RNA at two different concentrations were converted to target by random primed synthesis for each of the two different fingerprinting primer pairs. These radioactively labeled fingerprint targets were then used to probe by hybridizing to a set of identical arrays each containing 18,432 I.M.A.G.E. consortium cDNA clones. As controls, total genomic DNA and total poly(A) +  mRNA were also labeled by random priming, as described above, and used as targets on identical arrays.  
     [0235] The RAP-PCR fingerprint targets, the total mRNA target and the genomic target were hybridized individually against replicates of a Genome Systems colony array. Genomic DNA was used as a blocking agent and as a competitor for highly repetitive sequences. Washing at 68° C. in 0.1×SSC and 0.1% SDS removed virtually all hybridization to known Alu elements on the membrane, presumably because Alu elements are sufficiently diverged from each other at this wash stringency.  
     [0236] Shown in FIG. 3 are autoradiograms from the same half of each membrane. All images presented are autoradiograms of the bottom half of duplicates of the same filter (Genome Systems) probed by hybridization with radiolabeled DNA. Panels A and B show hybridization of two RAP-PCR reactions generated using the same primers (GP14 and GP16) and derived from untreated (Panel A) or EGF treated (Panel B) HaCaT cells. Three double-spotted clones that show differential hybridization signals are marked on each array. The GenBank Accession numbers of the clone and the corresponding genes are H10045 and H10098, corresponding to vav-3 and AF067817 (square)(Katzav et al.,  EMBO J.  8:2283-2290 (1989); H28735, gene unknown, similar to heparan sulfate 3-O-sulfotransferase-1, AF019386 (circle)(Shworak et al.,  J. Biol. Chem.  272:28008-28019 (1997); and R48633, gene unknown (diamond).  
     [0237]FIG. 3 shows the results of hybridization of targets from these fingerprints to the arrays. As shown in FIGS. 3A and 3B, arrayed clones corresponding to the 291 nucleotide (vav-3, marked by square) and 317 nucleotide (similar to heparin sulfate N-sulfotransferase (N-HSST), marked by circle) RAP-PCR fragments are indicated. The sequences of these RAP-PCR fragments were determined. Also indicated on this array is a differentially regulated gene that could not be visualized on the original fingerprint gel (marked by diamond).  
     [0238] Comparing FIGS. 3A and 3B, a more than 10-fold down-regulation was observed for vav-3 upon treatment with EGF. The gene corresponding to H28735 was up-regulated more than 10-fold with EGF treatment. The gene corresponding to R48633 was up-regulated about 3-fold with EGF treatment. These changes in gene expression in response to EGF were independently confirmed by RT-PCR.  
     [0239] These results indicate that RAP-PCR samples a population of mRNAs largely independently of message abundance. This is because the low abundance class of messages has much higher complexity than the abundant class, making it more likely that the arbitrary primers will find good matches. Unlike differential display, RAP-PCR demands two such arbitrary priming events, possibly biasing RAP-PCR toward the complex class. Overall, these data suggest that the majority of the mRNA population in a cell (&lt;20,000 mRNAs) can be found in as few as ten RAP-PCR fingerprints. This result indicates that differential gene regulation can be detected by the combined fingerprinting and array approach even when the event cannot be detected using the standard gel electrophoresis approach.  
     [0240]FIG. 3C shows an array hybridized with a RAP-PCR target using the same RNA as in panel A but with a different pair of primers, Nucl+ and OPN24. As shown in FIG. 3C, using a different set of primers yields an entirely different pattern of hybridizing genes. FIG. 3D shows an array hybridized with a cDNA generated by reverse transcription of 1 μg poly(A) + -selected mRNA. FIG. 3E shows an array hybridized with human genomic DNA labeled using random priming.  
     [0241] The data were analyzed in a number of ways. First, estimates were made of the overlap between the clones hybridized by each target. In all pairwise comparisons between all of the different types of targets, there was less than 5% overlap among the 500 clones that hybridized most intensely (compare FIGS. 3A, 3B,  3 D, and  3 E). Of the top 500 clones hybridized by the genomic target, which included nearly all clones known to contain the Alu repeats, less than 5% overlapped with the top 500 clones hybridized by the fingerprint targets or the total poly(A) +  mRNA target. This indicated that, except for the case of a genomic target, there was no significant hybridization to dispersed repeats. The overlap among the clones hybridized by the two RAP-PCR fingerprints generated with different primers was less than 3%, and the overlaps of either fingerprint with the poly(A) +  mRNA target were both less than 3%. Thus, most of the cDNAs detected using a target from the fingerprints could not be detected using the total mRNA target. These results indicate that RAP-PCR samples a population of mRNAs largely independently of message abundance. This is because the low abundance class of messages has much higher complexity than the abundant class, making it more likely that the arbitrary primers will find good matches. Unlike differential display, RAP-PCR demands two such arbitrary priming events, possibly biasing RAP-PCR toward the complex class. Overall, these data suggest that the majority of the mRNA population in a cell (&lt;20,000 mRNAs) can be found in as few as ten RAP-PCR fingerprints.  
     [0242] A total of 30 differentially hybridizing cDNA clones were detected among about 2000 hybridizing colonies using targets derived from both sets of arbitrary primers (FIG. 2) at a threshold of about three-fold differential hybridization. Twenty-two of these differentially hybridizing clones displayed differential hybridization at both RNA concentrations. These 22 were further characterized by RT-PCR. Differentially expressed genes exhibiting greater than a two-fold difference in expression in response to EGF treatment are shown in Table 1. For the results shown in Table 1, differential expression was confirmed by low stringency RT-PCR. The left column gives the accession numbers of the EST clones (5′ or 3′, or both when available). The right column gives the corresponding gene or the closest homolog. In cases of very low homologies, the gene is considered unknown. The cutoff for homology was p&lt;e-20 in tblastx.  
               TABLE 1                          Genes Regulated More than Two-fold After EGF       Treatment of HaCaT Keratinocytes.                             Accession number   Gene name                       Up-regulated               H11520 (3′)   unknown           H11161 (5′)/H11073 (3′)   TSC-22 (U35048)           R48633 (5′)   unknown           H28735 (3′)   similar to heparan sulfate 3-O-               sulfotransferase-1 precursor               (AF019386)           H25513 (5′)/H25514 (3′)   Fibronectin receptor α-subunit               (M13918)           H12999 (5′)/H05639 (3′)   similar to Focal adhesion kinase               (FAK2) (L49207)           H15184 (5′)/H15124 (3′)   ray gene (X79781)           H25195 (5′)/H24377 (3′)   X-box binding protein-1 (XBP-1)               (M31627)           H23972 (‘ ’)   unknown           H27350 (5′)   CPE-receptor (hCPE-R) (AB000712)           R75916 (5′)   similar to semaphorin C (X85992)           Down-regulated           R73021 (5′)/R73022 (3′)   epithelium-restricted Ets               protein ESX (U66894)           H10098 (5′)/H10045 (3′)   vav-3 (AF067817)                      
 
     [0243] The eight false-positive clones that appeared to be regulated at only one concentration were further characterized. Of these eight, five false-positive clones showed differential hybridization at one concentration but were present and not regulated on the membranes for the other concentration. The most likely source of this type of false-positive is the membranes. Although each clone is spotted twice, it is possible that occasionally one membrane received substantially more, or less, DNA in both spots than the other three membranes for these clones. However, this potential difference was easily detected and is rare, occurring only five times in over 2000 clones. The other three false-positive clones hybridized under only one treatment condition and at only one RNA concentration used for RAP-PCR. These three false-positive clones could be differentially expressed genes or could be false-positives from variable PCR products. However, the number of false positives was very low and were easily identified by comparing the results of two targets derived from PCR of different starting concentrations of RNA.  
     [0244] Differential expression was confirmed using low stringency RT-PCR. Only those hybridization events that indicated differential expression at both input RNA concentrations were further characterized. For confirmation of differential expression, RT-PCR was used with specific targets rather than Northern blots, which are much less sensitive than RT-PCR, because it was expected that many of the mRNAs would be rare and in low abundance. One of the advantages of using the arrays from the I.M.A.G.E. consortium is that more than 70% of the clones have single pass sequence reads from the 5′ or 3′ end, or both, deposited in the GenBank database.  
     [0245] Clones for which some sequence is available in the database were chosen for further characterization. Five of the 22 ESTs representing differentially regulated genes on the array had not been sequenced and two of the remaining 17 ESTs were from the same gene. The remaining 15 unique sequenced genes were aligned with other sequences in the database in order to derive a higher quality sequence from multiple reads and longer sequence from overlapping clones. The UniGene database clusters human and mouse ESTs that appear to be from the same gene (Schuler,  J. Mol. Med.  75:694-698 (1997)). This database greatly aids in the process of assembling a composite sequence from different clones of the same mRNA (http://www.ncbi.nlm.nih.gov/UniGene/index.html). These composite sequences were then used to choose primers for RT-PCR.  
     [0246] For each gene, two specific primers were used in RT-PCR under low stringency conditions similar to those used to generate RAP-PCR fingerprints. In addition to the product of interest, a pattern of arbitrary products was generated, which is largely invariant and behaves as an internal control for RNA quality and quantity, and for reverse transcription efficiency (Mathieu-Daude et al., supra, 1998). The number of PCR cycles was adjusted to between 14 to 25 cycles, according to the abundance of the product, in order to preserve the differences in starting template mRNA abundances. This is necessary because rehybridization of abundant products during the PCR inhibits their amplification, and the difference in product abundances diminishes as the number of PCR cycles increases (Mathieu-Daude et al.,  Nucleic Acids Res.  24:2080-2086 (1996)).  
     [0247] Low stringency RT-PCR experiments confirmed the differential expression of the two transcripts that were identified in the RAP-PCR fingerprints of FIG. 2A and showed differential hybridization to the cDNA array (compare FIG. 3A versus  3 B). One of these differentially expressed genes corresponds to a new family member of the vav protooncogene family (Katzav et al., supra, 1989; Katzav,  Crit. Rev. Oncog.  6:87-97 (1995); Bustelo,  Crit. Rev. Oncog.  7:65-88 (1996); Romero and Fischer,  Cell Signal.  8:545-553 (1996)). The other differentially expressed gene has homology to heparan sulfate 3-O-sulfotransferase-1 (Shworak et al., supra, 1997).  
     [0248] The other 13 differentially expressed were also tested and 11 were confirmed using low stringency RT-PCR. Some of the differentially expressed genes are shown in FIG. 4. Reverse transcription was performed at two RNA concentrations (500 ng, left column; 250 ng, right column). The reaction was diluted 4-fold in water and one fourth was used for low stringency RT-PCR at different cycle numbers. The RT-PCR products were resolved on polyacrylamide-urea gels. Shown are bands for the control (22 cycles); for GenBank accession number H11520 (22 cycles); for TSC-22, corresponding to GenBank accession numbers H11073 and H11161 (19 cycles) (Jay et al.,  Biochem. Biophys. Res. Commun.  222:821-826 (1996); Dmitrenko et al.,  Tsitol. Genet.  30:41-47 (1996); Ohta et al.,  Eur. J. Biochem.  242:460-466 (1996)); and for GenBank accession number R48633 (19 cycles). Genes corresponding to H11520 and TSC-22 are up-regulated about 8-10 fold with EGF treatment. The gene corresponding to R48633 is up-regulated about 3-fold with EGF treatment.  
     [0249] Of the two differentially expressed genes that were not confirmed, one proved unamplifiable. The other gene gave a product but appeared to not be differentially regulated when analyzed by RT-PCR.  
     [0250] RAP-PCR targets were very effective at detecting rare, low abundance mRNAs. Each fingerprint hybridized to a set of clones almost entirely different from the set hybridized by a target derived from poly(A) + -selected mRNA (see FIG. 3). In addition, numerous other primer pairs, membranes, and sources of RNA consistently showed less than a 5% overlap between clones hybridized by any two fingerprints, or between a fingerprint and a total poly(A) + -selected cDNA target. Detection of differentially expressed vav-3 mRNA, which is a new member of the vav oncogene family, was attempted using a Northern blot of poly(A) + -selected RNA. Despite being able to detect serially diluted vector down to the equivalent of a few copies per cell, vav-3 mRNA was undetectable on the Northern blot, whereas RT-PCR confirmed expression. A G3PDH control was used to confirm that the conditions used in the Northern blot could detect a control gene. Therefore, vav-3 appears to be a low abundance message that is represented in a RAP-PCR fingerprint as a prominent band.  
     [0251] The frequency of homologs of cDNAs detected by the RAP-PCR targets in the EST database was determined (&gt;98% identity). This was compared to the frequency of homologs for a random set of other cDNAs on the same membrane. If the RAP-PCR fingerprints were heavily biased towards common mRNAs, then many would occur often in the EST database because it is partly derived from cDNA libraries that are not normalized or incompletely normalized. However, the cDNAs detected by RAP-PCR had frequencies in the EST database comparable to the frequencies for randomly selected cDNAs, including cases where the clone was unique in the database. These results indicate that sampling by arbitrarily sampled targets generated by RAP-PCR is at least as good as random sampling of the partly normalized libraries used to construct the array, and very different from that obtained for a target such as total mRNA target.  
     [0252] These results demonstrate that an arbitrarily sampled target generated using RT-PCR and arbitrary primers can detect genes differentially expressed in response to EGF.  
     EXAMPLE III  
     An Arbitrarily Sampled Target Generated by Differential Display Detects Genes Differentially Expressed in Response to EGF  
     [0253] This example shows the use of differential display to generate an arbitrarily sampled target and detection of differentially expressed genes responsive to EGF.  
     [0254] RNA was prepared from the human keratinocyte cell line HaCaT as described in Example II. Briefly, cells were grown to confluence and maintained at confluence for 2 days. The medium was changed 1 day prior to the experiment. EGF (Gibco-BRL) was added at 20 ng/ml. Treated and untreated cells were harvested after 4 hrs and total RNA was prepared with the RNEASY total RNA purification kit (Qiagen) according to the manufacturer&#39;s protocol. To remove remaining genomic DNA, the extracted total RNA was treated with RNase-free DNase (Boehringer Mannheim) and cleaned again using the RNEASY kit. The purified RNA was adjusted to 400 ng/μl in water and checked for quality by agarose gel electrophoresis.  
     [0255] For standard differential display, differential display was performed using the materials supplied in the RNAIMAGE kit (GenHunter Corporation; Nashville Tenn.), AMPLITAQ DNA polymerase (Perkin-Elmer-ABI; Foster City Calif.) and α-( 32 P)-dCTP according to the manufacturer&#39;s protocol, except that each RNA template was used at four different concentrations, 800, 400, 200 and 100 ng per 20 μl reaction, with each anchored oligo(dT) primer (0.2 μM). The PCR reaction contained 2 μM dNTPs, for a total of 4 μM, including the carryover from the cDNA mix, 0.2 μM each primer, and one tenth of the newly synthesized cDNA, corresponding to 80, 40, 20 and 10 ng RNA. The anchored oligo(dT) primers were used in all possible combinations with four different arbitrary primers. The anchored oligo(dT) primers used were H-T 11 G (HTTTTTTTTTTTG; SEQ ID NO:); H-T 11 A (HTTTTTTTTTTTA; SEQ ID NO:); and H-T 11 C (HTTTTTTTTTTTC; SEQ ID NO:), where H is AAGC, which is an arbitrary sequence used as a clamp to ensure the primers stay in register and have a high Tm at subsequent PCR steps. The arbitrary primers used were H-AP1 (AAGCTTGATTGCC; SEQ ID NO:); H-AP2 (AAGCTTCGACTGT; SEQ ID NO:); H-AP3 (AAGCTTTGGTCAG; SEQ ID NO:); and H-AP4 (AAGCTTCTCAACG; SEQ ID NO:).  
     [0256] For modified differential display, reverse transcription was performed using four different concentrations of each RNA template, 1000, 500, 250 and 125 ng per 10 μl reaction. The reaction mix contained 1.5 μM oligo(dT) anchored primers AT 15 A, GT 15 G, and T 13 V, 50 mM Tris, pH 8.3, 75 mM KCl, 3 mM MgCl 2 , 20 mM DTT, 0.2 mM each dNTP, 8 U RNase inhibitor (Boehringer Mannheim) and 20 U MuLV reverse transcriptase (Promega). The anchored primers were AT 15 A (ATTTTTTTTTTTTTTTA; SEQ ID NO:); GT 15 G (GTTTTTTTTTTTTTTTG; SEQ ID NO:); and T 13 V (TTTTTTTTTTTTTV; SEQ ID NO:; where V is A, G or C)). The reaction mix was ramped for 5 min from 25° C. to 37° C., held at 37° C. for 1 hr, and finally the enzyme was inactivated at 94° C. for 5 min. The newly synthesized cDNA was diluted 4-fold in water.  
     [0257] The PCR was performed after adding 10 μl of reaction mix to 10 μl of the diluted cDNAs, corresponding to 250, 125, 62.5 and 31.25 ng of RNA, to yield a 20 μl final reaction volume containing 2 μM anchored oligo(dT) primer, 0.4 μM arbitrary primer, either KA2 (GGTGCCTTTGG; SEQ ID NO:) or OPN28 (GCACCAGGGG; SEQ ID NO:), 2.5 units AMPLITAQ DNA polymerase Stoffel fragment (Perkin Elmer-ABI), 2 μCi α-( 32 P)-dCTP, 175 μM each dNTP, 10 mM Tris, pH 8.3, 10 mM KCl, and 3.125 mM MgCl 2 . These concentrations do not include the carryover from the reverse transcription reaction. The reactions were thermocycled for 35 cycles of 94° C. for 40 sec, 40° C. for 1 min and 40 sec, and 72° C. for 40 sec.  
     [0258] An aliquot of the PCR products resulting from the four different concentrations of the same RNA template were displayed side by side on a 5% polyacrylamide gel and visualized by autoradiography as described in Example II.  
     [0259] For labeling of differential display products for use as targets to probe arrays, random primed labeling of the differential display products was performed as described in Example II. The differential display PCR reactions (14 pl) were purified using a QIAQUICK PCR Purification Kit (Qiagen) and the DNA was recovered in 50 μl 10 mM Tris, pH 8.3. Random primed synthesis was performed using a standard protocol. Briefly, 5 μl of the recovered differential display products were combined with 3 μg random hexamers, boiled for 3 min and placed on ice. The hexamer/DNA mix was combined with the reaction mix to yield a 25 μl reaction containing 0.05 mM three dNTPs (minus dCTP), 50 μCi of 3000 Ci/mmol α-( 32 P)-dCTP, 1× Klenow fragment buffer, and 4 U Klenow fragment (Gibco-BRL). The reaction was performed at room temperature for 4 hrs, chased for 15 min at room temperature by adding 1 μl of 1.25 mM dCTP, and incubated for an additional 15 min at 37° C. The unincorporated nucleotides and hexamers were removed with the Qiagen Nucleotide Removal Kit and the purified products were eluted using two aliquots of 140 μl 10 mM Tris, pH 8.3.  
     [0260] Hybridization to the array was performed essentially as described in Examples I and II. Briefly, the cDNA membranes (Genome Systems) were prewashed in three changes of prewash solution, containing 2×SSC and 0.1% SDS, in a horizontally shaking flat bottom container to reduce the residual bacterial debris. The first wash used 500 ml of prewash buffer for 10 min at room temperature. The second and third washes were each carried out in 1 liter of prewash solution, prewarmed to 55° C., for 10 min.  
     [0261] The membranes were transferred to large roller bottles and prehybridized in 60 ml prehybridization solution, prewarmed to 42° C., containing 6×SSC, 5× Denhardt&#39;s reagent, 0.5% SDS, 100 μg/ml fragmented, denatured salmon sperm DNA, and 50% formamide for 1-2 hrs at 42° C.  
     [0262] The prehybridization solution was removed, and 10 ml hybridization solution, prewarmed to 42° C. and containing 6×SSC, 0.5% SDS, 100 μg/ml fragmented, denatured salmon sperm DNA and 50% formamide, was added to the bottles. To decrease the background hybridization due to repeats such as Alu and Line elements, sheared human genomic DNA was denatured in a boiling water bath for 10 min and immediately added to the hybridization solution to a final concentration of 10 μg/ml. An aliquot of 10 ng/ml poly(dA) was added to block oligo (dT) stretches in the radiolabeled target. Simultaneously, the labeled target was denatured in a boiling water bath for 4 min and immediately added to the hybridization solution. The hybridizations-were carried out at 42° C. for 18-20 hrs.  
     [0263] Following hybridization, the hybridization solution was poured off and the membranes were thoroughly washed in six changes of wash solution, including a transfer of the membranes from the roller bottles to a horizontally shaking flat bottom container and back to the roller bottles, over 2-3 hrs. The stringency of the washes was increased stepwise from 2×SSC and 0.1% SDS at room temperature to 0.1×SSC and 0.1% SDS at 64° C. The separate washes were maintained at exactly the same indicated temperatures for all of the membranes. The last high stringency wash was at least 40 min to ensure exactly equilibrated temperatures in all bottles. The final wash solution was removed, and the membranes were briefly rinsed in 2×SSC at room temperature, blotted with 3MM paper, wrapped in SARAN wrap while moist, and placed against Kodak Biomax film (Eastman-Kodak; Rochester, N.Y.).  
     [0264] Differential expression was confirmed using low stringency RT-PCR. The first level of confirmation was the use of two RNA concentrations per sample. Only those hybridization events that indicated differential expression at both RNA concentrations in both RNA samples were further characterized.  
     [0265] Nucleotide sequences, which were available from Genome Systems, the commercial source of the array, or were sequenced, were used to derive PCR primers of 18 to 25 bases in length using MacVector 6.0 (Oxford Molecular Group). Generally, primers were chosen that generate PCR products of 100 to 250 base pairs, have melting temperatures of at least 60° C., and were preferably located close to the polyadenylation site of the mRNA so as to reduce the chance of sampling family members.  
     [0266] Reverse transcription was performed on total RNA using two RNA concentrations per sample and an oligo-(dT 15 ) primer (TTTTTTTTTTTTTTT; SEQ ID NO:; Genosys). The reactions contained 100 and 50 ng per liter total RNA, 0.5 μM oligo-(dT 15 ) primer (SEQ ID NO:), 50 mM Tris, pH 8.3, 75 mM KCl, 3 mM MgCl 2 , 20 mM DTT, 0.2 mM of each dNTP, 0.8 U/μl RNase inhibitor (Boehringer Mannheim) and 2 U/μl of MuLV-reverse transcriptase (Promega). The reactions were ramped for 5 min from 25° C. to 37° C. and held at 37° C. for 1 hr. The enzyme was inactivated by heating the reactions at 94° C. for 5 min and the newly synthesized cDNA was diluted 4-fold in water.  
     [0267] Diluted cDNAs (10 μl) were mixed with 2×PCR mixture containing 20 mM Tris, pH 8.3, 20 mM KCl, 6.25 mM MgCl 2 , 0.35 mM of each dNTP, 3 μM of each specific primer, 2 μCi α-( 32 P)-dCTP (ICN, Irvine, Calif.) and 2 U AMPLITAQ DNA polymerase Stoffel fragment (Perkin-Elmer-Cetus) for a 20 μl final reaction volume. A low stringency thermal profile was used: 94° C. for 40 sec, 40° C. for 40 sec, and 72° C. for 1 min, for 17 and 19 cycles in separate tubes. The reaction was carried out in two sets of tubes at different cycle numbers because the abundance of the transcripts, the performance of the primer pairs and the amplifiability of the PCR products can vary. PCR products were run under the same conditions as described above on a 5% polyacrylamide and 43% urea gel. The gel was dried and placed for 18 to 72 hours on a phosphoimager screen and read with a STORM phosphoimager (Molecular Dynamics; Sunnyvale Calif.). Invariance among the other arbitrary products in the fingerprint was used as an internal control to indicate the reliability of the relative quantitation. The gene-specific products from four sets of reactions per differentially regulated gene were quantitated using IMAGEQUANT Software (Molecular Dynamics).  
     [0268] Primer pairs were used to confirm differential expression.  
     [0269] For GenBank accession number R72714 (Egr-1)(155 nt product); primer A, CACGTCTTGGTGCCTTTTGTGTG (SEQ ID NO:); primer B, GAAGCTCAGCTCAGCCCTCTTCC (SEQ ID NO:).  
     [0270] For GenBank accession number H14529 (ACTB, β-actin)(174 nt product); primer A, CCAGGGAGACCAAAAGCCTTCATAC (SEQ ID NO:); primer B, CACAGGGGAGGTGATAGCATTGC (SEQ ID NO:).  
     [0271] For GenBank accession number H27389 (A+U-rich element RNA binding factor)(144 nt product); primer A, GTGCTTTTCAAAGATGCTGCTAGTG (SEQ ID NO:); primer B, GCTCAATCCACCCACAAAAACC (SEQ ID NO:).  
     [0272] For GenBank accession number H05545 (protein phosphatase 2A catalytic subunit)(141 nt product); primer A, TCCTCTCACTGCCTTGGTGGATG (SEQ ID NO:); primer B, CACAGCAAGTCACACATTGGACCC (SEQ ID NO:).  
     [0273] For GenBank accession number H27969 (103 nt product); primer A, CCAAAGACATTCAGAGGCATGG (SEQ ID NO:); primer B, GAGGTGGGGAAGGATACAGCAG (SEQ ID NO:).  
     [0274] For GenBank accession number R73247 (inositol tris phosphate kinase)(168 nt product); primer A, GAAAAGGGTTGGGGAGAAGCCTC (SEQ ID NO:); primer B, TCTCTAGCGTCCTCCATCTCACTGG (SEQ ID NO:).  
     [0275] For GenBank accession number H21777 (α-tubulin isoform 1) (155 nt product); primer A, ACAACTGCATCCTCACCACCCAC (SEQ ID NO:); primer B, GGACACAATCTGGCTAATAAGGCGG (SEQ ID NO:).  
     [0276] Total RNA was obtained from immortalized HaCaT keratinocytes, treated and untreated with EGF, as described in Example II (Boukamp et al., supra, 1997). The first differential display protocol tried was the RNAimage kit 1 (cut G50′; GenHunter. The anchor primers, oligo(dT)-G (H-T 11 G; SEQ ID NO:), oligo(dT)-C (H-T 11 C; SEQ ID NO:) or oligo(dT)-A (H-T 11 A; SEQ ID NO:), were used for reverse transcription, and then each cDNA was used for PCR in combination with four different arbitrary primers, H-AP1 (SEQ ID NO:), H-AP2 (SEQ ID NO:), H-AP3 (SEQ ID NO:) and H-AP4 (SEQ ID NO:).  
     [0277] As shown in FIG. 5, the fingerprints were resolved on a denaturing acrylamide gel to determine the quality of the reactions. Differential display reactions were performed using the RNAIMAGE kit protocol (GenHunter Corporation) according to the manufacturer&#39;s suggestion except that four different starting concentrations of 800, 400, 200 and 100 ng of total RNA were used. One tenth of this material was then used for PCR. The anchored oligo(dT) primer H-T 11 C (SEQ ID NO:) was used with two different arbitrary primers, H-AP3 (SEQ ID NO:) and H-AP4 (SEQ ID NO:), as indicated. The arbitrary primer H-AP4 (SEQ ID NO:) was used with two different anchored oligo(dT) primers, H-T 11 C (SEQ ID NO:) and H-T 11 A (SEQ ID NO:). The reactions that share either the arbitrary primer or the anchored oligo(dT) primer showed almost no visible overlap in the visible bands.  
     [0278]FIG. 5B shows differential display using a different set of primers. Differential display was performed using the arbitrary primer KA2 (SEQ ID NO:) with three different anchored oligo(dT) primers, T 13 V (SEQ ID NO:), AT 15 A (SEQ ID NO:), and GT 15 G (SEQ ID NO:), as indicated. The differential display protocol was adjusted to yield more mass and a higher complexity of the generated products. The starting concentrations of RNA were 1000, 500, 250 and 125 ng. One fourth of this material was then used for PCR. As observed in FIG. 5A, using different oligo(dT) anchored primers changes the pattern of the displayed bands almost entirely.  
     [0279] The fingerprints generated about 30 to 50 clearly visible products (see FIG. 5A). Fingerprints were generally reproducible in the range from 100 to 800 ng of total mRNA used in these experiments, with very few RNA concentration dependent products. Three of the most reproducible fingerprints that shared either a oligo(dT) anchored primer or an arbitrary primer (FIG. 5A) were radiolabeled by random priming in the presence of three unlabeled dNTPs and α-( 32 P)-dCTP, and each was used to probe identical arrays of 18,000 double spotted  E. coli  colonies carrying ESTs from the I.M.A.G.E. consortium. The arrays were hybridized and washed as described above.  
     [0280] The kit protocol used 0.2 μM of the arbitrary primer and 4 μM dNTPs compared to 1 μM primers and 200 μM dNTPs used in the RAP-PCR protocol described in Example II. The fingerprint reaction contained less than 40 ng of product in 20 μl, presumably because of limiting components. This was about five times less DNA than used in the method described in Example II. For this reason, it took about ten days with an intensifying screen in order to obtain an adequate exposure of X-ray film. Approximately 500 products were easily discernible with each target after a sufficient exposure. The number of reliably observable genes is usually increased by at least two-fold or more when using a phosphoimager screen, indicating the greater sensitivity of phosphoimaging compared to X-ray film. Furthermore, pooling of separate labeled fingerprints into the same target can increase throughput even further.  
     [0281] In order to reduce the exposure time for target hybridization to arrays, experiments were performed at the higher concentration of primer and dNTPs described in Example II using RAP-PCR protocols (FIG. 5B). These experiments yielded the expected increase in product mass and a corresponding reduction in exposure times for arrays.  
     [0282] The selectivity of oligo(dT) primers was determined using different anchor bases. As shown in FIG. 6, differential display reactions were hybridized to cDNA arrays. The differential display products generated as described in FIG. 5A, with the primers GT 15 G (SEQ ID NO:) and KA2 (SEQ ID NO:) from untreated (FIG. 6A) and EGF treated (FIG. 6B) HaCaT cells, were labeled by random priming and hybridized to cDNA arrays. A section representing less than 5% of a membrane is shown with a differentially regulated gene indicated by an arrow. FIG. 6C shows hybridization of differential display products generated with the primers AT 15 A (SEQ ID NO:) and KA2 (SEQ ID NO:) from untreated HaCaT cells. Comparing FIG. 6A versus  6 C, there is a significant overlap of hybridization signals that were not obvious from the polyacrylamide display (compare to FIG. 5B, lanes AT 15 A/KA2 versus GT 15 G/KA2).  
     [0283] When the arbitrary primer was changed while keeping the same anchor primer, the pattern of clones hybridized changed almost entirely, with typically less than 5% overlap between any two fingerprints. In contrast, targets containing the same arbitrary primer and different anchored primers shared about 30% of the clones to which they hybridized. FIGS. 6A and 6C show examples of such shared products from a small portion of an array.  
     [0284] Similar observations were made using fingerprints generated under a wide variety of conditions, including the protocols and primers from the GenHunter kit, modified protocols, and protocols using primers independent of those in the GenHunter kit. The possibility of this overlap being due to repeats was excluded by the use of genomic and total mRNA targets against the same membranes.  
     [0285] The overlap among targets that had different anchored primers but shared the same arbitrary primer was not reflected in any noticeable similarity in the fingerprint products when resolved on a denaturing polyacrylamide gel. For example, the targets used in FIGS. 6A and 6C are shown in FIG. 5B and show no easily discerned similarities, despite having 30% of the products in common. Many of the shared products were among the most intensely hybridizing clones on the array. Therefore, some of the products visible on the gel could share the arbitrary primer at one end but, during PCR, the products are preferentially primed at multiple different locations in the opposite direction by the different anchored primers. This would result in fingerprints that had little or no similarity in a polyacrylamide display while being compatible with the observation that targets with the same arbitrary primer but different anchored primers overlap by 30% in the clones to which they hybridize.  
     [0286] Shared products are a general phenomenon for anchored fingerprints that share an arbitrary primer under a fairly wide range of conditions. Overlap among fingerprints can be avoided by not using the same arbitrary primer with different anchored primers.  
     [0287] Comparison of the pattern of hybridizing clones with that generated by total genomic DNA indicated that the clones hybridizing to a target generated by the GenHunter fingerprint did not generally contain the Alu repetitive element that occurs in a few percent of mRNA 3′ untranslated regions (UTRs). The clones hybridized by the target did not overlap significantly with clones hybridized by a total cDNA target derived from reverse transcription of poly(A) +  mRNA, indicating that the genes sampled were not heavily biased towards the most abundant RNAs. These results are consistent with results obtained using only arbitrary primers for fingerprinting (see Example II) and indicate that arbitrary priming combined with anchored oligo(dT) priming can be used to monitor rare genes in cDNA arrays. These results also confirm that RAP-PCR and differential display are not heavily biased toward abundant transcripts.  
     [0288] Among over 2000 clones surveyed for differential gene expression between untreated and EGF treated HaCaT cells, there were 29 different clones that appeared to clearly reflect differential expression at one RNA concentration. The 12 clones having the highest signal to noise ratio and differential expression ratio were chosen and specific primers were designed for RT-PCR. An example of one of these differentially expressed genes is indicated by an arrow in FIG. 6A versus  6 B.  
     [0289] Differential expression of at least 1.5-fold was confirmed for seven genes, which are shown in FIG. 7. Reverse transcription was performed at twofold different RNA concentrations. The reactions were diluted 4 fold in water and low stringency PCR was performed at different cycle numbers. The amount of input RNA/cDNA for each PCR reaction was 125 ng, left column and 250 ng, right column. The reactions shown in FIG. 7 were carried out for 10 cycles and resolved on polyacrylamide-urea gels. Shown are products for the control (unregulated) and genes differing by at least 1.6-fold. The regulated genes shown correspond to GenBank accession numbers R72714, H14529, H27389, H05545, H27969, R73247, and H21777.  
     [0290] The regulation of the genes shown in FIG. 7 are summarized in Table 2. Identified genes regulated by four hr treatment with EGF, corresponding GenBank accession numbers, and the fold-increase in expression relative to untreated cells are shown.  
               TABLE 2                          EGF Regulated Genes.                                         Fold                   Up-regulation           Gene   Accession #   by EGF                       EGR1   R72714, X52541   8.3 ± 3.4           ACTB, beta-actin   H14529, M10277   2.0 ± 0.3           A + U-rich element RNA   H27389, D89092,   1.9 ± 0.3           binding factor   D89678           Protein phosphatase 2A   H05545, J03804   1.6 ± 0.4           catalytic subunit           Unknown   D31765, H27969   1.6 ± 0.4           Inositol tris phosphate   R73247, U51336   1.6 ± 0.3           kinase           Alpha-tubulin isoform 1   H21777, K00558   1.6 ± 0.3                      
 
     [0291] Egr-1 was previously known to be differentially regulated by EGF in other cell types (Iwami et al.,  Am. J. Physiol.  270:H2100-H2107 (1996); Kujubu et al.,  J. Neurosci. Res.  36:58-65 (1993); Cao et al.,  J. Biol. Chem.  267:1345-1349 (1992); Ito et al.,  Oncogene  5:1755-1760 (1990)). The observations of changes in β-actin and α-tubulin expression are likely associated with the dramatic change in morphology these cells undergo after EGF treatment. Regulation of β-actin and α-tubulin genes by EGF has been observed in other cell types (Torok et al.,  J. Cell Physiol.  167:422-433 (1996); Hazan and Norton,  J. Biol. Chem.  273:9078-9084 (1998); Shinji et al.,  Hepatogastroenterology  44:239-244 (1997); Ball et al.,  Cell Motil. Cytoskeleton  23:265-278 (1992)). These observations independently validate the treatments and the method used to detect differential expression. The regulation of protein phosphatase 2A mRNA has not previously been observed but is consistent with the role of this protein in transduction of the EGF signal (Chajry et al.,  Eur. J. Biochem.  235:97—102 (1996)). Similarly, the gene associated with the metabolism of inositol phosphates had not previously been shown to be regulated by EGF but such regulation is consistent with the previous observation of increases in the compounds generated by this enzyme after EGF treatment in another ectodermal cell type (Contreras,  J. Neurochem.  61:1035-1042 (1993)). Regulation of two other genes by EGF, an unknown gene, with GenBank accession number H27969, and an RNA binding protein, with GenBank accession number D89692, was not previously reported in any cell type. GenBank accesssion number D31765 corresponds to KIAA0061.  
     [0292] Five other genes were not confirmed to be regulated when RT-PCR was used. The number of false positives can vary from experiment to experiment and depends on the quality of the fingerprints and on the quality of the commercially available membranes. The number of false positives can be limited by using two RNA concentrations on arrays before confirmation by RT-PCR, as described in Example II. These experiments involved only a single concentration because the primary purpose was to determine the efficiency of coverage and overlap among targets made by the oligo(dT)-X anchored priming method. Nevertheless, over half of the differentially hybridizing clones observed at one concentration correspond to differentially expressed genes. When two array hybridizations were performed for each treatment at two different input template concentrations, the error rate was well below 10%.  
     [0293] These results demonstrate that an arbitrarily sampled target generated using differential display and arbitrary primers can detect genes differentially expressed in response to EGF.  
     [0294] Throughout this application various publications have been referenced. The disclosures of these publications in their entireties are hereby incorporated by reference in this application in order to more fully describe the state of the art to which this invention pertains.  
     [0295] Although the invention has been described with reference to the examples provided above, it should be understood that various modifications can be made without departing from the spirit of the invention. Accordingly, the invention is limited only by the claims.  
    
     
       
         1 
         
           
             85  
           
           
             1  
             395  
             DNA  
             Homo sapiens  
             
               unsure  
               (82)  
             
             
               unsure  
               (115)..(116)  
             
             
               unsure  
               (314)  
             
             
               unsure  
               (350)  
             
             
               unsure  
               (359)  
             
             
               unsure  
               (383)  
             
           
            1 

tttttttttt acaacaatgc agtcatttat ttattgagta tgtgcacatt atggtattat     60 

tactatactg attatattta anaagtgact tctaattaga aaatgtatcc aaaannaaaa    120 

cagcagatat acaaaattaa agagacagaa gatagacatt aacagataag gcaacttata    180 

cattgaggaa tccaaatcca atacatttaa acatttggga aatgaggggg acaaatggga    240 

agccagatca aatttgtgta aaactattca gtatgtttcc cttggcttca tgtctgagga    300 

agggctctcc cttncaatgg gggatggaca aactccaaat gccacacaan tgtttaacng    360 

gtatactagg tttcacactg ggnacggggg ttaaa                               395 

 
           
             2  
             389  
             DNA  
             Homo sapiens  
             
               unsure  
               (230)  
             
             
               unsure  
               (384)  
             
           
            2 

acacagcccc ccgcccagcc agcatcgcag ggcttcaggg accaaccgca tagctgccta     60 

tgcccccgca gaactggctg ctgcgtgtga actgaacaga cggagaagat gtgctaggga    120 

gaatctgcct ccacagtcac ccatttcatt gctcgctgcg aaagagacgt gagactgaca    180 

tatgccatta tctcttttcc agtattaaac actcatatgc ttatggcttn gagaaatttc    240 

ttagttgggt gaattaaagg ttaatccgag aattagcatg gatataccgg gtcctcatgc    300 

agcttggcag atatctgaga aatggtttaa ttcatgctca ggagctgtgt gccttttcca    360 

tcccttccgg gtcccttacc cctnacttt                                      389 

 
           
             3  
             465  
             DNA  
             Homo sapiens  
             
               unsure  
               (384)  
             
             
               unsure  
               (411)  
             
             
               unsure  
               (445)  
             
           
            3 

tttttttttt tatcaacatt tatatgcttt attgaaagtt gacaagtgca acagttaaat     60 

acagtgacac cttacaattg tgtagagaac atgcacagaa acatatgcat ataactacta    120 

tacaggtgat atgcagaaac ccctactggg aaatccattt cattagttag aactgagcat    180 

ttttcaaagt attcaaccag actcaattga aagacttcag tgaacaagga tttacttcag    240 

cgtattcagg caggctagga tttcaggatt acacaaagtg aggtaactgt gccaaattct    300 

taaaatttct ttagggtgtg ggtttttgtc atgtagcagt ttttatgtgg atctattata    360 

taaaagtcca cacctcctca gacngccaat ggaaacaact taaatttcca ntctgttaca    420 

acctaattgg taggttacag tcccnttttg ttacaaatgg ttaca                    465 

 
           
             4  
             1718  
             DNA  
             Homo sapiens  
           
            4 

ggcacgaggg gatccgcatc tgcctgggat catcaagccc tagaagctgg gtttctttaa     60 

attagggctg ccgttttctg tttctccctg ggctgcggaa agccagaaga ttttatctag    120 

cttatacaag gctgctggtg ttccctcttt ttttccacga gggtgttttt ggctgcaatt    180 

gcatgaaatc ccaatggtgt agaccagtgg cgatggatct aggagtttac caactgagac    240 

atttttcaat ttctttcttg tcatccttgc tggggactga aaacgcttct gtgagacttg    300 

ataatagctc ctctggtgca agtgtggtag ctattgacaa caaaatcgag caagctatgg    360 

atctagtgaa aagccatttg atgtatgcgg tcagagaaga agtggaggtc ctcaaagagc    420 

aaatcaaaga actaatagag aaaaattccc agctggagca ggagaacaat ctgctgaaga    480 

cactggccag tcctgagcag cttgcccagt ttcaggccca gctgcagact ggctcccccc    540 

ctgccaccac ccagccacag ggcaccacac agccccccgc ccagccagca tcgcagggct    600 

caggaccaac cgcatagctg cctatgcccc cgcagaactg gctgctgcgt gtgaactgaa    660 

cagacggaga agatgtgcta gggagaatct gcctccacag tcacccattt cattgctcgc    720 

tgcgaaagag acgtgagact gacatatgcc attatctctt ttccagtatt aaacactcat    780 

atgcttatgg cttggagaaa tttcttagtt gggtgaatta aaggttaatc cgagaattag    840 

catggatata ccgggacctc atgcagcttg gcagatatct gagaaatggt ttaattcatg    900 

ctcaggagct gtgtgccttt ccatcccttc cggctcccta cccctcactt ccaagggttc    960 

tctctcctgc ttgcgcttag tgtcctacat ggggttgtga agcgatggag ctcctcactg   1020 

gactcgcctc tctcctctcc tccccccagg aggaacttga aaggagggta aaaagactaa   1080 

aatgaggggg aacagagttc actgtacaaa tttgacaact gtcaccaaaa ttcataaaaa   1140 

acaatagtac tgtgcctctt tcttctcaaa caatggatga cacaaaacta tgagagtgac   1200 

aaaatggtga caggtagctg ggacctaggc tatcttacca tgaaggttgt tttgcttatt   1260 

gtatatttgt gtatgtagtg taactatttt gtacaataga ggactgtaac tactatttag   1320 

gttgtacaga ttgaaattta gttgtttcat tggctgtctg aggaggtgtg gacttttata   1380 

tatagatcta cataaaaact gctacatgac aaaaaccaca cctaaagaaa ttttaagaat   1440 

ttggcacagt tactcacttt gtgtaatctg aaatctagct gctgaatacg ctgaagtaaa   1500 

tccttgttca ctgaagtctt tcaattgagc tggttgaata ctttgaaaaa tgctcagttc   1560 

taactaatga aatggatttc ccagtagggg tttctgcata tcacctgtat agtagttata   1620 

tgcatatgtt tctgtgcatg ttctctacac aattgtaagg tgtcactgta tttaactgtt   1680 

gcacttgtca actttcaata aagcatataa atgttgat                           1718 

 
           
             5  
             392  
             DNA  
             Homo sapiens  
             
               unsure  
               (342)  
             
             
               unsure  
               (362)  
             
           
            5 

gctcctacca cccagacacc caaacagccg tggccccaga ggtcctggcc aaatatgggg     60 

gcctgcctag gttggtggaa cagtgctcct tatgtaaact gagccctttg tttagaaaac    120 

aattccaaat gtgaaactag aatgagaggg aagagatagc atggcatgca gcacacacgg    180 

ctgctccagt tcatggcctc ccaggggtgc tggggatgca tccaaagtgg ttgtctgaga    240 

cagagttggg aaaccctcac caactgggcc tctttcacct tccacattat cccgctgcca    300 

ccggttgccc tgttttcatt gcaggtttca gggaccagct tngggttgcg tgcgtttttg    360 

cntttgccag ttcaggccga gggtgttagt tt                                  392 

 
           
             6  
             429  
             DNA  
             Homo sapiens  
           
            6 

ttttttttta aggacacgag agagccatat ttatttcaca tggacaagca tgattccatt     60 

gcatgctgaa catgaaagct cgtatgagca aagtacccgt aacagcagaa ttatgtgctt    120 

ttgtccacag ggagcaggga gaatcacaaa gttgttttca gagacagtgt ttttcaagca    180 

cagttgagac cataggctct ggaagtcact ggtttatttc atcaccaaag ggtctgtctc    240 

ccagggagtg gccggagtgc tttcagcttt gcaatctctc aatgaattga taaggtctga    300 

ggagggctga ggatggtctc ccatcccacc acccagagca tctttgaagg aaatgaagct    360 

cagaggggaa ggttacatgc cattgggaat ttaacaaggg ccattcctgg gttggacaat    420 

gacagggga                                                            429 

 
           
             7  
             1305  
             DNA  
             Homo sapiens  
           
            7 

cgcggctcag taattgaagg cctgaaacgc ccatgtgcca ctgactagga ggcttccctg     60 

ctgcggcact tcatgaccca gcggcgcgcg gcccagtgaa gccaccgtgg tgtccagcat    120 

ggccgcgctg ctcctgggcg cggtgctgct ggtggcccag ccccagctag tgccttcccg    180 

ccccgccgag ctaggccagc aggagcttct gcggaaagcg gggaccctcc aggatgacgt    240 

ccgcgatggc gtggccccaa acggctctgc ccagcagttg ccgcagacca tcatcatcgg    300 

cgtgcgcaag ggcggcacgc gcgcactgct ggagatgctc agcctgcacc ccgacgtggc    360 

ggccgcggag aacgaggtcc acttcttcga ctgggaggag cattacagcc acggcttggg    420 

ctggtacctc agccagatgc ccttctcctg gccacaccag ctcacagtgg agaagacccc    480 

cgcgtatttc acgtcgccca aagtgcctga gcgagtctac agcatgaacc cgtccatccg    540 

gctgctgctc atcctgcgag acccgtcgga gcgcgtgcta tctgactaca cccaagtgtt    600 

ctacaaccac atgcagaagc acaagcccta cccgtccatc gaggagttcc tggtgcgcga    660 

tggcaggctc aatgtggact acaaggccct caaccgcagc ctctaccacg tgcacatgca    720 

gaactggctg cgctttttcc cgctgcgcca catccacatt gtggacggcg accgcctcat    780 

cagggacccc ttccctgaga tccaaaaggt cgagaggttc ctaaagctgt cgccgcagat    840 

caatgcttcg aacttctact ttaacaaaac caagggcttt tactgcctgc gggacagcgg    900 

ccgggaccgc tgcttacatg agtccaaagg ccgggcgcac ccccaagtcg atcccaaact    960 

actcaataaa ctgcacgaat attttcatga gccaaataag aagttcttcg agcttgttgg   1020 

cagaacattt gactggcact gatttgcaat aagctaagct cagaaacttt cctactgtaa   1080 

gttctggtgt acatctgagg ggaaaaagaa ttttaaaaaa gcatttaagg tataatttat   1140 

ttgtaaaatc cataaagtac ttctgtacag tattagattc acaattgcca tatatactag   1200 

ttatattttt ctacttgtta aatggagggc attttgtatt gtttttcatg gttgttaaca   1260 

ttgtgtaata tgtctctata tgaaggaact aaactatttc actga                   1305 

 
           
             8  
             331  
             DNA  
             Homo sapiens  
             
               unsure  
               (80)  
             
             
               unsure  
               (104)  
             
             
               unsure  
               (115)  
             
             
               unsure  
               (135)  
             
             
               unsure  
               (186)  
             
             
               unsure  
               (308)  
             
             
               unsure  
               (271)  
             
           
            8 

gctcaggaca gatgccacac aaggatagat gctggcccag ggccaagagc ccagctccaa     60 

ggggaatcag aactcaaatn gggccagatc cagcctgggg tctngagttg atctngaacc    120 

cagactcaga cattngcacc taatccaggc agatccagga ctatatttgg gcctgctcca    180 

gacctngatc ctggaggccc agttcaccct gatttaggag aagccaggaa tttcccagga    240 

ccctgaaggg gccatgatgg caacagatct ngaacctcag cctggccaga cacaggccct    300 

ccctgttncc cagagaaagg ggagcccact g                                   331 

 
           
             9  
             346  
             DNA  
             Homo sapiens  
             
               unsure  
               (40)  
             
             
               unsure  
               (286)  
             
             
               unsure  
               (320)  
             
           
            9 

tttattgcac ttgcaacaga gtttaaataa gtcctgggtn tctggtgcca aggtgaggga     60 

agggttgggc agagagatga ggggcagcat cagtgcagct ggcaggcaga acccaaattc    120 

tgcaggccca ggacagtggg ctcccctttc tctggggaac agggagggcc tgtgtctggc    180 

caggctgagg ttccagatct gttgccatca tggccccttc agggtcctgg ggaaattcct    240 

gggcttctcc taaatcaggg tgaactgggc ctccagggat caggtntggg agcaggccca    300 

aatataagtc ctgggatctn cctgggatta gggtgccaat gtctga                   346 

 
           
             10  
             4132  
             DNA  
             Homo sapiens  
           
            10 

cgctggggcc cccggcgccg acccccgctg ctgccgctgc tgttgctgct gctgccgccg     60 

ccacccaggg tcgggggctt caacttagac gcggaggccc cagcagtact ctcggggccc    120 

ccgggctcct tcttcggatt ctcagtggag ttttaccggc cgggaacaga cggggtcagt    180 

gtgctggtgg gagcacccaa ggctaatacc agccagccag gagtgctgca gggtggtgct    240 

gtctacctct gtccttgggg tgccagcccc acacagtgca cccccattga atttgacagc    300 

aaaggctctc ggctcctgga gtcctcactg tccagctcag agggagagga gcctgtggag    360 

tacaagtcct tgcagtggtt cggggcaaca gttcgagccc atggctcctc catcttggca    420 

tgcgctccac tgtacagctg gcgcacagag aaggagccac tgagcgaccc cgtgggcacc    480 

tgctacctct ccacagataa cttcacccga attctggagt atgcaccctg ccgctcagat    540 

ttcagctggg cagcaggaca gggttactgc caaggaggct tcagtgccga gttcaccaag    600 

actggccgtg tggttttagg tggaccagga agctatttct ggcaaggcca gatcctgtct    660 

gccactcagg agcagattgc agaatcttat taccccgagt acctgatcaa cctggttcag    720 

gggcagctgc agactcgcca ggccagttcc atctatgatg acagctacct aggatactct    780 

gtggctgttg gtgaattcag tggtgatgac acagaagact ttgttgctgg tgtgcccaaa    840 

gggaacctca cttacggcta tgtcaccatc cttaatggct cagacattcg atccctctac    900 

aacttctcag gggaacagat ggcctcctac tttggctatg cagtggccgc cacagacgtc    960 

aatggggacg ggctggatga cttgctggtg ggggcacccc tgctcatgga tcggacccct   1020 

gacgggcggc ctcaggaggt gggcagggtc tacgtctacc tgcagcaccc agccggcata   1080 

gagcccacgc ccacccttac cctcactggc catgatgagt ttggccgatt tggcagctcc   1140 

ttgacccccc tgggggacct ggaccaggat ggctacaatg atgtggccat cggggctccc   1200 

tttggtgggg agacccagca gggagtagtg tttgtatttc ctgggggccc aggagggctg   1260 

ggctctaagc cttcccaggt tctgcagccc ctgtgggcag ccagccacac cccagacttc   1320 

tttggctctg cccttcgagg aggccgagac ctggatggca atggatatcc tgatctgatt   1380 

gtggggtcct ttggtgtgga caaggctgtg gtatacaggg gccgccccat cgtgtccgct   1440 

agtgcctccc tcaccatctt ccccgccatg ttcaacccag aggagcggag ctgcagctta   1500 

gaggggaacc ctgtggcctg catcaacctt agcttctgcc tcaatgcttc tggaaaacac   1560 

gttgctgact ccattggttt cacagtggaa cttcagctgg actggcagaa gcagaaggga   1620 

ggggtacggc gggcactgtt cctggcctcc acgcaggcaa ccctgaccca gaccctgctc   1680 

atccagaatg gggctcgaga ggattgcaga gagatgaaga tctacctcag gaacgagtca   1740 

gaatttcgag acaaactctc gccgattcac atcgctctca acttctcctt ggacccccaa   1800 

gccccagtgg acagccacgg cctcaggcca gccctacatt atcagagcaa gagccggata   1860 

gaggacaagg ctcagatctt gctggactgt ggagaagaca acatctgtgt gcctgacctg   1920 

cagctggaag tgtttgggga gcagaaccat gtgtacctgg gtgacaagaa tgccctgaac   1980 

ctcactttcc atgcccagaa tgtgggtgag ggtggcgcct atgaggctga gcttcgggtc   2040 

accgcccctc cagaggctga gtactcagga ctcgtcagac acccagggaa cttctccagc   2100 

ctgagctgtg actactttgc cgtgaaccag agccgcctgc tggtgtgtga cctgggcaac   2160 

cccatgaagg caggagccag tctgtggggt ggccttcggt ttacagtccc tcatctccgg   2220 

gacactaaga aaaccatcca gtttgacttc cagatcctca gcaagaatct caacaactcg   2280 

caaagcgacg tggtttcctt tcggctctcc gtggaggctc aggcccaggt caccctgaac   2340 

ggtgtctcca agcctgaggc agtgctattc ccagtaagcg actggcatcc ccgagaccag   2400 

cctcagaagg aggaggacct gggacctgct gtccaccatg tctatgagct catcaaccaa   2460 

ggccccagct ccattagcca gggtgtgctg gaactcagct gtccccaggc tctggaaggt   2520 

cagcagctcc tatatgtgac cagagttacg ggactcaact gcaccaccaa tcaccccatt   2580 

aacccaaagg gcctggagtt ggatcccgag ggttccctgc accaccagca aaaacgggaa   2640 

gctccaagcc gcagctctgc ttcctcggga cctcagatcc tgaaatgccc ggaggctgag   2700 

tgtttcaggc tgcgctgtga gctcgggccc ctgcaccaac aagagagcca aagtctgcag   2760 

ttgcatttcc gagtctgggc caagactttc ttgcagcggg agcaccagcc atttagcctg   2820 

cagtgtgagg ctgtgtacaa agccctgaag atgccctacc gaatcctgcc tcggcagctg   2880 

ccccaaaaag agcgtcaggt ggccacagct gtgcaatgga ccaaggcaga aggcagctat   2940 

ggcgtcccac tgtggatcat catcctagcc atcctgtttg gcctcctgct cctaggtcta   3000 

ctcatctaca tcctctacaa gcttggattc ttcaaacgct ccctcccata tggcaccgcc   3060 

atggaaaaag ctcagctcaa gcctccagcc acctctgatg cctgagtcct cccaatttca   3120 

gactcccatt cctgaagaac cagtcccccc accctcattc tactgaaaag gaggggtctg   3180 

ggtacttctt gaaggtgctg acggccaggg agaagctcct ctccccagcc cagagacata   3240 

cttgaagggc cagagccagg ggggtgagga gctggggatc cctccccccc atgcactgtg   3300 

aaggaccctt gtttacacat accctcttca tggatggggg aactcagatc cagggacaga   3360 

ggcccagcct ccctgaagcc tttgcatttt ggagagtttc ctgaaacaac ttggaaagat   3420 

aactaggaaa tccattcaca gttctttggg ccagacatgc cacaaggact tcctgtccag   3480 

ctccaacctg caaagatctg tcctcagcct tgccagagat ccaaaagaag cccccagcta   3540 

agaacctgga acttggggag ttaagacctg gcagctctgg acagccccac cctggtgggc   3600 

caacaaagaa cactaactat gcatggtgcc ccaggaccag ctcaggacag atgccacaca   3660 

aggatagatg ctggcccagg gccagagccc agctccaagg ggaatcagaa ctcaaatggg   3720 

gccagatcca gcctggggtc tggagttgat ctggaaccca gactcagaca ttggcaccta   3780 

atccaggcag atccaggact atatttgggc ctgctccaga cctgatcctg gaggcccagt   3840 

tcaccctgat ttaggagaag ccaggaattt cccaggacct gaaggggcca tgatggcaac   3900 

agatctggaa cctcagcctg gccagacaca ggccctccct gttccccaga gaaaggggag   3960 

cccactgtcc tgggcctgca gaatttccct tctgcctgcc agctgcactg atgctgcccc   4020 

tcatctctct gcccaaccct tccctcacct tggcaccaga cacccaggac ttatttaaac   4080 

tctgttgcaa gtgcaataaa tctgacccag tgcccccact gaccagaact ag           4132 

 
           
             11  
             486  
             DNA  
             Homo sapiens  
             
               unsure  
               (376)  
             
             
               unsure  
               (395)  
             
             
               unsure  
               (467)  
             
           
            11 

agcctgatct ctgtccaccg gtcctttata ccctcatgac ccgctgctgg gactacgacc     60 

ccagtgaccg gccccgcttc accgagctgg tgtgcagcct cagtgacgtt tatcagatgg    120 

agaaggacat tgccatggag caagagagga atgctcgcta ccgaaccccc aaaatcttgg    180 

agcccacagc cttccaggaa cccccaccca agcccagccg acctaagtac agaccccctc    240 

cgcaaaccaa cctcctgggc tccaaagctg cagttccagg ttcctgaggg tctgtgtgcc    300 

agctctcctg acggcttcac cagccctatg ggagtattcc attcttcccg ttaaattcac    360 

tggcacaccc cacctnttcc accgggcaca atgtntttca aaacggccac aggatggggg    420 

ggagggaggg attttcattc caacccaggc aggccgagga agagggncca gcagttgttg    480 

gggagg                                                               486 

 
           
             12  
             393  
             DNA  
             Homo sapiens  
             
               unsure  
               (349)  
             
             
               unsure  
               (360)  
             
             
               unsure  
               (377)  
             
             
               unsure  
               (384)  
             
           
            12 

tttttttttt ttttgcaaat gggacaattt taattcaacc acaagtcaaa tagaaagaag     60 

ttaaaagaat gtttatgcaa acacatgaga aaagaagggt gcagatgaga atgggggttg    120 

gggagagaaa gaggaggagt aagaaaagag ggaaaagcaa gggaaagtaa aggaagaaag    180 

agaaagaggg gcaggaagag agcggatttg gcccaaggtc ctatcttggc cgcatctctc    240 

tgcttcttcc ccctgatgct tggtttgttg acaacacagc atcctgtgcc tgggactccc    300 

aattagcttg ttcctgggac tgtgccccag ggtcctccct caggagggnc acatgctgtn    360 

cagtccagac caaactncac attnaaataa ttt                                 393 

 
           
             13  
             4089  
             DNA  
             Homo sapiens  
           
            13 

gaattccgtc agccctttta ctcagccaca gcctccggag ccgttgcaca cctacctgcc     60 

cggccgactt acctgtactt gccgccgtcc cggctcacct ggcggtgccc gaggagtagt    120 

cgctggagtc cgcgcctccc tgggactgca atgtgccgat cttagctgct gcctgagagg    180 

atgtctgggg tgtccgagcc cctgagtcga gtaaagttgg gcacgttacg ccggcctgaa    240 

ggccctgcag agcccatggt ggtggtacca gtagatgtgg aaaaggagga cgtgcgtatc    300 

ctcaaggtct gcttctatag caacagcttc aatcctggga aaaacttcaa actggtcaaa    360 

tgcactgtcc agacggagat ccgggagatc atcacctcca tcctgctgag cgggcggatc    420 

gggcccaaca tccggttggc tgagtgctat gggctgaggc tgaagcacat gaagtccgat    480 

gagatccact ggctgcaccc acagatgacg gtgggtgagg tgcaggacaa gtatgagtgt    540 

ctgcacgtgg aagccgagtg gaggtatgac cttcaaatcc gctacttgcc agaagacttc    600 

atggagagcc tgaaggagga caggaccacg ctgctctatt tttaccaaca gctccggaac    660 

gactacatgc agcgctacgc cagcaaggtc agcgagggca tggccctgca gctgggctgc    720 

ctggagctca ggcggttctt caaggatatg ccccacaatg cacttgacaa gaagtccaac    780 

ttcgagctcc tagaaaagga agtggggctg gacttgtttt tcccaaagca gatgcaggag    840 

aacttaaagc ccaaacagtt ccggaagatg atccagcaga ccttccagca gtacgcctcg    900 

ctcagggagg aggagtgcgt catgaagttc ttcaacactc tcgccccgtt cgccaacatc    960 

gaccaggaga cctaccgctg tgaactcatt caaggatgga acattactgt ggacctggtc   1020 

attggcccta aagggatccg ccagctgact agtcaggacg caaagcccac ctgcctggcc   1080 

gagttcaagc agatcaggtc catcaggtgc ctcccgctgg aggagggcca ggcagtactt   1140 

cagctgggca ttgaaggtgc cccccaggcc ttgtccatca aaacctcatc cctagcagag   1200 

gctgagaaca tggctgacct catagacggc tactgccggc tgcagggtga gcaccaaggc   1260 

tctctcatca tccatcctag gaaagatggt gagaagcgga acagcctgcc ccagatcccc   1320 

atgctaaacc tggaggcccg gcggtcccac ctctcagaga gctgcagcat agagtcagac   1380 

atctacgcag agattcccga cgaaaccctg cgaaggcccg gaggtccaca gtatggcatt   1440 

gcccgtgaag atgtggtcct gaatcgtatt cttggggaag gcttttttgg ggaggtctat   1500 

gaaggtgtct acacaaatca taaaggggag aaaatcaatg tagctgtcaa gacctgcaag   1560 

aaagactgca ctctggacaa caaggagaag ttcatgagcg aggcagtgat catgaagaac   1620 

ctcgaccacc cgcacatcgt gaagctgatc ggcatcattg aagaggagcc cacctggatc   1680 

atcatggaat tgtatcccta tggggagctg ggccactacc tggagcggaa caagaactcc   1740 

ctgaaggtgc tcaccctcgt gctgtactca ctgcagatat gcaaagccat ggcctacctg   1800 

gagagcatca actgcgtgca cagggacatt gctgtccgga acatcctggt ggcctcccct   1860 

gagtgtgtga agctggggga ctttggtctt tcccggtaca ttgaggacga ggactattac   1920 

aaagcctctg tgactcgtct ccccatcaaa tggatgtccc cagagtccat taacttccga   1980 

cgcttcacga cagccagtga cgtctggatg ttcgccgtgt gcatgtggga gatcctgagc   2040 

tttgggaagc agcccttctt ctggctggag aacaaggatg tcatcggggt gctggagaaa   2100 

ggagaccggc tgcccaagcc tgatctctgt ccaccggtcc tttataccct catgacccgc   2160 

tgctgggact acgaccccag tgaccggccc cgcttcaccg agctggtgtg cagcctcagt   2220 

gacgtttatc agatggagaa ggacattgcc atggagcaag agaggaatgc tcgctaccga   2280 

acccccaaaa tcttggagcc cacagccttc caggaacccc cacccaagcc cagccgacct   2340 

aagtacagac cccctccgca aaccaacctc ctggctccaa agctgcagtt ccaggttcct   2400 

gagggtctgt gtgccagctc tcctacgctc accagcccta tggagtatcc atctcccgtt   2460 

aactcactgc acaccccacc tctccaccgg cacaatgtct tcaaacgcca cagcatgggg   2520 

gaggaggact tcatccaacc cagcagccga gaagaggccc agcagctgtg ggaggctgaa   2580 

aaggtcaaaa tgcggcaaat cctggacaaa cagcagaagc agatggtgga ggactaccag   2640 

tggctcaggc aggaggagaa gtccctggac cccatggttt atatgaatga taagtcccca   2700 

ttgacgccag agaaggaggt cggctacctg gagttcacag ggcccccaca gaagcccccg   2760 

aggctgggcg cacagtccat ccagcccaca gctaacctgg accggaccga tgacctggtg   2820 

tacctcaatg tcatggagct ggtgcgggcc gtgctggagc tcaagaatga gctctgtcag   2880 

ctgccccccg agggctacgt ggtggtggtg aagaatgtgg ggctgaccct gcggaagctc   2940 

atcgggagcg tggatgatct cctgccttcc ttgccgtcat cttcacggac agagatcgag   3000 

ggcacccaga aactgctcaa caaagacctg gcagagctca tcaacaagat gcggctggcg   3060 

cagcagaacg ccgtgacctc cctgagtgag gagtgcaaga ggcagatgct gacggcttca   3120 

cacaccctgg ctgtggacgc caagaacctg ctcgacgctg tggaccaggc caaggttctg   3180 

gccaatctgg cccacccacc tgcagagtga cggagggtgg gggccacctg cctgcgtctt   3240 

ccgcccctgc ctgccatgta cctcccctgc cttgctgttg gtcatgtggg tcttccaggg   3300 

agaaggccaa ggggagtcac cttcccttgc cactttgcac gacgccctct ccccacccct   3360 

acccctggct gtactgctca ggctgcagct ggacagaggg gactctgggc tatggacaca   3420 

gggtgacggt gacaaagatg gctcagaggg ggactgctgc tgcctggcca ctgctcccta   3480 

agccagcctg gtccatgcag ggggctcctg ggggtgggga ggtgtcacat ggtgccccta   3540 

gctttatata tggacatggc aggccgattt gggaaccaag ctattccttt cccttcctct   3600 

tctcccctca gatgtccctt gatgcacaga gaagctgggg aggagctttg ttttcggggg   3660 

tcaggcagcc agtgagatga gggatgggcc tggcattctt gtacagtgta tattgaaatt   3720 

tatttaatgt gaggtttggt ctggactgac agcatgtgcc ctcctgaggg aggaccaggg   3780 

cacagtccag gaacaagcta attgggagtc caggcacagg atgctgtgtt gtcaacaaac   3840 

caagcatcag ggggaagaag cagagagatg cggccaagat aggaccttgg gccaaatccg   3900 

ctctcttcct gcccctcttt ctctttcttc ctttactttc ccttgctttt ccctcttttc   3960 

ttactcctcc tctttctctc ccccaccccc attctcatct gcacccttct tttctcatgt   4020 

gtttgcataa acattctttt aacttctttc tatttgactt gtggttgaat taaaattgtc   4080 

ccatttgca                                                           4089 

 
           
             14  
             464  
             DNA  
             Homo sapiens  
             
               unsure  
               (146)  
             
             
               unsure  
               (448)  
             
           
            14 

gacctggaga tcaacgggga gaaggtgaag ctgcagatct gggacacagc ggggcaggag     60 

cgcttccgca ccatcacctc cacgtattat cgggggaccc acggggtcat ttgtggttta    120 

cgacgtcacc agtgccgagt cctttntcaa cgtcaagcgg tggcttcacg aaatcaacca    180 

gaactgtgat gatgtgtgcc gaatattagt gggtaataag aatgacgacc ctgagcggaa    240 

ggtggtggag acggaagatg cctacaaatt cgccgggcag atgggcatcc agttgttcga    300 

gaccagcgcc aaggagaatg tcaacgtggg aagagatgtt tcaactgcat tcacggagct    360 

ggtcctccga gcaaagaaag acaaccttgg gcaaaacagc agcagcaaca acagaacgat    420 

gttggttgaa gtttacgaag gaacattnaa cgaaagaaac gttt                     464 

 
           
             15  
             373  
             DNA  
             Homo sapiens  
           
            15 

tttttttttt tttttttttt taattgtgag gaatttaatt cacttgattt ggcttcattt     60 

tcttgatctg ttaaaataat cctcccatag cccccctgcc agccccatct ctgcacgaac    120 

ctaccccgac ctttctgttg gaactgaaac ctgttggtgt aaatgagaag ccatggctgc    180 

cctgggtttg gagctcagag gcatctagaa ggcaggacaa gaaatctgtt ggccaaaggg    240 

caagacctgc cacctctgtg gaactgcagg gcctgccttg agaccaggtt ccccagctcc    300 

cagaatggct gtggggacag gacaacgggg agggaaggga gctggcacag gccccggaga    360 

aggggcaaga ccc                                                       373 

 
           
             16  
             730  
             DNA  
             Homo sapiens  
           
            16 

gctgccggag cagcccgaag agctgcggat cgcgaggcca gtaccgaccc cgcccgcccg     60 

cgcgctccgc ccccgcccgc catggcccgg gactacgacc acctcttcaa gctgctcatc    120 

atcggcgaca gcggtgtggg caagagcagt ttactgttgc gttttgcaga caacactttc    180 

tcaggcagct acatcaccac gatcggagtg gatttcaaga tccggaccgt ggagatcaac    240 

ggggagaagg tgaagctgca gatctgggac acagcggggc aggagcgctt ccgcaccatc    300 

acctccacgt attatcgggg gacccacggg gtcattgtgg tttacgacgt caccagtgcc    360 

gagtcctttg tcaacgtcaa gcggtggctt cacgaaatca accagaactg tgatgatgtg    420 

tgccgaatat tagtgggtaa taagaatgac gaccctgagc ggaaggtggt ggagacggaa    480 

gatgcctaca aattcgccgg gcagatgggc atccagttgt tcgagaccag cgccaaggag    540 

aatgtcaacg tggaagagat gttcaactgc atcacggagc tggtcctccg agcaaagaaa    600 

gacaacctgg caaaacagca gcagcaacaa cagaacgatg tggtgaagct cacgaagaac    660 

agtaaacgaa agaaacgctg ctgctaatgg cacccagtcc actgcagaga ctgcactgcg    720 

gtccctcccc                                                           730 

 
           
             17  
             334  
             DNA  
             Homo sapiens  
             
               unsure  
               (61)  
             
             
               unsure  
               (223)  
             
             
               unsure  
               (230)  
             
             
               unsure  
               (304)  
             
           
            17 

acagagtagc agctcagatg ccagagatcg aaagaaggct cgaatgagtg agctggaaca     60 

naagtggtag atttagaaga agagaaccaa aaacttttgc tagaaaatca gcttttacga    120 

gagaaaactc atggccttgt agttgagaac caggagttaa gacagcgctt ggggatggat    180 

gccctggttg ctgaagagga ggcggagcaa ggggaatgaa gtnaggccan tgcgggtctg    240 

ctgagtccgc agcactcaga ctacgtgcac ctctgcagca ggtgcaggcc cagttgtcac    300 

cctncagaac atctccccat ggattctggc ggta                                334 

 
           
             18  
             412  
             DNA  
             Homo sapiens  
             
               unsure  
               (120)  
             
             
               unsure  
               (153)  
             
             
               unsure  
               (210)  
             
             
               unsure  
               (372)  
             
             
               unsure  
               (381)  
             
             
               unsure  
               (411)  
             
           
            18 

tttttttttg ctgcattgta ccttttaatt gcatgggtag ttttaaataa atggagaaag     60 

cacctttcag aagctacact agcaggaaaa aattccatca agcatttaca tagtaaattn    120 

ctataatttc acaaaagatt cttgatctta ctngaagtat acatgaggga aagagccccc    180 

tcagcaggtg ttcccgttgc ttacagaagn aaactaaagg acctaaaact ggaggcaagc    240 

cagggtgcca aaaaggggga agagaaatga taaagaacca ttcataaatt ccatgtctac    300 

ttcaaggaca tttgtctaat gacccttaca taataagtat tttaggggaa aactaccacc    360 

ctttttaagg tnaaagtaca nttcttaaaa ggctggtagg tttctcaatt nt            412 

 
           
             19  
             1818  
             DNA  
             Homo sapiens  
           
            19 

tagtctggag ctatggtggt ggtggcagcc gcgccgaacc cggccgacgg gacccctaaa     60 

gttctgcttc tgtcggggca gcccgcctcc gccgccggag ccccggcggc caggctgccg    120 

ctcatggtgc cagcccagag aggggccagc ccggaggcag cgagcggggg gctgccccag    180 

gcgcgcaagc gacagcgcct cacgcacctg agccccgagg agaaggcgct gaggaggaaa    240 

ctgaaaaaca gagtagcagc tcagactgcc agagatcgaa agaaggctcg aatgagtgag    300 

ctggaacagc aagtggtaga tttagaagaa gagaaccaaa aacttttgct agaaaatcag    360 

cttttacgag agaaaactca tggccttgta gttgagaacc aggagttaag acagcgcttg    420 

gggatggatg ccctggttgc tgaagaggag gcggaagcca aggggaatga agtgaggcca    480 

gtggccgggt ctgctgagtc cgcagcactc agactacgtg cacctctgca gcaggtgcag    540 

gcccagttgt cacccctcca gaacatctcc ccatggattc tggcggtatt gactcttcag    600 

attcagagtc tgatatcctg ttgggcattc tggacaactt ggacccagtc atgttcttca    660 

aatgcccttc cccagagcct gccagcctgg aggagctccc agaggtctac ccagaaggac    720 

ccagttcctt accagcctcc ctttctctgt cagtggggac gtcatcagcc aagctggaag    780 

ccattaatga actaattcgt tttgaccaca tatataccaa gcccctagtc ttagagatac    840 

cctctgagac agagagccaa gctaatgtgg tagtgaaaat cgaggaagca cctctcagcc    900 

cctcagagaa tgatcaccct gaattcattg tctcagtgaa ggaagaacct gtagaagatg    960 

acctcgttcc ggagctgggt atctcaaatc tgctttcatc cagccactgc ccaaagccat   1020 

cttcctgcct actggatgct acagtgactg tggatacggg ggttcccttt ccccattcag   1080 

tgacatgtcc tctctgcttg gtgtaaacat tcttgggagg acacttttgc caatgaactc   1140 

tttccccagc tgattagtgt ctaaggaatg atccaatact gttgcccttt tccttgacta   1200 

ttacactgcc tggaggatag cagagaagcc tgtctgtact tcattcaaaa agccaaaata   1260 

gagagtatac agtcctagag aatccctcta tttgttcaga tctcatagat gacccccagg   1320 

tattgccttt tgacatccag cagtccaagg tattgagaca tattactgga agtaagaaat   1380 

attactataa ttgagaacta cagcttttaa gattgtactt ttaagattgt acttttatct   1440 

taaaagggtg gtagttttcc ctaaaatact tattatgtaa gggtcattag acaaatgtct   1500 

tgaagtagac atggaattta tgaatggtct ttatcatttc tcttccccct ttttggcatc   1560 

ctggcttgcc tccagtttta ggtcctttag tttgcttctg caagcaacgg gaacacctgc   1620 

tgagggggct ctttccctca tgtatacttc aagtaagatc aagaatcttt tgtgaaatta   1680 

tagaaattta ctatgtaaat gcttgatgga attttttcct gctagtgtag cttctgaaag   1740 

gtgctttctc catttattta aaaactaccc atgcaattaa aaggtacaat gcaaaaaaaa   1800 

aaaaaaaaaa attttttt                                                 1818 

 
           
             20  
             350  
             DNA  
             Homo sapiens  
             
               unsure  
               (68)  
             
             
               unsure  
               (86)  
             
             
               unsure  
               (188)  
             
             
               unsure  
               (253)  
             
           
            20 

aaacagtaat tctttagact ttattaaaaa atgacataaa gtgcatctta ttaaaaaatg     60 

tataaaancc acataaattc cagggncccc tgtgcctggg cagtgttgat atcccttaga    120 

gtggaggaag gtgagggatg gagggtgaac tggggactgg ggagaggacc agggtgcagt    180 

tagttccncg tgtttgagtt caaagatgga gcgagggtgg atatggtggg aaggggcaca    240 

cgggttctca cgncaacaac ggaggaaggc aggcgacagt ctcttccctg aattctgagg    300 

gaaaggcgta cattgtcacg aaatctctcc tgagctcgcg ctgtcctctc               350 

 
           
             21  
             394  
             DNA  
             Homo sapiens  
             
               unsure  
               (208)  
             
             
               unsure  
               (345)  
             
             
               unsure  
               (361)  
             
             
               unsure  
               (373)  
             
             
               unsure  
               (378)  
             
           
            21 

gaaggaactg gtctgctcac acttgctggc ttgcgcatca ggactggctt tatctcctga     60 

ctcacggtgc aaaggtgcac tctgcgaacg ttaagtccgt ccccagcgct tggaatccta    120 

cggcccccac agccggatcc cctcagcctt ccaggtcctc aactcccgtg gacgctgaac    180 

aatggcctcc atggggctac aggtaatngg catcgcgctg gccgtcctgg gctggctggc    240 

cgtcatgctg tgctgcgcgc tgcccatgtg gcgcgtgacg gcctttcatc ggcagcaaca    300 

ttgtcaactt gcagaccatc tgggaagggc ctattggatg aactncgtgg ttcaaaagcc    360 

ngtccaagat tgnatttnaa aggttttaac gatt                                394 

 
           
             22  
             1665  
             DNA  
             Homo sapiens  
           
            22 

gaaggaactg gttctgctca cacttgctgg cttgcgcatc aggactggct ttatctcctg     60 

actcacggtg caaaggtgca ctctgcgaac gttaagtccg tccccagcgc ttggaatcct    120 

acggccccca cagccggatc ccctcagcct tccaggtcct caactcccgt ggacgctgaa    180 

caatggcctc catggggcta caggtaatgg gcatcgcgct ggccgtcctg ggctggctgg    240 

ccgtcatgct gtgctgcgcg ctgcccatgt ggcgcgtgac ggccttcatc ggcagcaaca    300 

ttgtcacctc gcagaccatc tgggagggcc tatggatgaa ctgcgtggtg cagagcaccg    360 

gccagatgca gtgcaaggtg tacgactcgc tgctggcact gccgcaggac ctgcaggcgg    420 

cccgcgccct cgtcatcatc agcatcatcg tggctgctct gggcgtgctg ctgtccgtgg    480 

tggggggcaa gtgtaccaac tgcctggagg atgaaagcgc caaggccaag accatgatcg    540 

tggcgggcgt ggtgttcctg ttggccggcc ttatggtgat agtgccggtg tcctggacgg    600 

cccacaacat catccaagac ttctacaatc cgctggtggc ctccgggcag aagcgggaga    660 

tgggtgcctc gctctacgtc ggctgggccg cctccggcct gctgctcctt ggcggggggc    720 

tgctttgctg caactgtcca ccccgcacag acaagcctta ctccgccaag tattctgctg    780 

cccgctctgc tgctgccagc aactacgtgt aaggtgccac ggctccactc tgttcctctc    840 

tgctttgttc ttccctggac tgagctcagc gcaggctgtg accccaggag ggccctgcca    900 

cgggccactg gctgctgggg actggggact gggcagagac tgagccaggc aggaaggcag    960 

cagccttcag cctctctggc ccactcggac aacttcccaa ggccgcctcc tgctagcaag   1020 

aacagagtcc accctcctct ggatattggg gagggacgga agtgacaggg tgtggtggtg   1080 

gagtggggag ctggcttctg ctggccagga tagcttaacc ctgactttgg gatctgcctg   1140 

catcggcgtt ggccactgtc cccatttaca ttttccccac tctgtctgcc tgcatctcct   1200 

ctgttccggg taggccttga tatcacctct gggactgtgc cttgctcacc gaaacccgcg   1260 

cccaggagta tggctgaggc cttgcccacc cacctgcctg ggaagtgcag agtggatgga   1320 

cgggtttaga ggggaggggc gaaggtgctg taaacaggtt tgggcagtgg tgggggaggg   1380 

ggccagagag gcggctcagg ttgcccagct ctgtggcctc aggactctct gcctcacccg   1440 

cttcagccca gggcccctgg agactgatcc cctctgagtc ctctgcccct tccaaggaca   1500 

ctaatgagcc tgggagggtg gcagggagga ggggacagct tcacccttgg aagtcctggg   1560 

gtttttcctc ttccttcttt gtggtttctg ttttgtaatt taagaagagc tattcatcac   1620 

tgtaattatt attattttct acaataaatg ggacctgtgc acagg                   1665 

 
           
             23  
             345  
             DNA  
             Homo sapiens  
             
               unsure  
               (291)  
             
           
            23 

aggtcctact ggaaggagtt cctggtgatg tgcacgctct ttgtgctggc cgtgctgctc     60 

ccagttttat tcttgctcta ccggcaccgg aacagcatga aagtcttcct gaagcagggg    120 

gaatgtgcca gcgtgcaccc caagacctgc cctgtggtgc tgccccctga gacccgccca    180 

ctcaacggcc tagggcccct agcaccccgc tcgatcaccg agggtaccag tccctgtcag    240 

acagcccccc ggggttcccg agtcttcact gagtcagaga agaggccact nagcatccaa    300 

gacagcttcg tgggaggtat ccccagtgtg cccccggccc cgggg                    345 

 
           
             24  
             2433  
             DNA  
             Homo sapiens  
           
            24 

gaagaaaggc tgattagaaa atttgaagct gaaaacatct ccaactacac ggcccttctg     60 

ctgagccagg atggaaagac gctgtatgtg ggggcccgag aggccctctt tgcacttaac    120 

agcaacctca gcttcttgcc aggcggggag taccaagagc tactgtggag tgcagatgct    180 

gacaggaagc agcagtgcag cttcaagggc aaggacccaa agcgtgactg tcaaaactac    240 

atcaagatcc tcctgccact caacagcagc cacctgctca cctgtggcac ggccgccttc    300 

agccccctgt gtgcttacat tcacatagcg agctttactt tagcccaaga tgaggccggt    360 

aatgtcattc tggaggatgg caagggtcat tgtccctttg accccaactt caagtccacg    420 

gctctggtgg ttgatggtga gctgtacact ggaacagtca gtagcttcca gggaaacgac    480 

ccagccattt cccggagcca gagttcccgc cccaccaaga ctgagagctc cctcaactgg    540 

ctacaagacc ctgcctttgt ggcctcggct acgtcccccg agagcctggg cagccccata    600 

ggtgatgatg ataagatcta cttcttcttc agcgagacgg gccaggagtt tgagttcttt    660 

gagaacacca tcgtgtcccg agttgcccga gtctgtaagg gcgatgaggg tggagagcgg    720 

gtgttgcagc aacgctggac ctcctttctc aaggctcagc tcctgtgctc ccggcctgat    780 

gatggctttc cctttaacgt gctacaagat gtcttcaccc tgaaccccaa ccctcaggat    840 

tggcgcaaga ccctttctat cggggtcttt acctcccagt ggcacagagg gaccacagaa    900 

ggctctgcca tctgcgtctt caccatgaat gatgtgcaga aggcctttga cggcctgtac    960 

aagaaagtaa acagagagac acagcagtgg tataccgaga cccaccaggt gcccacaccg   1020 

cggccgggag cgtgcattac caacagtgcc cgggaacgga agatcaactc gtccctgcag   1080 

ctcccagacc gagtgctgaa cttcctcaag gatcacttct tgatggatgg gcaggtccgc   1140 

agtcgcctgc tgctgctgca gcccagagcc cgctaccagc gtgtggctgt gcaccgtgtg   1200 

cctggcctgc acagcactta tgatgtccta tttctgggca ctggtgatgg ccgcctgcac   1260 

aaagcagtga ccctgagctc cagagtccac atcattgagg agctgcagat cttccctcaa   1320 

ggacagcctg tgcagaacct gctcttggac agccatgggg gactgttgta tgcctcctcc   1380 

cattccgggg tggtgcaagt gcccgtagcc aactgcagcc tgtacccaac ctgtggagac   1440 

tgcctcctgg ctcgagaccc ctactgcgcc tggactggct ctgcctgcag gctcgctagc   1500 

ctctaccagc ctgatctggc ctccaggcca tggacccagg acattgaggg tgccagtgtc   1560 

aaggaactct gcaagaattc ctcatacaag gcccggtttc ttgtgccagg taagccatgt   1620 

aaacaagtcc agatccaacc aaacacagtg aacaccctgg cctgcccact cctctcaaac   1680 

ctggccactc ggctctgggt gcacaatgga gccccagtca atgcctctgc ctcctgccgc   1740 

gtgttaccca ccggggacct gctgctggtg ggcagccagc agggtttggg ggtgttccag   1800 

tgttggtcga tagaagaagg attccagcag cttgtggcca gctactgccc agaggtgatg   1860 

gaggaggggg taatggacca aaagaaccag cgtgatggta ccccagtcat tatcaacaca   1920 

tcacgagtga gtgcaccggc tggtggcagg gacagctggg gtgcggacaa gtcctactgg   1980 

aatgaattcc tggtgatgtg tactctgttt gtgtttgcta tggtgctttt gtttctgttc   2040 

tttctctacc gacatcggga tggcatgaaa ctcttcctaa agcagggcga gtgtgccagt   2100 

gtgcacccca agactcgccc tatagtgcta ccacctgaga cccgaccgct gaatggtgtc   2160 

ggccctccta gcaccccact tgaccaccga ggctaccagg ctctgtcgga tagctcccca   2220 

gggcccagag tcttcactga atcagagaag aggccactga gcatccagga cagctttgta   2280 

gaggtgtctc ccgtgtgtcc ccggccccga gttcgactgg gctctgagat ccgagactct   2340 

gtggtatgag agctgacttt agatgtggtc accctgacct cagggttgtg agtgtcagtg   2400 

gaagtcagct acctctgctc tcacagaaca cag                                2433 

 
           
             25  
             463  
             DNA  
             Homo sapiens  
             
               unsure  
               (368)  
             
             
               unsure  
               (402)  
             
             
               unsure  
               (458)  
             
           
            25 

gtttggcaaa aactcaagcg gctggaagga ggaagaggtt ctccagagtc ggaactgagg     60 

gttggaacta tacccgggac caaactcacg gaccactcga ggcctgcaaa ccttcctggg    120 

aggacaggca ggccagatgg ccgctccact ggggaatgct cccagctgtg ctgtggagag    180 

aagctgatgt tttggtgtat tgtcagccat cgtccttgga ctcggagact atggcctcgc    240 

tccccaccct cctcttggaa ttacaagccc tggggtttga agctgacttt atagctgcaa    300 

gtgtatctcc ttttatctgg tgcctcctca aacccagtct cagacactta aatgcagaca    360 

acaccttnct cctgcagaca cctgggactg agccaaggag gncttgggga aggcccttag    420 

ggggagcacc ctgatgggag aggacagagc aggggttnca gca                      463 

 
           
             26  
             331  
             DNA  
             Homo sapiens  
             
               unsure  
               (13)  
             
             
               unsure  
               (15)  
             
             
               unsure  
               (322)  
             
           
            26 

agaaaaagcc cantnttcac tttattggag gtctctgcct ccattcacag gagaaaggag     60 

ctgggagccc catcctaagg gtcccagcat cagcccactg gagggcctgg aacagtccag    120 

cactctgtgg gagaggagtg gggaggggaa tgttttagaa aaaatagatc tctatgtaca    180 

tctgacatat ttatatagca cataaattag ggagtgctct gacccctgcc cgtggagccc    240 

aagcactgag cagggaggtg aacgccagtc cagaaagaag gtgctgggag cccctgctct    300 

gtcctctcca tccacggtgc tncccctagg g                                   331 

 
           
             27  
             1907  
             DNA  
             Homo sapiens  
           
            27 

cggccagata cctcagcgct acctggcgga actggatttc tctcccgcct gccggcctgc     60 

ctgccacagc cggactccgc cactccggta gcctcatggc tgcaacctgt gagattagca    120 

acatttttag caactacttc agtgcgatgt acagctcgga ggactccacc ctggcctctg    180 

ttccccctgc tgccaccttt ggggccgatg acttggtact gaccctgagc aacccccaga    240 

tgtcattgga gggtacagag aaggccagct ggttggggga acagccccag ttctggtcga    300 

agacgcaggt tctggactgg atcagctacc aagtggagaa gaacaagtac gacgcaagcg    360 

ccattgactt ctcacgatgt gacatggatg gcgccaccct ctgcaattgt gcccttgagg    420 

agctgcgtct ggtctttggg cctctggggg accaactcca tgcccagctg cgagacctca    480 

cttccagctc ttctgatgag ctcagttgga tcattgagct gctggagaag gatggcatgg    540 

ccttccagga ggccctagac ccagggccct ttgaccaggg cagccccttt gcccaggagc    600 

tgctggacga cggtcagcaa gccagcccct accaccccgg cagctgtggc gcaggagccc    660 

cctcccctgg cagctctgac gtctccaccg cagggactgg tgcttctcgg agctcccact    720 

cctcagactc cggtggaagt gacgtggacc tggatcccac tgatggcaag ctcttcccca    780 

gcgatggttt tcgtgactgc aagaaggggg atcccaagca cgggaagcgg aaacgaggcc    840 

ggccccgaaa gctgagcaaa gagtactggg actgtctcga gggcaagaag agcaagcacg    900 

cgcccagagg cacccacctg tgggagttca tccgggacat cctcatccac ccggagctca    960 

acgagggcct catgaagtgg gagaatcggc atgaaggcgt cttcaagttc ctgcgctccg   1020 

aggctgtggc ccaactatgg ggccaaaaga aaaagaacag caacatgacc tacgagaagc   1080 

tgagccgggc catgaggtac tactacaaac gggagatcct ggaacgggtg gatggccggc   1140 

gactcgtcta caagtttggc aaaaactcaa gcggctggaa ggaggaagag gttctccaga   1200 

gtcggaactg agggttggaa ctatacccgg gaccaaactc acggaccact cgaggcctgc   1260 

aaaccttcct gggaggacag gcaggccaga tggcccctcc actggggaat gctcccagct   1320 

gtgctgtgga gagaagctga tgttttggtg tattgtcagc catcgtcctt ggactcggag   1380 

actatggcct cgcctcccca ccctcctctt ggaattacaa gccctggggt ttgaagctga   1440 

ctttatagct gcaagtgtat ctccttttat ctggtgcctc ctcaaaccca gtctcagaca   1500 

cttaaatgca gacaacacct tcttcctgca gacacttgga ctgagccaag gaggcttggg   1560 

aggccctagg gagcaccgtg atggagagga cagagcaggg gctccagcac ttctttctgg   1620 

actggcgttc acctccctgc tcagtgcttg ggctccacgg gcaggggtca gagcactccc   1680 

taatttatgt gctatataaa tatgtcagat gtacatagag atctattttt tctaaaacat   1740 

tcccctcccc actcctctcc cacagagtgc tggactgttc caggccctcc agtgggctga   1800 

tgctgggacc cttaggatgg ggctcccagc tcctttctcc tgtgaatgga ggcagagacc   1860 

tccaataaag tgccttctgg gctttttcta aaaaaaaaaa aaaaaaa                 1907 

 
           
             28  
             467  
             DNA  
             Homo sapiens  
             
               unsure  
               (428)  
             
             
               unsure  
               (462)  
             
           
            28 

agtactacaa gcatcattct ctcaaggaag ggttcagaac cttagataca actctgcagt     60 

ttccatacaa ggagccagaa cattcagctg gacagagggg taatagagca ggcaacagct    120 

tgttaagtcc aaaagtgctg ggcattgcat cgctcggtat gacttctgtg caagagatat    180 

gagagagttg tccttgttga aaggagatgt ggtgaagatt tacacaaaga tgagtgcaaa    240 

tggctggtgg agaggagaag taaatggcag ggtgggctgg tttccatcca catatgtggg    300 

aaggaggatg aataaattca aatcccgtgt tgcaccctgc accaaaattt tcagaggaag    360 

gggataatta ggaagcctgc acagcttcgt ggatttaact tgaagtgttt ttaaaaagct    420 

ggcttttntg ggctgtttca acatcctccc tccttaggcc cntccta                  467 

 
           
             29  
             453  
             DNA  
             Homo sapiens  
             
               unsure  
               (240)  
             
             
               unsure  
               (387)  
             
             
               unsure  
               (438)  
             
           
            29 

ttttttttcc caacatgtaa ctctctcagt cttgtcagaa cacaacttct gctatggagg     60 

aaatatttcc atcaggaaag ggccaagtta gtgtcttaac ttgactgcct tgaatgggga    120 

ctctggaccc caggaagaat gtatttaggc tcctcacaaa aaagagtgat ggctgggcaa    180 

aacaaatgta ctgcaagacc catcttccct ccagttaata cactcccagg gatgggnctg    240 

cagaggggga gactctgaga gaagctggag gcccacaaaa gtccactgac cctctttctg    300 

tcccagaaat gaataaagga cccagttgtg ctttccttcc aaaatcctca acaaagttgt    360 

ttgtgctcca aggaaaatgt gggggantta aaaaaatcat gttcccgggt catctttgtg    420 

tgtgttgcgg gggaggtngg tggggaggga aaa                                 453 

 
           
             30  
             4762  
             DNA  
             Homo sapiens  
           
            30 

cccgccccgg cccagccgcg tcccggagcc gtcgggcatg gagccgtgga agcagtgcgc     60 

gcagtggctc atccattgca aggtgctgcc caccaaccac cgggtgacct gggactcggc    120 

tcaggtgttc gaccttgcgc agaccctccg cgatggagtc ctgctctgcc agctgcttaa    180 

caacctccgg gcgcactcca tcaacctgaa ggagatcaac ctgaggccgc agatgtccca    240 

gtttctctgt ttgaagaaca taaggacatt tctcacggcc tgttgtgaga cgtttggaat    300 

gaggaaaagt gaacttttcg aggcatttga cttgtttgat gttcgtgact ttggagaggt    360 

tatagaaaca ttatcacgac tttctcgaac acctatagca ttggccacag gaatcaggcc    420 

cttcccaaca gaagaaagca ttaatgatga agacatctac aaaggccttc ctgatttaat    480 

agatgaaacc cttgtggaag atgaagaaga tctctatgac tgtgtttatg gggaagatga    540 

aggtggagaa gtctatgagg acttaatgaa ggcagaggaa gcacatcagc ccaaatgtcc    600 

agaaaatgat atacgaagtt gttgtctagc agaaattaag cagacagaag aaaaatatac    660 

agaaactttg gagtcaatag aaaaatattt catggcacca ctaaaaagat ttctgacagc    720 

agcagaattt gattcagtat tcatcaacat tcctgaactt gtaaaacttc atcggaacct    780 

aatgcaagag attcatgatt ccattgtaaa taaaaatgac cagaacttgt accaagtttt    840 

tattaactac aaggaaagat tggttattta cgggcagtac tgcagtggag tggagtcagc    900 

catctctagt ttagactaca tttctaagtc aaaagaagat gtcaaactga aattagagga    960 

atgttccaaa agagcaaata atgggaaatt tactcttcga gacttgcttg tggttcctat   1020 

gcaacgtgtt ttaaagtacc accttctcct ccaggaactg gtcaaacata ccactgatcc   1080 

gactgagaag gcaaatctga aactggctct tgatgccatg aaggacttgg cacaatatgt   1140 

gaatgaagtg aaaagagata atgagaccct tcgtgaaatt aaacagtttc agctatctat   1200 

agagaatttg aaccaaccag ttttgctttt tggacgacct cagggagatg gtgaaattcg   1260 

aataaccact ctagacaagc ataccaaaca agaaaggcat atcttcttat ttgatttggc   1320 

agtgatcgta tgtaagagaa aaggtgataa ctatgaaatg aaggaaataa tagatcttca   1380 

gcagtacaag atagccaata atcctacaac cgataaagaa aacaaaaagt ggtcttatgg   1440 

cttctacctc atccataccc aaggacaaaa tgggttagaa ttttattgca aaacaaaaga   1500 

tttaaagaag aaatggctag aacagtttga aatggctttg tctaacataa gaccagacta   1560 

tgcagactcc aatttccacg acttcaagat gcataccttc actcgagtca catcctgcaa   1620 

agtctgccag atgctcctga ggggaacatt ttatcaaggc tatttatgtt ttaagtgtgg   1680 

agcgagagca cacaaagaat gtttgggaag agtagacaat tgtggcagag ttaattctgg   1740 

tgaacaaggg acactcaaac taccagagaa acggaccaat ggactgcgaa gaactcctaa   1800 

acaggtggat ccaggtttac caaagatgca ggtcattagg aactattctg gaacaccacc   1860 

cccagctctg catgaaggac cccctttaca gctccaggcc ggggataccg ttgaacttct   1920 

gaaaggagat gcacacagtc tgttttggca gggcagaaat ttagcatctg gagaggttgg   1980 

attttttcca agtgatgcag tcaagccttg cccatgtgtg cccaaaccag tagattattc   2040 

ttgccaaccc tggtatgctg gagcaatgga aagattgcaa gcagagaccg aacttattaa   2100 

tagggtaaat agtacttacc ttgtgaggca caggaccaaa gagtcaggag aatatgcaat   2160 

tagcattaag tacaataatg aagcaaagca catcaagatt ttaacaagag atggcttttt   2220 

tcacattgca gaaaatagaa aatttaaaag tttaatggaa cttgtggagt actacaagca   2280 

tcattctctc aaggaagggt tcagaacctt agatacaact ctgcagtttc catacaagga   2340 

gccagaacat tcagctggac agaggggtaa tagagcaggc aacagcttgt taagtccaaa   2400 

agtgctgggc attgccatcg ctcggtatga cttctgtgca agagatatga gagagttgtc   2460 

cttgttgaaa ggagatgtgg tgaagattta cacaaagatg agtgcaaatg gctggtggag   2520 

aggagaagta aatggcaggg tgggctggtt tccatccaca tatgtggaag aggatgaata   2580 

aattcaaatc ccgtgttgca ccctgcacca aaaatttcag agaagggata aatagaagcc   2640 

tgcacagcat cgtgaattaa ctgaagtgtt taaaaagctg catttctggc tgttcaacat   2700 

cctccctcct tagcccctcc taagtcttaa tgctgagatt tctaaagatg ctggtactga   2760 

cagattaatg gcttgcctag agctgtgcaa gaaacagcct gccagtctgt cattgtcagg   2820 

gaccagggca aaaccaagag ctgttcttcc cagaagagcc ctgcaaacac attggttcgt   2880 

gcttcccttt acttcttctg gtcagatacc atgaatgcca gtcatcagta aatcttaata   2940 

cacttttgct ttattctcac atgccattca ccagattatt tgatggtaca aagaagcaga   3000 

agtgtaattt tccttttccc agcatgacga aaaattggag ttctgccatt tgagcagctt   3060 

actggagaga tccagcctta cttgtcttaa attgtccaac aaggtgactc attgcccggc   3120 

aaacactttt accctcagat gttactcatg atattataaa atatgaggcc agtgctcagg   3180 

tttgcatcat aagtgagcta tccctgaagg gttttaatta cttatttggt gtcctgatta   3240 

tatttgcaaa cttctttata aaaggtgaaa aaagcacaca aaagagaggg tgtcttcata   3300 

ttaaaccttc acaaccttca tgatttcata ggattatttt ggaaatatag cacttgactt   3360 

tatgaaagga tctgggctag gtatattagg ggtagttgcc aataacctga agaagctggc   3420 

attgtttaca gaaacagatc aagggctata atttatgtca ttttatagca gcagtatcta   3480 

ttaatacatg ccttttcctc ccatccacct cccccgcaca cacacaaaga tgacctggga   3540 

catgattttt ttattcccac attttcttgg agcacaaaca actttgttga ggattttgga   3600 

aggaaagcac aactgggtcc tttattcatt tctgggacag aaagagggtc agtggacttt   3660 

tgtgggcctc cagcttctct cagagtctcc ccctctgcag cccatcctgg gagtgtatta   3720 

actggaggga agatgggtct tgcagtacat ttgttttgcc cagccatcac tcttttttgt   3780 

gaggagccta aatacattct tcctggggtc cagagtcccc attcaaggca gtcaagttaa   3840 

gacactaact tggccctttc ctgatggaaa tatttcctcc atagcagaag ttgtgttctg   3900 

acaagactga gagagttaca tgttgggaaa aaaaagaagc attaacttag tagaactgaa   3960 

ccaggagcat taagttctga aattttgaat catctctgaa atgaagcagg tgtagcctgc   4020 

cctctcatca atccgtccgt ctgggtgcca gaactcaagg ttcagtggac acatccccct   4080 

gttagagacc ctcatgggct aggacttttc atctaggata gattcaagac ctttacctca   4140 

gaattatgta aactgtgatt gtgttttaga aaaattatta tttgctaaaa ccatttaagt   4200 

ctttgtatat gtgtaaatga tcacaaaaat gtattttata aaatgttctg tacaataaag   4260 

ttacacctca aagtgtactc ttggaatgga ttctttcctg taaagtctta tctgcgactc   4320 

tgtctcggga atgttttgtc tgttgccgtc agccgaactt tgttatggag ggagcagcct   4380 

cacacaagca gaaacactcc tgtggatggt attgtagcat gtattgttta ttttagtcaa   4440 

tagaccctct ccttataaat ggtgtttagt cttcctgttg catttcatgg gcctgggggt   4500 

ttcctrgcag aggatattgg agcccctttt tgtgacatta ccaattacat ctttgtccac   4560 

gtttaatact ttgttttgga aaatttaaat gctgcagatt tgtgtagagt tctaatacca   4620 

aagacagaag taaatgtttt ccatatactt tgtcttgcct gtatgcagcc cttgtgtaat   4680 

atggtgaatt agagtggtat ttcactttgt attattttgt aaatatgtca atataataaa   4740 

tagtgactaa aaaaaaaaaa aa                                            4762 

 
           
             31  
             422  
             DNA  
             Homo sapiens  
             
               unsure  
               (80)  
             
             
               unsure  
               (217)  
             
             
               unsure  
               (254)  
             
             
               unsure  
               (311)  
             
             
               unsure  
               (321)  
             
             
               unsure  
               (375)  
             
             
               unsure  
               (381)  
             
             
               unsure  
               (386)  
             
             
               unsure  
               (394)  
             
             
               unsure  
               (416)  
             
           
            31 

ttttttactt tattttcgtt ttaatttttt ggaaggatat acaccacata tcccatgggc     60 

aataaagcgc attcaatgtn tttataagcc aaacagtcac tttgtttaag caaacacaag    120 

tacaaagtaa aatagaacca caaaataatg aactgcatgt tcataacata caaaaatcgc    180 

cgcctactca gtaggtaact acaacattcc aactccngaa tatatttata aatttacatt    240 

ttcagttaaa aaantagact tttgagagtt cagattttgt tttagatttt gttttcttac    300 

attctggaga ncccgaagct ncagctcagc ccctcttccc ttattttgct ccccaaagcc    360 

ttccccccaa atcancactg ncctgncccc cctntaaggg cttagaggtg agcatntccc    420 

ct                                                                   422 

 
           
             32  
             3132  
             DNA  
             Homo sapiens  
           
            32 

ccgcagaact tggggagccg ccgccgccat ccgccgccgc agccagcttc cgccgccgca     60 

ggaccggccc ctgccccagc ctccgcagcc gcggcgcgtc cacgcccgcc cgcgcccagg    120 

gcgagtcggg gtcgccgcct gcacgcttct cagtgttccc cgcgccccgc atgtaacccg    180 

gccaggcccc cgcaacggtg tcccctgcag ctccagcccc gggctgcacc cccccgcccc    240 

gacaccagct ctccagcctg ctcgtccagg atggccgcgg ccaaggccga gatgcagctg    300 

atgtccccgc tgcagatctc tgacccgttc ggatcctttc ctcactcgcc caccatggac    360 

aactacccta agctggagga gatgatgctg ctgagcaacg gggctcccca gttcctcggc    420 

gccgccgggg ccccagaggg cagcggcagc aacagcagca gcagcagcag cgggggcggt    480 

ggaggcggcg ggggcggcag caacagcagc agcagcagca gcaccttcaa ccctcaggcg    540 

gacacgggcg agcagcccta cgagcacctg accgcagagt cttttcctga catctctctg    600 

aacaacgaga aggtgctggt ggagaccagt taccccagcc aaaccactcg actgcccccc    660 

atcacctata ctggccgctt ttccctggag cctgcaccca acagtggcaa caccttgtgg    720 

cccgagcccc tcttcagctt ggtcagtggc ctagtgagca tgaccaaccc accggcctcc    780 

tcgtcctcag caccatctcc agcggcctcc tccgcctccg cctcccagag cccacccctg    840 

agctgcgcag tgccatccaa cgacagcagt cccatttact cagcggcacc caccttcccc    900 

acgccgaaca ctgacatttt ccctgagcca caaagccagg ccttcccggg ctcggcaggg    960 

acagcgctcc agtacccgcc tcctgcctac cctgccgcca agggtggctt ccaggttccc   1020 

atgatccccg actacctgtt tccacagcag cagggggatc tgggcctggg caccccagac   1080 

cagaagccct tccagggcct ggagagccgc acccagcagc cttcgctaac ccctctgtct   1140 

actattaagg cctttgccac tcagtcgggc tcccaggacc tgaaggccct caataccagc   1200 

taccagtccc agctcatcaa acccagccgc atgcgcaagt atcccaaccg gcccagcaag   1260 

acgccccccc acgaacgccc ttacgcttgc ccagtggagt cctgtgatcg ccgcttctcc   1320 

cgctccgacg agctcacccg ccacatccgc atccacacag gccagaagcc cttccagtgc   1380 

cgcatctgca tgcgcaactt cagccgcagc gaccacctca ccacccacat ccgcacccac   1440 

acaggcgaaa agcccttcgc ctgcgacatc tgtggaagaa agtttgccag gagcgatgaa   1500 

cgcaagaggc ataccaagat ccacttgcgg cagaaggaca agaaagcaga caaaagtgtt   1560 

gtggcctctt cggccacctc ctctctctct tcctacccgt ccccggttgc tacctcttac   1620 

ccgtccccgg ttactacctc ttatccatcc ccggccacca cctcataccc atcccctgtg   1680 

cccacctcct tctcctctcc cggctcctcg acctacccat cccctgtgca cagtggcttc   1740 

ccctccccgt cggtggccac cacgtactcc tctgttcccc ctgctttccc ggcccaggtc   1800 

agcagcttcc cttcctcagc tgtcaccaac tccttcagcg cctccacagg gctttcggac   1860 

atgacagcaa ccttttctcc caggacaatt gaaatttgct aaagggaaag gggaaagaaa   1920 

gggaaaaggg agaaaaagaa acacaagaga cttaaaggac aggaggagga gatggccata   1980 

ggagaggagg gttcctctta ggtcagatgg aggttctcag agccaagtcc tccctctcta   2040 

ctggagtgga aggtctattg gccaacaatc ctttctgccc acttcccctt ccccaattac   2100 

tattcccttt gacttcagct gcctgaaaca gccatgtcca agttcttcac ctctatccaa   2160 

agaacttgat ttgcatggat tttggataaa tcatttcagt atcatctcca tcatatgcct   2220 

gaccccttgc tcccttcaat gctagaaaat cgagttggca aaatggggtt tgggcccctc   2280 

agagccctgc cctgcaccct tgtacagtgt ctgtgccatg gatttcgttt ttcttggggt   2340 

actcttgatg tgaagataat ttgcatattc tattgtatta tttggagtta ggtcctcact   2400 

tgggggaaaa aaaaaaaaaa aagccaagca aaccaatggt gatcctctat tttgtgatga   2460 

tgctgtgaca ataagtttga accttttttt ttgaaacagc agtcccagta ttctcagagc   2520 

atgtgtcaga gtgttgttcc gttaaccttt ttgtaaatac tgcttgaccg tactctcaca   2580 

tgtggcaaaa tatggtttgg tttttctttt ttttttttga aagtgttttt tcttcgtcct   2640 

tttggtttaa aaagtttcac gtcttggtgc cttttgtgtg atgccccttg ctgatggctt   2700 

gacatgtgca attgtgaggg acatgctcac ctctagcctt aaggggggca gggagtgatg   2760 

atttggggga ggctttggga gcaaaataag gaagagggct gagctgagct tcggttctcc   2820 

agaatgtaag aaaacaaaat ctaaaacaaa atctgaactc tcaaaagtct atttttttaa   2880 

ctgaaaatgt aaatttataa atatattcag gagttggaat gttgtagtta cctactgagt   2940 

aggcggcgat ttttgtatgt tatgaacatg cagttcatta ttttgtggtt ctattttact   3000 

ttgtacttgt gtttgcttaa acaaagtgac tgtttggctt ataaacacat tgaatgcgct   3060 

ttattgccca tgggatatgt ggtgtatatc cttccaaaaa attaaaacga aaataaagta   3120 

gctgcgattg gg                                                       3132 

 
           
             33  
             464  
             DNA  
             Homo sapiens  
             
               unsure  
               (364)  
             
             
               unsure  
               (401)  
             
             
               unsure  
               (425)  
             
             
               unsure  
               (439)  
             
           
            33 

ttaaggtata cacttttatt caactggtct caagtcagtg tacaggtaag ccctggctgc     60 

ctccacccac tcccagggag accaaaagcc ttcatacatc tcaagttggg ggacaaaaaa    120 

gggggaaggg ggggcacgaa ggctcatcat tcaaaataaa acaaaataaa aaagtattaa    180 

ggcgaagatt aaaaaaattt tgcattacat aatttacacg aaagcaatgc tatcacctcc    240 

cctgtgtgga cttgggagag gactgggcca ttctccttag gagagaagtg ggggtgggct    300 

tttagggatg ggcaagggga ctttcctgtt aacaacggca tcttcatatt ttgggaattg    360 

actntttaaa aaaaaccaac aatgtggcaa ttcaaagtcc ntcgggccac atttgtggaa    420 

ctttnggggg gttgctcgnt cccacccgac tgttgttcac cttt                     464 

 
           
             34  
             3646  
             DNA  
             Homo sapiens  
           
            34 

gcccagcacc ccaaggcggc caacgccaaa actctccctc ctcctcttcc tcaatctcgc     60 

tctcgctctt tttttttttc gcaaaaggag gggagagggg gtaaaaaaat gctgcactgt    120 

gcggcgaagc cggtgagtga gcggcgcggg gccaatcagc gtgcgccgtt ccgaaagttg    180 

ccttttatgg ctcgagcggc cgcggcggcg ccctataaaa cccagcggcg cgacgcgcca    240 

ccaccgccga gaccgcgtcc gcccgcgagc acagagcctc gcctttgccg atccgccgcc    300 

cgtccacacc cgccgccagg taagcccggc cagccgaccg gggcatgcgg ccgcggccct    360 

tcgcccgtgc agagccgccg tctgggccgc agcggggggc gcatggggcg gaaccggacc    420 

gccgtggggg gcgcgggaga agcccctggg cctccggaga tgggggacac cccacgccag    480 

ttcgcaggcg cgaggccgcg ctcgggcggg cgcgctccgg gggtgccgct ctcggggcgg    540 

gggcaaccgg cggggtcttt gtctgagccg ggctcttgcc aatggggatc gcacggtggg    600 

cgcggcgtag cccccgtcag gcccggtggg ggctggggcg ccatgcgcgt gcgcgctggt    660 

cctttgggcg ctaactgcgt gcgcgctggg aattggcgct aattgcgcgt gcgcgctggg    720 

actcaatggc gctaatcgcg cgtgcgttct ggggcccggg cgcttgcgcc acttcctgcc    780 

cgagccgctg gcgcccgagg gtgtggccgc tgcgtgcgcg cgcgcgaccc ggtcgctgtt    840 

tgaaccgggc ggaggcgggg ctggcgcccg gttgggaggg ggttggggcc tggcttcctg    900 

ccgcgcgccg cggggacgcc tccgaccagt gtttgccttt tatggtaata acgcggccgg    960 

cccggcttcc tttgtcccca atctgggcgc gcgccggcgc cccctggcgg cctaaggact   1020 

cggcgcgccg gaagtggcca gggcgggggc gacttcggct cacagcgcgc ccggctattc   1080 

tcgcagctca ccatggatga tgatatcgcc gcgctcgtcg tcgacaacgg ctccggcatg   1140 

tgcaaggccg gcttcgcggg cgacgatgcc ccccgggccg tcttcccctc catcgtgggg   1200 

cgccccaggc accaggtagg ggagctggct gggtggggca gccccgggag cgggcgggag   1260 

gcaagggcgc tttctctgca caggagcctc ccggtttccg gggtgggctg cgcccgtgct   1320 

cagggcttct tgtcctttcc ttcccagggc gtgatggtgg gcatgggtca gaaggattcc   1380 

tatgtgggcg acgaggccca gagcaagaga ggcatcctca ccctgaagta ccccatcgag   1440 

cacggcatcg tcaccaactg ggacgacatg gagaaaatct ggcaccacac cttctacaat   1500 

gagctgcgtg tggctcccga ggagcacccc gtgctgctga ccgaggcccc cctgaacccc   1560 

aaggccaacc gcgagaagat gacccaggtg agtggcccgc tacctcttct ggtggccgcc   1620 

tccctccttc ctggcctccc ggagctgcgc cctttctcac tggttctctc ttctgccgtt   1680 

ttccgtagga ctctcttctc tgacctgagt ctcctttgga actctgcagg ttctatttgc   1740 

tttttcccag atgagctctt tttctggtgt ttgtctctct gactaggtgt ctgagacagt   1800 

gttgtgggtg taggtactaa cactggctcg tgtgacaagg ccatgaggct ggtgtaaagc   1860 

ggccttggag tgtgtattaa gtaggcgcac agtaggtctg aacagactcc ccatcccaag   1920 

accccagcac acttagccgt gttctttgca ctttctgcat gtcccccgtc tggcctggct   1980 

gtccccagtg gcttccccag tgtgacatgg tgcatctctg ccttacagat catgtttgag   2040 

accttcaaca ccccagccat gtacgttgct atccaggctg tgctatccct gtacgcctct   2100 

ggccgtacca ctggcatcgt gatggactcc ggtgacgggg tcacccacac tgtgcccatc   2160 

tacgaggggt atgccctccc ccatgccatc ctgcgtctgg acctggctgg ccgggacctg   2220 

actgactacc tcatgaagat cctcaccgag cgcggctaca gcttcaccac cacggccgag   2280 

cgggaaatcg tgcgtgacat taaggagaag ctgtgctacg tcgccctgga cttcgagcaa   2340 

gagatggcca cggctgcttc cagctcctcc ctggagaaga gctacgagct gcctgacggc   2400 

caggtcatca ccattggcaa tgagcggttc cgctgccctg aggcactctt ccagccttcc   2460 

ttcctgggtg agtggagact gtctcccggc tctgcctgac atgagggtta cccctcgggg   2520 

ctgtgctgtg gaagctaagt cctgccctca tttccctctc aggcatggag tcctgtggca   2580 

tccacgaaac taccttcaac tccatcatga agtgtgacgt ggacatccgc aaagacctgt   2640 

acgccaacac agtgctgtct ggcggcacca ccatgtaccc tggcattgcc gacaggatgc   2700 

agaaggagat cactgccctg gcacccagca caatgaagat caaggtgggt gtctttcctg   2760 

cctgagctga cctgggcagg tcagctgtgg ggtcctgtgg tgtgtgggga gctgtcacat   2820 

ccagggtcct cactgcctgt ccccttccct cctcagatca ttgctcctcc tgagcgcaag   2880 

tactccgtgt ggatcggcgg ctccatcctg gcctcgctgt ccaccttcca gcagatgtgg   2940 

atcagcaagc aggagtatga cgagtccggc ccctccatcg tccaccgcaa atgcttctag   3000 

gcggactatg acttagttgc gttacaccct ttcttgacaa aacctaactt gcgcagaaaa   3060 

caagatgaga ttggcatggc tttatttgtt ttttttgttt tgttttggtt tttttttttt   3120 

ttttggcttg actcaggatt taaaaactgg aacggtgaag gtgacagcag tcggttggag   3180 

cgagcatccc ccaaagttca caatgtggcc gaggactttg attgcattgt tgttttttta   3240 

atagtcattc caaatatgag atgcattgtt acaggaagtc ccttgccatc ctaaaagcca   3300 

ccccacttct ctctaaggag aatggcccag tcctctccca agtccacaca ggggaggtga   3360 

tagcattgct ttcgtgtaaa ttatgtaatg caaaattttt ttaatcttcg ccttaatact   3420 

tttttatttt gttttatttt gaatgatgag ccttcgtgcc cccccttccc cctttttgtc   3480 

ccccaacttg agatgtatga aggcttttgg tctccctggg agtgggtgga ggcagccagg   3540 

gcttacctgt acactgactt gagaccagtt gaataaaagt gcacacctta aaaatgaggc   3600 

caagtgtgac tttgtggtgt ggctgggttg ggggcagcag agggtg                  3646 

 
           
             35  
             318  
             DNA  
             Homo sapiens  
             
               unsure  
               (9)  
             
             
               unsure  
               (51)  
             
             
               unsure  
               (119)  
             
             
               unsure  
               (247)  
             
             
               unsure  
               (301)  
             
             
               unsure  
               (313)  
             
           
            35 

ctcgatttng ggaagttgta gactgcacaa ttaaaacaga tccagtcact nggagatcaa     60 

gaggatttgg atttgtgctt ttcaaagatg ctgctagtgt tgataaggtt ttggaactna    120 

aagaacacaa actggatggc aaattgatag atcccaaaag ggccaaagct ttaaaaggga    180 

aagaacctcc caaaaaggtt tttgtgggtg gattgagccc ggatacttct gaagaacaaa    240 

ttaaagnata ttttggagcc tttggagaga ttgaaaatat tgaacttccc atggatacaa    300 

naacaaattg aanggaag                                                  318 

 
           
             36  
             1291  
             DNA  
             Homo sapiens  
           
            36 

gatctcttcc gccgccattt taaatccagc tccatacaac gctccgccgc cgctgctgcc     60 

gcgacccgga ctgcgcgcca gcacccccct gccgacagct ccgtcactat ggaggatatg    120 

aacgagtaca gcaatataga ggaattcgca gagggatcca agatcaacgc gagcaagaat    180 

cagcaggatg acggtaaaat gtttattgga ggcttgagct gggatacaag caaaaaagat    240 

ctgacagagt acttgtctcg atttggggaa gttgtagact gcacaattaa aacagatcca    300 

gtcactggga gatcaagagg atttggattt gtgcttttca aagatgctgc tagtgttgat    360 

aaggttttgg aactgaaaga acacaaactg gatggcaaat tgatagatcc caaaagggcc    420 

aaagctttaa aagggaaaga acctcccaaa aaggtttttg tgggtggatt gagcccggat    480 

acttctgaag aacaaattaa agaatatttt ggagcctttg gagagattga aaatattgaa    540 

cttcccatgg atacaaaaac aaatgaaaga agaggatttt gttttatcac atatactgat    600 

gaagagccag taaaaaaatt gttagaaagc agataccatc aaattggttc tgggaagtgt    660 

gaaatcaaag ttgcacaacc caaagaggta tataggcagc aacagcaaca acaaaaaggt    720 

ggaagaggtg ctgcagctgg tggacgaggt ggtacgaggg gtcgtggccg aggtcagggc    780 

caaaactgga accaaggatt taataactat tatgatcaag gatatggaaa ttacaatagt    840 

gcctatggtg gtgatcaaaa ctatagtggc tatggcggat atgattatac tgggtataac    900 

tatgggaact atggatatgg acagggatat gcagactaca gtggccaaca gagcacttat    960 

ggcaaggcat ctcgaggggg tggcaatcac caaaacaatt accagccata ctaaaggaga   1020 

acattggaga aaacaggagg agatgttaaa gtaacccatc ttgcaggacg acattgaaga   1080 

ttggtcttct gttgatctaa gatgattatt ttgtaaaaga ctttctagtg tacaagacac   1140 

cattgtgtcc aactgtatat agctgccaat tagttttctt tgtttttact ttgtcctttg   1200 

ctatctgtgt tatgactcaa tgtggatttg tttatacaca ttttatttgt atcatttcat   1260 

gttaaacctc aaataaatgc ttccttatgt g                                  1291 

 
           
             37  
             2439  
             DNA  
             Homo sapiens  
           
            37 

gaattcgcag agggatccaa gatcaacgcg agcaagaatc agcaggatga cggtaaaatg     60 

tttattggag gcttgagctg ggatacaagc aaaaaagatc tgacagagta cttgtctcga    120 

tttggggaag ttgtagactg cacaattaaa acagatccag tcactgggag atcaagagga    180 

tttggatttg tgcttttcaa agatgctgct agtgttgata aggttttgga actgaaagaa    240 

cacaaactgg atggcaaatt gatagatccc aaaagggcca aagctttaaa agggaaagaa    300 

cctcccaaaa aggtttttgt gggtggattg agcccggata cttctgaaga acaaattaaa    360 

gaatattttg gagcctttgg agagattgaa aatattgaac ttcccatgga tacaaaaaca    420 

aatgaaagaa gaggattttg ttttatcaca tatactgatg aagagccagt aaaaaaattg    480 

ttagaaagca gataccatca aattggttct gggaagtgtg aaatcaaagt tgcacaaccc    540 

aaagaggtat ataggcagca acagcaacaa caaaaaggtg gaagaggtgc tgcagctggt    600 

ggacgaggtg gtacgagggg tcgtggccga ggtcagggcc aaaactggaa ccaaggattt    660 

aataactatt atgatcaagg atatggaaat tacaatagtg cctatggtgg tgatcaaaac    720 

tatagtggct atggcggata tgattatact gggtataact atgggaacta tggatatgga    780 

cagggatatg cagactacag tggccaacag agcacttatg gcaaggcatc tcgagggggt    840 

ggcaatcacc aaaacaatta ccagccatac taaaggagaa cattggagaa aacaggagga    900 

gatgttaaag taacccatct tgcaggacga cattgaagat tggtcttctg ttgatctaag    960 

atgattattt tgtaaaagac tttctagtgt acaagacacc attgtgtcca actgtatata   1020 

gctgccaatt agttttcttt gtttttactt tgtcctttgc tatctgtgtt atgactcaat   1080 

gtggatttgt ttatacacat tttatttgta tcatttcatg ttaaacctca aataaatgct   1140 

tccttatgtg attgcttttc tgcgtcaggt actacatagc tctgtaaaaa atgtaattta   1200 

aaataagcaa taattaaggc acagttgatt ttgtagagta ttggtccata cagagaaact   1260 

gtggtccttt ataaatagcc agccagcgtc accctcttct ccaatttgta ggtgtatttt   1320 

atgctcttaa ggcttcatct tctccctgta actgagattt ctaccacacc tttgaacaat   1380 

gttctttccc ttctggttat ctgaagactg tcctgaaagg aagacataag tgttgtgatt   1440 

agtagaagct ttgtaatcat aacacaatga gtaattcttg tataaaagtt cagatacaaa   1500 

aggagcactg taaaactggt aggagctatg gtttaagagc attggaagta gttacaactc   1560 

aaggattttg gtagaaaggt atgagtttgg tcgaaaaatt aaaatagtgg caaaataaga   1620 

tttagttgtg ttttctcaga gccgccacaa gattgaacaa aatgttttct gtttgggcat   1680 

cctgaggaag ttgtattagc tgttaatgct ctgtgagttt agagaaaagt cttgatagta   1740 

aatctagttt ttgacacagt gcatgaacta agtagttaaa tatttacata ttcagaaagg   1800 

aatagtggaa aaggtatctt ggttatgaca aagtcattac aaatgtgact aagtcattac   1860 

aaatgtgact gagtcattac agtggaccct ctgggtgcat tgaaaagaat ccgttttata   1920 

tccaggtttc agaggacctg gaataataat aagctttgga ttttgcattc agtgtagttg   1980 

gattttggga ccttggcctc agtgttattt actgggattg gcatacgtgt tcacaggcag   2040 

agtagttgat ctcacacaac gggtgatctc acaaaactgg taagtttctt atgctcatga   2100 

gccctccctt ttttttttta atttggtgcc tgcaactttc ttaacaatga ttctacttcc   2160 

tgggctatca cattataatg ctcttggcct cttttttgct gctgttttgc tattcttaaa   2220 

cttaggccaa gtaccaatgt tggctgttag aagggattct gttcattcaa catgcaactt   2280 

tagggaatgg aagtaagttc atttttaagt tgtgtggtca gtaggtgcgg tgtctagggt   2340 

agtgaatcct gtaagttcaa atttatgatt aggtgacgag ttgacattga gattgtcctt   2400 

ttcccctgat caaaaaaatg aataaagcct ttttaaacg                          2439 

 
           
             38  
             459  
             DNA  
             Homo sapiens  
             
               unsure  
               (426)  
             
             
               unsure  
               (445)  
             
           
            38 

ttttacagat ctttttgact atcttcctct cactgccttg gtggatgggc agatcttctg     60 

tctacatggt ggtctctcgc catctataga tacactggat catatcagag cacttgatcg    120 

cctacaagaa gttccccatg agggtccaat gtgtgacttg ctgtggtcag atccagatga    180 

ccgtggtggt tggggtatat ctcctcgagg agctggttac acctttgggc aagatatttc    240 

tgagacattt aatcatgcca atggcctcac gttggtgtct agagctcacc agctagtgat    300 

ggagggatat aactggtgcc atgaccggaa tgtagtaacg attttcagtg ctccaaacta    360 

ttgttatcgt tgtggtaacc aagctgcaat catgggaact tgacgatact ctaaaatact    420 

ctttcntgca gttttgaccc agcanctcgt agggccgag                           459 

 
           
             39  
             1787  
             DNA  
             Homo sapiens  
           
            39 

gagagctcgg ctctcggagg aggaggcgca cggccagcgg cagtactgcg gtgagagcca     60 

gcggccagcg ccacgctcaa cagccgccag aagtacacga ggaaccggcg gcggcgtgtg    120 

cgtgtaagcc ggcggcggcg cgggaggagc cggagcggca gccggctggg gcgggtggca    180 

tcatggacga gaaggtgttc accaaggagc tggaccagtg gatcgagcag ctgaacgagt    240 

gcaagcagct gtccgagtcc caggtcaaga gcctctgcga gaaggctaaa gaaatcctga    300 

caaaagaatc caacgtgcaa gaggttcgat gtccagttac tgtctgtgga gatgtgcatg    360 

ggcaatttca tgatctcatg gaactgttta gaattggtgg caaatcacca gatacaaatt    420 

acttgtttat gggagattat gttgacagag gatattattc agttgaaaca gttacactgc    480 

ttgtagctct taaggttcgt taccgtgaac gcatcaccat tcttcgaggg aatcatgaga    540 

gcagacagat cacacaagtt tatggtttct atgatgaatg tttaagaaaa tatggaaatg    600 

caaatgtttg gaaatatttt acagatcttt ttgactatct tcctctcact gccttggtgg    660 

atgggcagat cttctgtcta catggtggtc tctcgccatc tatagataca ctggatcata    720 

tcagagcact tgatcgccta caagaagttc cccatgaggg tccaatgtgt gacttgctgt    780 

ggtcagatcc agatgaccgt ggtggttggg gtatatctcc tcgaggagct ggttacacct    840 

ttgggcaaga tatttctgag acatttaatc atgccaatgg cctcacgttg gtgtctagag    900 

ctcaccagct agtgatggag ggatataact ggtgccatga ccggaatgta gtaacgattt    960 

tcagtgctcc aaactattgt tatcgttgtg gtaaccaagc tgcaatcatg gaacttgacg   1020 

atactctaaa atactctttc ttgcagtttg acccagcacc tcgtagaggc gagccacatg   1080 

ttactcgtcg taccccagac tacttcctgt aatgaaattt taaacttgta cagtattgcc   1140 

atgaaccata tatcgaccta atggaaatgg gaagagcaac agtaactcca aagtgtcaga   1200 

aaatagttaa cattcaaaaa acttgttttc acatggacca aaagatgtgc catataaaaa   1260 

tacaaagcct cttgtcatca acagccgtga ccactttaga atgaaccagt tcattgcatg   1320 

ctgaagcgac attgttggtc aagaaaccag tttctggcat agcgctattt gtagttactt   1380 

ttgtttctct gagagactgc agataataag atgtaaacat taacacctcg tgaatacaat   1440 

ttaacttcca tttagctata gctttactca gcatgactgt agataaggat agcagcaaac   1500 

aatcattgga gcttaatgaa catttttaaa aataattacc aaggcctccc ttctacttgt   1560 

gagttttgaa attgttcttt ttattttcag ggataccgtt taatttaatt atatgatttg   1620 

tctgcactca gtttattccc tactcaaatc tcagccccat gttgttcttt gttattgtca   1680 

gaacctggtg agttgttttg aacagaactg ttttttcccc ttcctgtaag acgatgtgac   1740 

tgcacaagag cactgcagtg tttttcataa taaacttgtg aactaac                 1787 

 
           
             40  
             452  
             DNA  
             Homo sapiens  
             
               unsure  
               (33)..(34)  
             
             
               unsure  
               (59)  
             
             
               unsure  
               (82)  
             
             
               unsure  
               (112)  
             
             
               unsure  
               (126)  
             
             
               unsure  
               (164)  
             
             
               unsure  
               (184)  
             
             
               unsure  
               (225)  
             
             
               unsure  
               (244)  
             
             
               unsure  
               (253)  
             
             
               unsure  
               (272)  
             
             
               unsure  
               (307)  
             
             
               unsure  
               (316)  
             
             
               unsure  
               (329)  
             
             
               unsure  
               (335)  
             
             
               unsure  
               (381)  
             
             
               unsure  
               (396)  
             
             
               unsure  
               (417)  
             
             
               unsure  
               (422)  
             
             
               unsure  
               (429)  
             
             
               unsure  
               (448)  
             
           
            40 

gtttacagat gccacttagt tacactggtt ttnntttttc agtctcatct gggttgganc     60 

caaagacatt cagaggcatg gnaagaggca aagcatcaga catctcattg gnggcaggta    120 

cttccngact actgtaccac ctgctgtatc cttccccacc tcancacccc caaagccatt    180 

tagngccaaa tgctacagta aaaacccaat gcatttacat aaaanaatgc ctaactgcat    240 

attnacattt ttnagaaaaa aaatcccatt angctcttct agaaagttat ggcaggaaag    300 

gtaaggncca aggctntgag caagccatnt gtggnaactt aaagtagatg agcactgagt    360 

ttctccatag ttggaaaaaa ngccacactg agcccncttt tcccgtggag ggcaagntga    420 

gnccctccnt ttataccccg ttgagatntc ag                                  452 

 
           
             41  
             263  
             DNA  
             Homo sapiens  
             
               unsure  
               (214)  
             
             
               unsure  
               (231)  
             
             
               unsure  
               (238)  
             
           
            41 

gagaaaaggg ttggggagaa gcctctgcag tcctggaaga tgtggggttc tgggtgagag     60 

gcatcagccc cacaagtatg tttttgtgtc ttaagatagc agtttacttt gaaaaagtga    120 

aaaaggcttc cgggctgtcc tctgcccagt gagatggagg acgctagaga aagtgctgag    180 

tgtcccgaga gaggcccccg agccagtgca tggnaggtcc ttcggcctgg ntcagctngg    240 

ctgcaggatg cccactttga gga                                            263 

 
           
             42  
             3049  
             DNA  
             Homo sapiens  
           
            42 

cccgcgggca ggggcggcga gtgcgcgggc cgccgccctt ctcggcgggc agcgcgcgag     60 

gaccaggccg aggaggaagt ggcggcggcg gcggcgggct ccccgcccga ggaggaagat    120 

gcagaccttt ctgaaaggga agagagttgg ctactggctg agcgagaaga aaatcaagaa    180 

gctgaatttc caggctttcg ccgagctgtg caggaagcga gggatggagg ttgtgcagct    240 

gaaccttagc cggccgatcg aggagcaggg ccccctggac gtcatcatcc acaagctgac    300 

tgacgtcatc cttgaagccg accagaatga tagccagtcc ctggagctgg tgcacaggtt    360 

ccaggagtac atcgatgccc accctgagac catcgtcctg gacccgctcc ctgccatcag    420 

aaccctgctt gaccgctcca agtcctatga gctcatccgg aagattgagg cctacatgga    480 

agacgacagg atctgctcgc cacccttcat ggagctcacg agcctgtgcg gggatgacac    540 

catgcggctg ctggagaaga acggcttgac tttcccattc atttgcaaaa ccagagtggc    600 

tcatggcacc aactctcacg agatggctat cgtgttcaac caggagggcc tgaacgccat    660 

ccagccaccc tgcgtggtcc agaatttcat caaccacaac gccgtcctgt acaaggtgtt    720 

cgtggttggc gagtcctaca ccgtggtcca gaggccctca ctcaagaact tctccgcagg    780 

cacatcagac cgtgagtcca tcttcttcaa cagccacaac gtgtcaaagc cggagtcgtc    840 

atcggtcctg acggagctgg acaagatcga gggcgtgttc gagcggccga gcgacgaggt    900 

catccgggag ctctcccggg ccctgcggca ggcactgggc gtgtcactct tcggcatcga    960 

catcatcatc aacaaccaga cagggcagca cgccgtcatt gacatcaatg ccttcccagg   1020 

ctacgagggc gtgagcgagt tcttcacaga cctcctgaac cacatcgcca ctgtcctgca   1080 

gggccagagc acagccatgg cagccacagg ggacgtggcc ctgctgaggc acagcaagct   1140 

tctggccgag ccggcgggcg gcctggtggg cgagcggaca tgcaacgcca gccccggctg   1200 

ctgcggcagc atgatgggcc aggacgcgcc ctggaaagct gaggccgacg cgggcggcac   1260 

cgccaagctg ccgcaccaga gactcggctg caacgccggc gtgtctccca gcttccagca   1320 

gcattgtgtg gcctccctgg ccaccaaggc ctcctcccag tagccacgga gccgggaccc   1380 

agagggcagc gcaggcgcag gagcacaccc gctgggccag cagctcccaa cggcgatgct   1440 

actactaaga atccccagtg atctgattct tctgtttttt aatttttaac ctgattttct   1500 

gatgtcatga tctaaatgag gggtagaaga gagtaccagg tggtccaccg ttggggagcg   1560 

gggccgtccg cctgctctct actgtgcaga cctcctaact gagtttacac acgcttgtgt   1620 

tgcaacacta ggtctggatg ggaggtgagg ggggtgcgta tactgccatg ccagtgtctg   1680 

tgcacatccc tgtctgttgt ctccatggcc actgtggact gggacccttg aagcctgccc   1740 

atgtgggtgt gggaggctga tcagtgcgtg tgagagtggc ttcccttctg cctgactccc   1800 

cactccctga cctgcccctt ccttgttttt cctcctactg gtctccacca aggctttgtt   1860 

agcccccacc ctgcctggtg tgcagctaac ccctccctcc ccacagccag aggaggccac   1920 

agacccctca gggagttccg cgctggggtc tgggctgtgc tccctcacta aagggaagga   1980 

aaggaagctg ggcgtcctcc gggcccccca acacacgtcc catttagccc tgcacagcgg   2040 

tctccttccc ctaagccagc actgctgctc cctggagccg ggaaggaggc tgcctggctg   2100 

gaggccgagc cgatgggcct gtgctgagga tttgtgctgt gatttgggca aatcattcca   2160 

ggtctttggg cctccacccc ctcgtctcta gtggacattt gagatcagag agcaccacag   2220 

ggctggcttt gtgccctaac ccctgggatg cagcctgcct ttccataaag tcacctaggt   2280 

gaggataggc gcgggagcct cggcatgaca ccatggagat cggggccctc ttcccagtgg   2340 

gttcactcct tttcacacct gctgggtccc tcctcgccca gcaggcctgg tccacctctc   2400 

attgcaagcc cgcaagcact gagccgagta aggtgcttag tgtgagccac ccgcccccca   2460 

tagcttctgc acacctcaga ctcaccccat caccttggca gcaaagcact gctctgccgt   2520 

ctgacccctg atccaggcag cagccccctc cgcagagaaa agggttgggg agaagcctct   2580 

gcagtcctgg aagatgtggg gtgctgggtg agaggcatca gcccccacaa gtatgttttt   2640 

gtgtcttaag atagcagttt actttgaaaa agtgaaaaag gcttccgggc tgtcctctgc   2700 

ccagtgagat ggaggacgct agagaaagtg ctgagtgtcc cgagagaggc ccccgagcca   2760 

gtgcatggag gtcttcggcc tggctcagct gggctgcagg atgcccactt tgaggaggga   2820 

ggcacagggc ttgggcgagg ggcagaggcc atcagaactg cccggctttt ttggaaactg   2880 

aggacccaac aactaaccac gtttacacga cttgagtttt gaaccccgat taatgtctgt   2940 

acgtcacctt tcctagttct gaccctgagc cctggggaac aggaaagcgt ggctggcctc   3000 

ttgcactgct ttgtctccaa aataaactac tgaaatcaaa ccgcatttc               3049 

 
           
             43  
             417  
             DNA  
             Homo sapiens  
             
               unsure  
               (198)  
             
             
               unsure  
               (260)  
             
             
               unsure  
               (299)  
             
             
               unsure  
               (344)  
             
             
               unsure  
               (373)  
             
             
               unsure  
               (378)  
             
             
               unsure  
               (384)  
             
             
               unsure  
               (410)  
             
           
            43 

ggttgagccc tacaactgca tcctcaccac ccacaccacc ctggagcact ctgattgtgc     60 

cttcatggta gacaatgagg ccatctatga catctgtcgt agaaacctcg atatcgagcg    120 

cccaacctac accaacctta accgccttat tagccagatt gtgtcctcca tcactgcttc    180 

cctgagattt gatggagncc tgaatgttga cctgacagaa ttccagacca acctgggtgc    240 

cctacccccg catccacttn cctctggcca catatgcccc tgtcatctct gctgagaang    300 

cctaccacga acagcttact gtagtagaga tcaccaatgc ttgntttgag ccagccaacc    360 

agatggtgaa atntggancc ttgncattgg taaattacat ggggtttgcn gtctgtt       417 

 
           
             44  
             1596  
             DNA  
             Homo sapiens  
           
            44 

tgtcggggac ggtaaccggg acccgtgctc tgctcctgtc gccttcgcct cctgaatccc     60 

tagccatatg cgtgagtgca tctccatcca cgttggccag gctggtgtcc agattggcaa    120 

tgcctgctgg gagctctact gcctggaaca cggcatccag cccgatggcc agatgccaag    180 

tgacaagacc attgggggag gagatgactc cttcaacacc ttcttcagtg agacgggcgc    240 

tggcaagcac gtgccccggg ctgtgtttgt agacttggaa cccacagtca ttgatgaagt    300 

tcgcactggc acctaccgcc agctcttcca ccctgagcag ctcatcacag gcaaggaaga    360 

tgctgccaat aactatgccc gagggcacta caccattggc aaggagatca ttgaccttgt    420 

gttggaccga attcgcaagc tggctgacca gtgcacccgt cttcagggct tcttggtttt    480 

ccacagcttt ggtgggggaa ctggttctgg gttcacctcc ctgctcatgg aacgcctgtc    540 

agttgattat ggcaagaaat ccaagctgga gttctccatt tacccggcac cccaggtttc    600 

cacagctgta gttgagccct acaactccat cctcaccacc cacaccaccc tggagcactc    660 

tgattgtgcc ttcatggtag acaatgaggc catctatgac atctgtcgta gaaacctcga    720 

tatcgagcgc ccaacctaca ctaaccttaa ccgccttatt agccagattg tgtcctccat    780 

cactgcttcc ctgagatttg atggagccct gaatgttgac ctgacagaat tccagaccaa    840 

cctggtcccc tacccccgca tccacttccc tctggccaca tatgcccctg tcatctctgc    900 

tgagaaagcc taccatgaac agctttctgt agcagacatc accaatgctt gctttgagcc    960 

agccaaccag atggtgaaat gtgaccctgg ccatggtaaa tacatggctt gctgcctgtt   1020 

gtaccgtggt gacgtggttc ccaaagatgt caatgctgcc attgccacca tcaaaaccaa   1080 

gcgcacgatc cagtttgtgg attggtgccc cactggcttc aaggttggca tcaactacca   1140 

gcctcccact gtggtgcctg gtggagacct ggccaaggta cagagagctg tgtgcatgct   1200 

gagcaacacc acagccattg ctgaggcctg ggctcgcctg gaccacaagt ttgacctgat   1260 

gtatgccaag cgtgcctttg ttcactggta cgtgggtgag gggatggagg aaggcgagtt   1320 

ttcagaggcc cgtgaagata tggctgccct tgagaaggat tatgaggagg ttggtgtgga   1380 

ttctgttgaa ggagagggtg aggaagaagg agaggaatac taattatcca ttccttttgg   1440 

ccctgcagca tgtcatgctc ccagaatttc agcttcagct taactgacag atgttaaagc   1500 

tttctggtta gattgttttc acttggtgat catgtctttt ccatgtgtac ctgtaatatt   1560 

tttccatcat atctcaaagt aaagtcatta acatca                             1596 

 
           
             45  
             4276  
             DNA  
             Homo sapiens  
           
            45 

ctgtgaccca gaagtcttcg aattcactgg tttttcagac tctgccacgg cacatgcgac     60 

gaagagccat gagccacaac gtcaaacgcc ttcccagacg gttacaggag attgcccaga    120 

aagaggcgga gaaagccgta catcagaaaa aagaacattc aaaaaataaa tgccataaag    180 

ctcgaagatg tcacatgaac cggacgctag aatttaaccg tagacaaaag aagaacattt    240 

ggttagaaac tcacatctgg cacgccaagc ggtttcatat ggtcaagaag tggggctact    300 

gccttgggga gaggccaaca gtcaagagcc acagagcctg ctatcgagcc atgacgaacc    360 

ggtgcctcct gcaggattta tcctattact gttgtttgga gttgaaaggc aaagaggaag    420 

aaatactaaa ggcgctttct ggaatgtgta acatagacac agggctgacg tttgcagcag    480 

ttcactgctt gtctggaaag cgccaaggga gccttgtgct ttatcgggtg aataaatatc    540 

ccagagaaat gcttgggcct gttacgttta tctggaagtc ccagaggacc ccgggtgacc    600 

cttctgagag caggcagctg tggatctggc tgcatccaac ccttaaacag gatatcttag    660 

aggaaataaa agcagcgtgc cagtgtgtgg aacccatcaa atcagctgtc tgcatcgctg    720 

acccacttcc aacaccatcc caagaaaaaa gccaaactga attgcctgac gagaaaattg    780 

gcaagaaaag aaaaaggaaa gatgatggag aaaatgctaa accaattaaa aaaattatcg    840 

gtgatggaac tagagatcca tgtctaccat actcttggat ctctccaacc acaggcatta    900 

taatcagcga tttgacgatg gagatgaaca gattccggct gattgggcca ctttcccact    960 

ccatcctaac tgaagcaata aaagctgctt ctgtccacac tgtgggagag gacacagagg   1020 

agacacctca ccgctggtgg atagaaacct gtaagaaacc tgacagcgtt tcccttcatt   1080 

gcagacaaga agccattttc gagttgttgg gaggaataac atcaccagca gaaattccgg   1140 

caggtactat tctgggactg acagttgggg atcctcgaat aaatttgccc caaaagaagt   1200 

ccaaagcttt gcccaatcca gaaaaatgcc aagataatga gaaagttaga cagctgcttc   1260 

tggagggtgt gcctgtggaa tgtacgcata gctttatctg gaaccaagat atctgtaaga   1320 

gtgtcacaga gaataaaatc tcggatcagg atttaaaccg gatgaggagt gaattgctgg   1380 

tgcctgggtc acagcttatt ttaggtcccc atgaatccaa gatacctata cttttgattc   1440 

agcagccagg aaaagtgact ggtgaagatc gactaggctg gggaagtggc tgggatgtcc   1500 

tactcccaaa gggctggggc atggctttct ggattccatt tatttatcga ggtgtgagag   1560 

tcggagggtt gaaagagtct gcagtgcatt ctcagtataa gaggtcgcct aatgtcccag   1620 

gcgattttcc agactgccct gccgggatgc tgtttgcgga agagcaagct aagaatcttc   1680 

ttgaaaagta caaaagacgc cctcctgcaa aacggcccaa ctacgttaag cttggcactc   1740 

tggcaccttt ctgctgtccc tgggagcagt taactcaaga ctgggagtca agagtccagg   1800 

cttacgaaga accttctgta gcttcatctc caaatggtaa ggagagtgac ctaagaagat   1860 

ctgaggtgcc ttgtgctccc atgcctaaaa aaactcatca gccatctgat gaagtgggca   1920 

catccataga gcaccccagg gaggcagagg aggtaatgga tgcagggtgt caagaatcgg   1980 

cagggcctga gaggatcaca gaccaggagg ccagtgaaaa ccatgttgct gccacaggga   2040 

gtcacctctg cgttctcagg agtagaaaat tactgaagca actgtcagcc tggtgtgggc   2100 

ccagttctga ggatagtcgg ggaggccggc gagctcccgg cagaggccag caaggattga   2160 

ccagagaggc ttgcctgtcc atcttgggcc acttccccag ggccctggtt tgggtcagcc   2220 

tgtccctgct cagcaagggc agccccgagc ctcacaccat gatctgtgtc ccagccaagg   2280 

aggacttcct ccagctccat gaggactggc attactgtgg gccccaggaa tccaaacaca   2340 

gtgacccatt caggagcaag atcctgaaac agaaagagaa gaagaaaagg gagaagaggc   2400 

agaagccagg acgtgcctct tctgatggcc cggcggggga agagcccgtg gctgggcagg   2460 

aagctctgac tctagggctg tggtcaggcc ctctgccgcg tgtgacgttg cactgctcca   2520 

gaactctcct aggctttgtg actcagggag atttttccat ggctgttggc tgtggagaag   2580 

ccctggggtt tgttagcttg acaggcttgc tggatatgct gtccagccag cctgcagcgc   2640 

agaggggctt agtgctactg aggcctcccg cctctctgca gtatcgattt gcgaggattg   2700 

ctattgaggt gtgaatgcgt gcttgtatcc cagcagggca tagataatac gttattattg   2760 

tctgccaagt tctacatgtg gagaatctgc ttctgcttta aaatatcatg tgaaactccc   2820 

tggaaacaag aataaaaaat tatgtattat gcagatgatg aaatgtttac atcattccag   2880 

taatgtcatt gattttcatc tttccctgtc cttgctgtaa tacttttaaa ttatttggcc   2940 

aaaagctttg tattatgatc tcttggtctg tgtagttgtg gctgaaaata atgagaagct   3000 

ctacgagtta tcatcccctt tttttgttag aaacaaaggg cttgtcaggt ctatttgaaa   3060 

aacctcatag tcatgtgata agcaacaata gatgtttaat gatttcactg ttatagcaga   3120 

agacaagaga agacgcttgg cctctgtaca tgaaatatgg gctcctgatg gacctcattc   3180 

aattctgtac tgtgatttcc atgccgaaca actcaagcct taaagagaga aatcatggac   3240 

aactgatttc tgcctgtttt caggcaggca cagtttatgg cgtcagtgct aggctggaat   3300 

tagaaagtgg gggtctatga cgtggacttc ctgactcttt gatctctttg ttgttgacca   3360 

acacttgatc ctactagtta cttaattttt ttaagtaaaa aattattatt attttgtttc   3420 

tgcaaagatt ttctcaaagc catagaggag catttctcag aatatgttct atgatatgtg   3480 

tcacctaaaa aagtaagaga ttccaaggtc aggttgatat ggaaactcta ggttaaataa   3540 

agttaagcat ttctttatga aagaacttct ggaaacttcc atgtgataat gtgcattgcg   3600 

gatctctagg aaggaaatga tagtgtatag tattttctaa atacttgtga ttcctaaagt   3660 

tctcttacaa ggagcccttt gtaggaccag tgttcttagt agcgcgcttt gggcagtgtg   3720 

gctgtgtagt gcatagctac ctctgcaagg tgataactaa gccggcaagc tgcctttcaa   3780 

cactcatgca gtcacgttgt ccacctgaga ttctcaacag ggtataaaag gaaggtctca   3840 

tcttgcctca caggaagagt gggctcagtg tggctttttt ccaactatgg agaaactcag   3900 

tgctcatcta ctttaagttt ccacatatgg cttgctcata gccttggtcc ttacctttcc   3960 

tgccataact ttctagaaga gcttaatggg atttttttct aaaaaatgta aatatgcagt   4020 

taggcattat tttatgtaaa tgcattgggt ttttactgta gcatttggca ctaaatggct   4080 

ttgggggtga tgaggtgggg aaggatacag caggtggtac agtagtcagg aagtacctgc   4140 

caccaatgag atgtctgatg ctttgcctct taccatgcct ctgaatgtct ttggatccaa   4200 

cccagatgag actgaaaaaa aaaaaacagt gtaactaagt ggcatctgta aacagaataa   4260 

atgaaaatgt cacctg                                                   4276 

 
           
             46  
             10  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            46 

gtagcccagc                                                            10 

 
           
             47  
             10  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            47 

gccacccaga                                                            10 

 
           
             48  
             14  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            48 

acgaagaaga agag                                                       14 

 
           
             49  
             10  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            49 

aggggcacca                                                            10 

 
           
             50  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            50 

aatgaggggg acaaatggga agc                                             23 

 
           
             51  
             25  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            51 

ggagagccct tcctcagaca tgaag                                           25 

 
           
             52  
             25  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            52 

tgacaaaatg gtgacaggta gctgg                                           25 

 
           
             53  
             24  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            53 

aagtccacac ctcctcagac agcc                                            24 

 
           
             54  
             22  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            54 

cccagacacc caaacagccg tg                                              22 

 
           
             55  
             20  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            55 

tggagcagcc gtgtgtgctg                                                 20 

 
           
             56  
             16  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            56 

aagctttttt tttttg                                                     16 

 
           
             57  
             16  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            57 

aagctttttt ttttta                                                     16 

 
           
             58  
             16  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            58 

aagctttttt tttttc                                                     16 

 
           
             59  
             13  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            59 

aagcttgatt gcc                                                        13 

 
           
             60  
             13  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            60 

aagcttcgac tgt                                                        13 

 
           
             61  
             13  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            61 

aagctttggt cag                                                        13 

 
           
             62  
             13  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            62 

aagcttctca acg                                                        13 

 
           
             63  
             17  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            63 

attttttttt tttttta                                                    17 

 
           
             64  
             17  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            64 

gttttttttt ttttttg                                                    17 

 
           
             65  
             14  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            65 

tttttttttt tttv                                                       14 

 
           
             66  
             11  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            66 

ggtgcctttg g                                                          11 

 
           
             67  
             10  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            67 

gcaccagggg                                                            10 

 
           
             68  
             15  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            68 

tttttttttt ttttt                                                      15 

 
           
             69  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            69 

cacgtcttgg tgccttttgt gtg                                             23 

 
           
             70  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            70 

gaagctcagc tcagccctct tcc                                             23 

 
           
             71  
             25  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            71 

ccagggagac caaaagcctt catac                                           25 

 
           
             72  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            72 

cacaggggag gtgatagcat tgc                                             23 

 
           
             73  
             25  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            73 

gtgcttttca aagatgctgc tagtg                                           25 

 
           
             74  
             22  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            74 

gctcaatcca cccacaaaaa cc                                              22 

 
           
             75  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            75 

tcctctcact gccttggtgg atg                                             23 

 
           
             76  
             24  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            76 

cacagcaagt cacacattgg accc                                            24 

 
           
             77  
             22  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            77 

ccaaagacat tcagaggcat gg                                              22 

 
           
             78  
             22  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            78 

gaggtgggga aggatacagc ag                                              22 

 
           
             79  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            79 

gaaaagggtt ggggagaagc ctc                                             23 

 
           
             80  
             25  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            80 

tctctagcgt cctccatctc actgg                                           25 

 
           
             81  
             23  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            81 

acaactgcat cctcaccacc cac                                             23 

 
           
             82  
             25  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            82 

ggacacaatc tggctaataa ggcgg                                           25 

 
           
             83  
             16  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            83 

aagctttttt tttttg                                                     16 

 
           
             84  
             16  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            84 

aagctttttt tttttc                                                     16 

 
           
             85  
             16  
             DNA  
             Artificial Sequence  
             
               Description of Artificial Sequence Primer  
             
           
            85 

aagctttttt ttttta                                                     16