Abstract

Isolated GSTO2 nucleic acid molecules that include a nucleotide sequence variant and nucleotides flanking the sequence variant are described, as well as GSTO2 allozymes. Methods for determining if a subject contains a GSTO2 sequence variant also are described.

Description

CROSS-REFERENCE TO RELATED APPLICATIONS

This application is a divisional of U.S. patent application Ser. No. 11/482,525, filed on Jul. 7, 2006, now U.S. Pat. No. 7,390,894, which claims benefit of priority from U.S. Provisional Patent Application Ser. No. 60/697,128, filed on Jul. 7, 2005, the entire contents of which are hereby incorporated by reference.

STATEMENT AS TO FEDERALLY SPONSORED RESEARCH

The invention was made with government support under grant nos. GM28157, GM35720, and GM61388, awarded by the National Institutes of Health. The government has certain rights in the invention.

The Omega class of glutathione transferases (GSTs) have been identified in organisms including human, mouse, rat, pig, Caenorhabditis elegans, and Drosophila melanogaster. These GSTs, which include glutathione S-transferase omega 2 (GSTO2) have poor activity with common GST substrates, but exhibit novel glutathione-dependent thioltransferase, dehydroascorbate reductase, and monomethylarsonate reductase activities, and also modulate calcium release by ryanodine receptors. The GSTO2 gene contains seven exons, the first of which does not contain GSTO2 coding sequence. The GSTO2 gene is separated from the GSTO1 gene by 7.5 kb on chromosome 10q24.3. GSTO2 is expressed at high levels in the testes, and also is expressed in liver, kidney, skeletal muscle, heart, and lung, with low level expression in the cervix, ovary, and prostate. Substrates for GSTO2 include monomethylarsonic acid (MMV), which is reduced from the +V form to the +III form, as well as dehydroascorbate (DHA) and glutathione (GSH), which can be reacted to give ascorbate and glutathione disulfate.

SUMMARY

This document is based on the discovery of sequence variants that occur in both coding and non-coding regions of GSTO2 nucleic acids. Certain GSTO2 nucleotide sequence variants encode GSTO2 enzymes that are associated with individual differences in enzymatic activity. Other GSTO2 sequence variants in non-coding regions of the GSTO2 nucleic acid may alter regulation of transcription and/or splicing of the GSTO2 nucleic acid. Discovery of these sequence variants allows individual differences in the metabolism of drugs and other xenobiotics in humans to be assessed such that particular treatment regimens can be tailored to an individual based on the presence or absence of one or more sequence variants. For example, the presence of absence of a particular variant may indicate an individual's ability to metabolize arsenic, which is a contaminant in drinking water, and which can cause skin, bladder, kidney, liver, and lung cancer in individuals subjected to chronic arsenic exposure. The presence of absence of a particular variant also may indicate an individual's ability to metabolize arsenic trioxide, an agent that can be used to treat acute promyelocytic leukemia. Identification of GSTO2 sequence variants also may allow predisposition to cancer or conditions such as Alzheimer's or Parkinson's disease to be assessed in individuals.

In one aspect, this document features an isolated nucleic acid molecule consisting essentially of a variant GSTO2 nucleic acid sequence, wherein the variant GSTO2 nucleic acid sequence is selected from the group consisting of: a) at least ten contiguous nucleotides of SEQ ID NO:8, wherein said sequence includes nucleotide 121, 389, or 472 of SEQ ID NO:5, with the proviso that the nucleotide at position 121 is adenine, the nucleotide at position 389 is adenine, or the nucleotide at position 472 is adenine; b) at least ten contiguous nucleotides of SEQ ID NO:1, wherein said sequence includes nucleotide 15569, 15811, 15918, 15984, 16006, 16016, 16041, 16403, 16477, 16801, 16851, 16945, or 16948 of SEQ ID NO:1, with the proviso that a guanine is inserted after the nucleotide at position 15569, the nucleotide at position 15811 is cytosine, the nucleotide at position 15918 is thymine, the nucleotide at position 15984 is guanine, the nucleotide at position 16006 is cytosine, the nucleotide at position 16016 is adenine, the nucleotide at position 16041 is thymine, the nucleotide at position 16403 is thymine, the nucleotide at position 16477 is adenine, the nucleotide at position 16801 is cytosine, a thymine is inserted after the nucleotide at position 16851, the nucleotide at position 16945 is thymine, or the nucleotide at position 16948 is guanine; c) at least ten contiguous nucleotides of SEQ ID NO:1, wherein said sequence includes nucleotide 16957, 17056, 17119, 17132, 17158, 17172, or 17348 of SEQ ID NO:1, with the proviso that the nucleotide at position 16957 is guanine, the nucleotide at position 17056 is adenine, the nucleotide at position 17119 is adenine, the nucleotide at position 17132 is cytosine, the nucleotide at position 17158 is adenine, the nucleotide at position 17172 is guanine, or the nucleotide at position 17348 is thymine; d) at least ten contiguous nucleotides of SEQ ID NO:1, wherein said sequence includes nucleotide 17402, 17502, 17659 and 17662, 17683, 17752, 17785, 17827, 19044, 19496, 21635, 22747, or 22754 of SEQ ID NO:1, with the proviso that the nucleotide at position 17402 is guanine, the nucleotide at position 17502 is thymine, the nucleotides at positions 17660 and 17661 are deleted, the nucleotide at position 17683 is cytosine, the nucleotide at position 17752 is thymine, the nucleotide at position 17785 is adenine, a thymine is inserted after the nucleotide at position 17827, the nucleotide at position 19044 is guanine, the nucleotide at position 19496 is guanine, the nucleotide at position 21635 is thymine, the nucleotide at position 22747 is adenine, or the nucleotide at position 22754 is adenine; e) at least ten contiguous nucleotides of SEQ ID NO:1, wherein said sequence includes nucleotide 23291, 26236, 27722, 47154, or 47188 of SEQ ID NO:1, with the proviso that the nucleotide at position 23291 is thymine, the nucleotide at position 26236 is thymine, the nucleotide at position 27722 is cytosine, the nucleotide at position 47154 is guanine, or the nucleotide at position 47188 is thymine; f) at least ten contiguous nucleotides of SEQ ID NO:1, wherein said sequence includes nucleotide 47713, 47876, 48248, 48282, or 48357 of SEQ ID NO:1, with the proviso that the nucleotide at position 47713 is thymine, the nucleotide at position 47876 is guanine, the nucleotide at position 48248 is thymine, the nucleotide at position 48282 is cytosine, or the nucleotide at position 48357 is thymine; and g) the complement of a), b), c), d), e), or f). The isolated nucleic acid can be from 10 to 100 nucleotides in length or from 20 to 50 nucleotides in length.

This document also features a vector comprising the nucleic acid molecule described herein. The nucleic acid molecule can be from 20 to 50 nucleotides in length.

In another aspect, this document features an isolated nucleic acid encoding a GSTO2 polypeptide, wherein the polypeptide contains a GSTO2 amino acid sequence variant relative to the amino acid sequence of SEQ ID NO:9, and wherein the amino acid sequence variant is at a residue selected from the group consisting of 41, 130, and 158. The amino acid sequence variant can be an isoleucine at residue 41, a tyrosine at residue 130, or an isoleucine at residue 158.

In another aspect, this document features an isolated GSTO2 polypeptide, wherein the polypeptide contains a GSTO2 amino acid sequence variant relative to the amino acid sequence of SEQ ID NO:9, and wherein the amino acid sequence variant is at a residue selected from the group consisting of 41, 130, and 158. The amino acid sequence variant can be an isoleucine at residue 41, a tyrosine at residue 130, or an isoleucine at residue 158.

In still another aspect, this document features an article of manufacture comprising a substrate, wherein the substrate comprises a population of isolated variant GSTO2 nucleic acid molecules as disclosed herein. The article of manufacture can comprises a plurality of discrete regions, wherein each region comprises a different population of isolated variant GSTO2 nucleic acid molecules, and wherein each population of molecules comprises a different GSTO2 nucleotide sequence variant.

This document also features a method for determining a GSTO2 genotype. The method can comprise: a) providing a nucleic acid sample from a subject, and b) screening the nucleic acid sample for one or more genetic markers in linkage disequilibrium with a GSTO2 allele. The GSTO2 allele can be selected from the group consisting of a G121A allele, a G389A allele, an A424G allele, or a C472A allele.

Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention pertains. Although methods and materials similar or equivalent to those described herein can be used to practice the invention, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.

Other features and advantages of the invention will be apparent from the following detailed description and drawings, and from the claims.

FIG. 2 shows the open reading frame (SEQ ID NO:8) of the reference GSTO2, and also shows the reference amino acid sequence (SEQ ID NO:9) of the encoded GSTO2. Positions of cSNPs are circled, as are the positions of amino acid changes that result from the cSNPs. Start and stop codons are double-underlined.

This document features GSTO2 nucleotide and amino acid sequence variants. GSTO2 is thought to be a cytoplasmic enzyme that has glutathione transferase and glutathione dehydrogenase activity. For example, GSTO2 can catalyze the transfer of glutathione to an aliphatic, aromatic, or heterocyclic molecule, wherein the linkage is via a thiol group. GSTO2 also can catalyze the dehydrogenation of DHA, which can be reacted with GSH to give ascorbate and glutathione disulfate. Genetically-based variations in GSTO2 activity that lead to altered levels of GSTO2 or altered GSTO2 activity may be important in particular diseases or clinical conditions.

Nucleic Acid Molecules

Provided herein are isolated nucleic acids that include a GSTO2 nucleic acid sequence. The GSTO2 nucleic acid sequence includes a nucleotide sequence variant and nucleotides flanking the sequence variant. As used herein, “isolated nucleic acid” refers to a nucleic acid that is separated from other nucleic acid molecules that are present in a mammalian genome, including nucleic acids that normally flank one or both sides of the nucleic acid in a mammalian genome (e.g., nucleic acids that encode non-GSTO2 proteins). The term “isolated” as used herein with respect to nucleic acids also includes any non-naturally-occurring nucleic acid sequence since such non-naturally-occurring sequences are not found in nature and do not have immediately contiguous sequences in a naturally-occurring genome.

An isolated nucleic acid can be, for example, a DNA molecule, provided one of the nucleic acid sequences normally found immediately flanking that DNA molecule in a naturally-occurring genome is removed or absent. Thus, an isolated nucleic acid includes, without limitation, a DNA molecule that exists as a separate molecule (e.g., a chemically synthesized nucleic acid, or a cDNA or genomic DNA fragment produced by PCR or restriction endonuclease treatment) independent of other sequences as well as DNA that is incorporated into a vector, an autonomously replicating plasmid, a virus (e.g., a retrovirus, lentivirus, adenovirus, or herpes virus), or into the genomic DNA of a prokaryote or eukaryote. In addition, an isolated nucleic acid can include an engineered nucleic acid such as a recombinant DNA molecule that is part of a hybrid or fusion nucleic acid. A nucleic acid existing among hundreds to millions of other nucleic acids within, for example, cDNA libraries or genomic libraries, or gel slices containing a genomic DNA restriction digest, is not to be considered an isolated nucleic acid.

Nucleic acids provided herein are at least about 8 nucleotides in length. For example, the nucleic acid can be about 8, 9, 10-20 (e.g., 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20 nucleotides in length), 20-50, 50-100 or greater than 100 nucleotides in length (e.g., greater than 150, 200, 250, 300, 350, 400, 450, 500, 750, or 1000 nucleotides in length). Nucleic acids can be in a sense or antisense orientation, can be complementary to the GSTO2 reference sequence, and can be DNA, RNA, or nucleic acid analogs. Nucleic acid analogs can be modified at the base moiety, sugar moiety, or phosphate backbone to improve, for example, stability, hybridization, or solubility of the nucleic acid. Modifications at the base moiety include deoxyuridine for deoxythymidine, and 5-methyl-2′-deoxycytidine or 5-bromo-2′-doxycytidine for deoxycytidine. Modifications of the sugar moiety include modification of the 2′ hydroxyl of the ribose sugar to form 2′-O-methyl or 2′-O-allyl sugars. The deoxyribose phosphate backbone can be modified to produce morpholino nucleic acids, in which each base moiety is linked to a six membered, morpholino ring, or peptide nucleic acids, in which the deoxyphosphate backbone is replaced by a pseudopeptide backbone and the four bases are retained. See, Summerton and Weller (1997) Antisense Nucleic Acid Drug Dev. 7:187-195; and Hyrup et al. (1996) Bioorgan. Med. Chem. 4:5-23. In addition, the deoxyphosphate backbone can be replaced with, for example, a phosphorothioate or phosphorodithioate backbone, a phosphoroamidite, or an alkyl phosphotriester backbone.

As used herein, “nucleotide sequence variant” refers to any alteration in a GSTO2 reference sequence, and includes variations that occur in coding and non-coding regions, including exons, introns, and untranslated sequences. Nucleotides are referred to herein by the standard one-letter designation (A, C, G, or T). Variations include single nucleotide substitutions, deletions of one or more nucleotides, and insertions of one or more nucleotides. The reference GSTO2 nucleic acid sequence is provided in FIG. 1 (SEQ ID NO:1) and in GenBank (Accession No. NT—030059). The reference GSTO2 ORF is provided in FIG. 2 (SEQ ID NO:8), as is the corresponding reference GSTO2 amino acid sequence (SEQ ID NO:9). The mRNA and amino acid reference sequences also are found in GenBank (Accession No. NM—183239). The nucleic acid and amino acid reference sequences also are referred to herein as “wild type.”

As used herein, “untranslated sequence” includes 5′ and 3′ flanking regions that are outside of the messenger RNA (mRNA) as well as 5′ and 3′ untranslated regions (5′-UTR or 3′-UTR) that are part of the mRNA, but are not translated. Positions of nucleotide sequence variants in 5′ untranslated sequences can be designated as “−X” relative to the “A” in the translation initiation codon; positions of nucleotide sequence variants in the coding sequence and 3′ untranslated sequence can be designated as “+X” or “X” relative to the “A” in the translation initiation codon. Nucleotide sequence variants that occur in introns can be designated as “+X” or “X” relative to the “G” in the splice donor site (GT) or as “−X” relative to the “G” in the splice acceptor site (AG).

In some embodiments, a GSTO2 nucleotide sequence variant encodes a GSTO2 polypeptide having an altered amino acid sequence. The term “polypeptide” refers to a chain of at least four amino acid residues (e.g., 4-8, 9-12, 13-15, 16-18, 19-21, 22-100, 100-150, 150-200, 200-225 residues, or a full-length GSTO2 polypeptide). GSTO2 polypeptides may or may not have GSTO2 catalytic activity, or may have altered activity relative to the reference GSTO2 polypeptide. Polypeptides that do not have activity or have altered activity can be useful for diagnostic purposes (e.g., for producing antibodies having specific binding affinity for variant GSTO2 polypeptides).

Corresponding GSTO2 polypeptides, irrespective of length, that differ in amino acid sequence are herein referred to as allozymes. For example, a GSTO2 nucleic acid sequence that includes an adenine at position 121 relative to the adenine in the translation initiation codon (i.e., nucleotide 23393 of SEQ ID NO:1 or nucleotide 121 of SEQ ID NO:8) encodes a GSTO2 polypeptide having an isoleucine at amino acid residue 41. This polypeptide (Val41Ile) would be considered an allozyme with respect to the reference GSTO2 polypeptide that contains a valine at amino acid residue 41. Additional non-limiting examples of GSTO2 sequence variants that alter amino acid sequence include variants at nucleotides 389, 424, and 472 relative to the adenine in the translation initiation codon (i.e., positions 27473, 27508, and 45662, respectively, of SEQ ID NO:1, or positions 389, 424, and 472, respectively, of SEQ ID NO:8). For example, a GSTO2 nucleic acid molecule can include an adenine at nucleotide 389 and encode a GSTO2 polypeptide having a tyrosine at amino acid residue 130 in place of a cysteine residue (Cys130Tyr); a guanine at nucleotide 424 and encode a GSTO2 polypeptide having an aspartic acid at amino acid residue 142 in place of an asparagine residue (Asn142Asp), or an adenine at nucleotide 472 and encode a GSTO2 polypeptide having an isoleucine at amino acid 158 in place of a leucine residue (Leu158Ile).

GSTO2 allozymes as described above are encoded by a series of GSTO2 alleles. These alleles represent nucleic acid sequences containing sequence variants, typically multiple sequence variants, within coding and non-coding sequences. Representative examples of single nucleotide variants are described herein. Table 2 sets out a series of GSTO2 alleles that encode GSTO2. Some alleles are commonly observed, i.e., have allele frequencies >1%, such as the allele having a guanine at nucleotide 424 place of an adenine. The relatively large number of alleles and allozymes for GSTO2 indicates the potential complexity of GSTO2 pharmacogenetics. Such complexity emphasizes the need for determining single nucleotide variants, (i.e., single nucleotide polymorphisms, SNPs) as well as complete GSTO2 haplotypes (i.e., the set of alleles on one chromosome or a part of a chromosome) of patients.

GSTO2 nucleotide sequence variants that do not change the amino acid sequence also can be within an exon or in 5′ or 3′ untranslated sequences. For example, a GSTO2 nucleotide sequence variant within exon 2 can be a guanine substitution at nucleotide −183 relative to the adenine in the translation initiation codon. Exon 7 sequence variants can include, for example, a thymine substitution at nucleotide 591 (i.e., nucleotide 47224 of SEQ ID NO:1 or nucleotide 591 of SEQ ID NO:8), or a thymine substitution at nucleotide 630 (i.e., nucleotide 47263 of SEQ ID NO:1 or nucleotide 630 of SEQ ID NO:8).

In some embodiments, a GSTO2 nucleic acid molecule can consist essentially of at least ten (e.g., at least 12, at least 15, at least 18, at least 20, or at least 25) contiguous nucleotides of an GSTO2 reference sequence (e.g., SEQ ID NO:1 or SEQ ID NO:8). Such nucleic acids can contain one or more variant positions, with the proviso that the nucleotides at those positions are variant nucleotides as disclosed herein in Table 2, for example. A GSTO2 nucleic acid “consisting essentially of” a particular sequence has the basic and novel characteristic that it can be used to distinguish, based upon hybridization, a nucleic acid having a sequence that contains a variant from a corresponding nucleic acid having a sequence that does not contain the variant (e.g., a wild type sequence). Such nucleic acid molecules can include additional sequences or labels (e.g., a tag or a fluorescent label as disclosed herein), provided that such additions do not affect the basic and novel characteristic of the nucleic acid molecules.

In some embodiments, a GSTO2 nucleic acid sequence can have at least 90% (e.g., 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 97.5%, 98%, 98.5%, 99.0%, 99.5%, 99.6%, 99.7%, 99.8%, 99.9%, or 100%) sequence identity with a region of SEQ ID NO:1 or SEQ ID NO:8. The region of SEQ ID NO:1 or SEQ ID NO:8 is at least ten nucleotides in length (e.g., 10, 15, 20, 50, 60, 70, 75, 100, 150 or more nucleotides in length). For example, a nucleic acid can include a GSTO2 nucleic acid sequence with at least 90% identity to nucleotides 75 to 175, 350 to 450, 375 to 475, 425 to 525, 550 to 650, or 600 to 700 of SEQ ID NO:8. A variant GSTO2 nucleotide sequence can have an adenine at nucleotide 121, an adenine at nucleotide 389, an adenine at nucleotide 424, an adenine at nucleotide 472, a thymine at nucleotide 591, or a thymine at nucleotide 630, and combinations thereof, wherein the numbering of the positions is with respect to the numbering of SEQ ID NO:8.

Percent sequence identity is calculated by determining the number of matched positions in aligned nucleic acid sequences, dividing the number of matched positions by the total number of aligned nucleotides, and multiplying by 100. A matched position refers to a position in which identical nucleotides occur at the same position in aligned nucleic acid sequences. Percent sequence identity also can be determined for any amino acid sequence. To determine percent sequence identity, a target nucleic acid or amino acid sequence is compared to the identified nucleic acid or amino acid sequence using the BLAST 2 Sequences (B12seq) program from the stand-alone version of BLASTZ containing BLASTN version 2.0.14 and BLASTP version 2.0.14. This stand-alone version of BLASTZ can be obtained on the World Wide Web from Fish & Richardson's web site (fr.com/blast) or the U.S. government's National Center for Biotechnology Information web site (ncbi.nlm.nih.gov). Instructions explaining how to use the B12seq program can be found in the readme file accompanying BLASTZ.

B12seq performs a comparison between two sequences using either the BLASTN or BLASTP algorithm. BLASTN is used to compare nucleic acid sequences, while BLASTP is used to compare amino acid sequences. To compare two nucleic acid sequences, the options are set as follows: −i is set to a file containing the first nucleic acid sequence to be compared (e.g., C:\seq1.txt); −j is set to a file containing the second nucleic acid sequence to be compared (e.g., C:\seq2.txt); −p is set to blastn; −o is set to any desired file name (e.g., C:\output.txt); −q is set to −1; −r is set to 2; and all other options are left at their default setting. The following command will generate an output file containing a comparison between two sequences: C:\B12seq −i c:\seq1.txt −j c:\seq2.txt −p blastn −o c:\output.txt −q −1−r 2. If the target sequence shares homology with any portion of the identified sequence, then the designated output file will present those regions of homology as aligned sequences. If the target sequence does not share homology with any portion of the identified sequence, then the designated output file will not present aligned sequences.

Once aligned, a length is determined by counting the number of consecutive nucleotides from the target sequence presented in alignment with sequence from the identified sequence starting with any matched position and ending with any other matched position. A matched position is any position where an identical nucleotide is presented in both the target and identified sequence. Gaps presented in the target sequence are not counted since gaps are not nucleotides. Likewise, gaps presented in the identified sequence are not counted since target sequence nucleotides are counted, not nucleotides from the identified sequence.

The percent identity over a particular length is determined by counting the number of matched positions over that length and dividing that number by the length followed by multiplying the resulting value by 100. For example, if (1) a 1000 nucleotide target sequence is compared to the sequence set forth in SEQ ID NO:1, (2) the B12seq program presents 969 nucleotides from the target sequence aligned with a region of the sequence set forth in SEQ ID NO:1 where the first and last nucleotides of that 969 nucleotide region are matches, and (3) the number of matches over those 969 aligned nucleotides is 900, then the 1000 nucleotide target sequence contains a length of 969 and a percent identity over that length of 93 (i.e., 900÷969×100=93).

It will be appreciated that different regions within a single nucleic acid target sequence that aligns with an identified sequence can each have their own percent identity. It is noted that the percent identity value is rounded to the nearest tenth. For example, 78.11, 78.12, 78.13, and 78.14 are rounded down to 78.1, while 78.15, 78.16, 78.17, 78.18, and 78.19 are rounded up to 78.2. It also is noted that the length value will always be an integer.

Isolated nucleic acid molecules can be produced using standard techniques, including, without limitation, common molecular cloning and chemical nucleic acid synthesis techniques. For example, polymerase chain reaction (PCR) techniques can be used to obtain an isolated nucleic acid containing a GSTO2 nucleotide sequence variant. PCR refers to a procedure or technique in which target nucleic acids are enzymatically amplified. Sequence information from the ends of the region of interest or beyond typically is employed to design oligonucleotide primers that are identical in sequence to opposite strands of the template to be amplified. PCR can be used to amplify specific sequences from DNA as well as RNA, including sequences from total genomic DNA or total cellular RNA. Primers are typically 14 to 40 nucleotides in length, but can range from 10 nucleotides to hundreds of nucleotides in length. General PCR techniques are described, for example in PCR Primer: A Laboratory Manual, ed. by Dieffenbach and Dveksler, Cold Spring Harbor Laboratory Press, 1995. When using RNA as a source of template, reverse transcriptase can be used to synthesize complementary DNA (cDNA) strands. Ligase chain reaction, strand displacement amplification, self-sustained sequence replication, or nucleic acid sequence-based amplification also can be used to obtain isolated nucleic acids. See, for example, Lewis (1992) Genetic Engineering News 12:1; Guatelli et al. (1990) Proc. Natl. Acad. Sci. USA 87:1874-1878; and Weiss (1991) Science 254:1292.

Isolated nucleic acids also can be chemically synthesized, either as a single nucleic acid molecule (e.g., using automated DNA synthesis in the 3′ to 5′ direction using phosphoramidite technology) or as a series of oligonucleotides. For example, one or more pairs of long oligonucleotides (e.g., >100 nucleotides) can be synthesized that contain the desired sequence, with each pair containing a short segment of complementarity (e.g., about 15 nucleotides) such that a duplex is formed when the oligonucleotide pair is annealed. DNA polymerase is used to extend the oligonucleotides, resulting in a single, double-stranded nucleic acid molecule per oligonucleotide pair, which then can be ligated into a vector.

Isolated nucleic acids also can be obtained by mutagenesis. For example, the reference sequences depicted in FIG. 1 or 2 can be mutated using standard techniques including oligonucleotide-directed mutagenesis and site-directed mutagenesis through PCR. See, Short Protocols in Molecular Biology, Chapter 8, Green Publishing Associates and John Wiley & Sons, edited by Ausubel et al., 1992. Examples of positions that can be modified include those described herein.

GSTO2 Polypeptides

Isolated GSTO2 polypeptides provided herein include an amino acid sequence variant relative to the reference GSTO2 (SEQ ID NO:9, FIG. 2; GenBank Accession No. NM—183239). The term “isolated” with respect to a GSTO2 polypeptide refers to a polypeptide that has been separated from cellular components by which it is naturally accompanied. Typically, the polypeptide is isolated when it is at least 60% (e.g., 70%, 80%, 90%, 95%, or 99%), by weight, free from proteins and naturally-occurring organic molecules with which it is naturally associated. In general, an isolated polypeptide will yield a single major band on a non-reducing polyacrylamide gel.

GSTO2 polypeptides can include variants at one or more of amino acid residues 41, 130, 142, and 158. In particular, an isoleucine residue can be substituted at position 41, a tyrosine residue at position 130, an aspartic acid residue at position 142, or an isoleucine residue at position 158. In some embodiments, activity of GSTO2 polypeptides is altered relative to the reference GSTO2. Certain GSTO2 allozymes can have reduced activity, while other allozymes can have activity that is comparable to the reference GSTO2. Other allozymes can have increased activity relative to the reference GSTO2. Activity of GSTO2 polypeptides can be assessed in vitro. For example, reduction of DHA catalyzed by bacterially expressed GSTO2 can be detected by monitoring an increase in absorbance at 265 nm. Alternatively, activity of GSTO2 polypeptides can be assessed in vivo (e.g., in COS cells).

Other biochemical properties of allozymes, such as apparent Km values, also can be altered relative to the reference GSTO2. Apparent Km values can be calculated using, for example, the method of Wilkinson with a computer program written by Cleland. Wilkinson (1961) Biochem. J. 80:324-332; and Cleland (1963) Nature 198:463-365.

Isolated polypeptides can be obtained, for example, by extraction from a natural source (e.g., testes, liver, kidney, or skeletal muscle), chemical synthesis, or by recombinant production in a host cell. To recombinantly produce GSTO2 polypeptides, a nucleic acid encoding a GSTO2 nucleotide sequence variant can be ligated into an expression vector and used to transform a prokaryotic (e.g., bacteria) or eukaryotic (e.g., insect, yeast, or mammal) host cell. In general, nucleic acid constructs include a regulatory sequence operably linked to a GSTO2 nucleic acid sequence. Regulatory sequences (e.g., promoters, enhancers, polyadenylation signals, or terminators) do not typically encode a gene product, but instead affect the expression of the nucleic acid sequence. In addition, a construct can include a tag sequence designed to facilitate subsequent manipulations of the expressed nucleic acid sequence (e.g., purification, localization). Tag sequences, such as green fluorescent protein (GFP), glutathione S-transferase (GST), six histidine (His6), c-myc, hemagglutinin, or Flag™ tag (Kodak) sequences are typically expressed as a fusion with the expressed nucleic acid sequence. Such tags can be inserted anywhere within the polypeptide including at either the carboxyl or amino termini. The type and combination of regulatory and tag sequences can vary with each particular host, cloning or expression system, and desired outcome. A variety of cloning and expression vectors containing combinations of regulatory and tag sequences are commercially available. Suitable cloning vectors include, without limitation, pUC18, pUC19, and pBR322 and derivatives thereof (New England Biolabs, Beverly, Mass.), and pGEN (Promega, Madison, Wis.). Additionally, representative prokaryotic expression vectors include pBAD (Invitrogen, Carlsbad, Calif.), the pTYB family of vectors (New England Biolabs), and pGEMEX vectors (Promega); representative mammalian expression vectors include pTet-On/pTet-Off (Clontech, Palo Alto, Calif.), pIND, pVAX1, pCR3.1, pcDNA3.1, pcDNA4, or pUni (Invitrogen), and pCI or pSI (Promega); representative insect expression vectors include pBacPAK8 or pBacPAK9 (Clontech), and p2Bac (Invitrogen); and representative yeast expression vectors include MATCHMAKER (Clontech) and pPICZ A, B, and C (Invitrogen).

In bacterial systems, a strain of Escherichia coli can be used to express GSTO2 variant polypeptides. For example, BL-21 cells can be transformed with a pGEX vector containing a GSTO2 nucleic acid sequence. The transformed bacteria can be grown exponentially and then stimulated with isopropylthiogalactopyranoside (IPTG) prior to harvesting. In general, the GSTO2-GST fusion proteins produced from the pGEX expression vector are soluble and can be purified easily from lysed cells by adsorption to glutathione-agarose beads followed by elution in the presence of free glutathione. The pGEX vectors are designed to include thrombin or factor Xa protease cleavage sites so that the expressed GSTO2 polypeptide can be released from the GST moiety.

In eukaryotic host cells, a number of viral-based expression systems can be utilized to express GSTO2 variant polypeptides. A nucleic acid encoding a polypeptide can be cloned into, for example, a baculoviral vector such as pBlueBac (Invitrogen) and then used to co-transfect insect cells such as Spodoptera frugiperda (Sf9) cells with wild type DNA from Autographa californica multinuclear polyhedrosis virus (AcMNPV). Recombinant viruses producing polypeptides can be identified by standard methodology. Alternatively, a nucleic acid encoding a polypeptide can be introduced into a SV40, retroviral, or vaccinia based viral vector and used to infect suitable host cells.

GSTO2 variant polypeptides can be purified by known chromatographic methods including ion exchange and gel filtration chromatography. See, for example, Caine et al. (1996) Protein Expr. Purif 8:159-166. GSTO2 polypeptides can be “engineered” to contain a tag sequence describe herein that allows the polypeptide to be purified (e.g., captured onto an affinity matrix). Immunoaffinity chromatography also can be used to purify GSTO2 polypeptides.

Non-Human Mammals

Also featured herein are non-human mammals that include GSTO2 nucleic acids as described herein, as well as progeny and cells of such non-human mammals. Non-human mammals include, for example, rodents such as rats, guinea pigs, and mice, and farm animals such as pigs, sheep, goats, horses, and cattle. Non-human mammals can express a GSTO2 variant nucleic acid in addition to an endogenous GSTO2 (e.g., a transgenic non-human that includes a GSTO2 nucleic acid randomly integrated into the genome of the non-human mammal). Alternatively, an endogenous GSTO2 nucleic acid can be replaced with a GSTO2 variant nucleic acid by homologous recombination. See, Shastry (1998) Mol. Cell. Biochem. 181:163-179, for a review of gene targeting technology.

In one embodiment, non-human mammals are produced that lack an endogenous GSTO2 nucleic acid (i.e., a knockout), and then a GSTO2 variant nucleic acid as provided herein is introduced into the knockout non-human mammal. Nucleic acid constructs used for producing knockout non-human mammals can include a nucleic acid sequence encoding a selectable marker, which is generally used to interrupt the targeted exon site by homologous recombination. Typically, the selectable marker is flanked by sequences homologous to the sequences flanking the desired insertion site. It is not necessary for the flanking sequences to be immediately adjacent to the desired insertion site. Suitable markers for positive drug selection include, for example, the aminoglycoside 3N phosphotransferase gene that imparts resistance to geneticin (G418, an aminoglycoside antibiotic), and other antibiotic resistance markers, such as the hygromycin-B-phosphotransferase gene that imparts hygromycin resistance. Other selection systems include negative-selection markers such as the thymidine kinase (TK) gene from herpes simplex virus. Constructs utilizing both positive and negative drug selection also can be used. For example, a construct can contain the aminoglycoside phosphotransferase gene and the TK gene. In this system, cells are selected that are resistant to G418 and sensitive to gancyclovir.

To create non-human mammals having a particular gene inactivated in all cells, it is necessary to introduce a knockout construct into the germ cells (sperm or eggs, i.e., the “germ line”) of the desired species. Genes or other DNA sequences can be introduced into the pronuclei of fertilized eggs by microinjection. Following pronuclear fusion, the developing embryo may carry the introduced gene in all its somatic and germ cells because the zygote is the mitotic progenitor of all cells in the embryo. Since targeted insertion of a knockout construct is a relatively rare event, it is desirable to generate and screen a large number of animals when employing such an approach. Because of this, it can be advantageous to work with the large cell populations and selection criteria that are characteristic of cultured cell systems. However, for production of knockout animals from an initial population of cultured cells, it is necessary that a cultured cell containing the desired knockout construct be capable of generating a whole animal. This is generally accomplished by placing the cell into a developing embryo environment of some sort.

Cells capable of giving rise to at least several differentiated cell types are “pluripotent.” Pluripotent cells capable of giving rise to all cell types of an embryo, including germ cells, are hereinafter termed “totipotent” cells. Totipotent murine cell lines (embryonic stem, or “ES” cells) have been isolated by culture of cells derived from very young embryos (blastocysts). Such cells are capable, upon incorporation into an embryo, of differentiating into all cell types, including germ cells, and can be employed to generate animals lacking an endogenous GSTO2 nucleic acid. That is, cultured ES cells can be transformed with a knockout construct and cells selected in which the GSTO2 gene is inactivated.

Nucleic acid constructs can be introduced into ES cells, for example, by electroporation or other standard technique. Selected cells can be screened for gene targeting events. For example, the polymerase chain reaction (PCR) can be used to confirm the presence of the transgene.

The ES cells further can be characterized to determine the number of targeting events. For example, genomic DNA can be harvested from ES cells and used for Southern analysis. See, for example, Section 9.37-9.52 of Sambrook et al., Molecular Cloning A Laboratory Manual, second edition, Cold Spring Harbor Press, Plainview; NY, 1989.

To generate a knockout animal, ES cells having at least one inactivated GSTO2 allele are incorporated into a developing embryo. This can be accomplished through injection into the blastocyst cavity of a murine blastocyst-stage embryo, by injection into a morula-stage embryo, by co-culture of ES cells with a morula-stage embryo, or through fusion of the ES cell with an enucleated zygote. The resulting embryo is raised to sexual maturity and bred in order to obtain animals, whose cells (including germ cells) carry the inactivated GSTO2 allele. If the original ES cell was heterozygous for the inactivated GSTO2 allele, several of these animals can be bred with each other in order to generate animals homozygous for the inactivated allele.

Alternatively, direct microinjection of DNA into eggs can be used to avoid the manipulations required to turn a cultured cell into an animal. Fertilized eggs are totipotent, i.e., capable of developing into an adult without further substantive manipulation other than implantation into a surrogate mother. To enhance the probability of homologous recombination when eggs are directly injected with knockout constructs, it is useful to incorporate at least about 8 kb of homologous DNA into the targeting construct. In addition, it is also useful to prepare the knockout constructs from isogenic DNA.

Embryos derived from microinjected eggs can be screened for homologous recombination events in several ways. For example, if the GSTO2 gene is interrupted by a coding region that produces a detectable (e.g., fluorescent) gene product, then the injected eggs are cultured to the blastocyst stage and analyzed for presence of the indicator polypeptide. Embryos with fluorescing cells, for example, are then implanted into a surrogate mother and allowed to develop to term. Alternatively, injected eggs are allowed to develop and DNA from the resulting pups analyzed by PCR or RT-PCR for evidence of homologous recombination.

Nuclear transplantation also can be used to generate non-human mammals. For example, fetal fibroblasts can be genetically modified such that they contain an inactivated endogenous GSTO2 gene and express a GSTO2 nucleic acid, and then fused with enucleated oocytes. After activation of the oocytes, the eggs are cultured to the blastocyst stage, and implanted into a recipient. See, Cibelli et al. (1998) Science 280:1256-1258. Adult somatic cells, including, for example, cumulus cells and mammary cells, can be used to produce animals such as mice and sheep, respectively. See, for example, Wakayama et al. (1998) Nature 394:369-374; and Wilmut et al. (1997) Nature 385:810-813. Nuclei can be removed from genetically modified adult somatic cells, and transplanted into enucleated oocytes. After activation, the eggs can be cultured to the 2-8 cell stage, or to the blastocyst stage, and implanted into a suitable recipient. Wakayama et al., supra.

Non-human mammals such as mice can be used, for example, to screen toxicity of compounds that are substrates for GSTO2, drugs that alter GSTO2 activity, or for carcinogenesis. For example, GSTO2 activity or toxicity can be assessed in a first group of such non-human mammals in the presence of a compound, and compared with GSTO2 activity or toxicity in a corresponding control group in the absence of the compound. As used herein, suitable compounds include biological macromolecules such as an oligonucleotide (RNA or DNA), or a polypeptide of any length, a chemical compound, a mixture of chemical compounds, or an extract isolated from bacterial, plant, fungal, or animal matter. The concentration of compound to be tested depends on the type of compound and in vitro test data.

Non-human mammals can be exposed to test compounds by any route of administration, including enterally (e.g., orally) and parenterally (e.g., subcutaneously, intravascularly, intramuscularly, or intranasally). Suitable formulations for oral administration can include tablets or capsules prepared by conventional means with pharmaceutically acceptable excipients such as binding agents (e.g., pregelatinized maize starch, polyvinylpyrrolidone or hydroxypropyl methylcellulose); fillers (e.g., lactose, microcrystalline cellulose or calcium hydrogen phosphate); lubricants (e.g., magnesium stearate, talc or silica); disintegrants (e.g., potato starch or sodium starch glycolate); or wetting agents (e.g., sodium lauryl sulfate). Tablets can be coated by methods known in the art. Preparations for oral administration can also be formulated to give controlled release of the compound.

Compounds can be prepared for parenteral administration in liquid form (e.g., solutions, solvents, suspensions, and emulsions) including sterile aqueous or non-aqueous carriers. Aqueous carriers include, without limitation, water, alcohol, saline, and buffered solutions. Examples of non-aqueous carriers include, without limitation, propylene glycol, polyethylene glycol, vegetable oils, and injectable organic esters. Preservatives and other additives such as, for example, antimicrobials, anti-oxidants, chelating agents, inert gases, and the like may also be present. Pharmaceutically acceptable carriers for intravenous administration include solutions containing pharmaceutically acceptable salts or sugars. Intranasal preparations can be presented in a liquid form (e.g., nasal drops or aerosols) or as a dry product (e.g., a powder). Both liquid and dry nasal preparations can be administered using a suitable inhalation device. Nebulised aqueous suspensions or solutions can also be prepared with or without a suitable pH and/or tonicity adjustment.

Genomic DNA generally is used in the analysis of GSTO2 nucleotide sequence variants, although mRNA also can be used. Genomic DNA is typically extracted from a biological sample such as a peripheral blood sample, but can be extracted from other biological samples, including tissues (e.g., mucosal scrapings of the lining of the mouth or from renal or hepatic tissue). Routine methods can be used to extract genomic DNA from a blood or tissue sample, including, for example, phenol extraction. Alternatively, genomic DNA can be extracted with kits such as the QIAamp® Tissue Kit (Qiagen, Chatsworth, Calif.), Wizard® Genomic DNA purification kit (Promega) and the A.S.A.P.™ Genomic DNA isolation kit (Boehringer Mannheim, Indianapolis, Ind.).

Typically, an amplification step is performed before proceeding with the detection method. For example, exons or introns of the GSTO2 gene can be amplified then directly sequenced. Dye primer sequencing can be used to increase the accuracy of detecting heterozygous samples.

Nucleic acid molecules provided herein can be used to detect variant GSTO2 sequences. For example, allele specific hybridization also can be used to detect sequence variants, including complete haplotypes of a subject (e.g., a mammal such as a human). See, Stoneking et al. (1991) Am. J. Hum. Genet. 48:370-382; and Prince et al. (2001) Genome Res. 11:152-162. In practice, samples of DNA or RNA from one or more mammals can be amplified using pairs of primers and the resulting amplification products can be immobilized on a substrate (e.g., in discrete regions). Hybridization conditions are selected such that a nucleic acid probe can specifically bind to the sequence of interest, e.g., the variant nucleic acid sequence. Such hybridizations typically are performed under high stringency as some sequence variants include only a single nucleotide difference. High stringency conditions can include the use of low ionic strength solutions and high temperatures for washing. For example, nucleic acid molecules can be hybridized at 42° C. in 2×SSC (0.3M NaCl/0.03 M sodium citrate/0.1% sodium dodecyl sulfate (SDS) and washed in 0.1×SSC (0.015M NaCl/0.0015 M sodium citrate), 0.1% SDS at 65° C. Hybridization conditions can be adjusted to account for unique features of the nucleic acid molecule, including length and sequence composition. Probes can be labeled (e.g., fluorescently) to facilitate detection. In some embodiments, one of the primers used in the amplification reaction is biotinylated (e.g., 5′ end of reverse primer) and the resulting biotinylated amplification product is immobilized on an avidin or streptavidin coated substrate.

Allele-specific restriction digests can be performed in the following manner. For nucleotide sequence variants that introduce a restriction site, restriction digest with the particular restriction enzyme can differentiate the alleles. For GSTO2 sequence variants that do not alter a common restriction site, mutagenic primers can be designed that introduce a restriction site when the variant allele is present or when the wild type allele is present. A portion of GSTO2 nucleic acid can be amplified using the mutagenic primer and a wild type primer, followed by digest with the appropriate restriction endonuclease.

Certain variants, such as insertions or deletions of one or more nucleotides, change the size of the DNA fragment encompassing the variant. The insertion or deletion of nucleotides can be assessed by amplifying the region encompassing the variant and determining the size of the amplified products in comparison with size standards. For example, a region of GSTO2 can be amplified using a primer set from either side of the variant. One of the primers is typically labeled, for example, with a fluorescent moiety, to facilitate sizing. The amplified products can be electrophoresed through acrylamide gels with a set of size standards that are labeled with a fluorescent moiety that differs from the primer.

PCR conditions and primers can be developed that amplify a product only when the variant allele is present or only when the wild type allele is present (MSPCR or allele-specific PCR). For example, patient DNA and a control can be amplified separately using either a wild type primer or a primer specific for the variant allele. Each set of reactions is then examined for the presence of amplification products using standard methods to visualize the DNA. For example, the reactions can be electrophoresed through an agarose gel and the DNA visualized by staining with ethidium bromide or other DNA intercalating dye. In DNA samples from heterozygous patients, reaction products would be detected in each reaction. Patient samples containing solely the wild type allele would have amplification products only in the reaction using the wild type primer. Similarly, patient samples containing solely the variant allele would have amplification products only in the reaction using the variant primer. Allele-specific PCR also can be performed using allele-specific primers that introduce priming sites for two universal energy-transfer-labeled primers (e.g., one primer labeled with a green dye such as fluorescein and one primer labeled with a red dye such as sulforhodamine). Amplification products can be analyzed for green and red fluorescence in a plate reader. See, Myakishev et al. (2001) Genome 11:163-169.

Mismatch cleavage methods also can be used to detect differing sequences by PCR amplification, followed by hybridization with the wild type sequence and cleavage at points of mismatch. Chemical reagents, such as carbodiimide or hydroxylamine and osmium tetroxide can be used to modify mismatched nucleotides to facilitate cleavage.

Alternatively, GSTO2 variants can be detected by antibodies that have specific binding affinity for variant GSTO2 polypeptides. Variant GSTO2 polypeptides can be produced in various ways, including recombinantly, as discussed above. Host animals such as rabbits, chickens, mice, guinea pigs, and rats can be immunized by injection of a GSTO2 variant polypeptide. Various adjuvants that can be used to increase the immunological response depend on the host species and include Freund's adjuvant (complete and incomplete), mineral gels such as aluminum hydroxide, surface active substances such as lysolecithin, pluronic polyols, polyanions, peptides, oil emulsions, keyhole limpet hemocyanin, and dinitrophenol. Polyclonal antibodies are heterogeneous populations of antibody molecules that are contained in the sera of the immunized animals. Monoclonal antibodies, which are homogeneous populations of antibodies to a particular antigen, can be prepared using a GSTO2 variant polypeptide and standard hybridoma technology. In particular, monoclonal antibodies can be obtained by any technique that provides for the production of antibody molecules by continuous cell lines in culture such as described by Kohler et al. (1975) Nature 256:495, the human B-cell hybridoma technique (Kosbor et al. (1983) Immunology Today 4:72; Cote et al. (1983) Proc. Natl. Acad. Sci. USA 80:2026), and the EBV-hybridoma technique (Cole et al., Monoclonal Antibodies and Cancer Therapy, Alan R. Liss, Inc., pp. 77-96 (1983). Such antibodies can be of any immunoglobulin class including IgG, IgM, IgE, IgA, IgD and any subclass thereof. A hybridoma producing monoclonal antibodies can be cultivated in vitro or in vivo.

Antibody fragments that have specific binding affinity for a GSTO2 variant polypeptide can be generated using known techniques. For example, such fragments include but are not limited to F(ab′)2 fragments that can be produced by pepsin digestion of the antibody molecule, and Fab fragments that can be generated by reducing the disulfide bridges of F(ab′)2 fragments. Alternatively, Fab expression libraries can be constructed. See, for example, Huse et al., Science, 246:1275 (1989). Once produced, antibodies or fragments thereof are tested for recognition of GSTO2 variant polypeptides by standard immunoassay methods including ELISA techniques, radioimmunoassays and Western blotting. See, Short Protocols in Molecular Biology, Chapter 11, Green Publishing Associates and John Wiley & Sons, edited by Ausubel et al., 1992.

Methods

As a result of the present disclosure, it is possible to determine GSH transferase status of a subject (e.g., a mammal such as a human). “GSH transferase status” refers to the ability of a subject to transfer a glutathione group to a substrate. GSH transferase status of a subject can be determined by, for example, measuring the level of GSTO2 activity in the subject using, for example, the methods described herein. Alternatively, GSH transferase status can be evaluated by determining whether a GSH transferase nucleic acid sequence (e.g., a GSTO2 nucleic acid sequence) of a subject contains one or more variants (e.g., one or more variants that are correlated with increased or decreased GSH transferase activity). A variant that results in decreased or increased GSTO2 activity, for example, can be said to result in “reduced” or “enhanced” GSH transferase status, respectively. In some embodiments, the variant profile of a subject can be used to determine the GSH transferase status of the subject.

“Variant profile” refers to the presence or absence of a plurality (e.g., two or more) of GSTO2 nucleotide sequence variants or GSTO2 amino acid sequence variants. For example, a variant profile can include the complete GSTO2 haplotype of the mammal (e.g., see Tables 5-9) or can include the presence or absence of a set of particular non-synonymous SNPs (e.g., single nucleotide substitutions that alter the amino acid sequence of a GSTO2 polypeptide). In one embodiment, the variant profile includes detecting the presence or absence of two or more non-synonymous SNPs (e.g., 2, 3, or 4 non-synonymous SNPs) described herein. There may be ethnic-specific pharmacogenetic variation, as certain of the nucleotide and amino acid sequence variants described herein were detected solely in African-American, Caucasian-American, Han Chinese-American, or Mexican-American subjects. In addition, the variant profile can include detecting the presence or absence of any type of GSTO2 SNP together with any other GSTO2 SNP (e.g., a polymorphism pair or a group of polymorphism pairs). Such polymorphism pairs include, without limitation, the pairs described in Table 4. Further, a variant profile can include detecting the presence or absence of any GSTO2 SNP together with any SNP from other GSH transferases.

GSH transferase activity of an enzyme such as GSTO2 can be measured using, for example, in vitro methods such as those described herein. As used herein, the term “reduced GSH transferase status” refers to a decrease (e.g., a 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 75%, 80%, 90%, 95%, or 100% decrease) in GSH transferase activity (e.g., GSTO2 activity) of a subject, as compared to a control level of GSH transferase activity. Similarly, the term “enhanced GSH transferase status” refers to an increase (e.g., a 5%, 10%, 15%, 20%, 25%, 30%, 40%, 50%, 60%, 70%, 75%, 80%, 90%, 95%, 100%, or more than 100% increase) in GSH transferase activity of a subject, as compared to a control level of GSH transferase activity. A control level of GSH transferase activity can be, for example, an average level of GSH transferase activity in a population of individuals. In one embodiment, the population includes individuals that do not contain particular GSTO2 nucleotide sequence variants or particular GSTO2 amino acid sequence variants (e.g., particular variants that affect GSH transferase status). Alternatively, a control level of GSH transferase activity can refer to the level of GSH transferase activity in a control subject (e.g., a subject that does not contain a GSTO2 nucleic acid containing a variant).

In some embodiments, evaluation of GSH transferase status can be used in diagnostic assays to determine whether a particular therapy may be useful in an individual, or to tailor particular treatment regimens to an individual. For example, the presence of absence of one or more variants may indicate an individual's ability to metabolize arsenic, a drinking water contaminant, or arsenic trioxide, an agent that can be used to treat acute promyelocytic leukemia. In further embodiments, GSH transferase status can be linked to predisposition to a particular disease or condition. For example, identification of GSTO2 sequence variants also may allow predisposition to conditions such as cancer, Alzheimer's, or Parkinson's disease to be assessed in individuals. Additional risk factors including, for example, family history and other genetic factors can be considered when determining risk. Predisposition to such diseases can be determined based on the presence or absence of a single GSTO2 sequence variant or based on a variant profile.

Determination of GSH transferase status and predisposition to conditions such as cancer, Alzheimer's, or Parkinson's disease can include identification of genetic markers (e.g., polymorphisms) in linkage disequilibrium with particular GSTO2 alleles. Although such markers may not be relevant from a functional perspective (i.e., may not directly affect function of GSTO2), their presence can be predictive of functional/clinically relevant polymorphisms. Thus, this document also provides methods for detecting a genotype by screening for genetic markers in linkage disequilibrium with particular GSTO2 alleles. The methods can include providing a nucleic acid sample from a subject (e.g., a human subject), and screening the sample for one or more markers in linkage disequilibrium with a particular GSTO2 allele. For example, a method can include screening for markers in linkage disequilibrium with a SNP at nucleotide 121, 389, 424, or 472 relative to the adenine in the translation initiation codon (nucleotide 23393, 27473, 27508, or 45662, respectively, of SEQ ID NO:1). Methods also can include screening for genetic markers in linkage disequilibrium with any of the other SNPs shown in Table 2.

Articles of Manufacture

Articles of manufacture also are provided herein, and can include populations of isolated GSTO2 nucleic acid molecules or GSTO2 polypeptides immobilized on a substrate. Suitable substrates provide a base for the immobilization of the nucleic acids or polypeptides, and in some embodiments, allow immobilization of nucleic acids or polypeptides into discrete regions. In embodiments in which the substrate includes a plurality of discrete regions, different populations of isolated nucleic acids or polypeptides can be immobilized in each discrete region. Thus, each discrete region of the substrate can include a different GSTO2 nucleic acid or GSTO2 polypeptide sequence variant. Such articles of manufacture can include two or more sequence variants of GSTO2, or can include all of the sequence variants known for GSTO2. For example, the article of manufacture can include two or more of the sequence variants identified herein and one or more other GSTO2 sequence variants, such as nucleic acid variants that occur in the promoter region of the GSTO2 gene. Furthermore, nucleic acid molecules containing sequence variants for other GSH transferases can be included on the substrate.

Suitable substrates can be of any shape or form and can be constructed from, for example, glass, silicon, metal, plastic, cellulose, or a composite. For example, a suitable substrate can include a multiwell plate or membrane, a glass slide, a chip, or polystyrene or magnetic beads. Nucleic acid molecules or polypeptides can be synthesized in situ, immobilized directly on the substrate, or immobilized via a linker, including by covalent, ionic, or physical linkage. Linkers for immobilizing nucleic acids and polypeptides, including reversible or cleavable linkers, are known in the art. See, for example, U.S. Pat. No. 5,451,683 and WO98/20019. Immobilized nucleic acid molecules are typically about 20 nucleotides in length, but can vary from about 10 nucleotides to about 1000 nucleotides in length.

In practice, a sample of DNA or RNA from a subject can be amplified, the amplification product hybridized to an article of manufacture containing populations of isolated nucleic acid molecules in discrete regions, and hybridization can be detected. Typically, the amplified product is labeled to facilitate detection of hybridization. See, for example, Hacia et al. (1996) Nature Genet. 14:441-447; and U.S. Pat. Nos. 5,770,722 and 5,733,729.

The invention will be further described in the following examples, which do not limit the scope of the invention described in the claims.

EXAMPLESExample 1Methods and Materials

PCR Amplification and DNA Sequencing: DNA samples from 60 Caucasian-American, 60 African-American, 60 Han Chinese-American, and 60 Mexican-American subjects were obtained from the Coriell Institute Cell Repository (Camden, N.J.). These samples had been anonymized, and written informed consent had been obtained from all donors for the use of their DNA for this purpose. All experiments were reviewed and approved by the Mayo Clinic Institutional Review Board. Fourteen PCR reactions were performed with each DNA sample to amplify all GSTO2 exons and splice junctions, as well as a portion of the 5′ flanking region. Primers were designed to hybridize within introns, the 5′-FR, or the 3′-UTR. This approach was used to ensure that a GSTO2 processed pseudogene that maps to chromosome 3 (Whitbread et al. (2003) Pharmacogenet. 13:131-144) would not be amplified. Amplicons were sequenced using dye-primer sequencing chemistry to facilitate the identification of heterozygous bases (Chadwick et al. (1996) Biotechniques 20:676-683). To make that possible, universal M13 sequencing tags were added to the 5′-ends of each forward and reverse primer. Forward primers contained the M13 forward sequence (5′-TGTAAAACGACGGCCAGT-3′; SEQ ID NO:10), and reverse primers contained the M13 reverse sequence (5′-CAGGAAACAGCTATGACC-3′; SEQ ID NO:11). In addition, resequencing of the GSTO2 5′-FR and the intergenic area utilized dye terminator chemistry. The sequences and locations of each primer within the gene are listed in Table 1. “F” represents forward; “R,” reverse; “U,” upstream; “D,” downstream; “I,” intron; “FR,” flanking region; and “UTR,” untranslated region. The locations of primers within the gene were chosen to avoid repetitive sequence.

Amplicons were sequenced in the Mayo Molecular Biology Core Facility with an ABI 377 DNA sequencer using BigDye™ (Perkin Elmer) dye-primer sequencing chemistry. Both DNA strands were sequenced in all cases. To exclude PCR-induced artifacts, independent amplification followed by DNA sequencing was performed for all samples in which a SNP was only observed once among the samples resequenced, or for samples with ambiguous chromatograms. DNA sequence chromatograms were analyzed using the PolyPhred 3.0 (Nickerson et al. (1997) Nucl. Acids Res. 25:2745-2751) and Consed 8.0 (Gordon et al. (1998) Genome Res. 8:195-202) programs developed by the University of Washington (Seattle, Wash.). The University of Wisconsin GCG software package, Version 10, was also used to analyze nucleotide sequence. GenBank accession numbers for the GSTO2 reference sequences were NT—030059 and NM—183239.

WT and variant GSTO2 expression constructs were transfected into COS-1 cells, together with pSV-β-galactosidase DNA (Promega) to make it possible to use β-galactosidase activity to correct for transfection efficiency. Transfections were performed with Lipofectamine 2000 (Invitrogen), and the cells were cultured for 48 hours in Dulbecco's modified Eagles medium (DMEM; Bio Whittaker, Walkersville, Md.) with 10% fetal bovine serum (Clontech). Cells were homogenized with a Polytron homogenizer (Brinkmann Instruments, Westbury, N.Y.). The homogenates were centrifuged at 100,000×g for 1 hour; and supernatant preparations were stored at −80° C.

Western Blot Analysis: Rabbit polyclonal antibodies were generated to GSTO2 (Cocalico Biologicals, Reamstown, Pa.) using a polypeptide that corresponded to amino acids 10 to 30, with an additional cysteine residue at the amino terminus. This synthetic polypeptide was conjugated to keyhole limpet hemocyanin prior to immunization of rabbits. For quantitative Western blot analysis, cytosol preparations of COS-1 cells transfected with the expression constructs were loaded on a 12% polyacrylamide gel on the basis of β-galactosidase activity to correct for possible variation in transfection efficiency. After SDS-PAGE, the proteins were transferred to a nitrocellulose membrane that was probed with rabbit antiserum to GSTO2, diluted 1:2000 with “blocking buffer.” Goat anti-rabbit horseradish peroxidase, diluted 1:20,000, was used as the secondary antibody. Bound antibody was detected with the ECL Western blotting system (Amersham Pharmacia Biotech, Piscataway, N.J.). The level of immunoreactive protein was quantitated using IP Gel Lab (Biosystemetica, Plymouth, UK).

In vitro Translation and Degradation: Selected GSTO2 expression constructs were transcribed and translated in vitro in the presence of 35S-methionine (1000 Ci/mmol, 2.5 mCi total activity) (Amersham Pharmacia Biotech) using the TNT coupled rabbit reticulocyte lysate (RRL) system (Promega, Madison, Wis.). Reaction mixtures used to generate radioactively labeled protein were incubated at 30° C. for 90 minutes, and 5 μl aliquots were used to perform SDS-PAGE with 10% gels. Protein degradation studies were performed as described (Wang et al. (2003) Pharmacogenet. 13:555-564). Briefly, 50 μl of an adenosine 5′-triphosphate (ATP) generating system, 50 μl of “untreated” RRL, and 10 μl of 35S-methionine radioactively labeled GSTO2 allozyme were mixed. This mixture was incubated at 37° C., and 10 μl aliquots were removed at various time intervals, followed by SDS-PAGE. After electrophoresis, the gels were dried, exposed to X-ray film and the bands of radioactively labeled protein were quantified.

Data analysis: Statistical comparison of the data was performed by ANOVA using the StatView program, version 4.5 (Abacus Concepts, Inc., Berkeley, Calif.). Values for π, θ, and Tajima's D were calculated as described by Tajima (Tajima (1989) Genet. 123:585-595), followed by correction for length. Linkage analysis was performed after all DNA samples had been genotyped at each of the polymorphic sites observed, using the EH program developed by Terwilliger and Ott, Handbook of Human Genetic Linkage, The Johns Hopkins University Press, Baltimore, pp. 188-193 (1994). D′ values, a quantitative method for reporting linkage data that is independent of allele frequency (Hartl and Clark Principles of Population Genetics, 3rd edition, Sinauer Associates, Inc., (Sunderland, Mass.), pp. 96-106 (1997); and Hedrick Genetics of Populations, 2nd edition, Jones and Bartlett (Sudbury, Mass.), pp. 396-405 (2000)), were calculated. The genotype data also were used to assign inferred haplotypes using a program based on the E-M algorithm (Long et al. (1995) Am. J. Hum. Genet. 56:799-810; Excoffier and Slatkin (1995) Mol. Biol. Evol. 12:921-927; and Schaid et al. (2002) Am. J. Hum. Genet. 70:425-434). Unambiguous haplotype assignment was possible on the basis of genotype for samples that contained no more than one heterozygous polymorphism. Group mean values were compared using student's t-test. Degradation data were plotted using GraphPad Prism (GraphPad Software, San Diego, Calif.).

GSTO2 Enzyme Activity: GSTO2 activity is measured using a bacterially expressed GSTO2 polypeptide. DHA is combined with DHA and GSH, and GSTO2 activity is measured as increased absorbance at 265 nm, which results from reduction of the DHA. This method also can be used to estimate apparent Km values. Alternatively, GSTO2 enzyme activity is measured in vivo (e.g., in CHO cells).

Quantitative RT-PCR: mRNA is isolated from COS-1 cells cotransfected with a control β-Galactosidase reporter and a GSTO2 WT expression construct or GSTO2 variant expression construct using an RNeasy Mini Kit (Qiagen, Valencia, Calif.), according to the manufacturer's instructions. RT-PCR is performed with primers for both GSTO2 and β-Galactosidase as an internal control.

GSTO2 reporter gene constructs and luciferase assay: Luciferase reporter gene constructs are created for the most common human GSTO2 5′-FR haplotypes (e.g., haplotypes having a frequency ≦3% in any of the populations). Specifically, a segment containing about 1000 bp of the GSTO2 5′-FR is amplified from human genomic DNA samples that contained the desired haplotypes. The forward and reverse primers include ACC65I and XhoI restriction sites, respectively, to enable subcloning of the amplicons into pGL-3 Basic (Promega), upstream of the firefly luciferase gene open reading frame (ORF). The insert is sequenced in both directions to ensure that the correct sequence is present.

The luciferase reporter gene constructs are used to transiently transfect primary cells or cell lines (e.g., liver, kidney, skeletal muscle, heart, lung, cervix, ovary, or prostate cells). Cells are transfected with 5 μg purified plasmid DNA with 50 ng pRL-TK (Promega) DNA. The Renilla luciferase activity expressed by pRL-TK is used as a control for transfection efficiency. Cells also are transfected with pGL-3 Basic without insert as a control. Transfection is performed using the TransFast reagent (Promega). After 48 hours, cells are lysed and reporter gene activity is measured using a Promega dual-luciferase assay system. Results are reported as the ratio of firefly luciferase light units to Renilla luciferase light units, and all values are expressed as a percentage of the activity of the pGL3 WT 5′-FR construct. All assays in three cell lines are performed at least in duplicate (e.g., in triplicate).

TABLE 1

PCR primers used for resequencing and site directed mutagenesis of GSTO2

Region

SEQ ID

Primer Name

Amplified

Primer Sequence (5′ to 3′ direction)

NO

Resequencing primers

GSTO1 3′UTR F(730)

5′ FR

TGTAAAACGACGGCCAGTGGGCAGGAGTCAGCAATAAA*

12

UR(−1348)

5′ FR

CAGGAAACAGCTATGACCGGTCTTCAGGGCCAACAATA

13

UF(−1462)

5′ FR

TGTAAAACGACGGCCAGTTGCTCTCTTGTTGCTTGTCTTC

14

UR(−791)

5′ FR

CAGGAAACAGCTATGACCAGCCTGCTGAATCGGAAAT

15

UF(−888)

5′ FR

TGTAAAACGACGGCCAGTTTCGCCATTGAGAGAAACCT

16

I1R(48)

Exon 1 (5′ UTR)

CAGGAAACAGCTATGACCCTCAAACCCCTCTTCCCTTC

17

UF(−368)

Exon 1/Intron 1

TGTAAAACGACGGCCAGTTAGTTGGCGGGTAGGATCAC

18

I1R(539)

Exon 1/Intron 1

CAGGAAACAGCTATGACCGTGGTTTGCGAAGGTTTCAT

19

I1F(1608)

Intron 1

TGTAAAACGACGGCCAGTTCTTTGGTGGCATTATTTTCCTA

20

I1R(2265)

Intron 1

CAGGAAACAGCTATGACCAACCTTTTTGGATTTCACTTTCC

21

I1F(−1423)

Intron 1

TGTAAAACGACGGCCAGTAGCCGGGAACCAATATGTCT

22

I1R(−1042)

Intron 1

CAGGAAACAGCTATGACCCGTAATGGCTGCTCAACAAA

23

I1F(−242)

Exon 2

TGTAAAACGACGGCCAGTGCTGTACTTACTTACCAAAGAGT

24

I2R(21)

Exon 2

CAGGAAACAGCTATGACCGACCCCATGGAGAGCACTCACCT

25

E2F(−34)

Exons 2 and 3

TGTAAAACGACGGCCAGTGAGCTCCGGGAGCTGCGCAAACCA

26

I3R(54)

Exons 2 and 3

CAGGAAACAGCTATGACCGGGCTCCTGTGAGGCGCTGGGTTTGCT

27

I3F(−85)

Exon 4

TGTAAAACGACGGCCAGTCTGATTGCTTCTGCTTTCAAGAAGA

28

I4R(140)

Exon 4

CAGGAAACAGCTATGACCCATGCTAGCTAATCTCCTGAGGAT

29

I4F(−94)

Exon 5

TGTAAAACGACGGCCAGTGCTGGAGTTATAAAGCTTCGCTGCCT

30

I5R(227)

Exon 5

TGTAAAACGACGGCCAGTGAAAGGAAGCAATTCATGGCATGT

31

I5F(−148)

Exon 6

CAGGAAACAGCTATGACCGAGAAGAAGCATTCTCCTAACCAGA

32

I6R(115)

Exon 6

TGTAAAACGACGGCCAGTCTGATTCACTTCACAGACATCCTGA

33

I6F(−97)

Exon 7

CAGGAAACAGCTATGACCCTCATTCTTCATATCTTGCAAATCTA

34

3′UTR R(775)

Exon 7

CAGGAAACAGCTATGACCCAAGCTGGGGAATTCTGGACAGCGA

35

I6F(−15)

Exon 7/3′ FR

TGTAAAACGACGGCCAGTTTTCCTGTCTTGCAGCTGTG

36

3′FR R(1288)

Exon 7/3′ FR

CAGGAAACAGCTATGACCTCCCAAAATGTGCTCTGACA

37

3′FR R(1218)

3′ FR

TGTAAAACGACGGCCAGTCCTGCCTAGCCTGTAGCTTCT

38

3′FR R(1855)

3′ FR

CAGGAAACAGCTATGACCTCAACTTTATCCCCAGTGTGC

39

Mutagenesis primers

O2_V41I_F

AGGACCCGCCTCATCCTCAAG**

40

O2_V41I_R

CTTGAGGATGAGGCGGGTCCT

41

O2_C130Y_F

ACCAAGGAGTACCTGGTAGCGTTG

42

O2_C130Y_R

CAACGCTACCAGGTACTCCTTGGT

43

O2_N142D_F

GAATGCACTAATCTGAAGGCA

44

O2_N142D_R

TGCCTTCAGATTAGTGCATTC

45

O2_L158I_F

GAAGAGATTATTGAGTATCAG

46

O2_L158I_R

CTGATACTCAATAATCTCTTC

47

*Underlined nucleotides in resequencing primers indicate M13 tag

**Underlined nucleotides in mutagenesis primers indicate variant

Example 2GSTO2 Polymorphisms

Fourteen separate PCR amplifications were performed for each of the DNA samples studied. All PCR amplicons were sequenced on both strands, making it possible to verify the presence of polymorphisms using data from the complimentary strand. Approximately 3×106 pg of DNA was sequenced and analyzed, including the “intergenic” area located between the GSTO1 and GSTO2 genes (i.e., between the final exon of GSTO1 and the first exon of GSTO2). A total of 66 polymorphisms were observed (Table 2). Polymorphisms in exons and flanking regions (FR) are numbered relative to the adenine in the GSTO2 translation initiation codon (ATG, adenine is +1). Polymorphisms in introns are numbered separately, either as positive numbers relative to the guanine in the splice donor site (GT, guanine is +1), or as negative numbers relative to the guanine in the splice acceptor site (AG, guanine is −1).

Variant allele frequencies ranged from 0.8% to 78.3%, with differences between the African-American (AA), Caucasian-American (CA), Han Chinese-American (HCA), and Mexican-American (MA) subjects. Fifty-four SNPs were observed in AA subjects, 19 in CA subjects, 17 in HCA subjects, and 36 in MA subjects. Nineteen of these polymorphisms were specific to AA, 5 to CA, and 2 each were specific to HCA and MA subjects.

Six SNPs were observed within the coding-region (cSNPs), and four of those cSNPs—located in exons 2, 4, and 5—were nonsynonymous, resulting in the amino acid alterations Val41Ile, Cys130Tyr, Asn142Asp, and Leu158Ile. The Val41Ile polymorphism had a frequency of 0.9% in African-Americans but was not observed in DNA from Caucasian-American, Han Chinese-American, or Mexican-American subjects. The Cys130Tyr polymorphism had frequencies of 0.8% in African Americans and 1.7% in Han Chinese-Americans, but was not observed in Caucasian-Americans or Mexican-Americans. The Asn142Asp polymorphism had frequencies of 23.3% in African-Americans, 69.2% in Caucasian-Americans, 78.3% in Han Chinese-Americans, and 74.2% in Mexican-Americans. The Leu158Ile polymorphism had a frequency of 0.8% in Caucasian-Americans, but was not observed in African-Americans, Han Chinese-Americans, or Mexican-Americans. To exclude artifacts introduced by PCR-dependent misincorporation, independent amplifications were performed and the amplicons were sequenced in all cases in which a polymorphism was observed only once among the DNA samples studied. All polymorphisms were in Hardy-Weinberg equilibrium (P>0.05).

Nucleotide diversity, a measure of genetic variation that is adjusted for the number of alleles studied, was calculated for the GSTO2 gene and for the intergenic region. Two standard measures of nucleotide diversity are π, average heterozygosity per site, and θ, a population mutation measure that is theoretically equal to the neutral mutation parameter (Fullerton et al. (2000) Am. J. Hum. Genet. 67:881-900). Values for Tajima's D, a test for the neutral mutation hypothesis (Fullerton et al., supra), also were calculated (Table 3). Under conditions of neutrality, Tajima's D should equal 0. For GSTO2, none of the values for Tajima's D were statistically significant. For the intergenic region, Tajima's D value for AA subjects was statistically significant, raising the possibility that this area may be undergoing selection.

TABLE 2

Human GSTO2 polymorphisms and frequencies

Polymorphism

Location

Position in

Amino Acid

WT

Variant Sequence

Position

in Gene

SEQ ID 1

Change

Nucleotide

Nucleotide

AA

CA

HCA

MA

−2012-2013

5′FR

15569/15570

Insertion of G

0.292

0.000

0.000

0.008

−1785

5′FR

15797

T

C

0.050

0.000

0.000

0.000

−1771

5′FR

15811

G

C

0.008

0.000

0.000

0.000

−1664

5′FR

15918

C

T

0.300

0.000

0.000

0.008

−1598

5′FR

15984

C

G

0.300

0.000

0.000

0.008

−1576

5′FR

16006

G

C

0.300

0.000

0.000

0.008

−1566

5′FR

16016

G

A

0.300

0.000

0.000

0.008

−1541

5′FR

16041

C

T

0.008

0.000

0.000

0.000

−1357

5′FR

16225

C

A

0.300

0.000

0.000

0.008

−1189

5′FR

16393

T

C

0.300

0.000

0.000

0.008

−1179

5′FR

16403

G

T

0.008

0.000

0.000

0.000

−1105

5′FR

16477

G

A

0.117

0.258

0.167

0.242

−1102

5′FR

16480

T

G

0.217

0.683

0.767

0.692

−1009

5′FR

16573

C

T

0.300

0.000

0.000

0.008

−781

5′FR

16801

A

C

0.283

0.000

0.000

0.008

−730-731

5′FR

16851/16852

Insertion of T

0.283

0.000

0.000

0.008

−637

5′FR

16945

C

T

0.275

0.000

0.000

0.008

−634

5′FR

16948

C

G

0.275

0.000

0.000

0.008

−625

Exon 1

16957

A

G

0.000

0.000

0.017

0.000

−530

Exon 1

17052

C

G

0.271

0.000

0.000

0.008

−526

Exon 1

17056

C

A

0.000

0.008

0.000

0.000

−463

Exon 1

17119

G

A

0.271

0.000

0.000

0.008

−450

Exon 1

17132

A

C

0.271

0.000

0.000

0.008

−424

Exon 1

17158

G

A

0.017

0.000

0.000

0.000

−410

Exon 1

17172

T

G

0.008

0.000

0.000

0.000

−293

Exon 1

17289

A

T

0.050

0.000

0.000

0.000

−234

Exon 1

17348

C

T

0.008

0.000

0.000

0.000

52

Intron 1

17402

A

G

0.008

0.000

0.000

0.000

152

Intron 1

17502

G

T

0.033

0.000

0.000

0.000

310-311

Intron 1

17660/17661

Deletion of CG

0.008

0.000

0.000

0.000

311

Intron 1

17661

G

A

0.300

0.000

0.000

0.008

333

Intron 1

17683

T

C

0.008

0.000

0.000

0.000

402

Intron 1

17752

C

T

0.042

0.025

0.017

0.025

426

Intron 1

17776

T

C

0.050

0.000

0.000

0.000

435

Intron 1

17785

T

A

0.017

0.000

0.000

0.000

477-478

Intron 1

17827/17828

Insertion of T

0.008

0.033

0.033

0.000

1694

Intron 1

19044

A

G

0.050

0.000

0.000

0.008

2003

Intron 1

19353

A

G

0.300

0.000

0.000

0.008

2146

Intron 1

19496

A

G

0.000

0.017

0.000

0.008

2168

Intron 1

19518

G

A

0.200

0.000

0.017

0.017

−1131

Intron 1

21635

C

T

0.417

0.025

0.017

0.025

−185

Intron 1

22581

A

G

0.060

0.000

0.000

0.000

−19

Intron 1

22747

T

A

0.000

0.000

0.000

0.008

−12

Intron 1

22754

G

A

0.008

0.000

0.000

0.008

−183

Exon 2

22814

A

G

0.108

0.242

0.175

0.183

−16

Intron 2

23291

C

T

0.008

0.000

0.000

0.000

121

Exon 3

23393

Val41Ile

G

A

0.009

0.000

0.000

0.000

20

Intron 4

26217

T

C

0.200

0.560

0.642

0.667

39

Intron 4

26236

C

T

0.000

0.000

0.000

0.008

−38

Intron 4

27413

A

C

0.000

0.017

0.000

0.017

389

Exon 5

27473

Cys130Tyr

G

A

0.008

0.000

0.017

0.000

424

Exon 5

27508

Asn142Asp

A

G

0.233

0.692

0.783

0.742

74

Intron 5

27626

A

G

0.233

0.692

0.783

0.742

170

Intron 5

27722

T

C

0.000

0.008

0.000

0.000

472

Exon 6

45662

Leu158Ile

C

A

0.000

0.008

0.000

0.000

−55

Intron 6

47154

C

G

0.058

0.000

0.008

0.000

−21

Intron 6

47188

C

T

0.108

0.000

0.008

0.000

591

Exon 7

47224

G

T

0.167

0.033

0.133

0.017

630

Exon 7

47263

C

T

0.117

0.000

0.000

0.008

1038

3′FR

47671

C

T

0.175

0.033

0.133

0.017

1080

3′FR

47713

C

T

0.000

0.008

0.000

0.008

1243

3′FR

47876

A

G

0.000

0.017

0.000

0.000

1343

3′FR

47976

T

G

0.125

0.000

0.000

0.000

1615

3′FR

48248

C

T

0.000

0.008

0.000

0.000

1649

3′FR

48282

T

C

0.000

0.000

0.000

0.008

1724

3′FR

48357

C

T

0.008

0.000

0.000

0.000

TABLE 3

Nucleotide diversity

π × 10−3*

θ × 10−3

Tajima's D

P value**

GSTO2

AA

1.30 ± 0.70

1.40 ± 0.40

−0.37

0.72

CA

0.50 ± 0.30

0.40 ± 0.20

0.62

0.55

HCA

0.50 ± 0.30

0.60 ± 0.20

−0.67

0.52

MA

0.50 ± 0.30

0.80 ± 0.30

−1.30

0.20

Intergenic region

AA

3.80 ± 2.03

2.20 ± 0.73

1.94

0.05

CA

0.55 ± 0.44

0.25 ± 0.18

1.79

0.07

HCA

0.43 ± 0.37

0.25 ± 0.18

1.10

0.30

MA

0.66 ± 0.50

1.70 ± 0.60

−1.65

0.10

*Values are parameter estimates mean ± SE.

**P values refer to Tajima's D.

Example 3Linkage Disequilibrium Analysis and Haplotype Analysis

Linkage disequilibrium analysis was performed after all of the DNA samples had been genotyped at each of the polymorphic sites. Pairwise combinations of these polymorphisms were tested for linkage disequilibrium using the EH program developed by Terwilliger and Ott, supra. The output of this program was used to calculate d′ values, a method for reporting linkage data that is independent of sample size. Pairwise combinations with a statistically significant linkage disequilibrium (P value<0.001) are shown in Table 4. D′ values greater than 0 indicate a positive association, while d′ values less than 0 indicate a negative association.

The genotype data also were used for haplotype analysis. Haplotypes typically are defined as combinations of alleles on a single chromosome, although they are sometimes referred to in a more restricted definition as all polymorphisms present on a single allele (Altshuler et al. Nature 437:1299-1320 (2005)). Haplotypes can be determined unequivocally if not more than one polymorphism per allele is heterozygous, but haplotypes also can be inferred computationally (Schaid, supra). GSTO2 haplotypes, both observed and inferred, with frequencies ≧0.08%, are listed in Tables 5-8. As shown in the tables, the identified haplotypes accounted for 88.9% of all DNA samples from Caucasian-American subjects, 80.8% of all DNA samples from African-American subjects, 95.0% of all DNA samples from Han Chinese-American subjects, and 90.9% of all DNA samples from Mexican-American subjects. A summary of haplotypes having frequencies ≧1% is presented in Table 9. Haplotypes including one or more cSNPs are included even if their frequency is less than 1%.

TABLE 4

Linkage disequilibrium statistics for GSTO2

D′

Location 1

Location 2

value

P value

African American Population

5′FR(−2012)

5′FR(−1664)

1

0

5′FR(−2012)

5′FR(−1598)

1

0

5′FR(−2012)

5′FR(−1576)

1

0

5′FR(−2012)

5′FR(−1566)

1

0

5′FR(−2012)

5′FR(−1357)

0.96

0

5′FR(−2012)

5′FR(−1189)

1

0

5′FR(−2012)

5′FR(−1102)

1

0

5′FR(−2012)

5′FR(−1009)

1

0

5′FR(−2012)

5′FR(−781)

0.96

0

5′FR(−2012)

5′FR(−730)

0.96

0

5′FR(−2012)

5′FR(−637)

0.96

0

5′FR(−2012)

5′FR(−634)

0.96

0

5′FR(−2012)

5′FR(−530)

0.96

0

5′FR(−2012)

5′FR(−463)

0.96

0

5′FR(−2012)

5′FR(−450)

0.96

0

5′FR(−2012)

Intron 1(311)

1

0

5′FR(−2012)

Intron 6(−55)

1

0

5′FR(−2012)

Intron 6(−21)

1

0

5′FR(−2012)

3′FR(1343)

0.77

0.0001

5′FR(−1785)

Exon 1(−293)

1

0

5′FR(−1785)

Intron 1(426)

1

5′FR(−1785)

Intron 1(−185)

1

0

5′FR(−1785)

Exon 7(630)

1

0

5′FR(−1664)

5′FR(−1598)

1

0

5′FR(−1664)

5′FR(−1576)

1

0

5′FR(−1664)

5′FR(−1566)

1

0

5′FR(−1664)

5′FR(−1357)

1

0

5′FR(−1664)

5′FR(−1189)

1

0

5′FR(−1664)

5′FR(−1109)

1

0

5′FR(−1664)

5′FR(−781)

1

0

5′FR(−1664)

5′FR(−730)

1

0

5′FR(−1664)

5′FR(−637)

1

0

5′FR(−1664)

5′FR(−634)

1

0

5′FR(−1664)

5′FR(−530)

1

0

5′FR(−1664)

5′FR(−463)

1

0

5′FR(−1664)

5′FR(−450)

1

0

5′FR(−1664)

Intron 1(311)

1

0

5′FR(−1664)

Intron 1(2003)

1

0

5′FR(−1664)

Intron 6(−55)

1

0.0001

5′FR(−1664)

Intron 6(−21)

1

0

5′FR(−1664)

3′FR(1343)

0.76

0.0001

5′FR(−1598)

5′FR(−1576)

1

0

5′FR(−1598)

5′FR(−1566)

1

0

5′FR(−1598)

5′FR(−1357)

1

0

5′FR(−1598)

5′FR(−1189)

1

0

5′FR(−1598)

5′FR(−1009)

1

0

5′FR(−1598)

5′FR(−781)

1

0

5′FR(−1598)

5′FR(−730)

1

0

5′FR(−1598)

5′FR(−637)

1

0

5′FR(−1598)

5′FR(−634)

1

0

5′FR(−1598)

5′FR(−530)

1

0

5′FR(−1598)

5′FR(−463)

1

0

5′FR(−1598)

5′FR(−450)

1

0

5′FR(−1598)

Intron 1(311)

1

0

5′FR(−1598)

Intron 1(2003)

1

0

5′FR(−1598)

Exon 5(424)

1

0

5′FR(−1598)

Intron 5(74)

1

0

5′FR(−1598)

Intron 6(−21)

1

0

5′FR(−1598)

3′FR(1343)

0.76

0.0001

5′FR(−1576)

5′FR(1566)

1

0

5′FR(−1576)

5′FR(1357)

1

0

5′FR(−1576)

5′FR(1189)

1

0

5′FR(−1576)

5′FR(1009)

1

0

5′FR(−1576)

5′FR(−781)

1

0

5′FR(−1576)

5′FR(−730)

1

0

5′FR(−1576)

5′FR(−637)

1

0

5′FR(−1576)

5′FR(−634)

1

0

5′FR(−1576)

5′FR(−530)

1

0

5′FR(−1576)

5′FR(−463)

1

0

5′FR(−1576)

5′FR(−450)

1

0

5′FR(−1576)

Intron 1(311)

1

0

5′FR(−1576)

Intron 1(2003)

1

0

5′FR(−1576)

Intron 6(−55)

1

0.0001

5′FR(−1576)

Intron 6(−21)

1

0

5′FR(−1576)

3′FR(1343)

0.76

0.0001

5′FR(−1566)

5′FR(−1357)

1

0

5′FR(−1566)

5′FR(−1189)

1

0

5′FR(−1566)

5′FR(−1009)

1

0

5′FR(−1566)

5′FR(−781)

0.96

0

5′FR(−1566)

5′FR(−730)

0.96

0

5′FR(−1566)

5′FR(−637)

0.96

0

5′FR(−1566)

5′FR(−634)

0.96

0

5′FR(−1566)

Exon 1(−530)

0.96

0

5′FR(−1566)

Exon 1(−463)

0.96

0

5′FR(−1566)

Exon 1(−450)

0.96

0

5′FR(−1566)

Intron 1(311)

1

0

5′FR(−1566)

Intron 1(2003)

1

0

5′FR(−1566)

Intron 6(−21)

0.87

0

5′FR(−1566)

3′FR(1343)

0.65

0.0004

5′FR(−1541)

Exon 1(−410)

1

0

5′FR(−1357)

5′FR(−1189)

1

0

5′FR(−1357)

5′FR(−1009)

1

0

5′FR(−1357)

5′FR(−781)

1

0

5′FR(−1357)

5′FR(−730)

1

0

5′FR(−1357)

5′FR(−637)

1

0

5′FR(−1357)

5′FR(−634)

1

0

5′FR(−1357)

Exon 1(−530)

1

0

5′FR(−1357)

Exon 1(−463)

1

0

5′FR(−1357)

Intron 1(311)

1

0

5′FR(−1357)

Intron 1(2003)

1

0

5′FR(−1357)

Intron 6(−55)

1

0.0001

5′FR(−1357)

Intron 6(−21)

1

0

5′FR(−1357)

3′FR(1343)

0.76

0.0001

5′FR(−1189)

5′FR(−1009)

1

0

5′FR(−1189)

5′FR(−781)

1

0

5′FR(−1189)

5′FR(−730)

1

0

5′FR(−1189)

5′FR(−637)

1

0

5′FR(−1189)

5′FR(−634)

1

0

5′FR(−1189)

Exon 1(−530)

1

0

5′FR(−1189)

Exon 1(−463)

1

0

5′FR(−1189)

Exon 1(−450)

1

0

5′FR(−1189)

Intron 1(311)

1

0

5′FR(−1189)

Intron 1(2003)

1

0

5′FR(−1189)

Intron 6(−55)

1

0.0001

5′FR(−1189)

Intron 6(−21)

1

0

5′FR(−1189)

3′FR(1343)

0.76

0.0001

5′FR(−1105)

Exon 2(−183)

1

0

5′FR(−1102)

Intron 4(20)

1

0

5′FR(−1102)

Exon 5(424)

1

0

5′FR(−1102)

Intron 5(74)

1

0

5′FR(−1009)

5′FR(−781)

1

0

5′FR(−1009)

5′FR(−730)

1

0

5′FR(−1009)

5′FR(−637)

1

0

5′FR(−1009)

5′FR(−634)

1

0

5′FR(−1009)

Exon 1(530)

1

0

5′FR(−1009)

Exon 1(463)

1

0

5′FR(−1009)

Exon 1(450)

1

0

5′FR(−1009)

Intron 1(311)

1

0

5′FR(−1009)

Intron 1(2003)

1

0

5′FR(−1009)

Intron 6(−55)

1

0.0001

5′FR(−1009)

Intron 6(−21)

1

0

5′FR(−1009)

3′FR(1343)

0.76

0.0001

5′FR(−781)

5′FR(−730)

1

0

5′FR(−781)

5′FR(−637)

1

0

5′FR(−781)

5′FR(−634)

1

0

5′FR(−781)

Exon 1(−530)

1

0

5′FR(−781)

Exon 1(−463)

1

0

5′FR(−781)

Exon 1(−450)

1

0

5′FR(−781)

Intron 1(311)

1

0

5′FR(−781)

Intron 1(2003)

1

0

5′FR(−781)

Intron 6(−21)

0.87

0

5′FR(−730)

5′FR(−637)

1

0

5′FR(−730)

5′FR(−634)

1

0

5′FR(−730)

Exon 1(−530)

1

0

5′FR(−730)

Exon 1(−463)

1

0

5′FR(−730)

Exon 1(−450)

1

0

5′FR(−730)

Intron 1(311)

1

0

5′FR(−730)

Intron 1(2003)

1

0

5′FR(−730)

Intron 6(−21)

0.87

0

5′FR(637)

5′FR(634)

1

0

5′FR(637)

Exon 1(−530)

1

0

5′FR(637)

Exon 1(−463)

1

0

5′FR(637)

Exon 1(−450)

1

0

5′FR(637)

Intron 1(311)

1

0

5′FR(637)

Intron 1(2003)

1

0

5′FR(637)

Intron 6(−21)

0.87

0

5′FR(634)

Exon 1(−530)

1

0

5′FR(634)

Exon 1(−463)

1

0

5′FR(634)

Exon 1(−450)

1

0

5′FR(634)

Intron 1(311)

1

0

5′FR(634)

Intron 1(2003)

1

0

5′FR(634)

Intron 6(−21)

0.87

0

5′FR(530)

Exon 1(−463)

1

0

5′FR(530)

Exon 1(−450)

1

0

5′FR(530)

Intron 1(2003)

1

0

5′FR(530)

Intron 6(−21)

0.87

0

Exon 1(−463)

Exon 1(−450)

1

0

Exon 1(−463)

Intron 1(311)

1

0

Exon 1(−463)

Intron 1(2003)

1

0

Exon 1(−463)

Intron 6(−21)

0.87

0

Exon 1(−450)

Intron 1(311)

1

0

Exon 1(−450)

Intron 1(2003)

1

0

Exon 1(−450)

Intron 6(−21)

0.87

0

Exon 1(−293)

Intron 1(426)

1

0

Exon 1(−293)

Intron 1(−185)

1

0

Exon 1(−293)

Exon 7(630)

1

0

Intron 1(52)

Intron 1(310)

1

0

Intron 1(52)

Intron 1(−12)

1

0

Intron 1(152)

Intron 6(−55)

1

0

Intron 1(310)

Intron 1(−12)

1

0

Intron 1(311)

Intron 1(2003)

1

0

Intron 1(311)

Intron 6(−21)

1

0

Intron 1(402)

Intron 1(−1131)

1

0

Intron 1(426)

Intron 1(−185)

1

0

Intron 1(426)

Exon 7(630)

1

0

Intron 1(1694)

Exon 7(630)

1

0

Intron 1(2003)

Intron 6(−55)

1

0.0001

Intron 1(2003)

Intron 6(−21)

1

0

Intron 1(2003)

3′FR(1343)

0.76

0.0001

Intron 1(2168)

Exon 7(591)

0.59

0

Intron 1(2168)

3′FR(1038)

0.61

0

Intron 1(−185)

Exon 7(630)

0.82

0

Intron 4(20)

Exon 5(424)

1

0

Intron 4(20)

Intron 5(74)

1

0

Exon 5(424)

Intron 5(74)

1

0

Intron 6(−55)

Intron 6(−21)

1

0

Intron 6(−55)

3′FR(1343)

1

0

Intron 6(−21)

3′FR(1343)

1

0

Exon 7(591)

3′FR(1038)

1

0

Caucasian American Population

5′FR(−2012)

Intron 1(2146)

0.49

0.0002

5′FR(−1105)

5′FR(−1102)

1

0

5′FR(−1105)

Exon 2(−183)

0.9

0

5′FR(−1105)

Intron 4(20)

0.79

0

5′FR(−1105)

Exon 5(424)

0.8

0

5′FR(−1105)

Intron 5(74)

0.8

0

5′FR(−1102)

Intron 1(402)

0.85

0.0001

5′FR(−1102)

Intron 1(−1131)

0.85

0.0001

5′FR(−1102)

Exon 2(−183)

0.89

0

5′FR(−1102)

Intron 4(20)

0.77

0.0002

5′FR(−1102)

Intron 4(−38)

0.78

0

5′FR(−1102)

Exon 5(424)

0.82

0

5′FR(−1102)

Intron 5(74)

0.82

0

Exon 1(−526)

Exon 7(591)

1

0.0002

Intron 1(402)

Intron 1(−1131)

1

0

Intron 1(402)

Intron 4(−38)

1

0

Intron 1(477)

Exon 7(591)

0.74

0

Intron 1(477)

3′FR(1038)

1

0

Intron 1(477)

3′FR(1243)

1

0

Intron 1(−1131)

Intron 4(−38)

1

0

Exon 2(−183)

Intron 4(20)

1

0

Exon 2(−183)

Exon 5(424)

1

0

Exon 2(−183)

Intron 5(74)

1

0

Intron 4(20)

Exon 5(424)

1

0

Intron 4(20)

Intron 5(74)

1

0

Exon 5(424)

Intron 5(74)

1

0

Exon 7(591)

3′FR(1038)

0.74

0

Mexican American Population

5′FR(−2012)

5′FR(−1664)

1

0

5′FR(−2012)

5′FR(−1598)

1

0

5′FR(−2012)

5′FR(−1576)

1

0

5′FR(−2012)

5′FR(−1566)

1

0

5′FR(−2012)

5′FR(−1357)

1

0

5′FR(−2012)

5′FR(−1189)

1

0

5′FR(−2012)

5′FR(−1009)

1

0

5′FR(−2012)

5′FR(−781)

1

0

5′FR(−2012)

5′FR(−730)

1

0

5′FR(−2012)

5′FR(−637)

1

0

5′FR(−2012)

5′FR(−634)

1

0

5′FR(−2012)

Exon 1(−530)

1

0

5′FR(−2012)

Exon 1(−463)

1

0

5′FR(−2012)

Exon 1(−450)

1

0

5′FR(−2012)

Intron 1(311)

1

0

5′FR(−2012)

Intron 1(2003)

1

0

5′FR(−2012)

3′FR(1649)

1

0

5′FR(−1664)

5′FR(1598)

1

0

5′FR(−1664)

5′FR(1576)

1

0

5′FR(−1664)

5′FR(1566)

1

0

5′FR(−1664)

5′FR(1357)

1

0

5′FR(−1664)

5′FR(1189)

1

0

5′FR(−1664)

5′FR(1009)

1

0

5′FR(−1664)

5′FR(781)

1

0

5′FR(−1664)

5′FR(730)

1

0

5′FR(−1664)

5′FR(637)

1

0

5′FR(−1664)

5′FR(634)

1

0

5′FR(−1664)

Exon 1(−530)

1

0

5′FR(−1664)

Exon 1(−463)

1

0

5′FR(−1664)

Exon 1(−450)

1

0

5′FR(−1664)

Intron 1(311)

1

0

5′FR(−1664)

Intron 1(2003)

1

0

5′FR(−1664)

3′FR(1649)

1

0

5′FR(−1598)

5′FR(−1576)

1

0

5′FR(−1598)

5′FR(−1566)

1

0

5′FR(−1598)

5′FR(−1357)

1

0

5′FR(−1598)

5′FR(−1189)

1

0

5′FR(−1598)

5′FR(−1009)

1

0

5′FR(−1598)

5′FR(−781)

1

0

5′FR(−1598)

5′FR(−730)

1

0

5′FR(−1598)

5′FR(−637)

1

0

5′FR(−1598)

5′FR(−634)

1

0

5′FR(−1598)

Exon 1(−530)

1

0

5′FR(−1598)

Exon 1(−463)

1

0

5′FR(−1598)

Exon 1(−450)

1

0

5′FR(−1598)

Intron 1(311)

1

0

5′FR(−1598)

Intron 1(2003)

1

0

5′FR(−1598)

3′FR(1649)

1

0

5′FR(−1576)

5′FR(−1566)

1

0

5′FR(−1576)

5′FR(−1357)

1

0

5′FR(−1576)

5′FR(−1189)

1

0

5′FR(−1576)

5′FR(−1009)

1

0

5′FR(−1576)

5′FR(−781)

1

0

5′FR(−1576)

5′FR(−730)

1

0

5′FR(−1576)

5′FR(−637)

1

0

5′FR(−1576)

5′FR(−634)

1

0

5′FR(−1576)

Exon 1(−530)

1

0

5′FR(−1576)

Exon 1(−463)

1

0

5′FR(−1576)

Exon 1(−450)

1

0

5′FR(−1576)

Intron 1(311)

1

0

5′FR(−1576)

Intron 1(2003)

1

0

5′FR(−1576)

3′FR(1649)

1

0

5′FR(−1566)

5′FR(−1357)

1

0

5′FR(−1566)

5′FR(−1189)

1

0

5′FR(−1566)

5′FR(−1009)

1

0

5′FR(−1566)

5′FR(−781)

1

0

5′FR(−1566)

5′FR(−730)

1

0

5′FR(−1566)

5′FR(−637)

1

0

5′FR(−1566)

5′FR(−634)

1

0

5′FR(−1566)

Exon 1(−530)

1

0

5′FR(−1566)

Exon 1(−463)

1

0

5′FR(−1566)

Exon 1(−450)

1

0

5′FR(−1566)

Intron 1(311)

1

0

5′FR(−1566)

Intron 1(2003)

1

0

5′FR(−1566)

3′FR(1649)

1

0

5′FR(−1357)

5′FR(−1189)

1

0

5′FR(−1357)

5′FR(−1109)

1

0

5′FR(−1357)

5′FR(−781)

1

0

5′FR(−1357)

5′FR(−730)

1

0

5′FR(−1357)

5′FR(−637)

1

0

5′FR(−1357)

5′FR(−634)

1

0

5′FR(−1357)

Exon 1(−530)

1

0

5′FR(−1357)

Exon 1(−463)

1

0

5′FR(−1357)

Exon 1(−450)

1

0

5′FR(−1357)

Intron 1(311)

1

0

5′FR(−1357)

Intron 1(2003)

1

0

5′FR(−1357)

3′FR(1649)

1

0

5′FR(−1189)

5′FR(−1009)

1

0

5′FR(−1189)

5′FR(−781)

1

0

5′FR(−1189)

5′FR(−730)

1

0

5′FR(−1189)

5′FR(−637)

1

0

5′FR(−1189)

5′FR(−634)

1

0

5′FR(−1189)

Exon 1(−530)

1

0

5′FR(−1189)

Exon 1(−463)

1

0

5′FR(−1189)

Exon 1(−450)

1

0

5′FR(−1189)

Intron 1(311)

1

0

5′FR(−1189)

Intron 1(2003)

1

0

5′FR(−1189)

3′FR(1649)

1

0

5′FR(−1105)

5′FR(−1102)

1

0

5′FR(−1105)

Exon 2(−183)

1

0

5′FR(−1105)

Intron 4(20)

0.76

0

5′FR(−1105)

Exon 5(424)

0.75

0

5′FR(−1105)

Intron 5(74)

0.75

0

5′FR(−1102)

Exon 2(−183)

1

0

5′FR(−1102)

Intron 4(20)

0.78

0

5′FR(−1102)

Exon 5(424)

1

0

5′FR(−1102)

Intron 5(74)

1

0

5′FR(−1009)

5′FR(−781)

1

0

5′FR(−1009)

5′FR(−730)

1

0

5′FR(−1009)

5′FR(−637)

1

0

5′FR(−1009)

5′FR(−634)

1

0

5′FR(−1009)

Exon 1(−530)

1

0

5′FR(−1009)

Exon 1(−463)

1

0

5′FR(−1009)

Exon 1(−450)

1

0

5′FR(−1009)

Intron 1(311)

1

0

5′FR(−1009)

Intron 1(2003)

1

0

5′FR(−1009)

3′FR(1649)

1

0

5′FR(−781)

5′FR(−730)

1

0

5′FR(−781)

5′FR(−637)

1

0

5′FR(−781)

5′FR(−634)

1

0

5′FR(−781)

Exon 1(−530)

1

0

5′FR(−781)

Exon 1(−463)

1

0

5′FR(−781)

Exon 1(−450)

1

0

5′FR(−781)

Intron 1(311)

1

0

5′FR(−781)

Intron 1(2003)

1

0

5′FR(−781)

3′FR(1649)

1

0

5′FR(−730)

5′FR(−637)

1

0

5′FR(−730)

5′FR(−634)

1

0

5′FR(−730)

Exon 1(−530)

1

0

5′FR(−730)

Exon 1(−463)

1

0

5′FR(−730)

Exon 1(−450)

1

0

5′FR(−730)

Intron 1(311)

1

0

5′FR(−730)

Intron 1(2003)

1

0

5′FR(−730)

3′FR(1649)

1

0

5′FR(−637)

5′FR(−634)

1

0

5′FR(−637)

Exon 1(−530)

1

0

5′FR(−637)

Exon 1(−463)

1

0

5′FR(−637)

Exon 1(−450)

1

0

5′FR(−637)

Intron 1(311)

1

0

5′FR(−637)

Intron 1(2003)

1

0

5′FR(−637)

3′FR(1649)

1

0

5′FR(−634)

Exon 1(−530)

1

0

5′FR(−634)

Exon 1(−463)

1

0

5′FR(−634)

Exon 1(−450)

1

0

5′FR(−634)

Intron 1(311)

1

0

5′FR(−634)

Intron 1(2003)

1

0

5′FR(−634)

3′FR(1649)

1

0

Exon 1(−530)

Exon 1(−463)

1

0

Exon 1(−530)

Exon 1(−450)

1

0

Exon 1(−530)

Intron 1(311)

1

0

Exon 1(−530)

Intron 1(2003)

1

0

Exon 1(−530)

3′FR(1649)

1

0

Exon 1(−463)

Exon 1(−450)

1

0

Exon 1(−463)

Intron 1(311)

1

0

Exon 1(−463)

Intron 1(2003)

1

0

Exon 1(−463)

3′FR(1649)

1

0

Exon 1(−450)

Intron 1(311)

1

0

Exon 1(−450)

Intron 1(2003)

1

0

Exon 1(−450)

3′FR(1649)

1

0

Intron 1(311)

Intron 1(2003)

1

0

Intron 1(311)

3′FR(1649)

1

0

Intron 1(402)

Intron 1(−1131)

1

0

Intron 1(402)

Intron 1(−38)

1

0

Intron 1(1694)

Exon 7(630)

1

0

Intron 1(2003)

3′FR(1649)

1

0

Intron 1(2168)

Exon 7(591)

0.49

0.0002

Intron 1(2168)

3′FR(1038)

0.49

0.0002

Intron 1(1131)

Intron 4(−38)

1

0

Exon 2(−183)

Intron(20)

1

0

Exon 2(−183)

Exon 5(424)

1

0

Exon 2(−183)

Intron 5(74)

1

0

Intron 4(20)

Exon 5(424)

1

0

Intron 4(20)

Intron 5(74)

1

0

Exon 5(424)

Intron 5(74)

1

0

Exon 7(591)

3′FR(1038)

1

0

Han Chinese American Population

5′FR(−1105)

5′FR(−1102)

1

0

5′FR(−1105)

Exon 2(−183)

1

0

5′FR(−1105)

Exon 2(424)

0.86

0

5′FR(−1105)

Intron 5(74)

0.86

0

5′FR(−1102)

Exon 2(−183)

1

0

5′FR(−1102)

Intron 5(20)

0.8

0

5′FR(−1102)

Exon 5(424)

0.95

0

5′FR(−1102)

Intron 5(74)

0.95

0

Exon 1(−625)

Intron 1(477)

0.95

0

Intron 1(402)

Intron 1(−1131)

1

0

Intron 1(477)

3′FR(1038)

1

1E−04

Exon 2(−183)

Exon 5(424)

1

0

Exon 2(−183)

Intron 5(74)

1

0

Intron 4(20)

Exon 5(424)

1

0

Intron 4(20)

Intron 5(74)

1

0

Exon 5(424)

Intron 5(74)

1

0

Intron 6(−55)

Intron 6(−21)

1

0

Exon 7(591)

3′FR(1038)

0.93

0

TABLE 5

GSTO2 Haplotypes for Caucasian-American Population

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

Freq.

o/i*

(−2012)

(−1785)

(−1771)

(−1664)

(−1598)

(−1576)

(−1566)

(−1541)

(−1357)

(−1189)

(−1179)

(−1105)

(−1102)

50.06%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

21.60%

o

D

T

G

C

C

G

G

C

C

T

G

A

T

9.92%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

3.14%

i

D

T

G

C

C

G

G

C

C

T

G

A

T

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

G

G

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

0.83%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

5′FR

5′FR

5′FR

5′FR

5′FR

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Freq.

(−1009)

(−781)

(−730)

(−637)

(−634)

(−625)

(−530)

(−526)

(−463)

(−450)

(−424)

(−410)

(−293)

(−234)

50.06%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

21.60%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

9.92%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

3.14%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

0.83%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Freq.

(+52)

(+152)

(+310)

(+311)

(+333)

(+402)

(+426)

(+435)

(+477)

(+1694)

(+2003)

(+2146)

(+2168)

(−1131)

50.06%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

21.60%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

9.92%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

3.14%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

C

T

T

I

A

A

A

G

C

0.83%

A

G

I

G

T

T

T

T

D

A

A

A

G

T

Intron1

Intron1

Intron1

Exon2

Intron2

Exon3

Intron4

Intron4

Intron4

Exon5

Exon5

Intron5

Intron5

Exon6

Freq.

(−185)

(−19)

(−12)

(−183)

(−16)

(+121)

(+20)

(+39)

(−38)

(+389)

(+424)

(+74)

(+170)

(+472)

50.06%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

21.60%

A

T

G

G

C

G

C

C

A

G

G

G

T

C

9.92%

A

T

G

A

C

G

C

C

A

G

A

A

T

C

3.14%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

1.67%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

1.67%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

0.83%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

Intron6

Intron6

Exon7

Exon7

3′ FR

3′FR

3′ FR

3′FR

3′ FR

3′FR

3′FR

Freq.

(−55)

(−21)

(+591)

(+630)

(+1038)

(+1080)

(+1243)

(+1343)

(+1615)

(+1649)

(+1724)

50.06%

C

C

G

C

C

C

A

T

C

T

C

21.60%

C

C

G

C

C

C

A

T

C

T

C

9.92%

C

C

G

C

C

C

A

T

C

T

C

3.14%

C

C

G

C

C

C

A

T

C

T

C

1.67%

C

C

G

C

C

C

A

T

C

T

C

1.67%

C

C

T

C

T

C

G

T

C

T

C

0.83%

C

C

G

C

C

C

A

T

C

T

C

*o = observed; i = inferred

TABLE 6

GSTO2 Haplotypes for African-American Population

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

Freq.

o/i*

(−2012)

(−1785)

(−1771)

(−1664)

(−1598)

(−1576)

(−1566)

(−1541)

(−1357)

(−1189)

(−1179)

(−1105)

(−1102)

16.70%

i

I

T

G

T

G

C

A

C

A

C

G

G

T

14.30%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

14.20%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

6.53%

i

D

T

G

C

C

G

G

C

C

T

G

A

T

5.00%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

4.17%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

3.33%

i

D

C

G

C

C

G

G

C

C

T

G

G

T

3.33%

i

D

T

G

C

C

G

G

C

C

T

G

G

G

3.33%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

2.50%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

2.50%

i

I

T

G

T

G

C

A

C

A

C

G

G

T

1.67%

i

D

C

G

C

C

G

G

C

C

T

G

G

T

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

1.53%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

5′FR

5′FR

5′FR

5′FR

5′FR

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Freq.

(−1009)

(−781)

(−730)

(−637)

(−634)

(−625)

(−530)

(−526)

(−463)

(−450)

(−424)

(−410)

(−293)

(−234)

16.70%

T

C

I

T

G

A

G

C

A

C

G

T

A

C

14.30%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

14.20%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

6.53%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

5.00%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

4.17%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

3.33%

C

A

D

C

C

A

C

C

G

A

G

T

T

C

3.33%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

3.33%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

2.50%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

2.50%

T

C

I

T

G

A

G

C

A

C

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

T

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.53%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Freq.

(+52)

(+152)

(+310)

(+311)

(+333)

(+402)

(+426)

(+435)

(+477)

(+1694)

(+2003)

(+2146)

(+2168)

(−1131)

16.70%

A

G

I

A

T

C

T

T

D

A

G

A

G

C

14.30%

A

G

I

G

T

C

T

T

D

A

A

A

A

C

14.20%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

6.53%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

5.00%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

4.17%

A

G

I

G

T

C

T

T

D

G

A

A

G

C

3.33%

A

G

I

G

T

C

C

T

D

A

A

A

G

C

3.33%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

3.33%

A

G

I

G

T

T

T

T

D

A

A

A

G

T

2.50%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

2.50%

A

G

I

A

T

C

T

T

D

A

G

A

G

C

1.67%

A

G

I

G

T

C

C

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

C

T

A

D

A

A

A

A

C

1.53%

A

G

I

G

T

C

T

T

D

A

A

A

A

C

Intron1

Intron1

Intron1

Exon2

Intron2

Exon3

Intron4

Intron4

Intron4

Exon5

Exon5

Intron5

Intron5

Exon6

Freq.

(−185)

(−19)

(−12)

(−183)

(−16)

(+121)

(+20)

(+39)

(−38)

(+389)

(+424)

(+74)

(+170)

(+472)

16.70%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

14.30%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

14.20%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

6.53%

A

T

G

G

C

G

C

C

A

G

G

G

T

C

5.00%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

4.17%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

3.33%

G

T

G

A

C

G

C

C

A

G

G

G

T

C

3.33%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

3.33%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

2.50%

A

T

G

A

C

G

C

C

A

G

A

A

T

C

2.50%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

1.67%

G

T

G

A

C

G

C

C

A

G

G

G

T

C

1.67%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

1.53%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

Intron6

Intron6

Exon7

Exon7

3′ FR

3′FR

3′ FR

3′FR

3′ FR

3′FR

3′FR

Freq.

(−55)

(−21)

(+591)

(+630)

(+1038)

(+1080)

(+1243)

(+1343)

(+1615)

(+1649)

(+1724)

16.70%

C

C

G

C

C

C

A

T

C

T

C

14.30%

C

C

G

C

C

C

A

T

C

T

C

14.20%

C

C

G

C

C

C

A

T

C

T

C

6.53%

C

C

G

C

C

C

A

T

C

T

C

5.00%

C

C

G

C

C

C

A

T

C

T

C

4.17%

C

C

G

T

C

C

A

T

C

T

C

3.33%

C

C

G

T

C

C

A

T

C

T

C

3.33%

C

C

T

C

T

C

A

T

C

T

C

3.33%

C

C

G

C

C

C

A

T

C

T

C

2.50%

C

C

G

C

C

C

A

T

C

T

C

2.50%

C

T

G

C

C

C

A

G

C

T

C

1.67%

C

C

T

T

T

C

A

T

C

T

C

1.67%

C

C

T

C

T

C

A

T

C

T

C

1.53%

C

C

T

C

T

C

A

G

C

T

C

*o = observed; i = inferred

TABLE 7

GSTO2 Haplotypes for Han Chinese-American Population

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

Freq.

o/i*

(−2012)

(−1785)

(−1771)

(−1664)

(−1598)

(−1576)

(−1566)

(−1541)

(−1357)

(−1189)

(−1179)

(−1105)

(−1102)

55.20%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

13.30%

o

D

T

G

C

C

G

G

C

C

T

G

A

T

12.30%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

6.43%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

A

T

1.67%

o

D

T

G

C

C

G

G

C

C

T

G

A

T

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

1.07%

i

D

T

G

C

C

G

G

C

C

T

G

G

G

5′FR

5′FR

5′FR

5′FR

5′FR

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Freq.

(−1009)

(−781)

(−730)

(−637)

(−634)

(−625)

(−530)

(−526)

(−463)

(−450)

(−424)

(−410)

(−293)

(−234)

55.20%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

13.30%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

12.30%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

6.43%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.07%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Freq.

(+52)

(+152)

(+310)

(+311)

(+333)

(+402)

(+426)

(+435)

(+477)

(+1694)

(+2003)

(+2146)

(+2168)

(−1131)

55.20%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

13.30%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

12.30%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

6.43%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

C

T

T

D

A

A

A

A

C

1.67%

A

G

I

G

T

C

T

T

I

A

A

A

G

C

1.07%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

Intron1

Intron1

Intron1

Exon2

Intron2

Exon3

Intron4

Intron4

Intron4

Exon5

Exon5

Intron5

Intron5

Exon6

Freq.

(−185)

(−19)

(−12)

(−183)

(−16)

(+121)

(+20)

(+39)

(−38)

(+389)

(+424)

(+74)

(+170)

(+472)

55.20%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

13.30%

A

T

G

G

C

G

C

C

A

G

G

G

T

C

12.30%

A

T

G

A

C

G

C

C

A

G

A

A

T

C

6.43%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

1.67%

A

T

G

G

C

G

C

C

A

A

G

G

T

C

1.67%

A

T

G

G

C

G

T

C

A

G

A

A

T

C

1.67%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

1.67%

A

T

G

A

C

G

C

C

A

G

G

G

T

C

1.07%

A

T

G

A

C

G

C

C

A

G

A

A

T

C

Intron6

Intron6

Exon7

Exon7

3′ FR

3′FR

3′ FR

3′FR

3′ FR

3′FR

3′FR

Freq.

(−55)

(−21)

(+591)

(+630)

(+1038)

(+1080)

(+1243)

(+1343)

(+1615)

(+1649)

(+1724)

55.20%

C

C

G

C

C

C

A

T

C

T

C

13.30%

C

C

G

C

C

C

A

T

C

T

C

12.30%

C

C

G

C

C

C

A

T

C

T

C

6.43%

C

C

T

C

T

C

A

T

C

T

C

1.67%

C

C

G

C

C

C

A

T

C

T

C

1.67%

C

C

G

C

C

C

A

T

C

T

C

1.67%

C

C

T

C

T

C

A

T

C

T

C

1.67%

C

C

T

C

T

C

A

T

C

T

C

1.07%

C

C

T

C

T

C

A

T

C

T

C

*o = observed; i = inferred

TABLE 8

GSTO2 Haplotypes for Mexican-American Population

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

Freq.

o/i*

(−2012)

(−1785)

(−1771)

(−1664)

(−1598)

(−1576)

(−1566)

(−1541)

(−1357)

(−1189)

(−1179)

(−1105)

(−1102)

60.70%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

16.70%

o

D

T

G

C

C

G

G

C

C

T

G

A

T

6.84%

o

D

T

G

C

C

G

G

C

C

T

G

G

G

3.51%

i

D

T

G

C

C

G

G

C

C

T

G

A

T

1.67%

i

D

T

G

C

C

G

G

C

C

T

G

G

T

1.49%

i

D

T

G

C

C

G

G

C

C

T

G

A

T

5′FR

5′FR

5′FR

5′FR

5′FR

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Exon1

Freq.

(−1009)

(−781)

(−730)

(−637)

(−634)

(−625)

(−530)

(−526)

(−463)

(−450)

(−424)

(−410)

(−293)

(−234)

60.70%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

16.70%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

6.84%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

3.51%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.67%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

1.49%

C

A

D

C

C

A

C

C

G

A

G

T

A

C

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

Freq.

(+52)

(+152)

(+310)

(+311)

(+333)

(+402)

(+426)

(+435)

(+477)

(+1694)

(+2003)

(+2146)

(+2168)

(−1131)

60.70%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

16.70%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

6.84%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

3.51%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

1.67%

A

G

I

G

T

T

T

T

D

A

A

A

G

T

1.49%

A

G

I

G

T

C

T

T

D

A

A

A

G

C

Intron1

Intron1

Intron1

Exon2

Intron2

Exon3

Intron4

Intron4

Intron4

Exon5

Exon5

Intron5

Intron5

Exon6

Freq.

(−185)

(−19)

(−12)

(−183)

(−16)

(+121)

(+20)

(+39)

(−38)

(+389)

(+424)

(+74)

(+170)

(+472)

60.70%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

16.70%

A

T

G

G

C

G

C

C

A

G

G

G

T

C

6.84%

A

T

G

A

C

G

C

C

A

G

A

A

T

C

3.51%

A

T

G

A

C

G

T

C

A

G

A

A

T

C

1.67%

A

T

G

A

C

G

C

C

C

G

G

G

T

C

1.49%

A

T

G

A

C

G

C

C

A

G

A

A

T

C

Intron6

Intron6

Exon7

Exon7

3′ FR

3′FR

3′ FR

3′FR

3′ FR

3′FR

3′FR

Freq.

(−55)

(−21)

(+591)

(+630)

(+1038)

(+1080)

(+1243)

(+1343)

(+1615)

(+1649)

(+1724)

60.70%

C

C

G

C

C

C

A

T

C

T

C

16.70%

C

C

G

C

C

C

A

T

C

T

C

6.84%

C

C

G

C

C

C

A

T

C

T

C

3.51%

C

C

G

C

C

C

A

T

C

T

C

1.67%

C

C

G

C

C

C

A

T

C

T

C

1.49%

C

C

G

C

C

C

A

T

C

T

C

*o = observed; i = inferred

TABLE 9

Summary of GSTO2 haplotypes

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

5′FR

AA

CA

HCA

MA

(−2012)

(−1785)

(−1771)

(−1664)

(−1598)

(−1576)

(−1566)

(−1357)

(−1189)

(−1105)

(−1102)

1

0.17

—

—

—

I

T

G

T

G

C

A

A

C

G

T

2

0.14

—

—

0.01

D

T

G

C

C

G

G

C

T

G

T

3

0.07

0.22

0.13

0.16

D

T

G

C

C

G

G

C

T

A

T

4

0.05

—

—

—

D

T

G

C

C

G

G

C

T

G

T

5

0.04

—

—

0.01

D

T

G

C

C

G

G

C

T

G

T

6

0.03

—

—

—

D

C

G

C

C

G

G

C

T

G

T

7

0.03

0.01

0.01

0.01

D

T

G

C

C

G

G

C

T

G

T

8

0.03

0.01

0.01

0.01

D

T

G

C

C

G

G

C

T

G

T

9

—

0.03

—

0.01

D

T

G

C

C

G

G

C

T

A

T

10

0.03

—

—

—

I

T

G

T

G

C

A

A

C

G

T

11

0.02

—

—

—

D

C

G

C

C

G

G

C

T

G

T

12

0.02

—

—

—

D

T

G

C

C

G

G

C

T

G

T

13

—

—

0.02

0.01

D

T

G

C

C

G

G

C

T

G

T

14

0.02

—

—

—

D

T

G

C

C

G

G

C

T

G

T

15

0.01

—

—

—

D

T

C

C

C

G

G

C

T

A

T

16

—

—

0.02

—

D

T

G

C

C

G

G

C

T

A

T

17

0.14

0.5

0.55

0.61

D

T

G

C

C

G

G

C

T

G

G

18

0.03

—

0.06

—

D

T

G

C

C

G

G

C

T

G

G

19

0.03

0.1

0.12

0.07

D

T

G

C

C

G

G

C

T

G

G

20

—

0.01

—

—

D

T

G

C

C

G

G

C

T

G

G

5′FR

5′FR

5′FR

5′FR

5′FR

Exon1

Exon1

Exon1

Exon1

Intron1

Intron1

Intron1

Intron1

Intron1

Intron1

(−1009)

(−781)

(−730)

(−637)

(−634)

(−530)

(−463)

(−450)

(−293)

(+311)

(+402)

(+426)

(+435)

(+1694)

(+2003)

1

T

C

I

T

G

G

A

C

A

A

C

T

T

A

G

2

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

3

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

4

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

5

C

A

D

C

C

C

G

A

A

G

C

T

T

G

A

6

C

A

D

C

C

C

G

A

T

G

C

C

T

A

A

7

C

A

D

C

C

C

G

A

A

G

T

T

T

A

A

8

C

A

D

C

C

C

G

A

A

G

T

T

T

A

A

9

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

10

T

C

I

T

G

G

A

C

A

A

C

T

T

A

G

11

C

A

D

C

C

C

G

A

T

G

C

C

T

A

A

12

C

A

D

C

C

C

G

A

A

G

C

T

A

A

A

13

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

14

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

15

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

16

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

17

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

18

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

19

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

20

C

A

D

C

C

C

G

A

A

G

C

T

T

A

A

Intron1

Intron1

Intron1

Exon2

Exon3

Intron4

Exon5

Exon5

Intron5

Exon6

Intron6

Exon7

Exon7

3′ FR

3′FR

(+2168)

(−1131)

(−185)

(−183)

(+121)

(+20)

(+389)

(+424)

(+74)

(+472)

(−21)

(+591)

(+630)

(+1038)

(+1343)

1

G

C

A

A

G

C

G

G

G

C

C

G

C

C

T

2

A

C

A

A

G

C

G

G

G

C

C

G

C

C

T

3

G

C

A

G

G

C

G

G

G

C

C

G

C

C

T

4

G

C

A

A

G

C

G

G

G

C

C

G

C

C

T

5

G

C

A

A

G

C

G

G

G

C

C

G

T

C

T

6

G

C

G

A

G

C

G

G

G

C

C

G

T

C

T

7

G

T

A

A

G

C

G

G

G

C

C

G

C

C

T

8

G

T

A

A

G

C

G

G

G

C

C

G

C

C

T

9

G

C

A

A

G

T

G

A

A

C

C

G

C

C

T

10

G

C

A

A

G

C

G

G

G

C

T

G

C

C

G

11

G

C

G

A

G

C

G

G

G

C

C

T

T

T

T

12

A

C

A

A

G

C

G

G

G

C

C

T

C

T

T

13

A

C

A

A

G

C

G

G

G

C

C

T

C

T

T

14

A

C

A

A

G

C

G

G

G

C

C

T

C

T

G

15

G

C

A

G

A

C

G

G

G

C

C

T

C

T

T

16

G

C

A

G

G

C

A

G

G

C

C

G

C

C

T

17

G

C

A

A

G

T

G

A

A

C

C

G

C

C

T

18

G

C

A

A

G

T

G

A

A

C

C

T

C

T

T

19

G

C

A

A

G

C

G

A

A

C

C

G

C

C

T

20

G

C

A

A

G

C

G

A

A

A

C

G

C

C

T

Example 4Immunoreactivity of GSTO2 Allozymes

As an initial step in the study of the possible functional implications of nonsynonymous cSNPs in the GSTO2 gene, expression constructs were created for all variant allozymes, and those constructs were used to transfect COS-1 cells. This mammalian cell line was transfected to ensure the presence of mammalian post-translational modification and the mammalian machinery for protein degradation. The GSTO2 expression constructs were cotransfected with β-galactosidase to make it possible to correct for variations in transfection efficiency. The antibodies for GSTO2 were directed against peptides that did not include any of the genetically variant amino acids. When GSTO2 constructs were transiently expressed in COS-1 cells, levels of the Cys130Tyr13 and Leu158Ile allozymes, as well as the double variant Asn142Asp/Leu 158Ile allozyme, were strikingly reduced. The Val41Ile and Asn142Asp allozymes were expressed at approximately 80% of the level of the WT allozyme (Table 10). The GSTO2 antibody cross-reacted with a protein in the COS-1 cells with a molecular mass slightly higher than that for recombinant GSTO2.

Attempts were made to complement expression studies with parallel assays of level of enzyme activity, but the assays currently available for GSTO2 activity-assays that have been used to study purified, bacterially expressed protein (Schmuck et al. (2005) Pharmacogenet Genomics 15:493-501)—lacked the sensitivity necessary for application to the cultured mammalian cell cytosol preparations.

TABLE 10

GSTO2 immunoreactive protein levels

WT

100 ± 2.3

Val41Ile

82.8 ± 1.6

Cys130Tyr

49.4 ± 2.8

Asn142Asp

76.1 ± 3.0

Leu158Ile

17.3 ± 1.7

Asn142Asp/Leu158Ile

15.2 ± 2.5

Example 5GSTO2 Degradation Studies

To investigate the mechanisms responsible for the striking decrease in expression of the GSTO2 Cys130Tyr and Leu158Ile allozymes, protein degradation studies were performed using RRL. Radioactive WT and variant GSTO2 allozymes were synthesized using “treated” RRL, and the recombinant proteins were incubated for various periods of time with “untreated” RRL and an ATP generating system to follow the time course of protein degradation. The CYS130Tyr, Leu158Ile, and Asn142Asp/Leu158Ile allozymes were degraded much more rapidly than WT GSTO2, with Leu158Ile and Asn142Asp/Leu158Ile being degraded more rapidly than Cys130Tyr (FIG. 3). After 24 hours, 26% of Cys130Tyr, 6.6% of Leu158Ile, 5.6% of Asn142Asp/Leu158Ile, and 73% of the WT protein remained. Therefore, accelerated degradation appeared to be responsible, at least in part, for differences in levels of the GSTO2 variant allozymes after expression in COS-1 cells.

Example 6Activity and Substrate Kinetics of GSTO2 Allozymes

Cell homogenate preparations containing recombinant GSTO2 allozymes, prepared as described in Example 1, are used to assess catalytic activity. The resulting activities are adjusted to a percentage of the WT GSTO2 enzyme activity. To determine whether alterations in activity of the GSTO2 variants might be related to quantity of immunoreactive enzyme protein, quantitative Western blot analyses are conducted.

Alterations in amino acid sequence can alter enzyme substrate affinity and/or catalytic efficiency. Substrate kinetic studies are conducted to determine whether the allozymes differ from the WT GSTO2 protein in these aspects. A series of DHA concentrations are used to estimate apparent Km values for recombinant wild type GSTO2 and for the variant allozymes.

Example 7mRNA Levels of GSTO2 Allozymes

RT-PCR was used to determine whether differences in protein levels were due to differences in mRNA levels. To evaluate the mRNA level after transient expression of GSTO2 WT and one or more GSTO allozymes, RT-PCR is performed using mRNA isolated from transfected COS-1 cells. β-Galactosidase is used to correct for transfection efficiency and also as an internal control for RT-PCR.

Example 8Reporter Assays

Luciferase reporter assays are conducted to determine whether polymorphisms in the 5′-FR affect transcription of GSTO2 allozymes. Reporter gene constructs including commonly occurring haplotypes (e.g., haplotypes with an allele frequency greater than 3%) are created within the first 1 kb upstream of the translation initiation codon. These constructs are transfected into primary cells or cell lines such as liver, kidney, skeletal muscle, heart, lung, cervix, ovary, or prostate cells. After incubation for 48 hours, cells are harvested and lysed, and luciferase activity is measured using a dual-luciferase assay system.

Other Embodiments

It is to be understood that while the invention has been described in conjunction with the detailed description thereof, the foregoing description is intended to illustrate and not limit the scope of the invention, which is defined by the scope of the appended claims. Other aspects, advantages, and modifications are within the scope of the following claims.

Claims (2)

1. An isolated glutathione S-transferase omega 2 (GSTO2) polypeptide, wherein said polypeptide comprises an amino acid sequence which differs from SEQ ID NO:9 solely by a modification selected from the group consisting of (a) an amino acid substitution at position 41 of SEQ ID NO:9, (b) an amino acid substitution at position 130 of SEQ ID NO:9, (c) an amino acid substitution at position 158 of SEQ ID NO:9, and (d) any combination of (a)-(c).

2. The isolated polypeptide of claim 1, wherein (a) the amino acid at position 41 of SEQ ID NO:9 has been replaced with isoleucine, (b) the amino acid at position 130 of SEQ ID NO:9 has been replaced with tyrosine, (c) the amino acid at position 158 of SEQ ID NO:9 has been replaced with isoleucine, or (d) any combination of (a)-(c).