The sequence listing that is contained in the file named “UFFL056WO_ST25.txt,” which is 54.1 kilobytes as measured in Microsoft Windows operating system and was created on Aug. 25, 2015, is filed electronically herewith and incorporated herein by reference.
The present invention relates to the fields of molecular biology and genetics. In specific embodiments, the invention relates to identification and characterization of polymorphisms in a nucleic acid sample.
A sequence listing contained in the file named “UFFL056USP1_ST25.txt” which is 55 kilobytes (measured in MS-Windows®) and created on Sep. 5, 2014, and comprises 304 nucleotide sequences, is filed electronically herewith and incorporated by reference in its entirety.
The ability to select individuals for breeding based on a favorable genotype at a polymorphic locus is an important tool in plant and animal breeding technology. In addition, the ability to screen polymorphic markers for a parentage assay that is reliable and effective in any species depends on the accurate genotyping of a number of polymorphic loci. However, the use of polymorphic markers in breeding programs is greatly complicated by polygenic inheritance and epistasis, which can necessitate genotyping of a number of distinct polymorphic loci in an individual to gain useful information regarding a particular trait. The process of genotyping numerous polymorphic loci simultaneously is laborious and costly, and existing methods of doing so are frequently inaccurate. In the absence of new methods for efficiently and reliably detecting the genotype of an individual at a plurality of polymorphic genomic loci, the use of large numbers of polymorphic markers for selection in breeding programs may not be practical.
In one aspect, the invention provides a method for genotyping, that may be implemented using next-generation sequencing, one or more target loci in a nucleic acid sample, comprising the steps of: a) providing a nucleic acid sample; b) adding a first set of primers to the sample to form a first amplification mixture, wherein the primers in the first set comprise a primer tail sequence and are capable of hybridizing to the target sequence within or adjacent to one or more of the target loci; c) performing a first amplification reaction on the first amplification mixture to produce a first library of amplicons, wherein the amplicons comprise the primer tail sequence; d) adding a second set of primers to the first library to form a second amplification mixture, wherein the primers in the second set are capable of hybridizing to the primer tail sequence; and e) performing a second amplification reaction on the second amplification mixture to generate a second library of amplicons; wherein for at least 90% of the target loci, the number of amplicons in the second library derived from each of the target loci deviates from the average number of amplicons for all target loci by less than one order of magnitude (+ or −10×).
In one embodiment of a method of the invention, for at least 90% of the target loci, the number of amplicons in the second library derived from each of the target loci deviates from the average number of amplicons for all target loci by less than 5× or less than 2.5×. In another embodiment, the number of target loci is greater than 10, greater than 100, or greater than 1,000. In a method of the invention, the first amplification reaction and the second amplification reaction can be carried out simultaneously or consecutively. In a specific embodiment, the first amplification reaction is carried out before the second amplification reaction. The method may comprise purifying the first library after the first amplification reaction and before the second amplification reaction. In certain further embodiment, the first and/or second amplification reaction comprises at least 2 cycles, at least 5 cycles, at least 10 cycles, at least 25 cycles, at least 50 cycles, between 5 and 50 cycles, between 1 and 15 cycles, between 2 and 10 cycles, or between 4 and 6 cycles. The primers in the first primer set are in one embodiment present in varying concentrations. The concentrations of the primers may be calculated according to a regression equation. In one embodiment, the target loci are polymorphic genomic loci within a population. In yet another embodiment, one or more primers used with the invention contain a unique index/barcode sequence to distinguish the sequencing of a sample or multiple samples in parallel.
In another aspect, a method of the invention further comprising the steps of: f) obtaining sequence data from the first or the second library; and h) determining the genotype at one or more of the target loci from the sequence data.
In yet another aspect, the invention provides a method for identifying a novel polymorphic genomic locus in a sample, comprising the steps of: a) providing two or more samples from individuals in a population; b) subjecting each of the samples to the method of claim 17; and c) aligning sequences corresponding to one or more target loci from two or more samples to identify target loci having sequence variation between individuals. Still further provided by the invention is a kit for use in genotyping one or more target loci in a nucleic acid sample, comprising: a) a first set of primers, wherein each primer in the first set comprises a primer tail sequence and is capable of hybridizing to a target sequence; and b) a second set of primers, wherein each primer in the second set is capable of hybridizing to the primer tail sequence.
Selection of individuals in a breeding population based on genetic markers linked to traits of interest is an essential tool in animal and plant breeding. Genetic markers can be used to accurately identify and track a desired trait within a breeding population, and allow for detection of a desired trait without the need to grow large populations to maturity in order to observe individual phenotypes. However, existing technology such as whole-genome sequencing or sequence-capture for the detection or discovery of polymorphic genetic markers is labor-intensive, expensive, and often inaccurate. In addition, the ability to screen polymorphic markers for a parentage assay that is reliable and effective in any species has been hampered by limitations on accurately genotyping multiple polymorphic loci.
While polymerase chain reaction (PCR)-based methods have been used for the detection of known genetic markers, these methods are typically not suitable for analysis of large numbers of polymorphic markers in a sample, such as with implementation of next generation sequencing methods. This is due in part to the difficulty in producing a consistent PCR amplification of a large number of target sequences in a single sample. This difficulty arises from variation in the efficiency by which distinct target sequences are amplified, due in part to differences in primer hybridization kinetics between target sequences. Small differences in primer annealing properties result in a biased amplification that is propagated exponentially with additional amplification cycles. Although amplification bias can be minimized by reducing the number of cycles in an amplification reaction, this results in insufficient amplified product for subsequent analysis of genetic markers. On the other hand, if a large number of cycles are carried out, amplification bias results in the detection of high numbers of amplicons corresponding to one or a few loci, while other loci are under-represented or not detected at all.
The present invention solves this problem by providing methods for amplifying a large number of loci from genomic DNA in an unbiased manner. In some embodiments, the methods of the invention comprise a two-step amplification reaction. In the first amplification step, locus specific primers comprising a primer tail are used for amplification for only few cycles to prevent the development of significant amplification bias. In the second amplification step, universal primers specific to the primer tail introduced in the first amplification step are used for further amplification in an unbiased manner. Using this approach, the final number of reads obtained by sequencing of the amplification products is consistent across loci.
Embodiments of the present invention therefore advantageously provide methods for unbiased amplification of sequences from multiple loci within a sample. For example, the number of loci to be detected in a sample may be 10 or more, 100 or more, or 1,000 or more loci, including, for example, from a lower range of about 5, 10, 25, 50, 75, 100, 150, 200, 250 or 500 or more to about 50, 75, 100, 150, 200, 300, 400, 500, 750, 1,000, or 1500 or more, including all combinations thereof. In further embodiments, the invention provides methods for amplifying multiple loci within a sample such that the final number of amplicons derived from each locus is balanced across loci. In particular embodiments, the invention provides methods for amplifying multiple loci within a sample such that for at least 90% of the loci tested, the final number of amplicons derived from any particular locus deviates from the average number of amplicons for all loci by one order of magnitude (i.e. + or −10×). For instance, if on average the number of reads obtained for all amplicons is 1% of the total number of reads, then an expected maximum of 10% and a minimum of 0.1% of reads would be detected for at least 90% of the other loci. In other embodiments, the method for amplifying multiple loci within a sample is such that for at least 80%, 85%, 90% or 95% of the loci tested, the final number of amplicons derived from any particular locus deviates from the average number of amplicons for all loci by less than about 7.5×, 5×, 2.5×, 1.5×, 1×, or less than 0.5×.
The unbiased two-step amplification methods provided by the present invention may comprise a first amplification step which is carried out using the lowest number of cycles required to effectively create a first library of amplicons comprising primer tails. In certain embodiments, the first amplification reaction is carried out using between 1-15 cycles, 2-10 cycles, between 3-7 cycles, or between 4-6 cycles of amplification. In certain embodiments, the primers used in the first amplification step may be present in varying concentrations according to the specific loci to which they correspond. In other embodiments, the primers used in the first amplification step may be calculated or adjusted according to the following equation:
Wherein reads stands for the number of reads sequenced for primer i and dilution stands for the concentration level used in the experiment for the same primer i. The regression coefficients β0 and β1 are calculated based on the data as described in Example 1, and for the test set in the example consisted of 0.34678 and 1.42626, respectively. These values can change for other test sets. The methods of the invention may further provide a second amplification step which amplifies a first library of amplicons by using primers directed to a primer tail which was added to the amplicons in the first amplification step. In certain embodiments, the second amplification reaction can comprise at least one 1 cycle, between 5 and 50 cycles, or between 10 and 25 cycles.
In specific embodiments, the first and the second amplification steps may be carried out simultaneously or consecutively. For example, the first amplification step may be carried out using a first set of primers at a concentration such that the amount of primer remaining after the first amplification step will be negligible. The second amplification step can be carried out consecutively or sequentially after the addition of a second primer set without the need for removing residual first primer. In other embodiments, the first and second amplification reactions may be carried out consecutively, and residual primer may be removed after the first purification step is complete. For example, a purification step may be used between the first amplification reaction and the second amplification reaction. A purification step may include any means known in the art for separating amplicons from a reaction mixture. In yet further embodiments, the first primer set and the second primer set are designed such that they hybridize with their specific target sequences under different conditions. The first and second amplification steps can then be carried out using differing temperature cycling protocols without the need for removal of residual primer between the steps.
The invention further provides methods for obtaining sequence data at one or more loci in a nucleic acid sample. For example, sequence data can be obtained from a first or second library produced by the first or the second amplification step by methods known in the art. The invention further contemplates determining the genotype at one or more target loci, for example one or more polymorphic genetic loci within a genomic DNA sample, from the sequence data. Further envisaged applications of the methods of the present invention include screening enriched microsatellite libraries, performing transcript profiling cDNA-AFLP (digital Northern), sequencing of complex genomes, EST library sequencing (on whole cDNA or cDNA-AFLP), microRNA discovery (sequencing of small insert libraries), bacterial artificial chromosome (BAC) contig sequencing, bulked segregant analysis approach AFLP/cDNA-AFLP, and detection of AFLP fragments, e.g. for marker-assisted selection (MAS) or marker-assisted back-crossing (MABC).
The invention further provides for genotyping a sample at one or more known polymorphic loci using the unbiased amplification methods provided herein. In other embodiments, the methods provide for identification of new polymorphic genomic loci within a population. In an exemplary embodiment, two or more samples are obtained from individuals in a population, and each of the samples is processed according to the methods of the present invention to provide sequence information for one or more target loci within the samples. Sequence data from the one or more individuals is then aligned to detect variations in sequences between individuals in the population, and variations in sequence within the population are used as genetic markers for tracking or identifying traits.
In yet a further embodiment, the invention provides kits for use in genotyping one or more target loci in a nucleic acid sample using the unbiased amplification methods provided herein. An exemplary kit according to the present invention comprises a first set of primers wherein each primer in the first set comprises a primer tail sequence and is capable of hybridizing to a target sequence, and a second set of primers wherein each primer in the second set is capable of hybridizing to the primer tail sequence. The kits provided by the invention may further provide reagents for carrying out nucleic acid amplification reactions, such as DNA polymerase, dideoxyribonucleotides with or without detectable labels, and buffer solutions. The kits of the invention may further provide instructions for using the kit components according to the methods provided herein.
“Marker,” “genetic marker,” “molecular marker,” “marker nucleic acid,” and “marker locus” refer to a nucleotide sequence or encoded product thereof (e.g., a protein) used as a point of reference when identifying a DNA locus influencing a phenotype in an organism. A marker can be derived from genomic nucleotide sequence or from expressed nucleotide sequences (e.g., from a spliced RNA, a cDNA, etc.), or from an encoded polypeptide, and can be represented by one or more particular variant sequences, or by a consensus sequence. In another sense, a marker is an isolated variant or consensus of such a sequence. The term also refers to nucleic acid sequences complementary to or flanking the marker sequences, such as nucleic acids used as probes or primer pairs capable of amplifying the marker sequence. A “marker probe” is a nucleic acid sequence or molecule that can be used to identify the presence of a marker locus, e.g., a nucleic acid probe that is complementary to a marker locus sequence. Alternatively, in some aspects, a marker probe refers to a probe of any type that is able to distinguish (i.e., genotype) the particular allele that is present at a marker locus. A “marker locus” is a locus that can be used to track the presence of a second linked locus, e.g., a linked locus that encodes or contributes to expression of a phenotypic trait. For example, a marker locus can be used to monitor segregation of alleles at a locus, such as a quantitative trait locus (QTL), that are genetically or physically linked to the marker locus. Thus, a “marker allele,” alternatively an “allele of a marker locus,” is one of a plurality of nucleotide sequences found at a polymorphic marker locus in a population.
Markers that can be used in the practice of the present invention include, but are not limited to, unique expressed sequence tags (EST); restriction fragment length polymorphisms (RFLP), amplified fragment length polymorphisms (AFLP), simple sequence repeats (SSR), simple sequence length polymorphisms (SSLPs), single nucleotide polymorphisms (SNP), insertion/deletion polymorphisms (Indels), variable number tandem repeats (VNTR), and random amplified polymorphic DNA (RAPD), isozymes, and others known to those skilled in the art. Polymorphisms comprising as little as a single nucleotide change can be assayed in a number of ways. For example, detection can be made by electrophoretic techniques including a single strand conformational polymorphism (Orita et al. (1989) Genomics 8(2), 271-278), denaturing gradient gel electrophoresis (Myers (1985) EPO 0273085), or cleavage fragment length polymorphisms (Life Technologies, Inc., Gaithersburg, Md. 20877), or direct sequencing.
Once a polymorphism in a population is known, assays can be designed to detect alleles at the polymorphic locus in members of the population. Methods for detecting alleles at a polymorphic locus include, e.g., PCR-based sequence specific amplification methods, detection of restriction fragment length polymorphisms (RFLP), detection of isozyme markers, detection of polynucleotide polymorphisms by allele specific hybridization (ASH), detection of amplified variable sequences of the plant genome, detection of self-sustained sequence replication, detection of simple sequence repeats (SSRs), detection of single nucleotide polymorphisms (SNPs), or detection of amplified fragment length polymorphisms (AFLPs). Methods are also known for the detection of expressed sequence tags (ESTs) and SSR markers derived from EST sequences and randomly amplified polymorphic DNA (RAPD).
A marker sequence typically comprises two alleles at each polymorphic locus in a diploid organism. A diploid individual can therefore be either homozygous or heterozygous at a given locus. Homozygosity is a condition in which both alleles at a locus are characterized by the same nucleotide sequence. Heterozygosity refers to the presence of two different alleles at a given locus in a diploid organism.
A favorable allele of a marker is the allele of the marker that co-segregates with a desired phenotype. As used herein, a marker has a minimum of one favorable allele, although it is possible that the marker might have two or more favorable alleles found in the population. Any favorable allele of that marker can be used advantageously for the identification and tracking of favorable traits in a breeding program. Alternatively, a marker allele that co-segregates with an undesirable phenotype may be useful in the invention, since that allele can be used to identify and counter select an unfavorable genotype. Such an allele can be used for exclusionary purposes during breeding to identify individuals having genotypes that negatively correlate with a desired phenotype for elimination during subsequent rounds of breeding.
The more tightly linked a marker is with a polymorphic locus influencing a phenotype, the more reliable the marker is in marker-assisted selection (MAS), as the likelihood of a recombination event unlinking the marker and the locus decreases. Markers containing the causal mutation for a trait, or that are within the coding sequence of a causative gene, are ideal as no recombination is expected between them and the sequence of DNA responsible for the phenotype.
Genetic markers are distinguishable from one another (as well as from the plurality of alleles of any one particular marker) on the basis of polynucleotide length and/or sequence. Genetic markers are known in the art for many well-characterized organisms, and novel markers may also be developed by methods known in the art. In general, any differentially inherited polymorphic trait (including a nucleic acid polymorphism) that segregates among progeny is a potential genetic marker.
Methods for determining the genotype of an organism at a given marker locus include, but are not limited to, PCR-based detection methods, microarray methods, mass spectrometry-based methods and nucleic acid sequencing methods, including whole genome sequencing. In certain embodiments, the detection of alleles at polymorphic sites in a sample of DNA, RNA, or cDNA may be facilitated through the use of nucleic acid amplification methods. Such methods specifically increase the concentration of polynucleotides that span the polymorphic site, or include that site and sequences located either distal or proximal to it. Such amplified molecules can be readily detected by gel electrophoresis, fluorescence detection methods, or other means. One method of achieving such amplification employs the polymerase chain reaction (PCR) (Mullis et al., 1986 Cold Spring Harbor Symp Quant Biol 51:263-273; European Patent 50,424; European Patent 84,796; European Patent 258,017; European Patent 237,362; European Patent 201,184; U.S. Pat. No. 4,683,202; U.S. Pat. No. 4,582,788; and U.S. Pat. No. 4,683,194), using primer pairs that are capable of hybridizing to the proximal sequences that define a polymorphism in its double-stranded form.
In some aspects, methods of the invention utilize an amplification step to genotype a marker locus. Separate detection probes can also be omitted in amplification/detection methods, e.g., by performing a real time amplification reaction that detects product formation by modification of the relevant amplification primer upon incorporation into a product, incorporation of labeled nucleotides into an amplicon, or by monitoring changes in molecular rotation properties of amplicons as compared to unamplified precursors (e.g., by fluorescence polarization).
“Amplifying,” in the context of nucleic acid amplification, is any process whereby additional copies of a selected nucleic acid (or a transcribed form thereof) are produced. In some embodiments, an amplification-based marker technology is used wherein a primer or amplification primer pair is admixed with a nucleic acid sample from an organism, and wherein the primer or primer pair is complementary to or partially complementary to at least a portion of a marker locus, and is capable of initiating DNA polymerization by a DNA polymerase using the nucleic acid sample as a template. The primer or primer pair is extended in a DNA polymerization reaction having a DNA polymerase and a template genomic nucleic acid to generate at least one amplicon.
Typical amplification methods include various polymerase based replication methods, including the polymerase chain reaction (PCR), ligase mediated methods such as the ligase chain reaction (LCR) and RNA polymerase based amplification (e.g., by transcription) methods. An “amplicon” is an amplified nucleic acid, e.g., a nucleic acid that is produced by amplifying a template nucleic acid by any available amplification method (e.g., PCR, LCR, transcription, or the like). A “template nucleic acid” is a nucleic acid that serves as a template in an amplification reaction (e.g., a polymerase based amplification reaction such as PCR, a ligase mediated amplification reaction such as LCR, a transcription reaction, or the like). A template nucleic acid can be genomic in origin, or alternatively, can be derived from expressed sequences, e.g., a cDNA or an EST. Details regarding the use of these and other amplification methods are known in the art, and one of skill will appreciate that essentially any RNA can be converted into a double stranded DNA suitable for restriction digestion, PCR expansion, and sequencing using reverse transcriptase and a polymerase.
In some embodiments, the presence or absence of a molecular marker is determined through detection of a nucleic acid sequence at a polymorphic marker region. In alternative embodiments, in silico methods can be used to detect the marker loci of interest. For example, the sequence of a nucleic acid comprising a marker locus of interest can be stored in a computer. The desired marker locus sequence or its homolog can be identified using an appropriate nucleic acid search algorithm as provided by, for example, in such readily available programs as BLAST, or even simple word processors.
As used herein, the terms “nucleic acid” and “polynucleotide” refer to a deoxyribonucleotide, ribonucleotide, or a mixed deoxyribonucleotide and ribonucleotide polymer in either single- or double-stranded form, and unless otherwise limited, would encompass known analogs of natural nucleotides that can function in a similar manner as naturally-occurring nucleotides. Polynucleotide sequences include the DNA strand sequence that is transcribed into RNA and the strand sequence that is complementary to the DNA strand that is transcribed. Polynucleotide sequences also include both full-length sequences as well as shorter sequences derived from the full-length sequences. Allelic variations of the exemplified sequences also fall within the scope of the subject invention. Polynucleotide sequences include both the sense and antisense strands either as individual strands or in the duplex. The nomenclature used herein is that required by Title 37 of the United States Code of Federal Regulations §1.822 and set forth in the tables in WIPO Standard ST.25 (1998), Appendix 2, Tables 1 and 3.
As used herein, the term “recombinant nucleic acid,” “recombinant polynucleotide” or “recombinant DNA molecule” refers to a polynucleotide that has been altered from its native state, such as by linkage to one or more other polynucleotide sequences to which the recombinant polynucleotide molecule is not normally linked to in nature. Such molecules may or may not be present, for example, in a host genome or chromosome.
The subject invention also concerns oligonucleotide probes and primers, such as polymerase chain reaction (PCR) primers, that can hybridize to a coding or non-coding sequence of a polynucleotide of the present invention. Oligonucleotide probes of the invention can be used in methods for detecting and quantitating nucleic acid sequences. Oligonucleotide primers of the invention can be used in PCR methods and other methods involving nucleic acid amplification. In a preferred embodiment, a probe or primer of the invention can hybridize to a polynucleotide of the invention under stringent conditions. Probes and primers of the invention can optionally comprise a detectable label or reporter molecule, such as fluorescent molecules, enzymes, radioactive moiety (e.g., 3H, 35S, 125I, etc.), and the like. Probes and primers of the invention can be of any suitable length for the method or assay in which they are being employed. Typically, probes and primers of the invention will be 10 to 500 or more nucleotides in length. Probes and primers of the invention can have complete (100%) nucleotide sequence identity with the polynucleotide sequence, or the sequence identity can be less than 100%. For example, sequence identity between a probe or primer and a sequence can be 70% or greater, 75% or greater, 80% or greater, 85% or greater, 90% or greater, or 95% to 100%, or any other percentage sequence identity allowing the probe or primer to hybridize under stringent conditions to a nucleotide sequence of a polynucleotide of the invention. In one embodiment, a probe or primer of the invention has 70% or greater, 75% or greater, 80% or greater, 85% or greater, 90% or greater, or 95% to 100% sequence identity with a nucleotide sequence provided herein, including the complement thereof.
The subject invention also concerns variants of the polynucleotides of the present invention. Variant sequences include those sequences wherein one or more nucleotides of the sequence have been substituted, deleted, and/or inserted. The nucleotides that can be substituted for natural nucleotides of DNA have a base moiety that can include, but is not limited to, inosine, 5-fluorouracil, 5-bromouracil, hypoxanthine, 1-methylguanine, 5-methylcytosine, and tritylated bases. The sugar moiety of the nucleotide in a sequence can also be modified and includes, but is not limited to, arabinose, xylulose, and hexose. In addition, the adenine, cytosine, guanine, thymine, and uracil bases of the nucleotides can be modified with acetyl, methyl, and/or thio groups. Sequences containing nucleotide substitutions, deletions, and/or insertions can be prepared and tested using standard techniques known in the art.
As used herein, the term “percent sequence identity” or “percent identity” refers to the percentage of identical nucleotides in a linear polynucleotide sequence of a reference (“query”) polynucleotide molecule (or its complementary strand) as compared to a test (“subject”) polynucleotide molecule (or its complementary strand) when the two sequences are optimally aligned (with appropriate nucleotide insertions, deletions, or gaps totaling less than 20 percent of the reference sequence over the window of comparison). Optimal alignment of sequences for aligning a comparison window are well known to those skilled in the art and may be conducted by tools such as the local homology algorithm of Smith and Waterman, the homology alignment algorithm of Needleman and Wunsch, the search for similarity method of Pearson and Lipman, and preferably by computerized implementations of these algorithms such as GAP, BESTFIT, FASTA, and TFASTA available as part of the GCG® Wisconsin Package® (Accelrys Inc., Burlington, Mass.). Polynucleotides contemplated within the scope of the subject invention can also be defined in terms of identity and/or similarity ranges with those sequences of the invention specifically exemplified herein. In certain embodiments, the invention provides polynucleotide sequences having at least about 70, 80, 85, 90, 95, 99, or 99.5 percent identity to a polynucleotide sequence provided herein.
The invention also contemplates polynucleotide molecules having sequences which are sufficiently homologous with the polynucleotide sequences exemplified herein so as to permit hybridization with that sequence under standard stringent conditions and standard methods (Maniatis, et al., 1982). As used herein, “stringent” conditions for hybridization refers to conditions wherein hybridization is typically carried out overnight at 20-25 C below the melting temperature (Tm) of the DNA hybrid in 6×SSPE, 5×Denhardt's solution, 0.1% SDS, 0.1 mg/ml denatured DNA. The melting temperature, Tm, is described by the following formula (Beltz, et al., 1983):
T
m=81.5C+16.6 Log [Na+]+0.41(% G+C)−0.61(% formamide)−600/length of duplex in base pairs.
Washes are typically carried out as follows:
In general, synthetic methods for making oligonucleotides, including probes and primers, are known in the art. For example, oligonucleotides can be synthesized chemically according to the solid phase phosphoramidite triester method. Oligonucleotides, including modified oligonucleotides, can also be ordered from a variety of commercial sources.
Any suitable label can be used with a probe of the invention. Detectable labels suitable for use with nucleic acid probes include, for example, any composition detectable by spectroscopic, radioisotopic, photochemical, biochemical, immunochemical, electrical, optical or chemical means. Useful labels include biotin for staining with labeled streptavidin conjugate, magnetic beads, fluorescent dyes, radio labels, enzymes, and colorimetric labels. Other labels include ligands which bind to antibodies labeled with fluorophores, chemiluminescent agents, and enzymes. A probe can also constitute radio labeled PCR primers that are used to generate a radio labeled amplicon. It is not intended that the nucleic acid probes of the invention be limited to any particular size.
In some embodiments, the molecular markers of the invention are detected using a suitable PCR-based detection method, where the size or sequence of the PCR amplicon is indicative of the absence or presence of the marker (e.g., a particular marker allele). In these types of methods, PCR primers are hybridized to the conserved regions flanking the polymorphic marker region. As used in the art, PCR primers used to amplify a molecular marker are sometimes termed “PCR markers” or simply “markers.” It will be appreciated that, although specific examples of primers are provided herein, suitable primers to be used with the invention can be designed using any suitable method. It is not intended that the invention be limited to any particular primer or primer pair.
In some embodiments, the primers of the invention are radiolabelled, or labeled by any suitable means (e.g., using a non-radioactive fluorescent tag), to allow for rapid visualization of the different size amplicons following an amplification reaction without any additional labeling step or visualization step. In some embodiments, the primers are not labeled, and the amplicons are visualized following their size resolution, e.g., following agarose gel electrophoresis. In some embodiments, ethidium bromide staining of the PCR amplicons following size resolution allows visualization of the different size amplicons. It is not intended that the primers of the invention be limited to generating an amplicon of any particular size. For example, the primers used to amplify the marker loci and alleles herein are not limited to amplifying the entire region of the relevant locus. The primers can generate an amplicon of any suitable length that is longer or shorter than those disclosed herein. In some embodiments, marker amplification produces an amplicon at least 20 nucleotides in length, or alternatively, at least 50 nucleotides in length, or alternatively, at least 100 nucleotides in length, or alternatively, at least 200 nucleotides in length.
Marker discovery and development provides the initial framework for marker-assisted breeding programs. Marker-assisted selection (MAS) refers to the selection of individuals based on genetic markers linked to traits of interest during breeding. Individuals may be selected according to their genotype at one or a plurality of marker loci in MAS breeding programs.
In some embodiments of the invention, one or more marker alleles are selected for in a single organism or in a population. In these methods, individuals are selected that contain favorable alleles from more than one marker, or alternatively, favorable alleles from more than one marker are introgressed into a desired population. The determination of which marker alleles correlate with a favorable phenotype is determined for the particular germplasm under study. One of skill recognizes that methods for identifying the favorable alleles are routine and well known in the art, and furthermore, that the identification and use of such favorable alleles is well within the scope of this invention.
Methods of the present invention may evaluate traits including, but not limited to, complex/quantitative traits, monogenic traits, and/or polygenic traits. Such traits in plants may include, for example, reproductive health, plant height, yield, biomass, increased or decreased tolerance to stress, both biotic or abiotic, or to a chemical such as a pesticide or a herbicide, and the like. Such traits in animals may include, for example, weight, weaning weight, carcass composition such as marbling and back fat, hip structure, litter size, fertility, reproductive health, and the like. An “individual” or “subject” in accordance with the present invention may be a plant including, but not limited to an agricultural plant or tree. Agricultural plants or trees as used herein generally refer to plants and trees grown primarily for food or production purposes. Such plants and trees include but are not limited to rice, soybean, corn, canola, sorghum, sugarcane, cotton, coffee, tomato, pine, oak, maple, citrus, or the like. In addition, an “individual” or “subject” may be an animal including, but not limited to a livestock animal. Livestock animals as used herein generally refer to animals raised primarily for food. Such animals include, but are not limited to cattle, swine, horse, goat, sheep, dog, ostrich, chicken, turkey, and the like. As used herein, the term “plant” includes plant cells, plant protoplasts, plant cells of tissue culture from which plants can be regenerated, plant calli, plant clumps and plant cells that are intact in plants or parts of plants such as pollen, flowers, seeds, leaves, stems, and the like.
Certain embodiments of the invention provide early selection of an individual for breeding. Early selection may include selection of an individual for breeding before the individual fully exhibits a trait or phenotype, or before a trait is fully established in an individual.
Embodiments of the invention may provide a kit for determining the genotype of an individual. Such a kit may include means for detecting one or a plurality of genetic markers. In vitro test kits (e.g., reagent kits) for determining the genotype of an individual may include reagents, materials, and protocols for assessing one or more biomarkers (e.g., nucleic acids, proteins, or the like), instructions and, optionally, software for comparing the biomarker data between individuals. Useful reagents and materials for kits include, but are not limited to PCR primers, hybridization probes and primers (e.g., labeled probes or primers), allele-specific oligonucleotides, reagents for genotyping SNP markers, reagents for detection of labeled molecules, restriction enzymes (e.g., for RFLP analysis), DNA polymerases, RNA polymerases, DNA ligases, marker enzymes, microarrays, antibodies, means for amplification of nucleic acid fragments from one or more individuals, means for analyzing the nucleic acid sequence of one or more individuals or fragments thereof, or means for analyzing the sequence of one or more amino acid residues from one or more individuals to be selected for breeding.
The definitions and methods provided define the present invention and guide those of ordinary skill in the art in the practice of the present invention. Unless otherwise noted, terms are to be understood according to conventional usage by those of ordinary skill in the relevant art. Examples of resources describing many of the terms related to molecular biology used herein can be found in Alberts, et al., Molecular Biology of The Cell, 5th Edition, Garland Science Publishing, Inc.: New York, 2007; Rieger, et al., Glossary of Genetics: Classical and Molecular, 5th edition, Springer-Verlag: New York, 1991; King, et al., A Dictionary of Genetics, 6th ed, Oxford University Press: New York, 2002; and Lewin, Genes Icorn, Oxford University Press: New York, 2007. The nomenclature for DNA bases as set forth at 37 CFR §1.822 is used.
“Adjacent,” when used to describe a nucleic acid molecule that hybridizes to DNA containing a polymorphism, refers to DNA sequences that directly abut the polymorphic nucleotide base position. For example, a nucleic acid molecule that can be used in a single base extension assay is “adjacent” to the polymorphism.
“Allele” refers to an alternative nucleic acid sequence at a particular locus; the length of an allele can be as small as 1 nucleotide base, but is typically larger. For example, a first allele can occur on one chromosome, while a second allele occurs on a second homologous chromosome, e.g., as occurs for different chromosomes of a heterozygous individual, or between different homozygous or heterozygous individuals in a population. “Allele frequency” refers to the frequency (proportion or percentage) at which an allele is present at a locus within an individual, within a line, or within a population of lines. For example, for an allele “A,” diploid individuals of genotype “AA,” “Aa,” or “aa” have allele frequencies of 1.0, 0.5, or 0.0, respectively. One can estimate the allele frequency within a line by averaging the allele frequencies of a sample of individuals from that line. Similarly, one can calculate the allele frequency within a population of lines by averaging the allele frequencies of lines that make up the population. For a population with a finite number of individuals or lines, an allele frequency can be expressed as a count of individuals or lines (or any other specified grouping) containing the allele. An allele positively correlates with a trait when it is linked to that trait and when presence of the allele is an indictor that the trait will occur in an individual. An allele negatively correlates with a trait when it is linked to the trait and when presence of the allele is an indicator that the trait will not occur in an individual.
“Genetic element” or “gene” refers to a heritable sequence of DNA, i.e., a genomic sequence, with functional significance. The term “gene” can also be used to refer to, e.g., a cDNA and/or an mRNA encoded by a genomic sequence, as well as to that genomic sequence.
“Genotype” is the genetic constitution of an individual (or group of individuals) at one or more genetic loci, as contrasted with the observable trait (the phenotype). Genotype is defined by the allele(s) at one or more loci that the individual has inherited from its parents. The term genotype can be used to refer to an individual's genetic constitution at a single locus, at multiple loci, or, more generally, the term genotype can be used to refer to an individual's genetic make-up for all the genes in its genome. A “haplotype” is the genotype of an individual at a plurality of genetic loci. Typically, the genetic loci described by a haplotype are physically and genetically linked, i.e., on the same chromosome interval. The terms “phenotype,” or “phenotypic trait” or “trait” refers to one or more traits of an organism. The phenotype can be observable to the naked eye, or by any other means of evaluation known in the art, e.g., microscopy, biochemical analysis, genomic analysis, an assay for a particular disease resistance, etc. In some cases, a phenotype is directly controlled by a single gene or genetic locus, i.e., a “single gene trait.” In other cases, a phenotype is controlled by a plurality of genes or genetic loci.
“Germplasm” refers to genetic material of or from an individual (e.g., a plant), a group of individuals (e.g., a plant line, variety or family), or a clone derived from a line, variety, species, or culture. The germplasm can be part of an organism or cell, or can be separate from the organism or cell. In general, germplasm provides genetic material with a specific molecular makeup that provides a physical foundation for some or all of the hereditary qualities of an organism or cell culture. As used herein, germplasm includes cells, seed or tissues from which new plants may be grown, or plant parts, such as leaves, stems, pollen, or cells that can be cultured into a whole plant.
“Linkage disequilibrium” refers to a non-random segregation of genetic loci or traits (or both). In either case, linkage disequilibrium implies that the relevant loci are within sufficient physical proximity along a length of a chromosome so that they segregate together with greater than random (i.e., non-random) frequency (in the case of co-segregating traits, the loci that underlie the traits are in sufficient proximity to each other). Linked loci co-segregate more than 50% of the time, e.g., from about 51% to about 100% of the time. The term “physically linked” is sometimes used to indicate that two loci, e.g., two marker loci, are physically present on the same chromosome. Recombination between linked loci does not occur during meiosis with high frequency, e.g., linked loci cosegregate at least about 90% of the time, e.g., 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, 99.5%, 99.75%, or more of the time.
“Locus” a chromosome region where a polymorphic nucleic acid, trait determinant, gene, or marker is located. A “gene locus” is a specific chromosome location in the genome of a species where a specific gene can be found.
“Marker Assay” means a method for detecting a polymorphism at a particular locus using a particular method, e.g. measurement of at least one phenotype (such as seed color, flower color, or other visually detectable trait), restriction fragment length polymorphism (RFLP), single base extension, electrophoresis, sequence alignment, allelic specific oligonucleotide hybridization (ASO), random amplified polymorphic DNA (RAPD), microarray-based technologies, and nucleic acid sequencing technologies, etc.
“Marker Assisted Selection” (MAS) is a process by which phenotypes are selected based on marker genotypes.
“Molecular phenotype” is a phenotype detectable at the level of a population of one or more molecules. Such molecules can be nucleic acids, proteins, or metabolites. A molecular phenotype could be an expression profile for one or more gene products, e.g., at a specific stage of plant development, in response to an environmental condition or stress, etc.
“Nucleic acid” refers to any polymer or oligomer of pyrimidine and purine bases, preferably cytosine, thymine, and uracil, and adenine and guanine, respectively (See Albert L. Lehninger, Principles of Biochemistry, at 793-800 (Worth Pub. 1982) which is herein incorporated by reference in its entirety). The present invention contemplates any deoxyribonucleotide, ribonucleotide or peptide nucleic acid component, and any chemical variants thereof, such as methylated, hydroxymethylated or glycosylated forms of these bases, and the like. The polymers or oligomers may be heterogeneous or homogenous in composition, and may be isolated from naturally occurring sources or may be artificially or synthetically produced. In addition, the nucleic acids may be DNA or RNA, or a mixture thereof, and may exist permanently or transitionally in single-stranded or double-stranded form, including homoduplex, heteroduplex, and hybrid states.
“Percent identity” or “% identity” means the extent to which two optimally aligned polynucleotide segments are invariant throughout a window of alignment of components, for example nucleotide sequence or amino acid sequence. An “identity fraction” for aligned segments of a test sequence and a reference sequence is the number of identical components that are shared by sequences of the two aligned segments divided by the total number of sequence components in the reference segment over a window of alignment which is the smaller of the full test sequence or the full reference sequence.
“Phenotype” refers to the detectable characteristics of a cell or organism which can be influenced by genotype.
“Polymorphism” refers to the presence of one or more variations in a population. A polymorphism may manifest as a variation in the nucleotide sequence of a nucleic acid or as a variation in the amino acid sequence of a protein. Polymorphisms include the presence of one or more variations of a nucleic acid sequence or nucleic acid feature at one or more loci in a population of one or more individuals. The variation may comprise but is not limited to one or more nucleotide base changes, the insertion of one or more nucleotides or the deletion of one or more nucleotides. A polymorphism may arise from random processes in nucleic acid replication, through mutagenesis, as a result of mobile genomic elements, from copy number variation and during the process of meiosis, such as unequal crossing over, genome duplication and chromosome breaks and fusions. The variation can be commonly found or may exist at low frequency within a population, the former having greater utility in general breeding programs and the latter may be associated with rare but important phenotypic variation. Useful polymorphisms may include single nucleotide polymorphisms (SNPs), insertions or deletions in DNA sequence (Indels), simple sequence repeats of DNA sequence (SSRs), a restriction fragment length polymorphism, and a tag SNP. A genetic marker, a gene, a DNA-derived sequence, a RNA-derived sequence, a promoter, a 5′ untranslated region of a gene, a 3′ untranslated region of a gene, microRNA, siRNA, a resistance locus, a satellite marker, a transgene, mRNA, ds mRNA, a transcriptional profile, and a methylation pattern may also comprise polymorphisms. In addition, the presence, absence, or variation in copy number of the preceding may comprise polymorphisms. Variations in the DNA sequences of e.g. humans or plants can affect how they handle diseases, bacteria, viruses, chemicals, drugs, etc.
A “population” refers to a set comprising any number, including one, of individuals, objects, or data from which samples are taken for evaluation. Most commonly, the terms relate to a breeding population from which members are selected and crossed to produce progeny in a breeding program. A population can include the progeny of a single breeding cross or a plurality of breeding crosses. The population members need not be identical to the population members selected for use in subsequent cycles of analyses or those ultimately selected to obtain final progeny. Often, a population is derived from a single biparental cross, but may also derive from two or more crosses between the same or different parents. Although a population may comprise any number of individuals, those of skill in the art will recognize that breeders commonly use population sizes ranging from one or two hundred individuals to several thousand, and that the highest performing 5-20% of a population is what is commonly selected to be used in subsequent crosses in order to improve the performance of subsequent generations of the population.
“Primer” refers to an oligonucleotide capable of hybridizing to a target nucleotide sequence to prime the synthesis of DNA by a polymerase. Oligonucleotide primers of the invention can be used in PCR methods and other methods involving nucleic acid amplification. A primer may comprise a “primer tail” which refers to a portion of the primer oligonucleotide sequence which does not hybridize with the target nucleotide sequence.
“Tagging” refers to the addition of a detection label to a nucleic acid sample in order to distinguish it from a second or further nucleic acid sample. Tagging can be performed e.g. by the addition of a sequence identifier or by any other means known in the art. Such sequence identifier can be e.g. a unique base sequence of varying but defined length uniquely used for identifying a specific nucleic acid sample. Typical examples thereof are, for example, ZIP sequences. Using such tag, the origin of a sample can be determined upon further processing. In case of combining processed products originating from different nucleic acid samples, the different nucleic acid samples can be identified using different tags. A “tagged library” refers to a library of tagged nucleic acids.
“Target DNA region” refers to a segment of genomic DNA of one or more nucleotides in length that may or may not be polymorphic in a population.
“Target polymorphism” refers to a specific genomic locus that is known to exhibit one or more variations of a nucleic acid sequence in a population.
“Test sample nucleic acid” refers to a nucleic acid sample that is investigated for polymorphisms.
A set of 150 genomic regions of the cattle genome (Table 1) were amplified by PCR in a multiplex reaction comprising the following reagents.
The multiplex PCR mixture was amplified under the following conditions.
1. 98° C.—10 min
2. 98° C.—2 min
3. 72° C.—1 min
4. Go to step 2, 9 more times (10 cycles total)
5. 72° C.—10 min
6. Hold at 25° C.
Each primer pair used in the reaction comprises a sequence that binds specifically to a region upstream or downstream of a polymorphism of interest as shown in the Forward Primer and Reverse Primer columns of Table 1. Each forward primer sequence further comprises a tail having a sequence of 5′ ACACGACGCTCTTCCGATCT 3′ (SEQ ID NO: 301) at the 5′ end. Each reverse primer sequence further comprises a tail having a sequence of 5′ CTGAACCCTTGTCGCCATTC 3′ (SEQ ID NO: 302) on the 5′ end.
Concentrations of each primer pair were adjusted individually based on the following regression equation:
The equation was developed by counting the number of sequencing reads obtained for all primer pairs, at different concentrations. In the equation, reads stands for the number of reads sequenced for primer i and dilution stands for the concentration level used in the experiment for the same primer i. The equation was used to calculate a primer dilution for a given locus, which is equal to the 1/number of loci to be amplified in the multiplex PCR reaction. Therefore, if amplifying 100 loci, then the primer pairs used to amplify one locus are diluted 1/100. After the outcome of a sequencing run is obtained, the number of reads obtained for primer i, as well as the sum of reads obtained for all primer pairs was plugged into the left side of the formula. Based on that number a new dilution was generated (right side of the formula). That new dilution was then used in future experiments as the final dilution to be used for that pair of primers.
After the first amplification step, amplified DNA was separated from the reaction mixture using AMPure beads, using the following procedure:
After the separation step, the product of the first amplification step was further amplified in a second PCR step using a pair of universal primers that bind to the tail of each primer pair used in the first PCR step. Universal Primer 1 had SEQ ID NO: 303 (5′ AATGATACGGCGACCACCGAGATCTACACNNNNNNNNACACTCTTTCCCTACACGA CGCTCTTCCGATCT 3′) and Universal Primer 2 had SEQ ID NO: 304 (5′ CAAGCAGAAGACGGCATACGAGATNNNNNNNNCGGTCTCGGCATTCCTGCTGAACC CTTGTCGCCATTC 3′), where NNNNNNNN represents an optional index (e.g., bar code) that can be inserted into such primers or other primers prepared according to the invention representing the nucleotides of the index used to identify the sample being processed. The second PCR step included the following reagents and was carried out under the conditions described below. Reagent concentration and sources are the same as those described above for the first PCR step.
The PCR mixture for the second amplification step was amplified under the following conditions:
1. 98° C.—10 min
2. 98° C.—2 min
3. 57° C.—30 sec
4. 72° C.—30 sec
5. Go to step 2, 9 more times (10 cycles total)
6. 72° C.—10 min
7. Hold at 25° C.
After the second amplification step, amplified DNA was separated from the reaction mixture using a Macherey-Nagel NucleoSpin Gel and PCR Clean-up Kit (Clontech) according to the following procedure:
The resulting product corresponded to the DNA library of an individual, containing amplification products of each of 150 regions of the cattle genome. DNA libraries of 24 individuals were quantified, and pooled in equimolar amounts, for sequencing in a HiSeq2500 Illumina DNA sequencer.
The outcome of the sequencing and analysis of a DNA library from one individual according to Example 1 is presented in Table 2. The table shows the average percentage of reads obtained for each of 150 loci relative to the total number of reads produced in the sample, and the standard deviation of 24 replicates. These data demonstrated that the target region was successfully amplified at all target loci, and that the number of reads produced was well-balanced across loci and across replicates.
This application claims the benefit of U.S. provisional application No. 62/046,795, filed Sep. 5, 2014, and is herein incorporated by reference in its entirety.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US15/46899 | 8/26/2015 | WO | 00 |
Number | Date | Country | |
---|---|---|---|
62046795 | Sep 2014 | US |