WO2002006294A2 - Haplotypes of the mmp13 gene - Google Patents

Abstract

Novel genetic variants of the Matrix Metalloproteinase 13 (Collagenase 3) (MMP13) gene are described. Various genotypes, haplotypes, and haplotype pairs that exist in the general United States population are disclosed for the MMP13 gene. Compositions and methods for haplotyping and/or genotyping the MMP13 gene in an individual are also disclosed. Polynucleotides defined by the haplotypes disclosed herein are also described.

Description

HAPLOTYPES OF THE MMP13 GENE

RELATED APPLICATIONS

^•This application claims the benefit ofU.S. Provisional Application Serial No. 60/217,950 filed July 13, 2000, and PCT Application Serial No. PCT US00/22693 filed August 17, 2000.

FIELD OF THE INVENTION

This invention relates to variation in genes that encode pharmaceutically-important proteins. In particular, this invention provides genetic variants of the human matrix metalloproteinase 13 (collagenase 3) (MMP13) gene and methods for identifying which variant(s) of this gene is/are possessed by an individual.

BACKGROUND OF THE INVENTION

Current methods for identifying pharmaceuticals to treat disease often start by identifying, cloning, and expressing an important target protein related to the disease. A determination of whether an agonist or antagonist is needed to produce an effect that may benefit a patient with the disease is then made. Then, vast numbers of compounds are screened against the target protein to find new potential drugs. The desired outcome of this process is a lead compound that is specific for the target, thereby reducing the incidence of the undesired side effects usually caused by activity at non-intended targets. The lead compound identified in this screening process then undergoes further in vitro and in vivo testing to determine its absorption, disposition, metabolism and toxicological profiles. Typically, this testing involves use of cell lines and animal models with limited, if any, genetic diversity.

What this approach fails to consider, however, is that natural genetic variability exists between individuals in any and every population with respect to pharmaceutically-important proteins, including the protein targets of candidate drugs, the enzymes that metabolize these drugs and^" the proteins whose activity is modulated by such drug targets. Subtle alteration(s) in the primary nucleotide sequence of a gene encoding a pharmaceutically-important protein may be manifested as significant variation in expression, structure and/or function of the protein. Such alterations may explain the relatively high degree of uncertainty inherent in the treatment of individuals with a drug whose design is based upon a single representative example of the target or enzyme(s) involved in metabolizing the drug. For example, it is well-established that some drugs frequently have lower efficacy in some individuals than others, which means such individuals and their physicians must weigh the possible benefit of a larger dosage against a greater risk of side effects. Also, there is significant variation in how well people metabolize drugs and other exogenous chemicals, resulting in substantial interindividual variation in the toxicity and/or efficacy of such exogenous substances (Evans et al., 1999, Science 286:487-491). This variability in efficacy or toxicity of a drug in genetically-diverse patients makes many drugs ineffective or even dangerous in certain groups of the population, leading to the failure of such drugs in clinical trials or their early withdrawal from the market even though they could be highly beneficial for other groups in the population. This problem significantly increases the time and cost of drug discovery and development, which is a matter of great public concern.

It is well-recognized by pharmaceutical scientists that considering the impact of the genetic variability of pharmaceutically-important proteins in the early phases of drug discovery and development is likely to reduce the failure rate of candidate and approved drugs (Marshall A 1997 Nature Biotech 15:1249-52; Kleyn PW et al. 1998 Science 281: 1820-21; Kola 1 1999 Curr Opin Biotech 10:589-92; Hill AVS et al. 1999 in Evolution in Health and Disease Stearns SS (Ed.) Oxford University Press, New York, pp 62-76; Meyer U.A. 1999 in Evolution in Health and Disease Steams SS (Ed.) Oxford University Press, New York, pp 41-49; Kalow W et al. 1999 Clin. Pharm. Therap.

66:445-7; Marshall, E 1999 Science 284:406-7; Judson R et al. 2000 Pharmacogenomics 1:1-12; Roses AD 2000 Nature 405:857-65). However, in practice this has been difficult to do, in large part because of the time and cost required for discovering the amount of genetic variation that exists in the population (Chakravarti A 1998 Nature Genet 19:216-7; Wang DG et al 1998 Science 280:1077-82; Chakravarti A 1999 Nat Genet 21:56-60 (suppl); Stephens JC 1999 Mol. Diagnosis 4:309-317; Kwok PY and Gu S 1999 Mol. Med. Today 5:538-43; Davidson S 2000 Nature Biotech 18:1134-5).

The standard for measuring genetic variation among individuals is the haplotype, which is the ordered combination of polymorphisms in the sequence of each form of a gene that exists in the population. Because haplotypes represent the variation across each form of a gene, they provide a more accurate and reliable measurement of genetic variation than individual polymoφhisms. For example, while specific variations in gene sequences have been associated with a particular phenotype such as disease susceptibility (Roses AD supra; Ulbrecht M et al. 2000 Am JRespir Crit Care Med 161: 469- 74) and drug response (Wolfe CR et al. 20005/147320:987-90; Dahl BS 1997 Ada Psychiatr Scand 96 (Suppl 391): 14-21), in many other cases an individual polymoφhism may be found in a variety of genomic backgrounds, i.e., different haplotypes, and therefore shows no definitive coupling between the polymoφhism and the causative site for the phenotype (Clark AG et al. 1998 Am JHum Genet 63:595- 612; Ulbrecht M et al. 2000 supra; Drysdale et al. 2000 PNAS 97: 10483-10488). Thus, there is an unmet need in the pharmaceutical industry for information on what haplotypes exist in the population for pharmaceutically-important genes. Such haplotype information would be useful in improving the efficiency and output of several steps in the drug discovery and development process, including target validation, identifying lead compounds, and early phase clinical trials (Marshall et al., supra).

One pharmaceutically-important gene for the treatment of cancer and arthritis is the matrix metalloproteinase 13 (collagenase 3) (MMP13) gene or its encoded product. MMP13, also known as collagenase-3 (CLG3), is a metalloproteinase. Matrix metalloproteinases (MMPs) are a family of secreted or transmembrane proteins that can degrade all the proteins of the extracellular matrix and have been implicated in many abnormal physiological conditions, including arthritis and cancer metastasis (Sun et al. JBiol Chem 2000 Apr 14;275(15):11327-32). MMP13 expression has been found in breast tumors by northern blot analysis (Freije et al. JBiol Chem 1994 Jun 17;269(24):16766-73). Freije et al. (1994, supra) suggest a possible role for MMP13 in the tumoral process because of the increased expression of MMP13 in breast carcinomas and the absence of detectable expression in normal tissues. Mitchell et al. (J Clin Invest 1996 Feb l;97(3):761-8) found that MMP13 is expressed by chondrocytes in human osteoarthritic cartilage. They further conclude that the expression of MMP13 in osteoarthritic cartilage and its activity against type II collagen suggest that the enzyme plays a significant role in cartilage collagen degradation, and may consequently form part of a complex target for proposed therapeutic interventions based on collagenase inhibition.

The matrix metalloproteinase 13 (collagenase 3) gene is located on chromosome 1 lq22.3 and contains 10 exons that encode a 471 amino acid protein. Although a complete reference sequence for the MMP13 gene has not been published, a reference sequence for this gene comprises the partial sequences shown in Figure 1 (nucleotides 31767-43261 in GenBank Accession No: AP001018.1; SEQ ID NO:l), which includes the coding sequence for exons 1-7, Figure 2 (nucleotides 152006-155369 in GenBank Accession No: AP001018.1; SEQ ID NO:2), which includes the coding sequence for exons 8 and 9, and Figure 3, which includes a reference sequence for exon 10 (the reverse complment of nucleotides 69826-62706 in GenBank Accession No. AP001019.1; SEQ ID NO:3). These reference sequences for the MMP13 gene are based on Genaissance Reference No. 958944. Reference sequences for the complete coding sequence (GenBank Accession No: NM_002427.1 ; SEQ ID NO:4) and protein (GenBank Accessio No: NP_002418.1; SEQ ID NO:5) are shown in Figures 4 and 5, respectively. Because of the potential for variation in the MMP13 gene to affect the expression and function of the encoded protein, it would be useful to know whether polymoφhisms exist in the MMP13 gene, as well as how such polymoφhisms are combined in different copies of the gene. Such information could be applied for studying the biological function of MMP13 as well as in identifying drugs targeting this protein for the treatment of disorders related to its abnormal expression or function.

SUMMARY OF THE INVENTION

Accordingly, the inventors herein have discovered 18 novel polymoφhic sites in the MMP13 gene. These polymoφhic sites (PS) correspond to the following nucleotide positions in Figures 1, 2, and 3: 3915 (PS1), 4437 (PS2), 5008 (PS3), 5037 (PS4), 5102 (PS5), 5290 (PS6), 5363 (PS7), 5628 (PS8), 7721 (PS9), 10537 (PS10), 11663 (PSll), 11752 (PS12), 14059 (PS13), 14072 (PS14), 18295 (PS15), 18407 (PS16), 18609 (PS17) and 18866 (PS18). The polymoφhisms at these sites are adenine or guanine at PS1, guanine or thymine at PS2, cytosine or thymine at PS3, guanine or adenine at PS4, cytosine or guanine at PS5, cytosine or guanine at PS6, cytosine or adenine at PS7, guanine or adenine at PS8, guanine or adenine at PS9, cytosine or thymine at PS10, thymine or cytosine at PS11, adenine or guanine at PS12, cytosine or thymine at PS13, guanine or cytosine at PS14, thymine or cytosine at

PS15, adenine or guanine at PS16, adenine or guanine at PS17 and adenine or guanine at PS18. In addition, the inventors have determined the identity of the alleles at these sites in a human reference population of 79 unrelated individuals self-identified as belonging to one of four major population groups: African descent, Asian, Caucasian and Hispanic/Latino. From this information, the inventors deduced a set of haplotypes and haplotype pairs for PS 1 -PS 18 in the MMP13 gene, which are shown below in Tables 5 and 4, respectively. Each of these MMP13 haplotypes defines a naturally-occurring isoform (also referred to herein as an "isogene") of the MMP13 gene that exists in the human population. The frequency with which each haplotype and haplotype pair occurs within the total reference population and within each of the four major population groups included in the reference population was also determined.

Thus, in one embodiment, the invention provides a method, composition and kit for genotyping the MMP13 gene in an individual. The genotyping method comprises identifying the nucleotide pair that is present at one or more polymoφhic sites selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PS11, PS12, PS13, PS14, PS15, PS16, PS17 and PS18 in both copies of the MMP13 gene from the individual. A genotyping composition of the invention comprises an oligonucleotide probe or primer which is designed to specifically hybridize to a target region containing, or adjacent to, one of these novel MMP13 polymoφhic sites. A genotyping kit of the invention comprises a set of oligonucleotides designed to genotype each of these novel MMP13 polymoφhic sites. The genotyping method, composition, and kit are useful in determining whether an individual has one of the haplotypes in Table 5 below or has one of the haplotype pairs in Table 4 below. The invention also provides a method for haplotyping the MMP13 gene in an individual. In one embodiment, the haplotyping method comprises determining, for one copy of the MMP13 gene, the identity of the nucleotide at one or more polymoφhic sites selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS13, PS14, PS15, PS16, PS17 and PS18. In another embodiment, the haplotyping method comprises determining whether one copy of the individual's MMP13 gene is defined by one of the MMP13 haplotypes shown in Table 5, below, or a sub-haplotype thereof. In a preferred embodiment, the haplotyping method comprises determining whether both copies of the individual's MMP13 gene are defined by one of the MMP1 haplotype pairs shown in Table 4 below, or a sub-haplotype pair thereof. The method for establishing the MMP13 haplotype or haplotype pair of an individual is useful for improving the efficiency and reliability of several steps in the discovery and development of drugs for treating diseases associated with MMP13 activity, e.g., cancer and arthritis.

For example, the haplotyping method can be used by the pharmaceutical research scientist to validate MMP13 as a candidate target for treating a specific condition or disease predicted to be associated with MMP13 activity. Determining for a particular population the frequency of one or more of the individual MMP13 haplotypes or haplotype pairs described herein will facilitate a decision on whether to pursue MMP13 as a target for treating the specific disease of interest. In particular, if variable MMP13 activity is associated with the disease, then one or more MMP13 haplotypes or haplotype pairs will be found at a higher frequency in disease cohorts than in appropriately genetically matched controls. Conversely, if each of the observed MMP13 haplotypes are of similar frequencies in the disease and control groups, then it may be inferred that variable MMP13 activity has little, if any, involvement with that disease. In either case, the pharmaceutical research scientist can, without a priori knowledge as to the phenotypic effect of any MMP 13 haplotype or haplotype pair, apply the information derived from detecting MMP 13 haplotypes in an individual to decide whether modulating MMP 13 activity would be useful in treating the disease.

The claimed invention is also useful in screening for compounds targeting MMP 13 to treat a specific condition or disease predicted to be associated with MMP 13 activity. For example, detecting which of the MMP 13 haplotypes or haplotype pairs disclosed herein are present in individual members of a population with the specific disease of interest enables the pharmaceutical scientist to screen for a compound(s) that displays the highest desired agonist or antagonist activity for each of the most frequent MMP 13 isoforms present in the disease population. Thus, without requiring any a priori knowledge of the phenotypic effect of any particular MMP 13 haplotype or haplotype pair, the claimed haplotyping method provides the scientist with a tool to identify lead compounds that are more likely to show efficacy in clinical trials.

The method for haplotyping the MMP 13 gene in an individual is also useful in the design of clinical trials of candidate drugs for treating a specific condition or disease predicted to be associated with MMP13 activity. For example, instead of randomly assigning patients with the disease of interest to the treatment or control group as is typically done now, determining which of the MMP 13 haplotype(s) disclosed herein are present in individual patients enables the pharmaceutical scientist to distribute MMP 13 haplotypes and/or haplotype pairs evenly to treatment and control groups, thereby reducing the potential for bias in the results that could be introduced by a larger frequency of a MMP 13 haplotype or haplotype pair that had a previously unknown association with response to the drug being studied in the trial. Thus, by practicing the claimed invention, the scientist can more confidently rely on the information learned from the trial, without first determining the phenotypic effect of any MMP 13 haplotype or haplotype pair.

In another embodiment, the invention provides a method for identifying an association between a trait and a MMP 13 genotype, haplotype, or haplotype pair for one or more of the novel polymoφhic sites described herein. The method comprises comparing the frequency of the MMP 13 genotype, haplotype, or haplotype pair in a population exhibiting the trait with the frequency of the MMP 13 genotype or haplotype in a reference population. A higher frequency of the MMP 13 genotype, haplotype, or haplotype pair in the trait population than in the reference population indicates the trait is associated with the MMP 13 genotype, haplotype, or haplotype pair. In preferred embodiments, the trait is susceptibility to a disease, severity of a disease, the staging of a disease or response to a drug. In a particularly preferred embodiment, the MMP 13 haplotype is selected from the haplotypes shown in Table 5, or a sub-haplotype thereof. Such methods have applicability in developing diagnostic tests and therapeutic treatments for cancer and arthritis.

In yet another embodiment, the invention provides an isolated polynucleotide comprising a nucleotide sequence which is a polymoφhic variant of a reference sequence for the MMP 13 gene or a fragment thereof. The reference sequence comprises the partial sequences shown in Figures 1, 2, and 3 and the polymoφhic variant comprises at least one polymoφhism selected from the group consisting of guanine at PS1, thymine at PS2, thymine at PS3, adenine at PS4, guanine at PS5, guanine at PS6, adenine at PS7, adenine at PS8, adenine at PS9, thymine at PS 10, cytosme at PS 11, guanine at PS 12, thymine at PS13, cytosine at PS14, cytosine at PS15, guanine at PS16, guanine at PS17 and guanine at PS18. A particularly preferred polymoφhic variant is an isogene of the MMP 13 gene. A MMP 13 isogene of the invention comprises adenine or guanine at PS1, guanine or thymine at PS2, cytosine or thymine at PS3, guanine or adenine at PS4, cytosine or guanine at PS5, cytosine or guanine at PS6, cytosine or adenine at PS7, guanine or adenine at PS8, guanine or adenine at PS9, cytosine or thymine at PS10, thymine or cytosine at PS11, adenine or guanine at PS12, cytosine or thymine at PS13, guanine or cytosine at PS14, thymine or cytosine at PS15, adenine or guanine at PS16, adenme or guanine at PS17 and adenine or guanine at PS18. The invention also provides a collection of MMP13 isogenes, referred to herein as a MMP 13 genome anthology.

In another embodiment, the invention provides a polynucleotide comprising a polymoφhic variant of a reference sequence for a MMP 13 cDNA or a fragment thereof. The reference sequence comprises SEQ ID NO:4 (Fig.4) and the polymoφhic cDNA comprises at least one polymoφhism selected from the group consisting of thymine at a position corresponding to nucleotide 326, cytosine at a position corresponding to nucleotide 1080 and guanine at a position corresponding to nucleotide 1169. A particularly preferred polymoφhic cDNA variant comprises the coding sequence of a MMP 13 isogene defined by haplotypes 8c, 10c, and 12c. Polynucleotides complementary to these MMP 13 genomic and cDNA variants are also provided by the invention. It is believed that polymoφhic variants of the MMP 13 gene will be useful in studying the expression and function of MMP13, and in expressing MMP13 protein for use in screening for candidate drugs to treat diseases related to MMP 13 activity.

In other embodiments, the invention provides a recombinant expression vector comprising one of the polymoφhic genomic variants operably linked to expression regulatory elements as well as a recombinant host cell transformed or transfected with the expression vector. The recombinant vector and host cell may be used to express MMP 13 for protein structure analysis and drug binding studies. In yet another embodiment, the invention provides a polypeptide comprising a polymoφhic variant of a reference amino acid sequence for the MMP 13 protein. The reference amino acid sequence comprises SEQ ID NO:5 (Fig.5) and the polymoφhic variant comprises at least one variant amino acid selected from the group consisting of leucine at a position corresponding to amino acid position 109 and glycine at a position corresponding to amino acid position 390. A polymoφhic variant of MMP13 is useful in studying the effect of the variation on the biological activity of MMP 13 as well as on the binding affinity of candidate drugs targeting MMP 13 for the treatment of cancer and arthritis.

The present invention also provides antibodies that recognize and bind to the above polymoφhic MMP 13 protein variant. Such antibodies can be utilized in a variety of diagnostic and prognostic formats and therapeutic methods.

The present invention also provides nonhuman transgenic animals comprising one of the MMP 13 polymoφhic genomic variants described herein and methods for producing such animals. The transgenic animals are useful for studying expression of the MMP 13 isogenes in vivo, for in vivo screening and testing of drugs targeted against MMP 13 protein, and for testing the efficacy of therapeutic agents and compounds for cancer and arthritis in a biological system.

The present invention also provides a computer system for storing and displaying polymoφhism data determined for the MMP 13 gene. The computer system comprises a computer processing unit; a display; and a database containing the polymoφhism data. The polymoφhism data includes the polymoφhisms, the genotypes and the haplotypes identified for the MMP 13 gene in a reference population. In a preferred embodiment, the computer system is capable of producing a display showing MMP 13 haplotypes organized according to their evolutionary relationships.

BRIEF DESCRIPTION OF THE DRAWINGS

Figure 1 illustrates a partial reference sequence for the MMP 13 gene (contiguous lines; SEQ ID NO:l; corresponding to nucleotides 31767-43261 of GenBank Accession No. AP001018.1), with the start and stop positions of each region of coding sequence indicated with a bracket ([ or ]) and the numerical position below the sequence and the polymoφhic site(s) and polymoφhism(s) identified by Applicants in a reference population indicated by the variant nucleotide positioned below the polymoφhic site in the sequence. SEQ ID NO:l is equivalent to Figure 1, with the two alternative allelic variants of each polymoφhic site indicated by the appropriate nucleotide symbol (R= G or A, Y= T or C, M= A or C, K= G or T, S= G or C, and W= A or T; WIPO standard ST.25).

Figure 2 illustrates a partial reference sequence for the MMP 13 gene (contiguous lines; SEQ ID NO:2; corresponding to nucleotides 152006-155369 of GenBank Accession No. AP001018.1), with the start and stop positions of each region of coding sequence indicated with a bracket ([ or ]) and the numerical position below the sequence and the polymoφhic site(s) and polymoφhism(s) identified by Applicants in a reference population indicated by the variant nucleotide positioned below the polymoφhic site in the sequence. SEQ ID NO:2 is equivalent to Figure 2, with the two alternative allelic variants of each polymoφhic site indicated by the appropriate nucleotide symbol (R= G or A, Y= T or C, M= A or C, K= G or T, S= G or C, and W= A or T; WIPO standard ST.25). Figure 3 illustrates a partial reference sequence for the MMP 13 gene (contiguous lines; SEQ ID

NO:3; corresponding to the reverse complement of nucleotides 69826-62706 of GenBank Accession No. AP001018.1), with the start and stop positions of each region of coding sequence indicated with a bracket ([ or ]) and the numerical position below the sequence and the polymoφhic site(s) and polymoφhism(s) identified by Applicants in a reference population indicated by the variant nucleotide positioned below the polymoφhic site in the sequence. SEQ ID NO:3 is equivalent to Figure 3, with the two alternative allelic variants of each polymoφhic site indicated by the appropriate nucleotide symbol (R= G or A, Y= T or C, M= A or C, K= G or T, S= G or C, and W= A or T; WIPO standard ST.25).

SEQ ID NO:96 is a modified version of SEQ ID NOS: 1-3 that shows the context sequence of each polymoφhic site, PS 1 -PS 18, in a uniform format to facilitate electronic searching. For each polymoφhic site, SEQ ID NO:96 contains a block of 60 bases of the nucleotide sequence encompassing the centrally-located polymoφhic site at the 30^th position, followed by 60 bases of unspecified sequence to represent that each PS is separated by genomic sequence whose composition is defined elsewhere herein.

Figure 4 illustrates a reference sequence for the MMP 13 coding sequence (contiguous lines; SEQ ID NO:4), with the polymoφhic site(s) and polymoφhism(s) identified by Applicants in a reference population indicated by the variant nucleotide positioned below the polymoφhic site in the sequence.

Figure 5 illustrates a reference sequence for the MMP 13 protein (contiguous lines; SEQ ID NO:5), with the variant amino acid(s) caused by the polymoφhism(s) of Figure 4 positioned below the polymoφhic site in the sequence.

DESCRIPTION OF THE PREFERRED EMBODIMENTS

The present invention is based on the discovery of novel variants of the MMP 13 gene. As described in more detail below, the inventors herein discovered 15 isogenes of the MMP 13 gene by characterizing the MMP 13 gene found in genomic DNAs isolated from an Index Repository that contains immortalized cell lines from one chimpanzee and 93 human individuals. The human individuals included a reference population of 79 unrelated individuals self-identified as belonging to one of four major population groups: Caucasian (21 individuals), African descent (20 individuals), Asian (20 individuals), or Hispanic/Latino (18 individuals). To the extent possible, the members of this reference population were organized into population subgroups by their self-identified ethnogeographic origin as shown in Table 1 below.

In addition, the Index Repository contains three unrelated indigenous American Indians (one from each of North, Central and South America), one three-generation Caucasian family (from the CEPH Utah cohort) and one two-generation African- American family.

The MMP 13 isogenes present in the human reference population are defined by haplotypes for 18 polymoφhic sites in the MMP 13 gene, all of which are believed to be novel. The novel MMP 13 polymoφhic sites identified by the inventors are referred to as PS1-PS18 to designate the order in which they are located in the gene (see Table 3 below). Using the genotypes identified in the Index Repository for PS 1 -PS 18 and the methodology described in the Examples below, the inventors herein also determined the pair of haplotypes for the MMP 13 gene present in individual human members of this repository. The human genotypes and haplotypes found in the repository for the MMP 13 gene include those shown in Tables 4 and 5, respectively. The polymoφhism and haplotype data disclosed herein are useful for validating whether MMP 13 is a suitable target for drugs to treat cancer and arthritis, screening for such drugs and reducing bias in clinical trials of such drugs.

In the context of this disclosure, the following terms shall be defined as follows unless otherwise indicated:

Allele - A particular form of a genetic locus, distinguished from other forms by its particular nucleotide sequence. Candidate Gene - A gene which is hypothesized to be responsible for a disease, condition, or the response to a treatment, or to be correlated with one of these.

Gene - A segment of DNA that contains all the information for the regulated biosynthesis of an RNA product, including promoters, exons, introns, and other untranslated regions that control expression. Genotype - An unphased 5 ' to 3 ' sequence of nucleotide pair(s) found at one or more polymoφhic sites in a locus on a pair of homologous chromosomes in an individual. As used herein, genotype includes a full-genotype and/or a sub-genotype as described below.

Full-genotype - The unphased 5' to 3' sequence of nucleotide pairs found at all polymoφhic sites examined herein in a locus on a pair of homologous chromosomes in a single individual. Sub-genotype - The unphased 5' to 3' sequence of nucleotides seen at a subset of the polymoφhic sites examined herein in a locus on a pair of homologous chromosomes in a single individual.

Genotyping - A process for determining a genotype of an individual.

Haplotype - A 5' to 3' sequence of nucleotides found at one or more polymoφhic sites in a locus on a single chromosome from a single individual. As used herein, haplotype includes a full- haplotype and/or a sub-haplotype as described below.

Full-haplotype - The 5' to 3' sequence of nucleotides found at all polymoφhic sites examined herein in a locus on a single chromosome from a single individual.

Sub-haplotype - The 5' to 3' sequence of nucleotides seen at a subset of the polymoφhic sites examined herein in a locus on a single chromosome from a single individual.

Haplotype pair - The two haplotypes found for a locus in a single individual.

Haplotyping - A process for determining one or more haplotypes in an individual and includes use of family pedigrees, molecular techniques and or statistical inference.

Haplotype data - Information concerning one or more of the following for a specific gene: a listing of the haplotype pairs in each individual in a population; a listing of the different haplotypes in a population; frequency of each haplotype in that or other populations, and any known associations between one or more haplotypes and a trait.

Isoform - A particular form of a gene, mRNA, cDNA or the protein encoded thereby, distinguished from other forms by its particular sequence and/or structure. Isogene - One of the isoforms of a gene found in a population. An isogene contains all of the polymoφhisms present in the particular isoform of the gene.

Isolated - As applied to a biological molecule such as RNA, DNA, oligonucleotide, or protein, isolated means the molecule is substantially free of other biological molecules such as nucleic acids, proteins, lipids, carbohydrates, or other material such as cellular debris and growth media. Generally, the term "isolated" is not intended to refer to a complete absence of such material or to absence of water, buffers, or salts, unless they are present in amounts that substantially interfere with the methods of the present invention. Locus - A location on a chromosome or DNA molecule corresponding to a gene or a physical or phenotypic feature.

Naturally-occurring - A term used to designate that the object it is applied to, e.g., naturally- occurring polynucleotide or polypeptide, can be isolated from a source in nature and which has not been intentionally modified by man.

Nucleotide pair - The nucleotides found at a polymoφhic site on the two copies of a chromosome from an individual.

Phased - As applied to a sequence of nucleotide pairs for two or more polymoφhic sites in a locus, phased means the combination of nucleotides present at those polymoφhic sites on a single copy of the locus is known.

Polymorphic site (PS) - A position within a locus at which at least two alternative sequences are found in a population, the most frequent of which has a frequency of no more than 99%.

Polymorphic variant - A gene, mRNA, cDNA, polypeptide or peptide whose nucleotide or amino acid sequence varies from a reference sequence due to the presence of a polymoφhism in the gene.

Polymorphism - The sequence variation observed in an individual at a polymoφhic site. Polymoφhisms include nucleotide substitutions, insertions, deletions and microsatellites and may, but need not, result in detectable differences in gene expression or protein function.

Polymorphism data - Information concerning one or more of the following for a specific gene: location of polymoφhic sites; sequence variation at those sites; frequency of polymoφhisms in one or more populations; the different genotypes and/or haplotypes determined for the gene; frequency of one or more of these genotypes and/or haplotypes in one or more populations; any known association(s) between a trait and a genotype or a haplotype for the gene.

Polymorphism Database - A collection of polymoφhism data arranged in a systematic or methodical way and capable of being individually accessed by electronic or other means.

Polynucleotide - A nucleic acid molecule comprised of single-stranded RNA or DNA or comprised of complementary, double-stranded DNA.

Population Group - A group of individuals sharing a common ethnogeographic origin.

Reference Population - A group of subjects or individuals who are predicted to be representative of the genetic variation found in the general population. Typically, the reference population represents the genetic variation in the population at a certainty level of at least 85%, preferably at least 90%, more preferably at least 95% and even more preferably at least 99%.

Single Nucleotide Polymorphism (SNP) - Typically, the specific pair of nucleotides observed at a single polymoφhic site. In rare cases, three or four nucleotides may be found. Subject - A human individual whose genotypes or haplotypes or response to treatment or disease state are to be determined.

Treatment - A stimulus administered internally or externally to a subject. Unphased - As applied to a sequence of nucleotide pairs for two or more polymoφhic sites in a locus, unphased means the combination of nucleotides present at those polymoφhic sites on a single copy of the locus is not known.

As discussed above, information on the identity of genotypes and haplotypes for the MMP 13 gene of any particular individual as well as the frequency of such genotypes and haplotypes in any particular population of individuals is expected to be useful for a variety of drug discovery and development applications. Thus, the invention also provides compositions and methods for detecting the novel MMP 13 polymoφhisms and haplotypes identified herein.

The compositions comprise at least one MMP 13 genotyping oligonucleotide. In one embodiment, a MMP 13 genotyping oligonucleotide is a probe or primer capable of hybridizing to a target region that is located close to, or that contains, one of the novel polymoφhic sites described herein. As used herein, the term "oligonucleotide" refers to a polynucleotide molecule having less than about 100 nucleotides. A preferred oligonucleotide of the invention is 10 to 35 nucleotides long. More preferably, the oligonucleotide is between 15 and 30, and most preferably, between 20 and 25 nucleotides in length. The exact length of the oligonucleotide will depend on many factors that are routinely considered and practiced by the skilled artisan. The oligonucleotide may be comprised of any phosphorylation state of ribonucleotides, deoxyribonucleotides, and acyclic nucleotide derivatives, and other functionally equivalent derivatives. Alternatively, oligonucleotides may have a phosphate-free backbone, which may be comprised of linkages such as carboxymethyl, acetamidate, carbamate, polyamide (peptide nucleic acid (PNA)) and the like (Varma, R. in Molecular Biology and

Biotechnology, A Comprehensive Desk Reference, Ed. R. Meyers, VCH Publishers, Inc. (1995), pages 617-620). Oligonucleotides of the invention may be prepared by chemical synthesis using any suitable methodology known in the art, or may be derived from a biological sample, for example, by restriction digestion. The oligonucleotides may be labeled, according to any technique known in the art, including use of radiolabels, fluorescent labels, enzymatic labels, proteins, haptens, antibodies, sequence tags and the like.

Genotyping oligonucleotides of the invention must be capable of specifically hybridizing to a target region of a MMP13 polynucleotide, i.e., a MMP13 isogene. As used herein, specific hybridization means the oligonucleotide forms an anti-parallel double-stranded structure with the target region under certain hybridizing conditions, while failing to form such a structure when incubated with a non-target region or a non-MMP13 polynucleotide under the same hybridizing conditions. Preferably, the oligonucleotide specifically hybridizes to the target region under conventional high stringency conditions. The skilled artisan can readily design and test oligonucleotide probes and primers suitable for detecting polymoφhisms in the MMP 13 gene using the polymoφhism information provided herein in conjunction with the known sequence information for the MMP 13 gene and routine techniques.

A nucleic acid molecule such as an oligonucleotide or polynucleotide is said to be a "perfect" or

"complete" complement of another nucleic acid molecule if every nucleotide of one of the molecules is complementary to the nucleotide at the corresponding position of the other molecule. A nucleic acid molecule is "substantially complementary" to another molecule if it hybridizes to that molecule with sufficient stability to remain in a duplex form under conventional low-stringency conditions. Conventional hybridization conditions are described, for example, by Sambrook J. et al., in Molecular Cloning, A Laboratory Manual, 2^nd Edition, Cold Spring Harbor Press, Cold Spring Harbor, NY (1989) and by Haymes, B.D. et al. in Nucleic Acid Hybridization, A Practical Approach, IRL Press, Washington, D.C. (1985). While perfectly complementary oligonucleotides are preferred for detecting polymoφhisms, departures from complete complementarity are contemplated where such departures do not prevent the molecule from specifically hybridizing to the target region. For example, an oligonucleotide primer may have a non-complementary fragment at its 5' end, with the remainder of the primer being complementary to the target region. Alternatively, non-complementary nucleotides may be interspersed into the oligonucleotide probe or primer as long as the resulting probe or primer is still capable of specifically hybridizing to the target region.

Preferred genotyping oligonucleotides of the invention are allele-specific oligonucleotides. As used herein, the term allele-specific oligonucleotide (ASO) means an oligonucleotide that is able, under sufficiently stringent conditions, to hybridize specifically to one allele of a gene, or other locus, at a target region containing a polymoφhic site while not hybridizing to the corresponding region in another allele(s). As understood by the skilled artisan, allele-specificity will depend upon a variety of readily optimized stringency conditions, including salt and formamide concentrations, as well as temperatures for both the hybridization and washing steps. Examples of hybridization and washing conditions typically used for ASO probes are found in Kogan et al, "Genetic Prediction of Hemophilia A" in PCR Protocols, A Guide to Methods and Applications, Academic Press, 1990 and Ruano et al., 87 Proc. Natl. Acad. Sci. USA 6296-6300, 1990. Typically, an ASO will be perfectly complementary to one allele while containing a single mismatch for another allele. Allele-specific oligonucleotides of the invention include ASO probes and ASO primers. ASO probes which usually provide good discrimination between different alleles are those in which a central position of the oligonucleotide probe aligns with the polymoφhic site in the target region (e.g., approximately the 7^th or 8^th position in a 15mer, the 8^th or 9^th position in a 16mer, and the 10^th or 11^th position in a 20mer). An ASO primer of the invention has a 3 ' terminal nucleotide, or preferably a 3 ' penultimate nucleotide, that is complementary to only one nucleotide of a particular SNP, thereby acting as a primer for polymerase-mediated extension only if the allele containing that nucleotide is present. ASO probes and primers hybridizing to either the coding or noncoding strand are contemplated by the invention.

ASO probes and primers listed below use the appropriate nucleotide symbol (R= G or A, Y= T or C, M= A or C, K= G or T, S= G or C, and W= A or T; WIPO standard ST.25) at the position of the polymoφhic site to represent the two alternative allelic variants observed at that polymoφhic site.

A preferred ASO probe for detecting MMP 13 gene polymoφhisms comprises a nucleotide sequence, listed 5' to 3' selected from the group consisting of:

AGTGACTRGGAAGTG (SEQ ID NO 6) and its complement,

TTCCCTCKAACTCTT (SEQ ID NO 7) and its complement,

AGTAGAAYTTTCAAA (SEQ ID NO 8) and its complement,

TTGAATTRAAAAGCA (SEQ ID NO 9) and its complement,

ATTGAAASTTTCTCT (SEQ ID NO 10) and its complement,

TATTTTASAAACATT (SEQ ID NO 11) and its complement,

TTTAAAG ATTCAAC (SEQ ID NO 12) and its complement,

ACACACTRGTATGAA (SEQ ID NO 13) and its complement,

GTAACTCRGTTATAG (SEQ ID NO 14) and its complement,

AATGTGAYAGAAAAC (SEQ ID NO 15) and its_ complement,

ATGGTTAYGACATTC (SEQ ID NO 16) and its complement,

TTTGAGGRTACAGGC (SEQ ID NO 17) and its complement,

CATTGCAYTGACTTT (SEQ ID NO 18) and its complement,

TTAAAGASGAGCTTC (SEQ ID NO 19) and its complement,

TTAAAAAYATATTTA (SEQ ID NO 20) and its complement,

GAGAGCCRTACATAA (SEQ ID NO 21) and its complement,

TCTGTTTRATTTTCA (SEQ ID NO 22) and its complement, and

ATTGTCCRGTTTGTC (SEQ ID NO 23) and its complement .

A preferred ASO primer for detecting MMP 13 gene polymoφhisms comprises a nucleotide sequence, listed 5' to 3', selected from the group consisting of:

CCTTCAAGTGACTRG SEQ ID NO 24) ; GGTTTCCACTTCCYA (SEQ ID NO 25)

AATGTTTTCCCTCKA SEQ ID NO 26) ; CCATTTAAGAGTTMG (SEQ ID NO 27)

ACCTCTAGTAGAAYT SEQ ID NO 28) ; AATTTCTTTGAAART (SEQ ID NO 29)

TGCATTTTGAATTRA SEQ ID NO 30) ; GTATATTGCTTTTYA (SEQ ID NO 31)

AAGTTCATTGAAAST SEQ ID NO 32) ; TTTAACAGAGAAAST (SEQ ID NO 33)

CACATTTATTTTASA SEQ ID NO 34) ; TATAAAAATGTTTST (SEQ ID NO 35)

TTACTTTTTAAAGMA SEQ ID NO 36) ; ATTAGTGTTGAATKC (SEQ ID NO 37)

CAGTGAACACACTRG SEQ ID NO 38) ; TTTAACTTCATACYA (SEQ ID NO 39)

TCTATGGTAACTCRG SEQ ID NO 40) ; CTATTTCTATAACYG (SEQ ID NO 41)

TTTAAAAATGTGAYA SEQ ID NO 42) ; AATGGAGTTTTCTRT (SEQ ID NO 43)

CTCTTAATGGTTAYG SEQ ID NO 44) ; CTTCCAGAATGTCRT (SEQ ID NO 45)

GTTCACTTTGAGGRT SEQ ID NO 46) ; AGTCTTGCCTGTAYC (SEQ ID NO 47)

TAGTTACATTGCAYT SEQ ID NO 48) ; CTCTTTAAAGTCART (SEQ ID NO 49)

CTGACTTTAAAGASG SEQ ID NO 50) ; GCTACAGAAGCTCST (SEQ ID NO 51)

CATAATTTAAAAAYA SEQ ID NO 52) ; TCCTTATAAATATRT (SEQ ID NO 53)

ATTTGGGAGAGCCRT SEQ ID NO 54) ; TTTCTTTTATGTAYG (SEQ ID NO 55)

TTGAATTCTGTTTRA SEQ ID NO 56) ; CAAAAGTGAAAATYA (SEQ ID NO 57)

TATTTAATTGTCCRG SEQ ID NO 58) ; and ATCCAAGACAAACYG (SEQ ID NO 59 )

Other genotyping oligonucleotides of the invention hybridize to a target region located one to several nucleotides downstream of one of the novel polymoφhic sites identified herein. Such oligonucleotides are useful in polymerase-mediated primer extension methods for detecting one of the novel polymoφhisms described herein and therefore such genotyping oligonucleotides are referred to herein as "primer-extension oligonucleotides". In a preferred embodiment, the 3 '-terminus of a primer- extension oligonucleotide is a deoxynucleotide complementary to the nucleotide located immediately adjacent to the polymoφhic site. A particularly preferred oligonucleotide primer for detecting MMP 13 gene polymoφhisms by primer extension terminates in a nucleotide sequence, listed 5' to 3', selected from the group consisting of:

TCAAGTGACT (SEQ ID NO 60) TTCCACTTCC (SEQ ID NO 61) ;

GTTTTCCCTC (SEQ ID NO 62) TTTAAGAGTT (SEQ ID NO 63) ;

TCTAGTAGAA (SEQ ID NO 64) TTCTTTGAAA (SEQ ID NO 65) ;

ATTTTGAATT (SEQ ID NO 66 ) ^■ TATTGCTTTT (SEQ ID NO 67) ;

TTCATTGAAA (SEQ ID NO 68) ^• AACAGAGAAA (SEQ ID NO 69 ) ;

ATTTATTTTA (SEQ ID NO 70) ^• AAAAATGTTT (SEQ ID NO 71) ;

CTTTTTAAAG (SEQ ID NO 72) ^• AGTGTTGAAT (SEQ ID NO 73) ;

TGAACACACT (SEQ ID NO 74) ^■ AACTTCATAC (SEQ ID NO 75);

ATGGTAACTC (SEQ ID NO 76) TTTCTATAAC (SEQ ID NO 77) ;

AAAAATGTGA (SEQ ID NO 78) GGAGTTTTCT (SEQ ID NO 79) ;

TTAATGGTTA (SEQ ID NO 80) CCAGAATGTC (SEQ ID NO 81) ;

CACTTTGAGG (SEQ ID NO 82) CTTGCCTGTA (SEQ ID NO 83) ;

TTACATTGCA (SEQ ID NO 84) TTTAAAGTCA (SEQ ID NO 85) ;

ACTTTAAAGA (SEQ ID NO 86) ACAGAAGCTC (SEQ ID NO 87) ;

AATTTAAAAA (SEQ ID NO 88) TTATAAATAT (SEQ ID NO 89) ;

TGGGAGAGCC (SEQ ID NO 90) CTTTTATGTA (SEQ ID NO 91) ;

AATTCTGTTT (SEQ ID NO 92) AAGTGAAAAT (SEQ ID NO 93) ;

TTAATTGTCC (SEQ ID NO 94) and CAAGACAAAC (SEQ ID NO: 95

In some embodiments, a composition contains two or more differently labeled genotyping oligonucleotides for simultaneously probing the identity of nucleotides at two or more polymoφhic sites. It is also contemplated that primer compositions may contain two or more sets of allele-specific primer pairs to allow simultaneous targeting and amplification of two or more regions containing a polymoφhic site.

MMP 13 genotyping oligonucleotides of the invention may also be immobilized on or synthesized on a solid surface such as a microchip, bead, or glass slide (see, e.g., WO 98/20020 and WO 98/20019). Such immobilized genotyping oligonucleotides may be used in a variety of polymoφhism detection assays, including but not limited to probe hybridization and polymerase extension assays. Immobilized MMP 13 genotyping oligonucleotides of the invention may comprise an ordered array of oligonucleotides designed to rapidly screen a DNA sample for polymoφhisms in multiple genes at the same time.

In another embodiment, the invention provides a kit comprising at least two genotyping oligonucleotides packaged in separate containers. The kit may also contain other components such as hybridization buffer (where the oligonucleotides are to be used as a probe) packaged in a separate container. Alternatively, where the oligonucleotides are to be used to amplify a target region, the kit may contain, packaged in separate containers, a polymerase and a reaction buffer optimized for primer extension mediated by the polymerase, such as PCR.

The above described oligonucleotide compositions and kits are useful in methods for genotyping and/or haplotyping the MMP 13 gene in an individual. As used herein, the terms "MMP 13 genotype" and "MMP 13 haplotype" mean the genotype or haplotype contains the nucleotide pair or nucleotide, respectively, that is present at one or more of the novel polymoφhic sites described herein and may optionally also include the nucleotide pair or nucleotide present at one or more additional polymoφhic sites in the MMP 13 gene. The additional polymoφhic sites may be currently known polymoφhic sites or sites that are subsequently discovered.

One embodiment of the genotyping method involves isolating from the individual a nucleic acid sample comprising the two copies of the MMP 13 gene, or a fragment thereof, that are present in the individual, and determining the identity of the nucleotide pair at one or more polymoφhic sites selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS13, PS14, PS15, PS16, PS17 and PS18 in the two copies to assign a MMP13 genotype to the individual. As will be readily understood by the skilled artisan, the two "copies" of a gene in an individual may be the same allele or may be different alleles. In a particularly preferred embodiment, the genotyping method comprises determining the identity of the nucleotide pair at each of PS1-PS18.

Typically, the nucleic acid sample is isolated from a biological sample taken from the individual, such as a blood sample or tissue sample. Suitable tissue samples include whole blood, semen, saliva, tears, urine, fecal material, sweat, buccal, skin and hair. The nucleic acid sample may be comprised of genomic DNA, mRNA, or cDNA and, in the latter two cases, the biological sample must be obtained from a tissue in which the MMP 13 gene is expressed. Furthermore it will be understood by the skilled artisan that mRNA or cDNA preparations would not be used to detect polymoφhisms located in introns or in 5' and 3' untranslated regions. If a MMP 13 gene fragment is isolated, it must contain the polymoφhic site(s) to be genotyped. One embodiment of the haplotyping method comprises isolating from the individual a nucleic acid sample containing only one of the two copies of the MMP 13 gene, or a fragment thereof, that is present in the individual and determining in that copy the identity of the nucleotide at one or more polymoφhic sites selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS13, PS14, PS15, PS16, PS17 and PSlδ in that copy to assign a MMP13 haplotype to the individual. The nucleic acid may be isolated using any method capable of separating the two copies of the MMP 13 gene or fragment such as one of the methods described above for preparing MMP 13 isogenes, with targeted in vivo cloning being the preferred approach. As will be readily appreciated by those skilled in the art, any individual clone will only provide haplotype information on one of the two MMP 13 gene copies present in an individual. If haplotype information is desired for the individual's other copy, additional MMP13 clones will need to be examined. Typically, at least five clones should be examined to have more than a 90% probability of haplotyping both copies of the MMP 13 gene in an individual. In a particularly preferred embodiment, the nucleotide at each of PSl- PSlδ is identified.

In another embodiment, the haplotyping method comprises determining whether an individual has one or more of the MMP13 haplotypes shown in Table 5. This can be accomplished by identifying, for one or both copies of the individual's MMP 13 gene, the phased sequence of nucleotides present at each of PS 1 -PS 18. The present invention also contemplates that typically only a subset of PS 1 -PS 18 will need to be directly examined to assign to an individual one or more of the haplotypes shown in Table 5. This is because at least one polymoφhic site in a gene is frequently in strong linkage disequilibrium with one or more other polymoφhic sites in that gene (Drysdale, CM et al. 2000 PNAS 97:10483-10488; Rieder MJ et al. 1999 Nature Genetics 22:59-62). Two sites are said to be in linkage disequilibrium if the presence of a particular variant at one site enhances the predictability of another variant at the second site (Stephens, JC 1999, Mol. Diag. 4:309-317). Techniques for determining whether any two polymoφhic sites are in linkage disequilibrium are well-known in the art (Weir B.S. 1996 Genetic Data Analysis II, Sinauer Associates, Inc. Publishers, Sunderland, MA).

In a preferred embodiment, a MMP 13 haplotype pair is determined for an individual by identifying the phased sequence of nucleotides at one or more polymoφhic sites selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PS11, PS12, PS13, PS14, PS 15, PS 16, PS 17 and PS 18 in each copy of the MMP 13 gene that is present in the individual. In a particularly preferred embodiment, the haplotyping method comprises identifying the phased sequence of nucleotides at each of PS 1 -PS 18 in each copy of the MMP 13 gene. When haplotyping both copies of the gene, the identifying step is preferably performed with each copy of the gene being placed in separate containers. However, it is also envisioned that if the two copies are labeled with different tags, or are otherwise separately distinguishable or identifiable, it could be possible in some cases to perform the method in the same container. For example, if first and second copies of the gene are labeled with different first and second fluorescent dyes, respectively, and an allele-specific oligonucleotide labeled with yet a third different fluorescent dye is used to assay the polymoφhic site(s), then detecting a combination of the first and third dyes would identify the polymoφhism in the first gene copy while detecting a combination of the second and third dyes would identify the polymoφhism in the second gene copy.

In both the genotyping and haplotyping methods, the identity of a nucleotide (or nucleotide pair) at a polymoφhic site(s) may be determined by amplifying a target region(s) containing the polymoφhic site(s) directly from one or both copies of the MMP 13 gene, or a fragment thereof, and the sequence of the amplified region(s) determined by conventional methods. It will be readily appreciated by the skilled artisan that only one nucleotide will be detected at a polymoφhic site in individuals who are homozygous at that site, while two different nucleotides will be detected if the individual is heterozygous for that site. The polymoφhism may be identified directly, known as positive-type identification, or by inference, referred to as negative-type identification. For example, where a SNP is known to be guanine and cytosine in a reference population, a site may be positively determined to be either guanine or cytosine for an individual homozygous at that site, or both guanine and cytosine, if the individual is heterozygous at that site. Alternatively, the site may be negatively determined to be not guanine (and thus cytosine/cytosine) or not cytosine (and thus guanine/guanine).

The target region(s) may be amplified using any oligonucleotide-directed amplification method, including but not limited to polymerase chain reaction (PCR) (U.S. Patent No. 4,965,188), ligase chain reaction (LCR) (Barany et al., Proc. Natl. Acad. Sci. USA 88:189-193, 1991; WO90/01069), and oligonucleotide ligation assay (OLA) (Landegren et al., Science 241:1077-1080, 1988).

Other known nucleic acid amplification procedures may be used to amplify the target region including transcription-based amplification systems (U.S. Patent No. 5,130,238; EP 329,822; U.S. Patent No. 5, 169,766, WO89/06700) and isothermal methods (Walker et al., Proc. Natl. Acad. Sci. USA 89:392-396, 1992).

A polymoφhism in the target region may also be assayed before or after amplification using one of several hybridization-based methods known in the art. Typically, allele-specific oligonucleotides are utilized in performing such methods. The allele-specific oligonucleotides may be used as differently labeled probe pairs, with one member of the pair showing a perfect match to one variant of a target sequence and the other member showing a perfect match to a different variant. In some embodiments, more than one polymoφhic site may be detected at once using a set of allele- specific oligonucleotides or oligonucleotide pairs. Preferably, the members of the set have melting temperatures within 5°C, and more preferably within 2°C, of each other when hybridizing to each of the polymoφhic sites being detected.

Hybridization of an allele-specific oligonucleotide to a target polynucleotide may be performed with both entities in solution, or such hybridization may be performed when either the oligonucleotide or the target polynucleotide is covalently or noncovalently affixed to a solid support. Attachment may be mediated, for example, by antibody-antigen interactions, poly-L-Lys, streptavidin or avidin-biotin, salt bridges, hydrophobic interactions, chemical linkages, UV cross-linking baking, etc. Allele-specific oligonucleotides may be synthesized directly on the solid support or attached to the solid support subsequent to synthesis. Solid-supports suitable for use in detection methods of the invention include substrates made of silicon, glass, plastic, paper and the like, which may be formed, for example, into wells (as in 96-well plates), slides, sheets, membranes, fibers, chips, dishes, and beads. The solid support may be treated, coated or derivatized to facilitate the immobilization of the allele-specific oligonucleotide or target nucleic acid.

The genotype or haplotype for the MMP 13 gene of an individual may also be determined by hybridization of a nucleic acid sample containing one or both copies of the gene, or fragment(s) thereof, to nucleic acid arrays and subarrays such as described in WO 95/11995. The arrays would contain a battery of allele-specific oligonucleotides representing each of the polymoφhic sites to be included in the genotype or haplotype.

The identity of polymoφhisms may also be determined using a mismatch detection technique, including but not limited to the RNase protection method using riboprobes (Winter et al, Proc. Natl. Acad. Sci. USA 82:7575, 1985; Meyers et al., Science 230: 1242, 1985) and proteins which recognize nucleotide mismatches, such as the E. coli mutS protein (Modrich, P. Ann. Rev. Genet. 25:229-253,

1991). Alternatively, variant alleles can be identified by single strand conformation polymoφhism

(SSCP) analysis (Orita et al., Genomics 5:874-879, 1989; Humphries et al., in Molecular Diagnosis of Genetic Diseases, R. Elles, ed., pp. 321-340, 1996) or denaturing gradient gel electrophoresis (DGGE) (Wartell et al, Nucl. Acids Res. 18:2699-2706, 1990; Sheffield et al, Proc. Natl. Acad. Sci. USA 86:232-236, 1989).

A polymerase-mediated primer extension method may also be used to identify the polymoφhism(s). Several such methods have been described in the patent and scientific literature and include the "Genetic Bit Analysis" method (W092/15712) and the ligase/polymerase mediated genetic bit analysis (U.S. Patent 5,679,524. Related methods are disclosed in W091/02087, WO90/09455, W095/17676, U.S. Patent Nos. 5,302,509, and 5,945,283. Extended primers containing a polymoφhism may be detected by mass spectrometry as described in U.S. Patent No. 5,605,798. Another primer extension method is allele-specific PCR (Ruano et al, Nucl. Acids Res. 17:8392, 1989; Ruano et al, Nucl. Acids Res. 19, 6877-6882, 1991; WO 93/22456; Turki et al, J. Clin. Invest. 95:1635-1641, 1995). In addition, multiple polymoφhic sites may be investigated by simultaneously amplifying multiple regions of the nucleic acid using sets of allele-specific primers as described in Wallace et al. (WO89/10414). In addition, the identity of the allele(s) present at any of the novel polymoφhic sites described herein may be indirectly determined by genotyping another polymoφhic site that is in linkage disequilibrium with the polymoφhic site that is of interest. Polymoφhic sites in linkage disequilibrium with the presently disclosed polymoφhic sites may be located in regions of the gene or in other genomic regions not examined herein. Genotyping of a polymoφhic site in linkage disequilibrium with the novel polymoφhic sites described herein may be performed by, but is not limited to, any of the above-mentioned methods for detecting the identity of the allele at a polymoφhic site.

In another aspect of the invention, an individual's MMP13 haplotype pair is predicted from its MMP 13 genotype using information on haplotype pairs known to exist in a reference population. In its broadest embodiment, the haplotyping prediction method comprises identifying a MMP 13 genotype for the individual at two or more MMP 13 polymoφhic sites described herein, enumerating all possible haplotype pairs which are consistent with the genotype, accessing data containing MMP 13 haplotype pairs identified in a reference population, and assigning a haplotype pair to the individual that is consistent with the data. In one embodiment, the reference haplotype pairs include the MMP 13 haplotype pairs shown in Table 4. Generally, the reference population should be composed of randomly-selected individuals representing the major ethnogeographic groups of the world. A preferred reference population for use in the methods of the present invention comprises an approximately equal number of individuals from Caucasian, African-descent, Asian and Hispanic-Latino population groups with the minimum number of each group being chosen based on how rare a haplotype one wants to be guaranteed to see. For example, if one wants to have a q% chance of not missing a haplotype that exists in the population at a p% frequency of occurring in the reference population, the number of individuals (n) who must be sampled is given by 2n=log(l-q)/log(l-p) where p and q are expressed as fractions. A preferred reference population allows the detection of any haplotype whose frequency is at least 10% with about 99% certainty and comprises about 20 unrelated individuals from each of the four population groups named above. A particularly preferred reference population includes a 3-generation family representing one or more of the four population groups to serve as controls for checking quality of haplotyping procedures.

In a preferred embodiment, the haplotype frequency data for each ethnogeographic group is examined to determine whether it is consistent with Hardy- Weinberg equilibrium. Hardy- Weinberg equilibrium (D.L. Haiti et al. Principles of Population Genomics, Sinauer Associates (Sunderland, MA), 3^rd Ed, 1997) postulates that the frequency of finding the haplotype pair H_x I H₂is equal to p_H__w(H_x IH_) = 2p(H_l)p(H₂) if H ≠ H₂ and p_H__w(H_x IH_) - p(H_x)p(H₂) if H = H₂ . A statistically significant difference between the observed and expected haplotype frequencies could be due to one or more factors including significant inbreeding in the population group, strong selective pressure on the gene, sampling bias, and/or errors in the genotyping process. If large deviations from Hardy- Weinberg equilibrium are observed in an ethnogeographic group, the number of individuals in that group can be increased to see if the deviation is due to a sampling bias. If a larger sample size does not reduce the difference between observed and expected haplotype pair frequencies, then one may wish to consider haplotyping the individual using a direct haplotyping method such as, for example, CLASPER System™ technology (U.S. Patent No. 5,866,404), single molecule dilution, or allele-specific long-range PCR (Michalotos-Beloin et al. Nucleic Acids Res. 24:4841-4843, 1996). In one embodiment of this method for predicting a MMP 13 haplotype pair for an individual, the assigning step involves performing the following analysis. First, each of the possible haplotype pairs is compared to the haplotype pairs in the reference population. Generally, only one of the haplotype pairs in the reference population matches a possible haplotype pair and that pair is assigned to the individual. Occasionally, only one haplotype represented in the reference haplotype pairs is consistent with a possible haplotype pair for an individual, and in such cases the individual is assigned a haplotype pair containing this known haplotype and a new haplotype derived by subtracting the known haplotype from the possible haplotype pair. Alternatively, the haplotype pair in an individual may be predicted from the individual's genotype for that gene using reported methods (e.g., Clark et al. 1990 Mol Bio Evol 7:111-22) or through a commercial haplotyping service such as offered by Genaissance Pharmaceuticals, Inc. (New Haven, CT). In rare cases, either no haplotypes in the reference population are consistent with the possible haplotype pairs, or alternatively, multiple reference haplotype pairs are consistent with the possible haplotype pairs. In such cases, the individual is preferably haplotyped using a direct molecular haplotyping method such as, for example, CLASPER System™ technology (U.S. Patent No. 5,866,404), SMD, or allele-specific long-range PCR (Michalotos-Beloin et al, supra). A preferred process for predicting MMP13 haplotype pairs from MMP13 genotypes is described in U.S. Provisional Application Serial No. 60/198,340 and the corresponding International Application, The invention also provides a method for determining the frequency of a MMP 13. genotype, haplotype, or haplotype pair in a population. The method comprises, for each member of the population, deteπnining the genotype or the haplotype pair for the novel MMP 13 polymoφhic sites described herein, and calculating the frequency any particular genotype, haplotype, or haplotype pair is found in the population. The population may be a reference population, a family population, a same sex population, a population group, or a trait population (e.g., a group of individuals exhibiting a trait of interest such as a medical condition or response to a therapeutic treatment).

In another aspect of the invention, frequency data for MMP 13 genotypes, haplotypes, and/or haplotype pairs are determined in a reference population and used in a method for identifying an association between a trait and a MMP 13 genotype, haplotype, or haplotype pair. The trait may be any detectable phenotype, including but not limited to susceptibility to a disease or response to a treatment. The method involves obtaining data on the frequency of the genotype(s), haplotype(s), or haplotype pair(s) of interest in a reference population as well as in a population exhibiting the trait. Frequency data for one or both of the reference and trait populations may be obtained by genotyping or haplotyping each individual in the populations using one of the methods described above. The ' haplotypes for the trait population may be determined directly or, alternatively, by the predictive genotype to haplotype approach described above. In another embodiment, the frequency data for the reference and/or trait populations is obtained by accessing previously determined frequency data, which may be in written or electronic fonn. For example, the frequency data may be present in a database that is accessible by a computer. Once the frequency data is obtained, the frequencies of the genotype(s), haplotype(s), or haplotype pair(s) of interest in the reference and trait populations are compared. In a preferred embodiment, the frequencies of all genotypes, haplotypes, and/or haplotype pairs observed in the populations are compared. If a particular MMP 13 genotype, haplotype, or haplotype pair is more frequent in the trait population than in the reference population at a statistically significant amount, then the trait is predicted to be associated with that MMP 13 genotype, haplotype or haplotype pair. Preferably, the MMP 13 genotype, haplotype, or haplotype pair being compared in the trait and reference populations is selected from the full-genotypes and full-haplotypes shown in Tables 4 and 5, or from sub-genotypes and sub-haplotypes derived from these genotypes and haplotypes. In a preferred embodiment of the method, the trait of interest is a clinical response exhibited by a patient to some therapeutic treatment, for example, response to a drug targeting MMP 13 or response to a therapeutic treatment for a medical condition. As used herein, "medical condition" includes but is not limited to any condition or disease manifested as one or more physical and/or psychological symptoms for which treatment is desirable, and includes previously and newly identified diseases and other disorders. As used herein the term "clinical response" means any or all of the following: a quantitative measure of the response, no response, and adverse response (i.e., side effects).

In order to deduce a correlation between clinical response to a treatment and a MMP 13 genotype, haplotype, or haplotype pair, it is necessary to obtain data on the clinical responses exhibited by a population of individuals who received the treatment, hereinafter the "clinical population". This clinical data may be obtained by analyzing the results of a clinical trial that has already been run and/or the clinical data may be obtained by designing and carrying out one or more new clinical trials. As used herein, the term "clinical trial" means any research study designed to collect clinical data on responses to a particular treatment, and includes but is not limited to phase I, phase II and phase III clinical trials. Standard methods are used to define the patient population and to enroll subjects.

It is preferred that the individuals included in the clinical population have been graded for the existence of the medical condition of interest. This is important in cases where the symptom(s) being presented by the patients can be caused by more than one underlying condition, and where treatment of the underlying conditions are not the same. An example of this would be where patients experience breathing difficulties that are due to either asthma or respiratory infections. If both sets were treated with an asthma medication, there would be a spurious group of apparent non-responders that did not actually have asthma. These people would affect the ability to detect any correlation between haplotype and treatment outcome. This grading of potential patients could employ a standard physical exam or one or more lab tests. Alternatively, grading of patients could use haplotyping for situations where there is a strong correlation between haplotype pair and disease susceptibility or severity.

The therapeutic treatment of interest is administered to each individual in the trial population and each individual's response to the treatment is measured using one or more predetermined criteria. It is contemplated that in many cases, the trial population will exhibit a range of responses and that the investigator will choose the number of responder groups (e.g., low, medium, high) made up by the various responses. In addition, the MMP 13 gene for each individual in the trial population is genotyped and/or haplotyped, which may be done before or after administering the treatment.

After both the clinical and polymoφhism data have been obtained, correlations between individual response and MMP 13 genotype or haplotype content are created. Correlations may be produced in several ways. In one method, individuals are grouped by their MMP13 genotype or haplotype (or haplotype pair) (also referred to as a polymoφhism group), and then the averages and standard deviations of clinical responses exhibited by the members of each polymoφhism group are calculated. These results are then analyzed to determine if any observed variation in clinical response between polymoφhism groups is statistically significant. Statistical analysis methods which may be used are described in L.D. Fisher and G. vanBelle, "Biostatistics: A Methodology for the Health Sciences", Wiley-Interscience (New York) 1993. This analysis may also include a regression calculation of which polymoφhic sites in the MMP 13 gene give the most significant contribution to the differences in phenotype. One regression model useful in the invention is described in PCT Application

Serial No. PCT/US00/17540, entitled "Methods for Obtaining and Using Haplotype Data".

A second method for finding correlations between MMP 13 haplotype content and clinical responses uses predictive models based on error-minimizing optimization algorithms. One of many possible optimization algorithms is a genetic algorithm (R. Judson, "Genetic Algorithms and Their Uses in Chemistry" in Reviews in Computational Chemistry, Vol. 10, pp. 1-73, K. B. Lipkowitz and D. B. Boyd, eds. (VCH Publishers, New York, 1997). Simulated annealing (Press et al, "Numerical Recipes in C: The Art of Scientific Computing", Cambridge University Press (Cambridge) 1992, Ch. 10), neural networks (E. Rich and K. Knight, "Artificial Intelligence", 2^nd Edition (McGraw-Hill, New York, 1991, Ch. 18), standard gradient descent methods (Press et al, supra, Ch. 10), or other global or local optimization approaches (see discussion in Judson, supra) could also be used. Preferably, the correlation is found using a genetic algorithm approach as described in PCT Application Serial No. PCT/USOO/17540.

Correlations may also be analyzed using analysis of variation (ANOVA) techniques to determine how much of the variation in the clinical data is explained by different subsets of the polymoφhic sites in the MMP13 gene. As described in PCT Application Serial No. PCT/USOO/17540, ANOVA is used to test hypotheses about whether a response variable is caused by or correlated with one or more traits or variables that can be measured (Fisher and vanBelle, supra, Ch. 10).

From the analyses described above, a mathematical model may be readily constructed by the skilled artisan that predicts clinical response as a function of MMP 13 genotype or haplotype content. Preferably, the model is validated in one or more follow-up clinical trials designed to test the model.

The identification of an association between a clinical response and a genotype or haplotype (or haplotype pair) for the MMP 13 gene may be the basis for designing a diagnostic method to determine those individuals who will or will not respond to the treatment, or alternatively, will respond at a lower level and thus may require more treatment, i.e., a greater dose of a drug. The diagnostic method may take one of several forms: for example, a direct DNA test (i.e., genotyping or haplotyping one or more of the polymoφhic sites in the MMP 13 gene), a serological test, or a physical exam measurement. The only requirement is that there be a good correlation between the diagnostic test results and the underlying MMP 13 genotype or haplotype that is in turn correlated with the clinical response. In a preferred embodiment, this diagnostic method uses the predictive haplotyping method described above.

In another embodiment, the invention provides an isolated polynucleotide comprising a polymoφhic variant of the MMP 13 gene or a fragment of the gene which contains at least one of the novel polymoφhic sites described herein. The nucleotide sequence of a variant MMP 13 gene is identical to the reference genomic sequence for those portions of the gene examined, as described in the Examples below, except that it comprises a different nucleotide at one or more of the novel polymoφhic sites PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS13, PS14, PS15, PS16, PS17 and PS18. Similarly, the nucleotide sequence of a variant fragment of the MMP13 gene is identical to the corresponding portion of the reference sequence except for having a different nucleotide at one or more of the novel polymoφhic sites described herein. Thus, the invention specifically does not include polynucleotides comprising a nucleotide sequence identical to the reference sequence of the MMP 13 gene, which is defined by haplotype 3, (or other reported MMP 13 sequences) or to portions of the reference sequence (or other reported MMP 13 sequences), except for genotyping oligonucleotides as described above.

The location of a polymoφhism in a variant gene or fragment is identified by aligning its sequence against SEQ ID NOS: 1-3. The polymoφhism is selected from the group consisting of guanine at PS1, thymine at PS2, thymine at PS3, adenine at PS4, guanine at PS5, guanine at PS6, adenine at PS7, adenine at PS8, adenine at PS9, thymine at PS10, cytosine at PS11, guanine at PS12, thymine at PS13, cytosine at PS14, cytosine at PS15, guanine at PS16, guanine at PS17 and guanine at PS 18. In a preferred embodiment, the polymoφhic variant comprises a naturally-occurring isogene of the MMP 13 gene which is defined by any one of haplotypes 1- 2 and 4 - 15 shown in Table 5 below. Polymoφhic variants of the invention may be prepared by isolating a clone containing the MMP 13 gene from a human genomic library. The clone may be sequenced to determine the identity of the nucleotides at the novel polymoφhic sites described herein. Any particular variant claimed herein could be prepared from this clone by performing in vitro mutagenesis using procedures well-known in the art.

MMP 13 isogenes may be isolated using any method that allows separation of the two "copies" of the MMP 13 gene present in an individual, which, as readily understood by the skilled artisan, may be the same allele or different alleles. Separation methods include targeted in vivo cloning (TIVC) in yeast as described in WO 98/01573, U.S. Patent No. 5,866,404, and U.S. Patent No. 5,972,614. Another method, which is described in U.S. Patent No. 5,972,614, uses an allele specific oligonucleotide in combination with primer extension and exonuclease degradation to generate hemizygous DNA targets. Yet other methods are single molecule dilution (SMD) as described in Ruano et al, Proc. Natl. Acad. Sci. 87:6296-6300, 1990; and allele specific PCR (Ruano et al, 1989, supra; Ruano et al, 1991, supra; Michalatos-Beloin et al, supra). The invention also provides MMP 13 genome anthologies, which are collections of MMP 13 isogenes found in a given population. The population may be any group of at least two individuals, including but not limited to a reference population, a population group, a family population, a clinical population, and a same sex population. A MMP 13 genome anthology may comprise individual MMP 13 isogenes stored in separate containers such as microtest tubes, separate wells of a microtitre plate and the like. Alternatively, two or more groups of the MMP13 isogenes in the anthology may be stored in separate containers. Individual isogenes or groups of isogenes in a genome anthology may be stored in any convenient and stable form, including but not limited to in buffered solutions, as DNA precipitates, freeze-dried preparations and the like. A preferred MMP 13 genome anthology of the invention comprises a set of isogenes defined by the haplotypes shown in Table 5 below. An isolated polynucleotide containing a polymoφhic variant nucleotide sequence of the invention may be operably linked to one or more expression regulatory elements in a recombinant expression vector capable of being propagated and expressing the encoded MMP 13 protein in a prokaryotic or a eukaryotic host cell. Examples of expression regulatory elements which may be used include, but are not limited to, the lac system, operator and promoter regions of phage lambda, yeast promoters, and promoters derived from vaccinia virus, adenovirus, retroviruses, or SV40. Other regulatory elements include, but are not limited to, appropriate leader sequences, termination codons, polyadenylation signals, and other sequences required for the appropriate transcription and subsequent translation of the nucleic acid sequence in a given host cell. Of course, the correct combinations of expression regulatory elements will depend on the host system used. In addition, it is understood that the expression vector contains any additional elements necessary for its transfer to and subsequent replication in the host cell. Examples of such elements include, but are not limited to, origins of replication and selectable markers. Such expression vectors are commercially available or are readily constructed using methods known to those in the art (e.g., F. Ausubel et al, 1987, in "Current Protocols in Molecular Biology", John Wiley and Sons, New York, New York). Host cells which may be used to express the variant MMP 13 sequences of the invention include, but are not limited to, eukaryotic and mammalian cells, such as animal, plant, insect and yeast cells, and prokaryotic cells, such as E. coli, or algal cells as known in the art. The recombinant expression vector may be introduced into the host cell using any method known to those in the art including, but not limited to, microinjection, electroporation, particle bombardment, transduction, and transfection using DEAE-dextran, lipofection, or calcium phosphate (see e.g., Sambrook et al. (1989) in "Molecular Cloning. A Laboratory Manual", Cold Spring Harbor Press, Plainview, New York). In a preferred aspect, eukaryotic expression vectors that function in eukaryotic cells, and preferably mammalian cells, are used. Non-limiting examples of such vectors include vaccinia virus vectors, adenovirus vectors, heφes virus vectors, and baculovirus transfer vectors. Preferred eukaryotic cell lines include COS cells, CHO cells, HeLa cells, NIH/3T3 cells, and embryonic stem cells (Thomson, J. A. et al, 1998 Science 282:1145-1147). Particularly preferred host cells are mammalian cells. As will be readily recognized by the skilled artisan, expression of polymoφhic variants of the

MMP 13 gene will produce MMP 13 mRNAs varying from each other at any polymoφhic site retained in the spliced and processed mRNA molecules. These mRNAs can be used for the preparation of a MMP 13 cDNA comprising a nucleotide sequence which is a polymoφhic variant of the MMP 13 reference coding sequence shown in Figure 4. Thus, the invention also provides MMP 13 mRNAs and corresponding cDNAs which comprise a nucleotide sequence that is identical to SEQ ID NO:4 (Fig. 4), or its corresponding RNA sequence, except for having one or more polymoφhisms selected from the group consisting of thymine at a position corresponding to nucleotide 326, cytosine at a position corresponding to nucleotide 1080 and guanine at a position corresponding to nucleotide 1169. A particularly preferred polymoφhic cDNA variant comprises the coding sequence of a MMP 13 isogene defined by haplotypes 8c, 10c, and 12c. Fragments of these variant mRNAs and cDNAs are included in the scope of the invention, provided they contain the novel polymoφhisms described herein. The invention specifically excludes polynucleotides identical to previously identified and characterized MMP 13 cDNAs and fragments thereof. Polynucleotides comprising a variant RNA or DNA sequence may be isolated from a biological sample using well-known molecular biological procedures or may be chemically synthesized.

As used herein, a polymoφhic variant of a MMP 13 gene fragment comprises at least one novel polymoφhism identified herein and has a length of at least 10 nucleotides and may range up to the full length of the gene. Preferably, such fragments are between 100 and 3000 nucleotides in length, and more preferably between 200 and 2000 nucleotides in length, and most preferably between 500 and 1000 nucleotides in length.

In describing the MMP 13 polymoφhic sites identified herein, reference is made to the sense strand of the gene for convenience. However, as recognized by the skilled artisan, nucleic acid molecules containing the MMP 13 gene may be complementary double stranded molecules and thus reference to a particular site on the sense strand refers as well to the corresponding site on the complementary antisense strand. Thus, reference may be made to the same polymoφhic site on either strand and an oligonucleotide may be designed to hybridize specifically to either strand at a target region containing the polymoφhic site. Thus, the invention also includes single-stranded polynucleotides which are complementary to the sense strand of the MMP 13 genomic variants described herein.

Polynucleotides comprising a polymoφhic gene variant or fragment may be useful for therapeutic puφoses. For example, where a patient could benefit from expression, or increased expression, of a particular MMP 13 protein isoform, an expression vector encoding the isoform may be administered to the patient. The patient may be one who lacks the MMP 13 isogene encoding that isoform or may already have at least one copy of that isogene.

In other situations, it may be desirable to decrease or block expression of a particular MMP 13 isogene. Expression of a MMP 13 isogene may be turned off by transforming a targeted organ, tissue or cell population with an expression vector that expresses high levels of untranslatable mRNA for the isogene. Alternatively, oligonucleotides directed against the regulatory regions (e.g., promoter, introns, enhancers, 3 ' untranslated region) of the isogene may block transcription. Oligonucleotides targeting the transcription initiation site, e.g., between positions -10 and +10 from the start site are preferred. Similarly, inhibition of transcription can be achieved using oligonucleotides that base-pair with region(s) of the isogene DNA to form triplex DNA (see e.g., Gee et al. in Huber, B.E. and B.I. Carr, Molecular and Immunologic Approaches, Futura Publishing Co, Mt. Kisco, N.Y, 1994). Antisense oligonucleotides may also be designed to block translation of MMP 13 mRNA transcribed from a particular isogene. It is also contemplated that ribozymes may be designed that can catalyze the specific cleavage of MMP 13 mRNA transcribed from a particular isogene. The oligonucleotides may be delivered to a target cell or tissue by expression from a vector introduced into the cell or tissue in vivo or ex vivo. Alternatively, the oligonucleotides may be formulated as a pharmaceutical composition for administration to the patient. Oligoribonucleotides and/or oligodeoxynucleotides intended for use as antisense oligonucleotides may be modified to increase stability and half-life. Possible modifications include, but are not limited to phosphorothioate or 2' O-methyl linkages, and the inclusion of nontraditional bases such as inosine and queosine, as well as acetyl-, methyl-, thio-, and similarly modified forms of adenine, cytosine, guanine, thymine, and uracil which are not as easily recognized by endogenous nucleases.

The invention also provides an isolated polypeptide comprising a polymoφhic variant of the reference MMP 13 amino acid sequence shown in Figure 5. The location of a variant amino acid in a MMP 13 polypeptide or fragment of the invention is identified by aligning its sequence against SEQ ID NO:5 (Fig. 5). A MMP13 protein variant of the invention comprises an amino acid sequence identical to SEQ ID NO:5 except for having one or more variant amino acids selected from the group consisting ofleucine at a position corresponding to amino acid position 109 and glycine at a position corresponding to amino acid position 390. The invention specifically excludes amino acid sequences identical to those previously identified for MMP13, including SEQ ID NO:5, and previously described fragments thereof. MMP 13 protein variants included within the invention comprise all amino acid sequences based on SEQ ID NO: 5 and having the combination of amino acid variations described in Table 2 below. In preferred embodiments, a MMP 13 protein variant of the invention is encoded by an isogene defined by one of the observed haplotypes shown in Table 5.

Table 2. Novel Polymoφhic Variants of MMP 13

Polymoφhic Amino Acid Position and Identities Variant

Number 109 390

1 R G

2 L D

3 L G

The invention also includes MMP 13 peptide variants, which are any fragments of a MMP 13 protein variant that contain one or more of the amino acid variations shown in Table 2. A MMP 13 peptide variant is at least 6 amino acids in length and is preferably any number between 6 and 30 amino acids long, more preferably between 10 and 25, and most preferably between 15 and 20 amino acids long. Such MMP 13 peptide variants may be useful as antigens to generate antibodies specific for one of the above MMP 13 isoforms. In addition, the MMP 13 peptide variants may be useful in drug screening assays.

A MMP 13 variant protein or peptide of the invention may be prepared by chemical synthesis or by expressing one of the variant MMP 13 genomic and cDNA sequences as described above. Alternatively, the MMP 13 protein variant may be isolated from a biological sample of an individual having a MMP 13 isogene which encodes the variant protein. Where the sample contains two different MMP13 isoforms (i.e., the individual has different MMP13 isogenes), a particular MMP13 isoform of the invention can be isolated by immunoaffinity chromatography using an antibody which specifically binds to that particular MMP 13 isoform but does not bind to the other MMP 13 isoform. The expressed or isolated MMP 13 protein may be detected by methods known in the art, including Coomassie blue staining, silver staining, and Western blot analysis using antibodies specific for the isoform of the MMP13 protein as discussed further below. MMP13 variant proteins can be purified by standard protein purification procedures known in the art, including differential precipitation, molecular sieve chromatography, ion-exchange chromatography, isoelectric focusing, gel electrophoresis, affinity and immunoaffinity chromatography and the like. (Ausubel et. al, 1987, In Current Protocols in Molecular Biology John Wiley and Sons, New York, New York). In the case of immunoaffinity chromatography, antibodies specific for a particular polymoφhic variant may be used. A polymoφhic variant MMP 13 gene of the invention may also be fused in frame with a heterologous sequence to encode a chimeric MMP 13 protein. The non-MMP13 portion of the chimeric protein may be recognized by a commercially available antibody. In addition, the chimeric protein may also be engineered to contain a cleavage site located between the MMP13 and non-MMP13 portions so that the MMP 13 protein may be cleaved and purified away from the non-MMP13 portion.

An additional embodiment of the invention relates to using a novel MMP 13 protein isoform in any of a variety of drug screening assays. Such screening assays may be performed to identify agents that bind specifically to all known MMP 13 protein isoforms or to only a subset of one or more of these isoforms. The agents may be from chemical compound libraries, peptide libraries and the like. The MMP 13 protein or peptide variant may be free in solution or affixed to a solid support. In one embodiment, high throughput screening of compounds for binding to a MMP 13 variant may be accomplished using the method described in PCT application WO84/03565, in which large numbers of test compounds are synthesized on a solid substrate, such as plastic pins or some other surface, contacted with the MMP 13 protein(s) of interest and then washed. Bound MMP 13 protein(s) are then detected using methods well-known in the art.

In another embodiment, a novel MMP 13 protein isoform may be used in assays to measure the binding affinities of one or more candidate drugs targeting the MMP 13 protein.

In yet another embodiment, when a particular MMP 13 haplotype or group of MMP 13 haplotypes encodes a MMP 13 protein variant with an amino acid sequence distinct from that of MMP 13 protein isoforms encoded by other MMP 13 haplotypes, then detection of that particular MMP 13 haplotype or group of MMP 13 haplotypes may be accomplished by detecting expression of the encoded MMP 13 protein variant using any of the methods described herein or otherwise commonly known to the skilled artisan.

In another embodiment, the invention provides antibodies specific for and immunoreactive with one or more of the novel MMP 13 variant proteins described herein. The antibodies may be either monoclonal or polyclonal in origin. The MMP 13 protein or peptide variant used to generate the antibodies may be from natural or recombinant sources or produced by chemical synthesis using synthesis techniques known in the art. If the MMP 13 protein variant is of insufficient size to be antigenic, it may be conjugated, complexed, or otherwise covalently linked to a carrier molecule to enhance the antigenicity of the peptide. Examples of carrier molecules, include, but are not limited to, albumins (e.g., human, bovine, fish, ovine), and keyhole limpet hemocyanin (Basic and Clinical Immunology, 1991, Eds. D.P. Stites, and A.I. Terr, Appleton and Lange, Norwalk Connecticut, San Mateo, California). In one embodiment, an antibody specifically immunoreactive with one of the novel protein isoforms described herein is administered to an individual to neutralize activity of the MMP 13 isoform expressed by that individual. The antibody may be formulated as a pharmaceutical composition which includes a pharmaceutically acceptable carrier.

Antibodies specific for and immunoreactive with one of the novel protein isoforms described herein may be used to immunoprecipitate the MMP 13 protein variant from solution as well as react with MMP 13 protein isoforms on Western or immunoblots of polyacrylamide gels on membrane supports or substrates. In another preferred embodiment, the antibodies will detect MMP13 protein isoforms in paraffin or frozen tissue sections, or in cells which have been fixed or unfixed and prepared on slides, coverslips, or the like, for use in immunocytochemical, immunohistochemical, and immunofluorescence techniques.

In another embodiment, an antibody specifically immunoreactive with one of the novel MMP 13 protein variants described herein is used in immunoassays to detect this variant in biological samples. In this method, an antibody of the present invention is contacted with a biological sample and the formation of a complex between the MMP 13 protein variant and the antibody is detected. As described, suitable immunoassays include radioimmunoassay, Western blot assay, immunofluorescent assay, enzyme linked immunoassay (ELISA), chemiluminescent assay, immunohistochemical assay, immunocytochemical assay, and the like (see, e.g., Principles and Practice of Immunoassay, 1991, Eds. Christopher P. Price and David J. Neoman, Stockton Press, New York, New York; Current Protocols in Molecular Biology, 1987, Eds. Ausubel et al, John Wiley and Sons, New York, New York). Standard techniques known in the art for ELISA are described in Methods in Immunodiagnosis, 2nd Ed, Eds. Rose and Bigazzi, John Wiley and Sons, New York 1980; and Campbell et al, 1984, Methods in Immunology, W.A. Benjamin, Inc.). Such assays may be direct, indirect, competitive, or noncompetitive as described in the art (see, e.g., Principles and Practice of Immunoassay, 1991, Eds. Christopher P. Price and David J. Neoman, Stockton Pres, NY, NY; and Oellirich, M, 1984, J. Clin. Chem. Clin. Biochem, 22:895-904). Proteins may be isolated from test specimens and biological samples by conventional methods, as described in Current Protocols in Molecular Biology, supra. Exemplary antibody molecules for use in the detection and therapy methods of the present invention are intact immunoglobulin molecules, substantially intact immunoglobulin molecules, or those portions of immunoglobulin molecules that contain the antigen binding site. Polyclonal or monoclonal antibodies may be produced by methods conventionally known in the art (e.g., Kohler and

Milstein, 1975, Nature, 256:495-497; Campbell Monoclonal Antibody Technology, the Production and

Characterization of Rodent and Human Hybridomas, 1985, In: Laboratory Techniques in Biochemistry and Molecular Biology, Eds. Burdon et al. Volume 13, Elsevier Science Publishers, Amsterdam). The antibodies or antigen binding fragments thereof may also be produced by genetic engineering. The technology for expression of both heavy and light chain genes in E. coli is the subject of PCT patent applications, publication number WO 901443, WO 901443 and WO 9014424 and in Huse et al, 1989, Science, 246:1275-1281. The antibodies may also be humanized (e.g., Queen, C. et al. 1989 Proc. Natl. Acad. Sci.USA 86; 10029).

Effect(s) of the polymoφhisms identified herein on expression of MMP 13 may be investigated by preparing recombinant cells and/or nonhuman recombinant organisms, preferably recombinant animals, containing a polymoφhic variant of the MMP 13 gene. As used herein, "expression" includes but is not limited to one or more of the following: transcription of the gene into precursor mRNA; splicing and other processing of the precursor mRNA to produce mature mRNA; mRNA stability; translation of the mature mRNA into MMP 13 protein (including codon usage and tRNA availability); and glycosylation and/or other modifications of the translation product, if required for proper expression and function. To prepare a recombinant cell of the invention, the desired MMP 13 isogene may be introduced into the cell in a vector such that the isogene remains extrachromosomal. In such a situation, the gene will be expressed by the cell from the extrachromosomal location. In a preferred embodiment, the MMP 13 isogene is introduced into a cell in such a way that it recombines with the endogenous MMP 13 gene present in the cell. Such recombination requires the occurrence of a double recombination event, thereby resulting in the desired MMP 13 gene polymoφhism. Vectors for the introduction of genes both for recombination and for extrachromosomal maintenance are known in the art, and any suitable vector or vector construct may be used in the invention. Methods such as electroporation, particle bombardment, calcium phosphate co-precipitation and viral transduction for introducing DNA into cells are known in the art; therefore, the choice of method may lie with the competence and preference of the skilled practitioner. Examples of cells into which the MMP 13 isogene may be introduced include, but are not limited to, continuous culture cells, such as COS, NIH/3T3, and primary or culture cells of the relevant tissue type, i.e., they express the MMP13 isogene. Such recombinant cells can be used to compare the biological activities of the different protein variants.

Recombinant nonhuman organisms, i.e., transgenic animals, expressing a variant MMP13 gene are prepared using standard procedures known in the art. Preferably, a construct comprising the variant gene is introduced into a nonhuman animal or an ancestor of the animal at an embryonic stage, i.e., the one-cell stage, or generally not later than about the eight-cell stage. Transgenic animals carrying the constructs of the invention can be made by several methods known to those having skill in the art. One method involves transfecting into the embryo a retrovirus constructed to contain one or more insulator elements, a gene or genes of interest, and other components known to those skilled in the art to provide a complete shuttle vector harboring the insulated gene(s) as a transgene, see e.g., U.S. Patent No.

5,610,053. Another method involves directly injecting a transgene into the embryo. A third method involves the use of embryonic stem cells. Examples of animals into which the MMP 13 isogenes may be introduced include, but are not limited to, mice, rats, other rodents, and nonhuman primates (see "The Introduction of Foreign Genes into Mice" and the cited references therein, In: Recombinant DNA, Eds. J.D. Watson, M. Gilman, J. Witkowski, and M. Zoller; W.H. Freeman and Company, New York, pages 254-272). Transgenic animals stably expressing a human MMP13 isogene and producing human MMP 13 protein can be used as biological models for studying diseases related to abnormal MMP 13 expression and/or activity, and for screening and assaying various candidate drugs, compounds, and treatment regimens to reduce the symptoms or effects of these diseases.

An additional embodiment of the invention relates to pharmaceutical compositions for treating disorders affected by expression or function of a novel MMP 13 isogene described herein. The pharmaceutical composition may comprise any of the following active ingredients: a polynucleotide comprising one of these novel MMP 13 isogenes; an antisense oligonucleotide directed against one of the novel MMP 13 isogenes, a polynucleotide encoding such an antisense oligonucleotide, or another compound which inhibits expression of a novel MMP 13 isogene described herein. Preferably, the composition contains the active ingredient in a therapeuticaUy effective amount. By therapeuticaUy effective amount is meant that one or more of the symptoms relating to disorders affected by expression or function of a novel MMP 13 isogene is reduced and/or eliminated. The composition also comprises a pharmaceutically acceptable carrier, examples of which include, but are not limited to, saline, buffered saline, dextrose, and water. Those skilled in the art may employ a formulation most suitable for the active ingredient, whether it is a polynucleotide, oligonucleotide, protein, peptide or small molecule antagonist. The phannaceutical composition may be administered alone or in combination with at least one other agent, such as a stabilizing compound. Administration of the pharmaceutical composition may be by any number of routes including, but not limited to oral, intravenous, intramuscular, intra- arterial, intramedullary, intrathecal, intraventricular, intradermal, transdermal, subcutaneous, intraperitoneal, intranasal, enteral, topical, sublingual, or rectal. Further details on techniques for formulation and administration may be found in the latest edition of Remington's Pharmaceutical Sciences (Maack Publishing Co, Easton, PA).

For any composition, determination of the therapeuticaUy effective dose of active ingredient and/or the appropriate route of administration is well within the capability of those skilled in the art. For example, the dose can be estimated initially either in cell culture assays or in animal models. The animal model may also be used to determine the appropriate concentration range and route of administration. Such information can then be used to determine useful doses and routes for administration in humans. The exact dosage will be determined by the practitioner, in light of factors relating to the patient requiring treatment, including but not limited to severity of the disease state, general health, age, weight and gender of the patient, diet, time and frequency of administration, other drugs being taken by the patient, and tolerance/response to the treatment.

Any or all analytical and mathematical operations involved in practicing the methods of the present invention may be implemented by a computer. In addition, the computer may execute a program that generates views (or screens) displayed on a display device and with which the user can interact to view and analyze large amounts of information relating to the MMP 13 gene and its genomic variation, including chromosome location, gene structure, and gene family, gene expression data, polymoφhism data, genetic sequence data, and clinical data population data (e.g., data on ethnogeographic origin, clinical responses, genotypes, and haplotypes for one or more populations). The MMP13 polymoφhism data described herein may be stored as part of a relational database (e.g., an instance of an Oracle database or a set of ASCII flat files). These polymoφhism data may be stored on the computer's hard drive or may, for example, be stored on a CD-ROM or on one or more other storage devices accessible by the computer. For example, the data may be stored on one or more databases in communication with the computer via a network.

Preferred embodiments of the invention are described in the following examples. Other embodiments within the scope of the claims herein will be apparent to one skilled in the art from consideration of the specification or practice of the invention as disclosed herein. It is intended that the specification, together with the examples, be considered exemplary only, with the scope and spirit of the invention being indicated by the claims which follow the examples.

EXAMPLES The Examples herein are meant to exemplify the various aspects of carrying out the invention and are not intended to limit the scope of the invention in any way. The Examples do not include detailed descriptions for conventional methods employed, such as in the performance of genomic DNA isolation, PCR and sequencing procedures. Such methods are well-known to those skilled in the art and are described in numerous publications, for example, Sambrook, Fritsch, and Maniatis, "Molecular

Cloning: A Laboratory Manual", 2^nd Edition, Cold Spring Harbor Laboratory Press, USA, (1989).

EXAMPLE 1 This example illustrates examination of various regions of the MMP 13 gene for polymoφhic sites.

Amplification of Target Regions

The following target regions were amplified using the PCR primers represented below. The nucleotide positions of the first and last nucleotide of the forward and reverse primers for each region amplified are presented below and correspond to nucleotide positions in Figures 1, 2, and 3. PCR Primer Pairs

Fragment No. Forward Primer Reverse Primer PCR Product

Fragment 1 3395-3417 complement of 4064-4042 670 nt

Fragment 2 3763-3784 complement of 4425-4402 663 nt

Fragment 3 3799-3821 complement of 4346-4324 548 nt

Fragment 4 4052-4072 complement of 4783-4760 732 nt

Fragment 5 4753-4778 complement of 5542-5521 790 nt

Fragment 6 5182-5205 complement of 5890-5869 709 nt

Fragment 7 7172-7194 complement of 7833-7813 662 nt

Fragment 8 9294-9316 complement of 9956-9935 663 nt

Fragment 9 10330-10350 complement of 10864-10845 535 nt

Fragment 10 11533-11555 complement of 12089-12069 557 nt

Fragment 11 13732-13752 complement of 14363-14343 632 nt

Fragment 12 17573-17595 complement of 18252-18229 680 nt

Fragment 13 17923-17947 complement of 18646-18623 724 nt

Fragment 14 18236-18258 complement of 18881-18862 646 nt

Fragment 15 18508-18532 complement of 19146-19121 639 nt

These primer pairs were used in PCR reactions containing genomic DNA isolated from immortalized cell lines for each member of the Index Repository. The PCR reactions were carried out under the following conditions:

Reaction volume = 10 μl

10 x Advantage 2 Polymerase reaction buffer (Clontech) = l μl

100 ng of human genomic DNA = l μl lO mM dNTP = 0.4 μl

Advantage 2 Polymerase enzyme mix (Clontech) = 0.2 μl

Forward Primer (10 μM) = 0.4 μl

Reverse Primer (10 μM) = 0.4 μl

Water = 6.6μl

Amplification profile: 97°C - 2 min. 1 cycle

97°C - 15 sec. 70°C - 45 sec. 10 cycles 72°C - 45 sec.

97°C - 15 sec. 64°C - 45 sec. 35 cycles 72°C - 45 sec.

Sequencing of PCR Products

The PCR products were purified using a Whatman/Polyfiltronics 100 μl 384 well unifilter plate essentially according to the manufacturers protocol. The purified DNA was eluted in 50 μl of distilled water. Sequencing reactions were set up using Applied Biosystems Big Dye Terminator chemistry essentially according to the manufacturers protocol. The purified PCR products were sequenced in both directions using either the primer sets represented below with the positions of their first and last nucleotide corresponding to positions in Figures 1, 2, and 3. Reaction products were purified by isopropanol precipitation, and run on an Applied Biosystems 3700 DNA Analyzer.

Sequencing Primer Pairs

Fragment No. Forward Primer Reverse Primer Fragment 1 3478-3499 complement of 4032-4013 Fragment 2 3844-3864 complement of 4339-4320 Fragment 3 3855-3874 complement of 4272-4253 Fragment 4 4082-4101 complement of 4612-4592 Fragment 5 4915-4936 complement of 5454-5435 Fragment 6 5225-5242 complement of 5794-5773 Fragment 7 7235-7254 complement of 7799-7780 Fragment 8 9369-9388 complement of 9824-9804 Fragment 9 10439-10459 complement of 10850-10829 Fragment 10 11538-11559 complement of 11944-11924 Fragment 11 13749-13768 complement of 14203-14184 Fragment 12 17776-17793 complement of 18245-18224 Fragment 13 18044-18063 complement of 18526-18505 Fragment 14 18256-18275 complement of 18784-18765 Fragment 15 18553-18574 complement of 19094-19075

Analysis of Sequences for Polymorphic Sites Sequence information for a minimum of 80 humans was analyzed for the presence of polymoφhisms using the Polyphred program (Nickerson et al. Nucleic Acids Res. 14:2745-2751, 1997). The presence of a polymoφhism was confirmed on both strands. The polymoφhisms and their locations in the MMP 13 gene are listed in Table 3 below.

Table 3. Polymoφhic Sites Identified in the MMP 13 Gene

Polymoφhic Nucleotide Reference Variant CDS Variant AA

Site Number PolyId^a Position Allele Allele Position Variant

PS1 959081 3915 A G

PS2 959087 4437 G T 326 R109L

PS3 959093 5008 C T

PS4 959102 5037 G A

PS5 959104 5102 C G

PS6 959110 5290 C G

PS7 959112 5363 C A

PS8 959114 5628 G A

PS9 959116 7721 G A

PS10 959118 10537 C T

PSl l 959151 11663 T C 1080 Y360Y

PS12 959153 11752 A G 1169 D390G

PS13 959157 14059 C T

PS14 959159 14072 G C

PS15 959135 18295 T C

PS16 959133 18407 A G

PS17 959131 18609 A G

PS18 959129 18866 A G ^aPolyId is a unique identifier assigned to each PS by Genaissance Pharmaceuticals, Inc.

EXAMPLE 2 This example illustrates analysis of the MMP 13 polymoφhisms identified in the Index

Repository for human genotypes and haplotypes.

The different genotypes containing these polymoφhisms that were observed in the reference population are shown in Table 4 below, with the haplotype pair indicating the combination of haplotypes determined for the individual using the haplotype derivation protocol described below. In Table 4, homozygous positions are indicated by one nucleotide and heterozygous positions are indicated by two nucleotides. Missing nucleotides in any given genotype in Table 4 were inferred based on linkage disequilibrium and/or Mendelian inheritance.

Table 4 (Parti) . Genotypeis and Haplotype Pairs Observed for MMP13 Gem

Genotype Polymorphic Sites

Number PS1 PS2 PS3 PS4 PS5 PS6 PS7 PS8 PS9 PS10 HAP Pair

1 A G C G C C C G G C 3 3

2 G G C G C C A A G C 14 14

3 A G/T C G C C C G G C 3 8

4 A/G G C G C C C/A G/A G/A C 3 12

5 A G C G C C C G G C 3 4

6 G G C G/A C C A A G/A C 14 9

7 A G C G C C C G G/A C 3 2

8 A/G G C/T G C C C/A G/A G C 3 15

9 A G C G C C C G G C 3 5

10 A/G G C G C c C/A G/A G/A C 3 11

11 A/G G C G C c C/A G/A A C 1 11

12 A/G G C G C c C/A G/A G/A C 3 10

13 G/A G C G C c A/C A/G G C 14 4

14 G/A G c G C c A/C A/G G C 14 5

17 G/A G c G C c A/C A/G G/A C 14 1

18 A/G G c G C c C/A G/A G C 3 14

19 A G c G C c C G G/A C 3 1

20 A/G G c G/A C c C/A G/A G/A c 3 9

21 A/G G c G C c C/A G/A G/A C/T 3 13

22 A G c G C/G C/G C G G c 3 7

23 A G c G C c C G G c 3 6

Table 4(Part2) . Genotypes and Haplotype Pairs Observed for M P13 Gene Genotype Polymorphic Sites

Number PSll PS12 PS13 PS14 PS15 PS16 PS17 PS18 HAP Pair

1 T A C G T A A A 3 3

2 T A C G T A A G 14 14

3 T A C G T A A A 3 8

4 T A/G C G T A A A 3 12

5 T A C G T A A A/G 3 4

6 T A C G T/C A A G/A 14 9

7 T A C G T A A A/G 3 2

8 T A C G T A A A/G 3 15

9 T A C G T A A/G A 3 5

10 T A C G T A A A 3 11

11 T A C G T A A A 1 11

12 T/C A C G T A A A 3 10

13 T A C G T A A G 14 4

14 T A C G T A A/G G/A 14 5

15 T A C G T A A G/A 14 11

16 T A C G T A A G/A 14 13

17 T A C G T A A G/A 14 1

18 T A C G T A A A/G 3 14

19 T A C G T A A A 3 1

20 T A C G T/C A A A 3 9

21 T A C G T A A A 3 13

22 T A C/T G/C T A A A 3 7

23 G T A/G A

The haplotype pairs shown in Table 4 were estimated from the unphased genotypes using a computer-implemented extension of Clark's algorithm (Clark, A.G. 1990 Mol Bio Evol 7, 111-122) for assigning haplotypes to unrelated individuals in a population sample, as described in U.S. Provisional Application Serial No. 60/198,340 entitled "A Method and System for Determining Haplotypes from a Collection of Polymoφhisms" and the corresponding International Application, PCT/US01/12831. In this method, haplotypes are assigned directly from individuals who are homozygous at all sites or heterozygous at no more than one of the variable sites. This list of haplotypes is then used to deconvolute the unphased genotypes in the remaining (multiply heterozygous) individuals. In our analysis, the list of haplotypes was augmented with haplotypes obtained from two families (one three- generation Caucasian family and one two-generation African- American family).

By following this protocol, it was determined that the Index Repository examined herein and, by extension, the general population contains the 15 human MMP 13 haplotypes shown in Table 5 below.

A MMP 13 isogene defined by a full-haplotype shown in Table 5 below comprises the regions of the SEQ ID NOS indicated in Table 5, with their corresponding set of polymoφhic locations and identities, which are also set forth in Table 5. Table 5. Haplotypes Identified in the M P13 Gene

Haplotype i Number^a PS PS SEQ Region

1 2 3 4 5 6 7 8 9 10 No^b Pos° ID NO.^d Examined⁶

A A A A A A A A G G 1 3915/30 1/96 3395-5890

G G G G G G G T G G 2 4437/150 1/96 3395-5890

C C C C C C C C C C 3 5008/270 1/96 3395-5890

G G G G G G G G A G 4 5037/390 1/96 3395-5890

C C C C C C G C C C 5 5102/510 1/96 3395-5890

C C C C C C G C C C 6 5290/630 1/96 3395-5890

C C C C C C C C A A 7 5363/750 1/96 3395-5890

G G G G G G G G A A 8 5628/870 1/96 3395-5890

A A G G G G G G A A 9 7721/990 1/96 7172-7833

9294-9956

C C C C C C C C C C 10 10537/1110 1/96 10330-10864

T T T T T T T T T C 11 168/1230 2/96 38-594

A A A A A A A A A A 12 257/1350 2/96 38-594

C C C C C C T C C C 13 2564/1470 2/96 2237-2868

G G G G G G C G G G 14 2577/1590 2/96 2237-2868

T T T T T T T T C T 15 3436/1710 3/96 2714-4287

A A A A A G A A A A 16 3548/1830 3/96 2714-4287

A A A A G A A A A A 17 3750/1950 3/96 2714-4287

A G A G A A A A A A 18 4007/2070 3/96 2714-4287

Haplotype Number^a PS S SEQ Region

11 12 13 14 15 No^b Pos^c ID NO. . ^d Examined⁵

G G G G G 1 3915/30 1/96 3395- -5890

G G G G G 2 4437 /120 1/96 3395- -5890

C C C C T 3 5008/270 1/96 3395- -5890

G G G G G 4 5037/390 1/96 3395- -5890

C C C C C 5 5102/510 1/96 3395- -5890

C C C C C 6 5290/630 1/96 3395- -5890

A A A A A 7 5363/750 1/96 3395- -5890

A A A A A 8 5628/870 1/96 3395- -5890

A A A G G 9 7721/990 1/96 7172- -7833

9294- -9956

C C T C C 1C 1 10537/1110 1/96 10330-10864

T T T T T 11 168/ 1230 2/96 38-594

A G A A A 12 257/1350 2/96 38-594

C C C C C 13 2564/1470 2/96 2237- -2868

G G G G G 14 2577 /1590 2/96 2237- -2868

T T T T T IE 3436/1710 3/96 2714- -4287

A A A A A 16 3548/1830 3/96 2714- -4287

A A A A A 17 3750/1950 3/96 2714- -4287

A A A G G 18 4007 /2070 3/96 2714- -4287 ^aAlleles for MMP 13 haplotypes are presented 5' to 3' in each column ^bPS = polymoφhic site; ^cPosition of PS within the indicated SEQ ID NO, with the Imposition number referring to the first SEQ ID NO and the 2^nd position number referring to the 2^nd SEQ ID NO; ^dl^s SEQ ID NO refers to Figures 1, 2, or 3, with the two alternative allelic variants of each polymoφhic site indicated by the appropriate nucleotide symbol; 2^nd SEQ ID NO is a modified version of the 1^st SEQ ID NO that comprises the context sequence of each polymoφhic site, PS 1 -PS 18, to facilitate electronic searching of the haplotypes; ^eRegion examined represents the nucleotide positions defining the start and stop positions within the 1^st SEQ ID NO of the sequenced region. SEQ ID NOS: 1-3 refer to Figures 1-3, respectively, with the two alternative allelic variants of each polymoφhic site indicated by the appropriate nucleotide symbol. SEQ ID NO:96 is a modified version of SEQ ID NOS: 1-3 that shows the context sequence of each of PS 1 -PS 18 in a uniform format to facilitate electronic searching of the MMP13 haplotypes. For each polymoφhic site, SEQ ID NO:96 contains a block of 60 bases of the nucleotide sequence encompassing the centrally-located polymoφhic site at the 30^th position, followed by 60 bases of unspecified sequence to represent that each polymoφhic site is separated by genomic sequence whose composition is defined elsewhere herein. Table 6 below shows the percent of chromosomes characterized by a given MMP 13 haplotype for all unrelated individuals in the Index Repository for which haplotype data was obtained. The percent of these unrelated individuals who have a given MMP13 haplotype pair is shown in Table 7. In Tables 6 and 7, the "Total" column shows this frequency data for all of these unrelated individuals, while the other columns show the frequency data for these unrelated individuals categorized according to their self-identified ethnogeographic origin. Abbreviations used in Tables 6 and 7 are AF = African

Descent, AS = = Asian, CA = Caucasian, HL = Hispanic- Latino, and AM = Native American,

Table 6. Frequency of Observed MMP 13 Haplotypes In Unrelated Individuals

HAP No. HAP ID Total CA AF AS HL AM

1 962527 5.49 4.76 2.5 10.0 5.56 0.0

2 962539 0.61 2.38 0.0 0.0 0.0 0.0

3 962525 54.27 61.9 50.0 47.5 52.78 83.33

4 962531 1.83 0.0 7.5 0.0 0.0 0.0

5 962533 1.22 0.0 5.0 0.0 0.0 0.0

6 962537 0.61 0.0 2.5 0.0 0.0 0.0

7 962534 0.61 0.0 2.5 0.0 0.0 0.0

8 962538 0.61 0.0 0.0 0.0 2.78 0.0

9 962528 4.27 2.38 0.0 0.0 16.67 0.0

10 962536 0.61 0.0 2.5 0.0 0.0 0.0

11 962530 3.05 0.0 12.5 0.0 0.0 0.0

12 962535 0.61 0.0 2.5 0.0 0.0 0.0

13 962529 3.66 2.38 0.0 12.5 0.0 0.0

14 962526 21.34 23.81 12.5 30.0 19.44 16.67

15 962532 1.22 2.38 0.0 0.0 2.78 0.0

Table 7. Frequency of Observed MMP 13 Haplotype Pairs In U

HAP1 HAP2 Total CA AF AS HL AM

3 3 28.05 38.1 25.0 20.0 22.22 66.67

14 14 4.88 9.52 0.0 10.0 0.0 0.0

3 8 1.22 0.0 0.0 0.0 5.56 0.0

3 12 1.22 0.0 5.0 0.0 0.0 0.0

3 4 1.22 0.0 5.0 0.0 0.0 0.0

14 9 2.44 0.0 0.0 0.0 11.11 0.0

3 2 1.22 4.76 0.0 0.0 0.0 0.0

3 15 2.44 4.76 0.0 0.0 5.56 0.0

3 5 1.22 0.0 5.0 0.0 0.0 0.0

3 11 3.66 0.0 15.0 0.0 0.0 0.0

1 11 1.22 0.0 5.0 0.0 0.0 0.0

3 10 1.22 0.0 5.0 0.0 0.0 0.0

14 4 2.44 0.0 10.0 0.0 0.0 0.0

14 5 1.22 0.0 5.0 0.0 0.0 0.0

14 11 1.22 0.0 5.0 0.0 0.0 0.0

14 13 2.44 0.0 0.0 10.0 0.0 0.0

14 1 3.66 4.76 0.0 5.0 5.56 0.0

3 14 19.51 23.81 5.0 25.0 22.22 33.33

3 1 6.1 4.76 0.0 15.0 5.56 0.0

3 9 6.1 4.76 0.0 0.0 22.22 0.0

3 13 4.88 4.76 0.0 15.0 0.0 0.0

3 7 1.22 0.0 5.0 0.0 0.0 0.0

3 6 1.22 0.0 5.0 0.0 0.0 0.0

The size and composition of the Index Repository were chosen to represent the genetic diversity across and within four major population groups comprising the general United States population. For example, as described in Table 1 above, this repository contains approximately equal sample sizes of African-descent, Asian- American, European-American, and Hispanic-Latino population groups. Almost all individuals representing each group had all four grandparents with the same ethnogeographic background. The number of unrelated individuals in the Index Repository provides a sample size that is sufficient to detect SNPs and haplotypes that occur in the general population with high statistical certainty. For instance, a haplotype that occurs with a frequency of 5% in the general population has a probability higher than 99.9% of being observed in a sample of 80 individuals from the general population. Similarly, a haplotype that occurs with a frequency of 10% in a specific population group has a 99% probability of being observed in a sample of 20 individuals from that population group. In addition, the size and composition of the Index Repository means that the relative frequencies determined therein for the haplotypes and haplotype pairs of the MMP 13 gene are likely to be similar to the relative frequencies of these MMP13 haplotypes and haplotype pairs in the general U.S. population and in the four population groups represented in the Index Repository. The genetic diversity observed for the three Native Americans is presented because it is of scientific interest, but due to the small sample size it lacks statistical significance. In view of the above, it will be seen that the several advantages of the invention are achieved and other advantageous results attained.

As various changes could be made in the above methods and compositions without departing from the scope of the invention, it is intended that all matter contained in the above description and shown in the accompanying drawings shall be inteφreted as illustrative and not in a limiting sense.

All references cited in this specification, including patents and patent applications, are hereby incoφorated in their entirety by reference. The discussion of references herein is intended merely to summarize the assertions made by their authors and no admission is made that any reference constitutes prior art. Applicants reserve the right to challenge the accuracy and pertinency of the cited references.

Claims

What is Claimed is:

1. A method for haplotyping the matrix metalloproteinase 13 (collagenase 3) (MMP13) gene of an individual, which comprises determining which of the MMP 13 haplotypes shown in the table immediately below defines one copy of the individual's MMP13 gene, wherein each of the MMP 13 haplotypes comprises a set of polymoφhisms whose locations and identities are set forth in the table immediately below:

Haplot}φe Number .a PS PS SEQ ID

1 2 3 4 5 6 7 8 9 10 No.^b Pos.^c NO.

A A A A A A A A G G 1 3915 1

G G G G G G G T G G 2 4437 1

C C C C C C C C C C 3 5008 1

G G G G G G G G A G 4 503 1

C C C C C C G C C C 5 5102 1

C C C C C C G C C C 6 5290 1

C C C C C C C C A A 7 5363 1

G G G G G G G G A A 8 5628 1

A A G G G G G G A A 9 7721 1

C C C C C C C C C C 10 10537 1

T T T T T T T T T C 11 168 2

A A A A A A A A A A 12 257 2

C C C C C C T C C C 13 2564 2

G G G G G G C G G G 14 2577 2

T T T T T T T T C T 15 3436 3

A A A A A G A A A A 16 3548 3

A A A A G A A A A A 17 3750 3

A G A G A A A A A A 18 4007 3

Haplotype Numbei .a PS PS SEQ ID

11 12 13 14 15 No.^b Pos.^c NO.

G G G G G 1 3915 1

G G G G G 2 4437 1

C C C C T 3 5008 1

G G G G G 4 5037 1

C C C C C 5 5102 1

C C C C C 6 5290 1

A A A A A 7 5363 1

A A A A A 8 5628 1

A A A G G 9 7721 1

C C T C C 10 10537 1

T T T T T 11 168 2

A G A A A 12 257 2

C C C C C 13 2564 2

G G G G G 14 2577 2

T T T T T 15 3436 3

A A A A A 16 3548 3

A A A A A 17 3750 3

A A A G G 18 4007 3 ^aAlleles for haplotypes are presented 5 ' to 3 ' in each column ^bPS = polymoφhic site; "Position of PS within the indicated SEQ ID NO.

2. The method of claim 1, wherein the determining step comprises identifying the phased sequence of nucleotides present at each of PS 1 -PS 18 on the one copy of the individual's MMP 13 gene.

3. A method for haplotyping the matrix metalloproteinase 13 (collagenase 3) (MMP 13) gene of an individual, which comprises determining which of the MMP 13 haplotype pairs shown in the table immediately below defines both copies of the individual's MMP 13 gene, wherein each of the MMP 13 haplotype pairs consists of first and second haplotypes which comprise first and second sets of polymoφhisms whose locations and identities are set forth in the table immediately below:

Haplotype Pair^a PS PS SEQ ID

3/3 14/14 3/8 3/12 3/4 14/9 3/2 3/15 No.^b Pos.' NO.

A/A G/G A/A A/G A/A G/G A/A A/G 1 3915 1

G/G G/G G/T G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/T 3 5008 1

G/G G/G G/G G/G G/G G/A G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C A/A C/C C/A C/C A/A C/C C/A 7 5363 1

G/G A/A G/G G/A G/G A/A G/G G/A 8 5628 1

G/G G/G G/G G/A G/G G/A G/A G/G 9 7721 1

C/C C/C C/C C/C C/C C/C C/C C/C 10 10537 1

T/T T/T T/T T/T T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/G A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/C T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/A A/A A/A A/A A/A A/A A/A A/A 17 3750 3

A/A G/G A/A A/A A/G G/A A/G A/G 18 4007 3

Haplotype Pair^a PS PS SEQ ID

3/5 3/11 1/11 3/10 14/4 14/E 5 14/11 14/13 No.^b Pos.^c NO.

A/A A/G A/G A/G G/A G/A G/G G/G 1 3915 1

G/G G/G G/G G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/C 3 5008 1

G/G G/G G/G G/G G/G G/G G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C C/A C/A C/A A/C A/C A/A A/A 7 5363 1

G/G G/A G/A G/A A/G A/G A/A A/A 8 5628 1

G/G G/A A/A G/A G/G G/G G/A G/A 9 7721 2

C/C C/C C/C C/C C/C C/C C/C C/T 10 10537 2

T/T T/T T/T T/C T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/A A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/T T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/G A/A A/A A/A A/A A/G A/A A/A 17 3750 3

A/A A/A A/A A/A G/G G/A G/A G/A 18 4007 3

Haplotype Pair^a PS PS SEQ ID

14/1 3/14 3/1 3/9 3/13 3/7 3/6 No. Pos.^c NO.

G/A A/G A/A A/G A/G A/A A/A 1 3915 1

G/G G/G G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C 3 5008 1

G/G G/G G/G G/A G/G G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/G C/C 5 5102 1

C/C C/C C/C C/C C/C C/G C/C 6 5290 1

A/C C/A C/C C/A C/A C/C C/C 7 5363 1

A/G G/A G/G G/A G/A G/G G/G 8 5628 1

G/A G/G G/A G/A G/A G/G G/G 9 7721 1

C/C C/C C/C C/C C/T C/C C/C 10 loss¹ 1

T/T T/T T/T T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/T C/C 13 2564 2

G/G G/G G/G G/G G/G G/C G/G 14 2577 2

T/T T/T T/T T/C T/T T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/G 16 3548 3

A/A A/A A/A A/A A/A A/A A/A 17 3750 3

G/A A/G A/A A/A A/A A/A A/A 18 4007 3

Ηaplotype pairs are represented as 1^st haplotype/2^nd haplotype; with alleles of each haplotype shown 5 ' to 3 ' as 1^st polymoφhism/2ⁿ polymoφhism in each column; ^bPS = polymoφhic site;

'Position of PS within the indicated SEQ ID NO.

4. The method of claim 3, wherein the determining step comprises identifying the phased sequence of nucleotides present at each of PS 1 -PS 18 on both copies of the individual's MMP 13 gene.

5. A method for genotyping the matrix metalloproteinase 13 (collagenase 3) (MMP 13) gene of an individual, comprising determining for the two copies of the MMP 13 gene present in the individual the identity of the nucleotide pair at one or more polymoφhic sites (PS) selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS 13 , PS 14, PS 15 , PS 16, PS 17 and PS 18, wherein the one or more PS have the location and alternative alleles shown in SEQ ID NOS: 1-3.

6. The method of claim 5, wherein the determining step comprises:

(a) isolating from the individual a nucleic acid mixture comprising both copies of the MMP 13 gene, or a fragment thereof, that are present in the individual;

(b) amplifying from the nucleic acid mixture a target region containing the selected polymoφhic site;

(c) hybridizing a primer extension oligonucleotide to one allele of the amplified target region;

(d) performing a nucleic acid template-dependent, primer extension reaction on the hybridized genotyping oligonucleotide in the presence of at least two different terminators of the reaction, wherein said terminators are complementary to the alternative nucleotides present at the selected polymoφhic site; and

(e) detecting the presence and identity of the terminator in the extended genotyping oligonucleotide.

7. The method of claim 5, which comprises determining for the two copies of the MMP13 gene present in the individual the identity of the nucleotide pair at each of PS 1 -PS 18.

8. A method for haplotyping the matrix metalloproteinase 13 (collagenase 3) (MMP 13) gene of an individual which comprises determining, for one copy of the MMP 13 gene present in the individual, the identity of the nucleotide at two or more polymoφhic sites (PS) selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSll, PS12, PS13, PS14, PS15, PS16, PS17 and PS18, wherein the selected PS have the location and alternative alleles shown in SEQ ID NOS: 1-3.

9. The method of claim 8, wherein the determining step comprises:

(a) isolating from the individual a nucleic acid sample containing only one of the two copies of the MMP 13 gene, or a fragment thereof, that is present in the individual;

(b) amplifying from the nucleic acid sample a target region containing the selected polymoφhic site;

(c) hybridizing a primer extension oligonucleotide to one allele of the amplified target region;

(d) performing a nucleic acid template-dependent, primer extension reaction on the hybridized genotyping oligonucleotide in the presence of at least two different terminators of the reaction, wherein said terminators are complementary to the alternative nucleotides present at the selected polymoφhic site; and

(e) detecting the presence and identity of the terminator in the extended genotyping oligonucleotide.

10. A method for predicting a haplotype pair for the matrix metalloproteinase 13 (collagenase 3) (MMP13) gene of an individual comprising:

(a) identifying a MMP 13 genotype for the individual, wherein the genotype comprises the nucleotide pair at two or more polymoφhic sites (PS) selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSll, PS12, PS13, PS14, PS15,

PS 16, PS 17 and PS 18, wherein the selected PS have the location and alternative alleles shown in SEQ ID NO:96;

(b) enumerating all possible haplotype pairs which are consistent with the genotype;

(c) comparing the possible haplotype pairs to the haplotype pair data set forth in the table immediately below; and

(d) assigning a haplotype pair to the individual that is consistent with the data

Haplotype Pair³ PS PS SEQ ID

3/3 14/14 3/8 3/12 3/4 14/9 3/2 3/15 No.^b Pos.^c NO.

A/A G/G A/A A/G A/A G/G A/A A/G 1 3915 1

G/G G/G G/T G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/T 3 5008 1

G/G G/G G/G G/G G/G G/A G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C A/A C/C C/A C/C A/A C/C C/A 7 5363 1

G/G A/A G/G G/A G/G A/A G/G G/A 8 5628 1

G/G G/G G/G G/A G/G G/A G/A G/G 9 7721 1

C/C C/C C/C C/C C/C C/C C/C C/C 10 10537 1

T/T T/T T/T T/T T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/G A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/C T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/A A/A A/A A/A A/A A/A A/A A/A 17 3750 3

A/A G/G A/A A/A A/G G/A A/G A/G 18 4007 3

Haplotype Pair PS PS SEQ ID

3/5 3/11 1/11 3/10 14/4 14/E 5 14/11 14/13 No.^b Pos.° NO.

A/A A/G A/G A/G G/A G/A G/G G/G 1 3915 1

G/G G/G G/G G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/C 3 5008 1

G/G G/G G/G G/G G/G G/G G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C C/A C/A C/A A/C A/C A/A A/A 7 5363 1

G/G G/A G/A G/A A/G A/G A/A A/A 8 5628 1

G/G G/A A/A G/A G/G G/G G/A G/A 9 7721 2

C/C C/C C/C C/C C/C C/C C/C C/T 10 10537 2

T/T T/T T/T T/C T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/A A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/T T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/G A/A A/A A/A A/A A/G A/A A/A 17 3750 3

A/A A/A A/A A/A G/G G/A G/A G/A 18 4007 3

Haplotype Pair^a PS PS SEQ ID

14/1 3/14 3/1 3/9 3/13 3/7 3/6 No.^b Pos.^c NO.

G/A A/G A/A A/G A/G A/A A/A 1 3915 1

G/G G/G G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C 3 5008 1

G/G G/G G/G G/A G/G G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/G C/C 5 5102 1

C/C C/C C/C C/C C/C C/G C/C 6 5290 1

A/C C/A C/C C/A C/A C/C C/C 7 5363 1

A/G G/A G/G G/A G/A G/G G/G 8 5628 1

G/A G/G G/A G/A G/A G/G G/G 9 7721 1

C/C C/C C/C C/C C/T C/C C/C 10 loss¹/ 1

T/T T/T T/T T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/T C/C 13 2564 2

G/G G/G G/G G/G G/G G/C G/G 14 2577 2

T/T T/T T/T T/C T/T T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/G 16 3548 3

A/A A/A A/A A/A A/A A/A A/A 17 3750 3

G/A A/G A/A A/A A/A A/A A/A 18 4007 3

„ 1st t,„„ ,„: - „π„ shown 5' to 3' as 1^st polymoφhism/2ⁿ polymoφhism in each column; ^bPS = polymoφhic site; 'Position of PS within the indicated SEQ ID NO.

11. The method of claim 10, wherein the identified genotype of the individual comprises the nucleotide pair at each of PS 1 -PS 18, which have the location and alternative alleles shown in SEQ ID NOS: 1-3.

12. A method for identifying an association between a trait and at least one haplotype or haplotype pair of the matrix metalloproteinase 13 (collagenase 3) (MMP 13) gene which comprises comparing the frequency of the haplotype or haplotype pair in a population exhibiting the trait with the frequency of the haplotype or haplotype pair in a reference population, wherein the haplotype is selected from haplotypes 1-15 shown in the table presented immediately below, wherein each of the haplotypes comprises a set of polymoφhisms whose locations and identities are set forth in the table immediately below:

Haplotype Num bei .a PS PS SEQ ID

1 2 3 4 5 6 7 8 9 10 No . ^b Pos . ^c NO .

A A A A A A A A G G 1 3915 1

G G G G G G G T G G 2 4437 1

C C C C C C C C C C 3 5008 1

G G G G G G G G A G 4 503 1

C C C C C C G C C C 5 5102 1

C C C C C C G C C C 6 5290 1

C C C C C C C C A A 7 5363 1

G G G G G G G G A A 8 5628 1

A A G G G G G G A A 9 7721 1

C C C C C C C C C C 10 10537 1

T T T T T T T T T C 11 168 2

A A A A A A A A A A 12 257 2

C C C C C C T C C C 13 2564 2

G G G G G G C G G G 14 2577 2

T T T T T T T T C T 15 3436 3

A A A A A G A A A A 16 3548 3

A A A A G A A A A A 17 3750 3

A G A G A A A A A A 18 4007 3

Haplotype Number^a PS PS SEQ ID

11 12 13 14 15 No . ^b Pos . ^c NO .

G G G G G 1 3915 1

G G G G G 2 4437 1

C C C C T 3 5008 1

G G G G G 4 5037 1

C C C C C 5 5102 1

C C C C C 6 5290 1

A A A A A 7 5363 1

A A A A A 8 5628 1

A A A G G 9 7721 1

C C T C C 10 1053¹/ ' 1

T T T T T 11 168 2

A G A A A 12 257 2

C C C C C 13 2564 2

G G G G G 14 2577 2

T T T T T 15 3436 3

A A A A A 16 3548 3

A A A A A 17 3750 3

A A A G G 18 4007 3 ^aAlleles for haplotypes are presented 5' to 3' in each column ^bPS = polymoφhic site;

"Position of PS within the indicated SEQ ID NO; and wherein the haplotype pair is selected from the haplotype pairs shown in the table immediately below, wherein each of the MMP 13 haplotype pairs consists of first and second haplotypes which comprise first and second sets of polymoφhisms whose locations and identities are set forth in the table immediately below:

Haplotype Pair^a S PS SEQ ID

3/3 14/14 3/8 3/12 3/4 14/9 3/2 3/15 No.^b Pos.^c NO.

A/A G/G A/A A/G A/A G/G A/A A/G 1 3915 1

G/G G/G G/T G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/T 3 5008 1

G/G G/G G/G G/G G/G G/A G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C A/A C/C C/A C/C A/A C/C C/A 7 5363 1

G/G A/A G/G G/A G/G A/A G/G G/A 8 5628 1

G/G G/G G/G G/A G/G G/A G/A G/G 9 7721 1

C/C C/C C/C C/C C/C C/C C/C C/C 10 10537 1

T/T T/T T/T T/T T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/G A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/C T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/A A/A A/A A/A A/A A/A A/A A/A 17 3750 3

A/A G/G A/A A/A A/G G/A A/G A/G 18 4007 3

Haplotype Pa Lir^a PS PS SEQ ID

3/5 3/11 1/11 3/10 14/4 14/5 14/11 14/13 No.^{1 5} Pos.^c NO.

A/A A/G A/G A/G G/A G/A G/G G/G 1 3915 1

G/G G/G G/G G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/C 3 5008 1

G/G G/G G/G G/G G/G G/G G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C C/A C/A C/A A/C A/C A/A A/A 7 5363 1

G/G G/A G/A G/A A/G A/G A/A A/A 8 5628 1

G/G G/A A/A G/A G/G G/G G/A G/A 9 7721 2

C/C C/C C/C C/C C/C C/C C/C C/T 10 10537 2

T/T T/T T/T T/C T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/A A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2 τ/τ T/T T/T T/T T/T T/T T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/G A/A A/A A/A A/A A/G A/A A/A 17 3750 3

A/A A/A A/A A/A G/G G/A G/A G/A 18 4007 3 105 Haplotype Paii_rr'^a PS PS SEQ ID

14/1 3/14 3/1 3/9 3/13 3 3_//7 3/6 No.^b Pos.^c NO.

G/A A/G A/A A/G A/G A A_//A A/A 1 3915 1

G G//GG G G//GG G/G G/G G/G G G_//G G/G 2 4437 1

C C//CC C C//CC C/C C/C C/C C C_//C C/C 3 5008 1

110 G G//GG G G//GG G/G G/A G/G G G_//G G/G 4 5037 1

C C//CC C C//CC C/C C/C C/C C C_//G C/C 5 5102 1

C C//CC C C//CC C/C C/C C/C C C_//G C/C 6 5290 1

A A//CC C C//AA C/C C/A C/A C C_//C C/C 7 5363 1

A A//GG G G//AA G/G G/A G/A G G_//G G/G 8 5628 1

115 G G//AA G G//GG G/A G/A G/A G G,/G G/G 9 7721 1

C C//CC C C//CC C/C C/C C/T C C_//C C/C 10 10537 1

T T//TT T T//TT T/T T/T T/T T T_//T T/T 11 168 2

A A//AA A A//AA A/A A/A A/A A A_//A A/A 12 257 2

C C//CC C C//CC C/C C/C C/C C C_//T C/C 13 2564 2

120 G G//GG G G//GG G/G G/G G/G G G//C G/G 14 2577 2

T T//TT T T//TT T/T T/C T/T T T_//T T/T 15 3436 3

A A//AA A A//AA A/A A/A A/A A A_//A A/G 16 3548 3

A A//AA A A//AA A/A A/A A/A A A_//A A/A 17 3750 3

G G//AA A A//GG A/A A/A A/A A A_//A A/A 18 4007 3

125

Ηaplotype pairs are represented as ι Is^stt h,a„p. ιl„o+t,y„p,_Qe//<2->nⁿd haplotype; with alleles of each haplotype shown 5 ' to 3 ' as 1^st polymoφhism/2^nd polymoφhism in each column; ^bPS = polymoφhic site;

"Position of PS within the indicated SEQ ID NO;

130 wherein a higher frequency of the haplotype or haplotype pair in the trait population than in the reference population indicates the trait is associated with the haplotype or haplotype pair.

13. The method of claim 12, wherein the trait is a clinical response to a drug targeting MMP13.

14. An isolated genotyping oligonucleotide for detecting a polymoφhism in the matrix metalloproteinase 13 (collagenase 3) (MMP13) gene at a polymoφhic site (PS) selected from the group consisting of PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS13, PS14, PS15, PS16, PS17 and PS18, wherein the selected PS have the location and alternative alleles shown in SEQ ID NOS:l-3.

15. The isolated genotyping oligonucleotide of claim 14, which is an allele-specific oligonucleotide that specifically hybridizes to an allele of the MMP 13 gene at a region containing the polymoφhic site.

16. The allele-specific oligonucleotide of claim 15, which comprises a nucleotide sequence selected from the group consisting of SEQ ID NOS:6-23, the complements of SEQ ID NOS:6-23, and SEQ ID NOS:24-59.

17. The isolated genotyping oligonucleotide of claim 14, which is a primer-extension oligonucleotide.

18. The primer-extension oligonucleotide of claim 17,which comprises a nucleotide sequence selected from the group consisting of SEQ ID NOS:60-95.

19. A kit for genotyping the matrix metalloproteinase 13 (collagenase 3) (MMP13) gene of an ■ individual, which comprises a set of oligonucleotides designed to genotype each of polymoφhic sites (PS) PS1, PS2, PS3, PS4, PS5, PS6, PS7, PS8, PS9, PS10, PSl l, PS12, PS13, PS14, PS15, PS 16, PS 17 and PS 18, wherein the selected PS have the location and alternative alleles shown in SEQ ID NOS: 1-3. 20. An isolated polynucleotide comprising a nucleotide sequence selected from the group consisting of: (a) a first nucleotide sequence which comprises a matrix metalloproteinase 13 (collagenase 3) (MMP 13) isogene, wherein the MMP 13 isogene is selected from the group consisting of isogenes 1- 2 and 4 - 15 shown in the table immediately below and wherein each of the isogenes comprises the regions of SEQ ID NOS: 1-3 shown in the table immediately below and wherein each of the isogenes 1- 2 and 4 - 15 is further defined by the corresponding set of polymoφhisms whose locations and identities are set forth in the^able immediately below

Esogenei Number³ PS PS SEQ ID Region

1 2 3 4 5 6 7 8 9 10 No.^b Pos.^c No.^d Examined³

A A A A A A A A G G 1 3915 1 3395-5890

G G G G G G G T G G 2 4437 1 3395-5890

C C C C C C C C C C 3 5008 1 3395-5890

G G G G G G G G A G 4 5037 1 3395-5890

C C C C C C G C C C 5 5102 1 3395-5890

C C C C C C G C C C 6 5290 1 3395-5890

C C C C C C C C A A 7 5363 1 3395-5890

G G G G G G G G A A 8 5628 1 3395-5890

A A G G G G G G A A 9 7721 1 7172-7833 9294-9956

C C C C C C C C C C 10 10537 1 10330-10864

T T T T T T T T T C 11 168 2 38-594

A A A A A A A A A A 12 257 2 38-594

C C C C C C T C C C 13 2564 2 2237-2868

G G G G G G C G G G 14 2577 2 223-7-2868

T T T T T T T T C T 15 3436 3 2714-4287

A A A A A G A A A A 16 3548 3 2714-4287

A A A A G A A A A A 17 3750 3 2714-4287

A G A G A A A A A A 18 4007 3 2714-4287

Isogene s Number³ PS PS SEQ ID Region

11 12 13 14 15 No . ^b Pos . ^c No . ^{d 1} Examined⁶

G G G G G 1 3915 1 3395 -5890

G G G G G 2 4437 1 3395-5890

C C C C T 3 5008 1 3395 - 5890

G G G G G 4 5037 1 3395 - 5890

C C C C C 5 5102 1 3395 - 5890

C C C C C 6 5290 1 3395 - 5890

A A A A A 7 5363 1 3395-5890

A A A A A 8 5628 1 3395-5890

A A A G G 9 7721 1 7172 - 7833 9294 - 9956

C C T C C 10 10537 1 10330 -10864

T T T T T 11 168 2 38 -594

A G A A A 12 257 2 38 - 594

C C C C C 13 2564 2 2237-2868

G G G G G 14 2577 2 2237-2868

T T T T T 15 3436 3 2714 -4287

A A A A A 16 3548 3 2714 -4287

A A A A A 17 3750 3 2714 -4287

A A A G G 18 4007 3 2714 -4287 ^aAlleles for isogenes are presented 5 ' to 3 ' in each column ^bPS = polymoφhic site;

"Position of PS within the indicated SEQ ID NO; ^dSEQ ID NO refers to Figures 1, 2, or 3, with the two alternative allelic variants of each polymoφhic site indicated by the appropriate nucleotide symbol; "Region examined represents the nucleotide positions defining the start and stop positions within the indicated SEQ ID NO of the sequenced region.

(b) a second nucleotide sequence which comprises a fragment of the first nucleotide sequence, wherein the fragment comprises one or more polymoφhisms selected from the group consisting of guanine at PS1, thymine at PS2, thymine at PS3, adenine at PS4, guanine at PS5, guanine at PS6, adenine at PS7, adenine at PS8, adenine at PS9, thymine at PS10, cytosine at PS11, guanine at PS 12, thymine at PS 13, cytosine at PS 14, cytosine at PS 15, guanine at PS16, guanine at PS17 and guanine at PS18, wherein the selected polymoφhism has the location set forth in the table immediately above; and

(c) a third nucleotide sequence which is complementary to the first or second nucleotide sequence.

21. The isolated polynucleotide of claim 20, which is a DNA molecule and comprises both the first and third nucleotide sequences and further comprises expression regulatory elements operably linked to the first nucleotide sequence.

22. A recombinant nonhuman organism transformed or transfected with the isolated polynucleotide of claim 20, wherein the organism expresses a MMP 13 protein encoded by the first nucleotide sequence.

23. The recombinant nonhuman organism of claim 22, which is a transgenic animal.

24. The isolated polynucleotide of claim 20 which consists of the second nucleotide sequence.

25. An isolated polynucleotide comprising a nucleotide sequence selected from the group consisting of:

(a) a coding sequence for a matrix metalloproteinase 13 (collagenase 3) (MMP 13) isogene wherein the coding sequence is defined by a haplotype selected from the group consisting of 8c, 10c, and 12c shown in the table immediately below and wherein the coding sequence comprises SEQ ID NO:4 except at each of the polymoφhic sites which have the locations and polymoφhisms set forth in the table immediately below:

Coding Sequence Haplotype Number³ PS PS

8c 10c 12c No.^b Pos.^c

T G G 2 326

T C T 11 1080

A A G 12 1169

"Alleles for coding sequence haplotypes are presented 5 ' to 3 ' in each column; the numerical portion of the coding sequence haplotype number represents the number of the parent full MMP 13 haplotype; ^bPS = polymoφhic site; "Position of PS in SEQ ID NO:4; and

(b) a fragment of the coding sequence, wherein the fragment comprises at least one polymoφhism selected from the group consisting of thymine at a position corresponding to nucleotide 326, cytosine at a position corresponding to nucleotide 1080 and guanine at a position corresponding to nucleotide 1169, wherein said positions in the coding sequence and the fragment refer to SEQ ID NO:4.

26. A recombinant nonhuman organism transformed or transfected with the isolated polynucleotide of claim 25, wherein the organism expresses a matrix metalloproteinase 13 (collagenase 3) (MMP 13) protein encoded by the polymoφhic variant sequence.

27. The recombinant nonhuman organism of claim 26, which is a transgenic animal.

28. An isolated polypeptide comprising an amino acid sequence which is a polymoφhic variant of a reference sequence for the matrix metalloproteinase 13 (collagenase 3) (MMP 13) protein or a fragment thereof, wherein the reference sequence comprises SEQ ID NO:5 and the polymoφhic variant comprises one or more variant amino acids selected from the group consisting of leucine at a position corresponding to amino acid position 109 and glycine at a position corresponding to amino acid position 390.

29. An isolated monoclonal antibody specific for and immunoreactive with the isolated polypeptide of claim 28.

30. A method for screening for drugs targeting the isolated polypeptide of claim 28 which comprises contacting the MMP 13 polymoφhic variant with a candidate agent and assaying for binding activity.

1. A computer system for storing and analyzing polymoφhism data for the matrix metalloproteinase 13 (collagenase 3) gene, comprising:

(a) a central processing unit (CPU);

(b) a communication interface;

(c) a display device;

(d) an input device; and

(e) a database containing the polymoφhism data; wherein the polymoφhism data comprises the haplotypes set forth in the table immediately below:

Haplotype Number³ PS PS SEQ ID

1 2 3 4 5 6 7 8 9 10 No.^b Pos.^c NO.

A A A A A A A A G G 1 3915 1

G G G G G G G T G G 2 4437 1

C C C C C C C C C C 3 5008 1

G G G G G G G G A G 4 503 1

C C C C C C G C C C 5 5102 1

C C C C C C G C C C 6 5290 1

C C C C C C C C A A 7 5363 1

G G G G G G G G A A 8 5628 1

A A G G G G G G A A 9 7721 1

C C C C C C C C C C 10 10537 1

T T T T T T T T T C 11 168 2

A A A A A A A A A A 12 257 2

C C C C C C T C C C 13 2564 2

G G G G G G C G G G 14 2577 2

T T T T T T T T C T 15 3436 3

A A A A A G A A A A 16 3548 3

A^" A A A G A A A A A 17 3750 3

A G A G A A A A A A 18 4007 3

Haplotype Number³ PS PS SEQ ID

11 12 13 14 15 No.^b Pos.^c NO.

G G G G G 1 3915 1

G G G G G 2 4437 1

C C C C T 3 5008 1

G G G G G 4 5037 1

C C C C C 5 5102 1

C C C C C 6 5290 1

A A A A A 7 5363 1

A A A A A 8 5628 1

A A A G G 9 7721 1

C C T C C 10 10537 ' ^• 1

T T T T T 11 168 2

A G A A A 12 257 2

C C C C C 13 2564 2

G G G G G 14 2577 2

T T T T T 15 3436 3

A A A A A 16 3548 3

A A A A A 17 3750 3

A A A G G 18 4007 3 ^aAlleles for haplotypes are presented 5' to 3' in each column ^bPS = polymoφhic site;

"Position of PS within the indicated SEQ ID NO; and the haplotype pairs set forth in the table immediately below:

Haplotype Pair³ PS PS SEQ ID

3/3 14/14 3/8 3/12 3/4 14/9 3/2 3/1E No.^b Pos.^c NO.

A/A G/G A/A A/G A/A G/G A/A A/G 1 3915 1

G/G G/G G/T G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C ^' C/C C/T 3 5008 1

G/G G/G G/G G/G G/G G/A G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C A/A C/C C/A C/C A/A C/C C/A 7 5363 1

G/G A/A^' G/G G/A G/G A/A G/G G/A 8 5628 1

G/G G/G G/G G/A G/G G/A G/A G/G 9 7721 1

C/C C/C C/C C/C C/C C/C C/C C/C 10 10537 1

T/T T/T T/T T/T T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/G A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/C T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/A A/A A/A A/A A/A A/A A/A A/A 17 3750 3

A/A G/G A/A A/A A/G G/A A/G A/G 18 4007 3

Haplotype Pair³ PS PS SEQ ID

3/5 3/11 1/11 3/10 14/4 14/5 14/11 14/13 No.^{1 5} Pos.^c NO.

A/A A/G A/G A/G G/A G/A G/G G/G 1 3915 1

G/G G/G G/G G/G G/G G/G G/G G/G 2 4437 1

C/C C/C C/C C/C C/C C/C C/C C/C 3 5008 1

G/G G/G G/G G/G G/G G/G G/G G/G 4 5037 1

C/C C/C C/C C/C C/C C/C C/C C/C 5 5102 1

C/C C/C C/C C/C C/C C/C C/C C/C 6 5290 1

C/C C/A C/A C/A A/C A/C A/A A/A 7 5363 1

G/G G/A G/A G/A A/G A/G A/A A/A 8 5628 1

G/G G/A A/A G/A G/G G/G G/A G/A 9 7721 2

C/C C/C C/C C/C C/C C/C C/C C/T 10 10537 2

T/T T/T T/T T/C T/T T/T T/T T/T 11 168 2

A/A A/A A/A A/A A/A A/A A/A A/A 12 257 2

C/C C/C C/C C/C C/C C/C C/C C/C 13 2564 2

G/G G/G G/G G/G G/G G/G G/G G/G 14 2577 2

T/T T/T T/T T/T T/T T/T T/T T/T 15 3436 3

A/A A/A A/A A/A A/A A/A A/A A/A 16 3548 3

A/G A/A A/A A/A A/A A/G A/A A/A 17 3750 3

A/A A/A A/A A/A G/G G/A G/A G/A 18 4007 3 100 Haplotype PaiLrr¹³ PS PS SEQ ID

14/1 3/14 3/1 3/9 3/13 3 3_//7 3/6 No.^b Pos.^c NO.

G/A A/G' A/A A/G A/G A A_//A A/A 1 3915 1

G/G G/G G/G G/G G/G G G_//G G/G 2 4437 1

C/C C/C C/C C/C C/C C C_//C C/C 3 5008 1

105 G G//GG GG//GG G/G G/A G/G G G_//G G/G 4 5037 1

C/C C/C C/C C/C C/C C C_//G C/C 5 5102 1

C/C C/C C/C C/C C/C C C_//G C/C 6 5290 1

A/C C/A C/C C/A C/A C C_//C C/C 7 5363 1

A/G G/A G/G G/A G/A G G_//G G/G 8 5628 1

110 G G//AA GG//GG G/A G/A G/A G G_//G G/G 9 7721 1

C/C C/C C/C C/C C/T C C_//C C/C 10 10537 1

T/T T/T T/T T/T T/T T T_//T T/T 11 168 2

A/A A/A A/A A/A A/A A A_//A A/A 12 257 2

C/C C/C C/C C/C C/C C C_//T C/C 13 2564 2

115 G G//GG ^'' GG//GG G/G G/G G/G G G_//C G/G 14 2577 2

T/T T/T T/T T/C T/T T T_//T T/T 15 3436 3

A/A A/A A/A A/A A/A A A_//A A/G 16 3548 3

A/A A/A A/A A/A A/A A A_//A A/A 17 3750 3

G/A A/G A/A A/A A/A A A_//A A/A 18 4007 3

120

Haplotype pairs are represented as l 1Stτ Hτa„ „pTlo+t,y,p„_αe//2onⁿdl r Hraplotype; with alleles of each haplotype shown 5' to 3' as 1^st polymoφhism/2^nd polymoφhism in each column; ^bPS = polymoφhic site;

"Position of PS within the indicated SEQ ID NO.

32. A genome anthology for the matrix metalloproteinase 13 (collagenase 3) (MMP13) gene which comprises MMP 13 isogenes defined by any one of haplotypes 1-15 set forth in the table shown below:

Ilap] .otj e Nunibei -a PS PS SEQ ID

1 2 3 4 5 6 7 8 9 10 No.^b Pos.^c NO.

A A A A A A A A G G 1 3915 1

G G G G G G G T G G 2 4437 1

C C C C C C C C C C 3 5008 1

10 G G G G G G G G A G 4 503 1

C C C C C C G C C C 5 5102 1

C C C C C C G C C C 6 5290 1

C C C C C C C C A A 7 5363 1

G G G G G G G G A A 8 5628 1

15 A A G G G G G G A A 9 7721 1

C C C C C C C C C C 10 10537 1

T T T T T T T T T C 11 168 2

A A A A A A A A A A 12 257 2

C C C C C C T C C C 13 2564 2

20 G G G G G G C G G G 14 2577 2

T T T T T T T T C T 15 3436 3

A A A A A G A A A A 16 3548 3

A A A A G A A A A A 17 3750 3

A G A G A A A A A A 18 4007 3

25 Haplotype Number* PS PS SEQ ID

11 12 13 14 15 No.^b Pos.° NO.

G G G G G 1 3915 1

G G G G G 2 4437 1

C C C C T 3 5008 1

G G G G G 4 5037 1

C C C C C 5 5102 1

C C C C C 6 5290 1

A A A A A 7 5363 1

A A A A A 8 5628 1

A A A G G 9 7721 1

C C T C C 10 10537 1

T T T T T 11 168 2

A G A A A 12 257 - 2

C C C C C 13 2564 2

G G G G G 14 2577 2

T T T T T 15 3436 3

A A A A A 16 3548 3

A A A A A 17 3750 3

A A A G G 18 4007 3

"Alleles for haplotypes are presented 5' to 3' in each column ^bPS = polymoφhic site; "Position of PS within the indicated SEQ ID NO.