Abstract:
A method is disclosed for detecting galactosemia-causing mutations in the GALT gene, comprising amplifying a portion of the GALT gene from isolated DNA and allowing a pair of labeled probes to hybridize to the portion. One of the labeled probes is adapted to match to a sequence that includes the galactosemia-causing mutation, and another of the labeled probes hybridizes to an adjacent sequence, thereby forming a hybrid. Melting curves of each hybrid are then analyzed, wherein peaks of the curves are produced at an acquired fluorescence and melting temperature, T m ; and a genotype is assigned based on the T m  of the hybrid. Resulting melting peaks are compared to reference sample peaks derived from samples characterized to contain the mutations, wherein the reference sample curves indicate a temperature change, ΔT m , between mutant and wild type peaks.

Description:
SPECIFIC REFERENCE  
       [0001]    The present application claims benefit of priority date so established by provisional application serial No. 60/335,630, filed Oct. 31, 2001. 
     
    
     
       BACKGROUND  
         [0002]    1. Field of the invention  
           [0003]    The present invention relates generally to methods for detecting mutations in the galactose-1-phosphate uridyl transferase (GALT) gene. In particular, a set of methods are disclosed for genotyping specimens to determine their status (wild type or mutant) at five loci, each having a particular melting peak as they are subjected to melting curve analysis.  
           [0004]    2. Description of the Related Art  
           [0005]    Galactosemia is a standard component of most newborn metabolic screening programs. The classical form of galactosemia is caused by mutations in the galactose-1-phosphate uridyl-transferase (GALT) gene. Screening for galactosemia is achieved through analysis of total galactose (galactose and galactose-1-phosphate) and determining the activity of the GALT enzyme. This approach is effective, but environmental factors, specimen processing procedures, and the high frequency of the Duarte variant (N314D) necessitates further analysis to reduce false positive results.  
           [0006]    Prospective screening of newborns for galactosemia is a routine procedure in the United States and many foreign countries. Screening utilizes the universally collected Guthrie dried blood card specimen (DBS) to assay for total galactose (galactose plus glactose-1-phosphate), and the activity of the galactose-1-phosphate uridyl transferase (GALT) enzyme. Total galactose is typically assayed for using NAD reduction analysis while GALT activity is determined using the Beutler assay.  
           [0007]    Classical galactosemia results from mutations to the GALT gene, which cause severe perturbation in the activity of the corresponding enzyme. The GALT gene has been characterized and numerous mutations have been identified. Biochemical analysis for galactose and GALT activity are sound principle methods to prospectively assay newborns for galactosemia, however complicating factors can interfere with these results, necessitating further analysis.  
           [0008]    It is commonly observed in newborn screening laboratories that environmental factors and sample collection/handling procedures, practiced at the site of specimen collection, may have severe adverse effects upon GALT activity, thereby causing abnormally low results in the Beutler assay. The most notable environmental influences upon GALT activity are heat and humidity. Specimens collected during hot, humid summer months, or in climates where such conditions are persistent, often present with reduced GALT activity in the Beutler assay. The practice of batching, where dried blood spot (DBS) specimens are permitted to accumulate before being mailed to the screening lab, also adversely affects GALT activity. Enzyme activity deteriorates over time, which minimizes the optimum time period between specimen collection and analysis. In addition steps must be taken to avoid the specimen&#39;s exposure to heat and humidity, which are the best practices to ensure optimal performance in the Beutler assay.  
           [0009]    Additionally, reasons inherent to the screening process itself and the nature of mutations in the GALT gene necessitate analysis beyond the biochemical regimen. A false negative result may cause, at the minimum, serious medical consequences and, in a worse case scenario, death. Comparatively, a false positive result may lead to parental anxiety, mistrust of the screening process, and unnecessary medical procedures. To balance these concerns, screening labs require that critical result values fall within a range where false negative results are eliminated and false positive results are minimized.  
           [0010]    The Duarte variant (N314D), carried by upwards of 5% of the general population, causes a partial loss (˜25%) of GALT activity. The high frequency of the Duarte variant is yet another complicating factor-in galactosemia screening.  
           [0011]    Supplementing biochemical data using mutational analysis is a powerful method to reduce false positive results and may be used to provide unambiguous confirmation of true positive results. For example, a 2-tiered approach, biochemical analysis followed by gene-level analysis, is the standard employed by most newborn screening laboratories in their cystic fibrosis screening programs. Primary screening for cystic fibrosis is performed by analysis of circulating trypsinogen and those specimens having elevated trypsinogen are subsequently assayed for the CFTR Δ508 mutation. Delta F508 accounts for approximatelv 70% of CFTR mutations worldwide. The addition of Δ508 analysis has dramatically reduced the false positive rate in cystic fibrosis screening. A similar approach, as described herein, increases specificity in galactosemia screening.  
           [0012]    Barriers to gene-level analysis in the screening lab include complexity and turnover time. Traditional methods for genetic analysis such as DNA sequencing, allele specific cleavage, and allele specific oligonucleotide hybridization are time consuming and labor intensive, thus limiting their usefulness in a high throughput laboratory. A recently developed platform dubbed the Lightcycler® eliminates essentially all issues of complexity, turnover time, and labor intensity encountered with classical methods of mutation analysis. The Lightcycler® utilizes rapid air driven thermal cycling and in-line fluorescence analysis of hybridization probes to generate melting curves, which are subsequently used to generate melting peaks for genotype assignment.  
           [0013]    DNA melts at a defined temperature (T m ), wherein T m  is defined as the temperature at which half of the double helical structure is lost. This is generally the principle as known in the art supporting the Lightcycler®, since the melting temperature of a DNA molecule is dependent upon its nucleotide composition. DNA molecules rich in GC base pairs have a higher Tm than those having an abundance of AT base pairs. The Lightcycler® provides an innovative solution for identifying the base composition of a PCR amplified product, and in the present embodiment, allows for a qualitative method to assay for galactosemia-causing mutations.  
         SUMMARY  
         [0014]    Light Cycler technology and fluorescent-labeled hybridization probes are employed and a 5-mutation panel is described which includes the 4 most frequently encountered classical galactosemia alleles (Q188R, S135L, K285N, L195P) and the Duarte N314D variant. Five assays are performed simultaneously under a common set of conditions for both thermal cycling and melting curve analysis. Including DNA preparation, set-up, amplification, and analysis, the entire process requires less than 2 hours. These assays are useful to reduce false positive results, confirm classical galactosemia, and differentiate classical galactosemia from Duarte/Galactosemia (D/G) compound heterozygotes. These assays, owing to their speed and efficacy, are ideal for utilization in a high-throughput newborn screening laboratory.  
           [0015]    Using the presented methodologies for melting curve analysis of galactosemia-causing mutations, genotype data is generated in minutes, as opposed to hours or even days required when using traditional methods. A second advantage to this platform is that it is a “closed tube” assay where both amplification and analysis are performed in a common reaction vessel. Unified amplification and analysis allows for simplified sample tracking and greatly reduces the likelihood of amplicon contamination in the laboratory.  
           [0016]    Accordingly, what is provided is a method for detecting specific galactosemia-causing mutations in the GALT gene, comprising amplifying a portion of the GALT gene from isolated DNA, wherein said portion potentially contains the galactosemia-causing mutation, and thereby forming an amplification product; allowing a pair of labeled probes to hybridize to one strand of the amplification product, wherein one of the probes spans the allele of interest and can match to either the mutant or wild type allele and another of the labeled probes hybridizes to an adjacent sequence, thereby forming hybrids. The hybridized probes bring the donor fluorophore and acceptor fluorophore into close proximity allowing a fluorescent signal to be generated when the appropriate wavelength of light is provided. This fluorescent signal is used to generate melting curves. Genotype is assigned based on the T m  of the disassociated hybrid that forms the peak. The galactosemia-causing mutations include the four most frequently encountered classical galactosemia alleles (Q188R, S135L, K285N, L195P) and the Duarte N314D variant. Thus, the method further comprises the step of comparing the resulting melting curves to reference sample curves of samples characterized to contain the above mutations, wherein the reference sample curves indicate a temperature change, ΔT m , between mutant and wild type peaks. 
       
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0017]    [0017]FIG. 1 shows the portion of the sequence for each allele with the SNP, and the sequences of the detection, anchor probes used for hybridization, and the location of fluorescent moieties.  
         [0018]    [0018]FIG. 2 shows the results of the melting peaks for each allele after simultaneous analysis. Melting peaks are well-separated, thereby facilitating unambiguous genotype assignment for each loci.  
         [0019]    FIGS.  3 A-E display the resulting melting peaks and acquired fluorescence for the N314D, Q188R, S135L, K285N, and L195P assays respectively.  
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT  
       [0020]    The mutations represented in the panel described herein represent the 4 most frequently encountered classical galactosemia mutations observed in the general United States population. Q188R is the most frequently encountered mutation representing approximately 70% of classical galactosemia alleles. The S135L mutation is most frequently observed among African Americans and is the second most frequently encountered allele. K285N is common in those of eastern European descent while L195P represents approximately 2% of classical galactosemia alleles. The Duarte variant, N314D, is present in approximately 5% of the US population and is probably the individually most complicating factor in screening for galactosemia. The so-called D/G compound heterozygotes (where a 314D allele is paired with a classical galactosemia mutation) may display both elevated total galactose and reduced GALT activity effectively mimicking classical galactosemia.  
         [0021]    As used herein, “GALT” refers to the enzyme galactose-1-phosphate uridyltransferase, and “galactosemia” is the deficiency in the activity of the GALT enzyme. Screening for galactosemia is thus achieved through analysis of total galactose (galactose and galactose-1-phosphate) and determining the activity of the GALT enzyme.  
         [0022]    For patient and/or reference specimen preparation, DNA samples are collected from any traditional methods, such as from any tissue or organ from which DNA can be amplified, or by purification from a dried blood spot on filter paper. “Reference specimens” as used herein, are previously characterized DNA specimens containing the mutations of interest. “Patient specimens” are the routinely collected DNA specimens that are to be screened and compared to results obtained from the reference specimens to determine a normal, mutant, or heterozygote sample.  
         [0023]    The sequence of the human GALT gene (Genbank accession number M96264) is the basis for primer and probe design. A portion containing intron 4-exon 10 of the human GALT gene having the potential mutations or complementary positions is set forth in SEQ ID NO: 1 and is nucleotide region (1337-3193) of the above human GALT gene sequence.  
         [0024]    A “primer” as used herein is a short piece of artificially made DNA complementary to a given DNA sequence and which acts as the initiation point from which replication proceeds via polymerase chain reaction (PCR). A “probe” is a fluorescent-labeled synthetic strand of DNA that anneals or “hybridizes” to a complementary DNA sequence generated by PCR. A “detection probe” hybridizes with a sequence that includes the site of the mutation and an “anchor probe” is another half of the probe set that hybridizes to an adjacent sequence. When both are hybridized it brings their respective fluorescent moieties into close proximity. A signal is generated by providing a specific wavelength of light, and fluorescence is monitored during incremental temperature increase to produce a melting curve. Melting curves are utilized to produce melting peaks. As will be further described, these resulting melting peaks are analyzed to determine wild type and mutant indications, wherein a temperature change, ΔT m , separates the different peaks. “A” as used in the claims may mean one or more depending on the context of the claim. It is not necessary to use only two probes as in the present assay. There are a few other types of probes that can produce these melting peaks without using two probes. For example, there are single probe systems with a single labeled probe that will produce melting peaks.  
         [0025]    Primers, probes, the concentration at which each is used, and the associated SEQ ID NO are listed in table 1. 
                                               TABLE 1                           Forward       Reverse Primer       Anchor       Detection           Allele   Primer   [ ]   Sequence   [ ]   Probe   [ ]   Probe   [ ]           Sequence   μM       μM   Sequence   μM   Sequence   μM                   N314D   5′   0.5   5′   0.3   5′   0.3   5′-LC Red 640   0.14           ACTGTAAAAG       GCAAGCATTTCGT       CGCAGGAGCGGAG       CTGCCAATGGT                   GGCTCTCTCT       AGCCAA 3′       GGTAGTAATGAGC       CCCAGTTGG-PO4 3′                       CC 3′               GTGCA-FITC 3′               SEQ ID   2       3       4       5       NO       Q188R   5′   0.5   5′   0.5   5′   0.2   5′-LC Red 640   0.13           CTTTTGGCTA       TTCCCATGTCCACA       GCCAAGAAACCCA       ACACCCTTACC                   ACAGAGCTCC       GTGCTGG 3′       CTGGAGCCCCT-FITC 3′       CGGCAGTG-PO4 3′                   G 3′               SEQ ID   6       7       8       9       NO       S135L   5′   0.25   5′   0.125   5′-LC Red 640   0.1   5′   0.2           CACAGCCAAG       ACCTCACAAACCT       GAAGCACATGACC       CGTTACATCCA                   CCCTACCTCT       GCACCCAA 3′       TTACTGGGTGGTG       ACCAGGGGT-FITC 3′                   C3′               ACGG-PO4 3′               SEQ ID   10       11       12       13       NO       K285N   5′   0.125   5′   0.25   5′-LC Red 640   0.2   5′   0.1           GCTGAGAGTC       CCAGAAATGGTGT       CTTTGAGACGTCCT       GCTCTTGACCA                   AGGCTCTGAT       TGGGGCT 3′       TTCCCTACTCCATG-PO4 3′       ATTATGACAAC-FITC 3′               SEQ ID   14       15       16       17       NO       L195P   5′   0.5   5′   0.5   5′-LC Red 640   0.2   5′   0.1           GAGGCTTGGA       TCCATTAGCAGGG       TGCCCAGCGTGAG       CAGCAGTTTCC                   GGTAAAGGAC       GCTCTCC 3′       GAGCGATCTCAGC       CGCCAGATA-FITC 3′                       3′               AG-PO4 3′               SEQ ID   18       19       21       22       NO                  
 
         [0026]    The number of nucleotides in the primers and probes may vary slightly. For a melting analysis assay, adequate thermal stability is needed such that the melting peak (T m ) is in a useful range, which may be defined generally as 50-70° C. As such, the T m  of the mismatched and matched probes should be within this range. The anchor probes serve to hold the second fluorophore in proximity to the first throughout the melting transition of the detection probe. As such, the anchor probe must have a higher Tm than the detection probe to which it is paired. Generally, at least 15% is required. In this preferred embodiment, 20-25% is allowed. So the length of the probes may vary depending upon the G/C content of the DNA to which it is hybridizing. If it is A/T rich, the probe is longer, and if it is G/C rich the probe may be shorter. Generally, this method starts with 30 nucleotides and nucleotides are either added or subtracted therefrom until a probe having desirable qualities is obtained. Such qualities include an adequately high melting temperature, no serious self hybridization or cross hybridization interactions, no serious self-interactions (folding) and no “false hybridization” within the amplified DNA fragment.  
         [0027]    All primer pairs are designed to function under a common set of thermal cycling conditions. In this embodiment for example, all primer pairs are designed with Tm values between 59-64° C.  
         [0028]    The organization of the anchor and detection probes for each allele is shown in FIG. 1. In this embodiment, one probe is labeled 3′ with FITC while the other probe is labeled 5′ with LC-red 640 and 3′ phosphorylated. Obviously, these labels can change if analysis takes place at a different wavelength. For example LC-red 705 could be used, which would change the interpretive guidelines as would be known in the art. Fluorescent-labeled probes are analyzed by spectrophotometry for oligonucleotide and fluorophore concentration. Probes with fluorophore/oligonucleotide ratios of 0.8-1.2 are generally suitable.  
         [0029]    The amplicons, or PCR products, are preferably held to fewer than 200 base pairs. Though not necessary, this facilitates the optimum binding of hybridization probes. Table 2 shows the number of base pairs in the amplicons for the current assay.  
                           TABLE 2                                   Allele   PCR Product Length                           N314D   171 bp           Q188R   160 bp           S135L   190 bp           K285N   155 bp           L195P   149 bp                      
 
         [0030]    Following amplification, the cycling protocol proceeds seamlessly to melting analysis. Fluorescence is acquired continuously during melting curve analysis and melting curves are constructed from data acquired during the upward temperature ramp.  
         [0031]    The following example presents the recorded preparatory procedure and results obtained for specimen preparation, hybridization, and probe analysis.  
       EXAMPLE  
       [0032]    Specimens and DNA Preparation. Reference specimens, previously characterized to contain the mutations of interest, were utilized for assay development. Additional specimens were collected during routine newborn screening for galactosemia. Specimens whose reducing capacity were below 60 μM reduced NAD or whose total galactose was above 20 mg/dl were selected for mutation analysis. DNA was isolated from DBS specimens as previously described and 80-130 ηg is utilized as template in each reaction.  
         [0033]    Amplification and Hybridization Probe Analysis. The sequence of the human GALT gene (Genbank accession number M96264) was used as a basis for primer and probe design. Primers and probes were designed in silico using Primer Premier 5.0 software. All primers and probes were HPLC purified and obtained from Operon technology (Alameda, Calif.). PCR reaction buffers, 20 mM MgCl 2  for the K285N and S135L assays, 30 mM MgCl 2  for the N314D, Q188R, and L195P assays, are obtained from Idaho Technology (Salt Lake City, Utah). All reactions use 0.6 U Klen taq DNA polymerase (AB Peptides, St. Louis, Mo.) complexed with TaqStart antibody (Clontech, Palo Alto, Calif.) according to manufacturers instructions. Primers, probes, and the concentration at which each is used are listed in table 1. The number of base pairs in the amplicons is listed in Table 2. All primer pairs were designed with Tm values between 59-64° C. and as such function under a common set of thermal cycling conditions. Amplification was performed in a Roche Light Cycler (Manheim, Germany). Cycling conditions were 40 cycles of 94° C., for 0 seconds (20°/second ramp speed)&gt;60° C. for 20 seconds (20°/second ramp speed), &gt;72° C., 0 seconds (2°/second ramp speed). Fluorescence was acquired at the end of the 20-second primer-annealing segment of the amplification. In the cases of N314D, S135L, and K285N assays, amplification was performed in an asymmetric fashion favoring the strand to which the hybridization probes bind (sense strand for S135L and N314D, antisense strand for K285N) while the Q188R and L195P assays were amplified in a symmetric manner. Amplicons were held to fewer than 200 base pairs (See Table 2), which facilitated optimum binding of hybridization probes. Following amplification, the cycling protocol proceeds seamlessly to melting analysis. Melting curve analysis used the following conditions: 97° C., 0 seconds, ramping at 2°/second down to 40° C., and ramping back up to 76° C. at 0.1° C/second. Fluorescence was acquired continuously during melting curve analysis and melting curves were constructed from data acquired during the upward ramp from 40° C. to 76° C.  
         [0034]    The assays utilized “probe:probe” format for Light Cycler genotyping assays. The probe:probe format utilizes 2 oligonucleotide probes that hybridize to a selected strand of the amplicon. A detection probe was employed which matches a sequence that includes the site of the mutation and an anchor probe, which hybridizes to an adjacent sequence. In these assays there is a 1-nucleotide gap between the anchor and detection probes. One probe was labeled 3′ with FITC while the other probe was labeled 5′ with LC-red 640 and 3′ phosphorylated. Probes were designed to maximize destabilization of the mismatch hybrid and it was determined in all 5 instances that matching the mutant allele with subsequent mismatch to the wild type allele provided the most effective probe design.  
         [0035]    Results  
         [0036]    [0036]FIG. 2 display all 5 assays analyzed simultaneously as would be initially observed following routine analysis. The N314D, Q188R, and S135L assays display analysis of specimens that are homozygous wild type, homozygous mutant, heterozygous, and a no-amplification control. Assays for K285N and L195P show the analysis of specimens that are homozygous wild type, heterozygous, and no amplification control. Combined analysis as shown in FIG. 2 is complex, so individual assays are subsequently viewed by selecting individual specimens or groups of specimens (e.g. control and test specimens) as seen in FIGS.  3 A-E.  
         [0037]    FIGS.  3 A-E display melting peaks for the N314D, Q188R, S135L, K285N, and L195P assays respectively. In all cases the peak representing the mutant allele is high temperature melting peak (perfect match with the detection probe) while the wild type allele is the low-temperature melting peak (mismatch hybrid with the detection probe). Melting peaks are well separated facilitating unambiguous genotype assignment for each loci. The melting temperature of each peak and the ΔT m  separating the wild type and mutant peaks for individual assays are displayed in Table 3.  
                                                         TABLE 3                                   Allele   Wild Type Tm   Mutant Tm   ΔTm                                        N314D   57.83   66.16   8.33           Q188R   56.80   65.47   8.67           S135L   55.08   62.09   7.01           K285N   54.89   61.41   6.52           L195P   49.94   60.57   10.63                      
 
         [0038]    Peak separation ranges from 6.52°-10.630° C., enabling easy and unambiguous genotype assignment.  
     
       
       
         1 
         
           
             21  
           
           
             1  
             1857  
             DNA  
             Homo sapiens  
           
            1 

gtaactatgg atttcccctc ttacaacttt caaaccagag ttggagactc agcattgggg     60 

ttcgccctgc ccgtagcaca gccaagccct acctctcggt tatcttttct cccgtcacca    120 

cccagtaagg tcatgtgctt ccacccctgg tcggatgtaa cgctgccact catgtcggtc    180 

cctgagatcc gggctgttgt tgatgcatgg gcctcagtca cagaggagct gggtgcccag    240 

tacccttggg tgcaggtttg tgaggtcgcc ccttcccctg gatgggcagg gagggggtga    300 

tgaagctttg gttctgggga gtaacatttc tgtttccaca gggtgtggtc aggagggagt    360 

tgacttggtg tcttttggct aacagagctc cgtatcccta tctgatagat ctttgaaaac    420 

aaaggtgcca tgatgggctg ttctaacccc cacccccact gccaggtaag ggtgtcaggg    480 

gctccagtgg gtttcttggc tgagtctgag ccagcactgt ggacatggga acaggattaa    540 

tggatgggac agaggaaata tgccaatgat gtggaggctt ggaggtaaag gacctgcctg    600 

ttcttctctg cttttgcccc ttgacaggta tgggccagca gtttcctgcc agatattgcc    660 

cagcgtgagg agcgatctca gcaggcctat aagagtcagc atggagagcc cctgctaatg    720 

gagtacagcc gccaggagct actcaggaag gtgggagaga gccaagccct gtgtccccaa    780 

ggagtcccta actttcttat cccatgagag aggtgtgtaa aggagaaagc tagaggtgaa    840 

ctagtagaga gagacttgct aggaggcctt agcaataatc cagtaatcta aaggaaagat    900 

gatggtgact tagactcggg tggttagtgg tagaggtggt gagaagacat cagatcctgg    960 

gcacattctt ttcttctgct tcccttgcct atttgctgac cacactccgg ctcctatgtc   1020 

accttgatga cttcctatcc attctgtctt cctaggaacg tctggtccta accagtgagc   1080 

actggttagt actggtcccc ttctgggcaa catggcccta ccagacactg ctgctgcccc   1140 

gtcggcatgt gcggcggcta cctgagctga cccctgctga gcgtgatggt cagtctccca   1200 

agtaggatcc tggggctagg cactggatgg aggttgctcc cagtagggtc agcatctgga   1260 

ccccaggctg agagtcaggc tctgattcca gatctagcct ccatcatgaa gaagctcttg   1320 

accaagtatg acaacctctt tgagacgtcc tttccctact ccatgggctg gcatggtgag   1380 

gcttttcaag tacctatatt tagccccaac accatttctg ggctcctggg ctcagcctag   1440 

tgaactgcaa cctcaaagga gcaagccttg aaacagttgc tgggggaagt ggccagagta   1500 

gagatgctgg gactgagggt ggagcagcaa acttggtgaa actacatctc caatgtgctt   1560 

tctaatctcc tgccagctct tctcaagcag gggatcctgg gagatgtagt tttcagatac   1620 

ctggttgggt ttgggagtag gtgctaacct ggataactgt aaaagggctc tctctcccca   1680 

ctgtctctct tctttctgtc aggggctccc acaggatcag aggctggggc caactggaac   1740 

cattggcagc tgcacgctca ttactaccct ccgctcctgc gctctgccac tgtccggaaa   1800 

ttcatggttg gctacgaaat gcttgctcag gctcagaggg acctcacccc tgagcag      1857 

 
           
             2  
             22  
             DNA  
             Homo sapiens  
           
            2 

actgtaaaag ggctctctct cc                                              22 

 
           
             3  
             19  
             DNA  
             Homo sapiens  
           
            3 

gcaagcattt cgtagccaa                                                  19 

 
           
             4  
             31  
             DNA  
             Homo sapiens  
           
            4 

cgcaggagcg gagggtagta atgagcgtgc a                                    31 

 
           
             5  
             20  
             DNA  
             Homo sapiens  
           
            5 

ctgccaatgg tcccagttgg                                                 20 

 
           
             6  
             21  
             DNA  
             Homo sapiens  
           
            6 

cttttggcta acagagctcc g                                               21 

 
           
             7  
             21  
             DNA  
             Homo sapiens  
           
            7 

ttcccatgtc cacagtgctg g                                               21 

 
           
             8  
             24  
             DNA  
             Homo sapiens  
           
            8 

gccaagaaac ccactggagc ccct                                            24 

 
           
             9  
             19  
             DNA  
             Homo sapiens  
           
            9 

acacccttac ccggcagtg                                                  19 

 
           
             10  
             21  
             DNA  
             Homo sapiens  
           
            10 

cacagccaag ccctacctct c                                               21 

 
           
             11  
             21  
             DNA  
             Homo sapiens  
           
            11 

acctcacaaa cctgcaccca a                                               21 

 
           
             12  
             30  
             DNA  
             Homo sapiens  
           
            12 

gaagcacatg accttactgg gtggtgacgg                                      30 

 
           
             13  
             20  
             DNA  
             Homo sapiens  
           
            13 

cgttacatcc aaccaggggt                                                 20 

 
           
             14  
             23  
             DNA  
             Homo sapiens  
           
            14 

gctgagagtc aggctctgat tcc                                             23 

 
           
             15  
             20  
             DNA  
             Homo sapiens  
           
            15 

ccagaaatgg tgttggggct                                                 20 

 
           
             16  
             28  
             DNA  
             Homo sapiens  
           
            16 

ctttgagacg tcctttccct actccatg                                        28 

 
           
             17  
             22  
             DNA  
             Homo sapiens  
           
            17 

gctcttgacc aattatgaca ac                                              22 

 
           
             18  
             20  
             DNA  
             Homo sapiens  
           
            18 

gaggcttgga ggtaaaggac                                                 20 

 
           
             19  
             20  
             DNA  
             Homo sapiens  
           
            19 

tccattagca ggggctctcc                                                 20 

 
           
             20  
             28  
             DNA  
             Homo sapiens  
           
            20 

tgcccagcgt gaggagcgat ctcagcag                                        28 

 
           
             21  
             20  
             DNA  
             Homo sapiens  
           
            21 

cagcagtttc ccgccagata                                                 20