This disclosure relates to methods and kits for karyotyping in which at least one chromosome is interrogated by amplifying one or more loci that are not within copy number variable regions (CNVR) thereof.
The human haploid genome is a deoxyribonucleic acid (DNA) sequence and consists of approximately three billion base pairs grouped into 23 chromosomes. However, the human genome is diploid and consists of approximately six billion base pairs grouped into 46 chromosomes. Hence, two copies of each genomic segment and two copies of each chromosome are represented in most of the human genome. The exception is the male human, which has only one copy each of chromosome X and chromosome Y. Nevertheless, variable copy numbers (i.e. not two copies) of genomic segments and chromosomes are observed in individual genomic DNA (gDNA) samples. Such copy number variable regions (CNVR) that are typically greater than one kilobase in length and generally occur at a minor frequency of equal to or greater than 1% in the population are termed “copy number variants” (CNVs; see, e.g., Feuk, et al. Nat. Rev. Genet. 7:85-97 (2006)). Copy number variation (CNV) and its mechanisms of formation, associations with phenotype, and methods of analysis have been extensively reviewed (Feuk, supra; Freeman, et al. Genome Res. 16:949-961 (2006); Sharp, et al. Annu. Rev. Genomics Hum. Genet. 7:407-442 (2006); Nature Genetics 39:S1-S54 (2007), entire issue). Currently, approximately 8600 CNVs have been identified and cover about 5-10% of the human genome (Redon, et al. Nature 444:444-454 (2006); Conrad, et al. Nature 464:704-712 (2010)). Continuing studies toward finer mapping of the CNV map and interrogating more diverse gDNA samples are tracked at the Database of Genomic Variants (http://projects.tcag.ca/variation/). Aberrations in the normal complement of 46 human chromosomes have been identified by cytogenetic analysis. Cytogenetics and its history, linking of chromosomal defects and disease, and methods of analysis have been extensively reviewed (Speicher, et al. Nat. Rev. Genet. 6:782-792 (2005); Trask, Nat. Rev. Genet. 3:769-778 (2002)). More recently, polymerase chain reaction (PCR)-based assays have been used to identify such aberrations.
There is a need in the art for methods that provide fast, accurate, easy-to-use and reliable karyotype information. Karyotype includes an analysis of chromosome number, type, shape, and banding. Currently available methods for determining karyotype and chromosome number lack accuracy due to the presence of CNVRs, are labor-intensive and do not provide for simultaneous interrogation of multiple chromosomes without requiring multiple reporting labels. Unless chromosomes are interrogated outside of CNVRs, an inaccurate karyotype or chromosome number determination may result. Accordingly, this disclosure provides such methods by using as target sequences only those targets known to be outside of CNVRs. In addition, the methods described herein provide for simultaneous interrogation of multiple chromosomes using a single, or multiple reporting labels. Using these methods, karyotypes may be rapidly and accurately determined. These and other advantages may be drawn from the description provided below.
For
NA 10859—1, NA10859—2, NA 10859—3, NA 10859—4, NA 14474, NA 14476, NA 17102, NA 17103, NA 17104, NA 17105, NA 17106, NA 17107, NA 17108, NA 17109, NA 17110, NA 17111, NA 17112, NA 17113, NA 17114, NA 17115, NA 17116, NA 17117, NA 17118, NA 17119, NA 17120, NA 17121, NA 17122M NA 17123, NA 17124, NA 17125, NA 17126, NA 17127, NA 17128, NA 17129, NA 17130, NA 17131, NA 17132, NA 17134, NA 17136, NA 17137, NA 17139, NA 17140, NA 17144, NA 17147, NA 17148, NA 17149, NA 17155, NA 17188, NA 17194, NA 17201, NA 17202, NA 17203, NA 17204, NA 17205, NA 17206, NA 17207, NA 17208, NA 17209, NA 17210, NA 17211, NA 17212, NA 17213, NA 17214, NA 17215, NA 17216, NA 17217, NA 17220, NA 17221, NA 17223, NA 17225, NA 17226, NA 17227, NA 17228, NA 17230, NA 17231, NA 17232, NA 17235, NA 17237, NA 17239, NA 17240, NA 17241, NA 17242, NA 17245, NA 17247, NA 17251, NA 17253, NA 17254, NA 17255, NA 17258, NA 17259, NA 17260, NA 17261, NA 17262, NA 17263.
Disclosed herein are methods for karyotype analysis and for determining the copy number of a test locus on a chromosome in a test biological sample, the test locus being located outside of a copy number variable region (CNVR) of the chromosome, by amplifying a nucleic acid sequence corresponding to the test locus to produce an amplification product, or amplicon, and quantifying the relative copy number of the chromosome in the test biological sample relative to a control biological sample. Karyotyping methods are also provided here that do not require comparison to a reference assay or reference sample. Such methods include those using virtual reference assays or virtual calibrator samples.
Embodiments of the methods disclosed herein can be performed on one or multiple samples with or without a calibrator sample or reference assay. In other embodiments, the target copy number assays can serve as the copy number reference assay herein denoted as a “virtual copy number reference assay”.
Variations of these methods, as well as reagents, kits, arrays and devices for carrying out the same are also provided.
This disclosure relates to methods for measuring the copy number of one, more than one, or all chromosomes of a genome using real-time amplification methods (i.e., “karyotyping”). The term “karyotype” refers to the interrogation of one, more than one, or all chromosomes of a particular biological sample, such as a nucleic acid preparation including, but not limited to, deoxyribonucleic acid (DNA, including but not limited to, genomic (gDNA) and mitochondrial DNA (mtDNA)), ribonucleic acid (RNA), proteins, cells, tissues, and organisms.
To “interrogate” one or more chromosomes means to perform an assay, including the methods described herein, to obtain information about the one or more chromosomes contained in the sample. Typically, information relating to “copy number” (CN) is desired. The term copy number generally refers to the number of amplicons resulting from polymerase chain reaction (PCR) amplification, the number of genes or regions on a chromosome, the number of chromosomes in a genome, and the like. The term “amplicon” refers to a piece of DNA formed as the product of an amplification reaction, such as PCR, ligase chain reaction (LCR) or natural gene duplication. The skilled artisan will understand that, within this disclosure, any such molecules can be referred to as target or reference loci depending on the particular assay system. For example, a genome can be queried using as the biological sample a gDNA that preferably, but not necessarily, contains at least one of each chromosome that the user desires to interrogate. In some embodiments, the gDNA sample can contain a pair of each of chromosomes 1-22 and either two X chromosomes, or an X and a Y chromosome. However, it is understood that variations and abnormalities occur such that the number or composition of chromosomes among gDNA samples can vary. In fact, the methods described herein are directed to detecting such variations and abnormalities. In addition, the user may desire to interrogate only a subset of chromosomes (i.e., one or more) present in a genome and it is therefore necessary that only those chromosomes of interest be present. Furthermore, the interrogation typically includes a reference locus (i.e., a “reference locus” or “endogenous control”) which is a locus different from the test locus and is not found within a CNVR. The copy number of the reference locus is typically, but not necessarily, known. A reference locus can reside on a chromosome of interest or on another chromosome and can be interrogated to normalize sample input or to determine relative copy number (i.e., an “endogenous control”).
The methods and assays described herein, which include the determination of the copy number of target loci, are referred to herein as “copy number assays”. The methods and assays described herein for interrogating reference loci are referred to herein as “copy number reference assays”. Typically, the reference locus is present in all samples and is used to normalize for sample input variation and determine relative copy number. However, the test locus can be present or absent. For example, a Y chromosome test assay may interrogate as “present” in a normal male sample, but “absent” in a normal female sample. Existing assay systems, such as the TaqMan® Copy Number Assays (Applied Biosystems), have been used in some instances to interrogate target and reference loci. However, interrogation of test loci within copy number variable regions (CNVRs) has been found to lead to inaccurate karyotype determinations.
An important feature of the methods described herein is that the test loci are outside of (i.e., not within) CNVRs. As described herein, a CNVR includes, but is not limited to, any type of structural variation of a target locus that can affect copy number determination including, but not limited to, duplications, deletions, insertions, repeats, substitutions, and the like. Certain variations can imitate or mimic changes in chromosomal copy number or detrimentally affect the assay by making it less reliable or accurate. Typically, a “universal” set of karyotyping assays (i.e., applicable to all samples or persons) will not target CNVRs. However, it is important to note that CNVRs may vary between individuals or populations of individuals. For instance, a duplication observed in an individual or population of individuals may not be observed in another individual or population of individuals. As such, the methods described herein may target loci that are within a CNVR in one individual or population of individuals but not within a CNVR or another individual or population of individuals.
The difference in copy number between samples is typically determined by “relative quantitation”. Relative quantitation is used to report the copy number of a target locus or chromosome on which the target locus is found in a test sample (i.e., “test locus” or “unknown locus”) relative to one or more reference loci or “calibrator samples”. With reference to the use of a calibrator sample, the term “relative to” means compared to the copy number of a target locus or chromosome on which the target locus is found in a calibrator sample. A “calibrator sample” is a biological sample, often of known karyotype (i.e., a “control gDNA”), that is interrogated before, after or simultaneously with the test sample. The number of copies in the test sample can then be determined relative to the calibrator sample. One or more reference loci or calibrator samples can be used as desired by the user. The one or more reference loci can be the same or different and can be present in the test sample(s) or calibrator sample(s). Typically, the measure of relative quantitation is reported using the term “fold change”, which refers to the amount of amplified product (which relates to the copy number) in a test sample relative to a reference or calibrator sample. Fold change can be quantified using any of several available methods, including but not limited to those described by Livak, et al. (Methods, 25:402-408 (2001)) or commercially available products such as CopyCaller™ (Applied Biosystems, see below). The methods described herein can be performed on one or more than one sample with or without a calibrator sample or reference assay. Various embodiments of such methods are described herein.
The real-time polymerase chain reaction (RT-PCR) is a conventional tool which can be used to amplify and quantify a target nucleic acid molecule, including but not limited to, DNA and RNA, in a sample. The amount of target nucleic acid can be determined as an absolute copy number or as a relative amount. Specifically, the use of RT-PCR to quantify gene expression using the comparative CT method is known to one of skill in the art (Livak, et al. (supra)). In general, the threshold cycle (CT) for a given genetic locus can be determined by arbitrarily setting a signal intensity threshold that falls within the linear range of amplification of real-time PCR data. Previous application of this calculation has been used, for example, to normalize the amplification product of a target gene to the amplification product of an endogenous control gene and then to normalize the data to a calibrator sample, such as an untreated reference sample, to determine the expression of the target gene.
Certain currently available assay systems, such as PCR-based TaqMan® systems, can be useful for relative quantification of nucleic acids amplified from a biological sample to determine gene or chromosome copy number (see, e.g., the TaqMan® assays described in the following Applied Biosystems publications: “Product Bulletin: TaqMan® Copy Number Assays and Custom TaqMan® Copy Number Assays”, “Protocol: TaqMan® Copy Number Assays” (PN 4397425 Rev. C), “Quick Reference Card: TaqMan® Copy Number Assays” (PN 4397424 Rev. B)). The TaqMan® Copy Number Assay and the TaqMan® Copy Number Reference Assay have been used in tandem to interrogate chromosomes to amplify target and control loci, respectively. Usually, the assays are carried out as a single step, “two-plex” real-time PCR assay (i.e., two separate assays within the same reaction well or container). In this type of assay, target and control loci are detected using differently labeled probes. Often, the sample interrogated by such PCR-based assays is human genomic DNA (gDNA) and, typically a minimum of two human gDNA samples are interrogated, including a test and a calibrator sample. Currently, the TaqMan® Copy Number Assay interrogates test loci, which are located throughout the genome and may reside within a copy number variable region (CNVR) and provides copy number variation (CNV) information. The TaqMan® Copy Number Reference Assay interrogates one or more reference loci which, as defined above, are found outside of CNVRs (e.g., Ribonuclease P RNA component H1 (H1RNA) gene (RPPH1) on chromosome 14, cytoband 14q11.2, also known as “RNase P”, and the telomerase reverse transcriptase (TERT) gene located on chromosome 5, cytoband 5p15.33). The TaqMan® Copy Number Reference Assay thereby serves as the endogenous control by interrogating a non-CNVR sequence of interest to normalize for sample input variation by relative quantification (RQ). For example, the TaqMan® Copy Number Assay can utilize an FAM™ dye-labeled minor groove binder (MGB) probe and unlabeled PCR primers directed at the one or more test loci in the test sample or calibrator sample, which are typically assayed separately. The TaqMan® Copy Number Reference Assay can use VIC dye-labeled TAMRA™ probes directed at the one or more reference loci. By comparing the products of these reactions, the user normalizes for sample amount and determines the relative number of copies of the target loci in the sample. However, as described above, interrogation of test loci within copy number variable regions (CNVRs) has been found to lead to inaccurate karyotype determinations. Thus, it is important that the target loci are found outside of (e.g., not within) CNVRs.
The methods described herein can utilize a copy number reference assay. In certain embodiments, the target loci in the test and calibrator sample(s) of, for example, gDNA, can be the same or different. In certain embodimetns, the reference assay is used to normalize for sample input variation between samples (i.e., between multiple test samples or between the test and calibrator samples). Thus, in such assays, one or more reference loci can be amplified from the test or control samples and the copy number (i.e., fold change) is calculated as described herein as a copy number reference assay. Variations of such embodiments are also contemplated as would be understood by one of skill in the art.
The methods described herein can also be conducted using a calibrator sample that, as described above, is typically a biological sample of known karyotype, while the karyotype of the test sample is unknown. The relative difference in copy number between the calibrator sample and the test sample can be calculated as is known in the art or described herein. Typically, the “test assay” and “calibrator assay” (i.e., reference or endogenous control assay) are both used for relative quantification of all samples in order to normalize for sample input variation between samples prior to comparing the different samples themselves. The copy number of the target loci of the test samples is calculated by dividing the fold change of the target copy number in the test sample by that of the calibrator sample and then multiplying by the target copy number of the calibrator sample. This calculation is used in the analysis of the currently and commercially available TaqMan® assay systems and similar systems. The test and calibrator assays are, at least in part, directed to the same target loci, but the target loci can also be different. In some embodiments, the methods can include amplifying a target locus from outside (i.e., not within) a CNVR from a test biological sample (such as a gDNA sample in which the number of copies of a particular chromosome is not known), amplifying the same target loci from a calibrator biological sample, and calculating the copy number of the target loci in the test sample relative to the control sample. Variations of such embodiments, including those that do not use a reference assay or calibrator sample, are also contemplated as would be understood by one of skill in the art.
As discussed above, many embodiments of the methods described herein are contemplated as would be understood by one of skill in the art. Exemplary assay formats using reference or calibrator samples are illustrated in Table 1 and described below.
In certain formats outlined in Table 1, the detectable label can be the same (e.g., formats 2, 4, 5, 7, 8, 10, 11, 13, 14, and 16). In other formats, the detectable labels can be different (e.g., Dye #1, Dye #2, and DNA binding dye as in formats 1, 3, 6, 9, 12, or 15) and are conjugated to probes or primers to provide a differential value as the reaction generates amplification products. One or more DNA binding dyes can be used to detect amplification products generated during a reaction (such as real-time PCR) in a multiplex reaction or in separate 1-plex reactions. Typically, the DNA binding dyes are not conjugated to a probe or primer. The skilled artisan will also note the following from Table 1:
Formats 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, and 16 can also permute and combine (e.g., via average, median, etc.) ΔCT values.
The ΔCT value is defined herein as the CT of the test copy number assay minus the CT of the copy number reference assay. The ΔΔCT value is defined herein as: [the CT of the test assay minus the CT of the reference assay of each test sample] minus [the CT of the test assay minus the CT of the reference assay of the calibrator sample]. In some embodiments, the ΔΔCT values can also be permuted and combined. However, because the ΔCT of the calibrator sample is subtracted from each ΔCT of the test sample, the values being permuted and combined now (ΔΔCT vs. ΔCT) are of a different scale, and track in a “parallel” manner (e.g., ΔCTs of 2, 4, 5 vs. ΔΔCTs of 1, 3, 4 where a value of 1 is subtracted from each). Formats 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, and 16 can also be multiplexed even though the multiplexed assays may not provide a ΔCT value. Modified permutations of different formats are also possible (for example, formats 5 and 14 can be combined for analysis of multiple copy number assays per chromosome).
In one embodiment, karyotyping assay format 1 of Table 1 uses a 2-plex real-time PCR reaction comprised of a copy number assay (or test assay) and a copy number reference assay wherein target loci outside of CNVRs are interrogated. Generally, the detectable labels used by each component of this karyotyping assay are different because the amplified products are contained within the same reaction. Subsequent copy number analysis by the 2−ΔΔCT method is the same as that used in the TaqMan® Copy Number Assay. In this calculation, ΔΔCT determined for each test copy number assay serves as the negative exponent of 2 in order to calculate the fold change for each test copy number assay relative to its corresponding copy number reference assay. This assay format provides for multiplexing of the copy number assay and the copy number reference assay, which minimizes the error of the ΔCT calculation and minimizes the use of reagents and consumables. This karyotyping assay format can be used to interrogate a single chromosome or multiple chromosomes. In certain embodiments, the test copy number assay provides the copy number of the target locus of the test sample, as well as the copy number of the chromosome containing the target locus. Different samples having a particular target locus known to be outside of a CNVR can be used for the reference assay. To properly normalize for sample input variation, the target locus or chromosomal copy number for the reference assay should be the same for all samples. In such embodiments, the copy number assay and copy number reference assay can target different chromosomes. Alternatively, or in addition, each chromosome can be interrogated with both a copy number assay and a copy number reference assay. Typically, but not necessarily, each chromosome is assayed separately from one another.
Karyotyping assay formats 2, 3, and 4 of Table 1 each use two single-plex real-time PCR reactions and the same, two or three different, detectable label(s). One real-time PCR reaction is the copy number assay and the other real-time PCR reaction is the copy number reference assay. The detectable label used in each assay can be the same or different. Subsequent copy number analysis by the 2−ΔΔCT method is performed as previously described hereinabove, with the exception that the replicate CT values of the copy number assay and copy number reference assay are typically, but not necessarily, averaged separately prior to calculating the average ΔCT value. These karyotyping assay designs provide the following advantages: 1) the use of copy number assay and copy number reference assay utilizing the same detectable label conjugated to the probe or primer; 2) the use of previously incompatible assays that cannot be multiplexed due to assay interaction issues or PCR bias of one assay over another; and 3) the use of a DNA binding dye in order to minimize reagent costs. It is possible to use this assay format, or modified versions thereof, to interrogate a single chromosome or multiple chromosomes. Each chromosome is interrogated with both a copy number assay and a copy number reference assay and typically, but not necessarily, each chromosome is assayed separately from one another.
In some embodiments, the karyotyping assay uses one or more copy number assays specific for the chromosome of interest, and one or more copy number reference assays. In some embodiments, the average or median CT values can be used. An average average (or “double average”) CT value is calculated by averaging the multiple average CT values derived from the replicates of each copy number assay when more than one copy number assay is used. One or more copy number assays against a chromosome or chromosomes different from the test sample can be used as controls. These control copy number assays fill the role served by copy number reference assay in karyotyping assay formats 1-4 described above. A “double” average CT value is calculated by averaging the multiple average CT values derived from the replicates of each copy number assay. The “double” average ΔCT value can then be calculated from the two “double” average CT values. Subsequent copy number analysis by the 2−ΔΔCT method is carried out as previously described herein. Karyotyping assay formats 5, 6, 7, 8, 9, 10, 11, 12, and 13 provide the following advantages: 1) allow one to use previously incompatible assays that cannot be multiplexed due to assay interaction issues or PCR bias of one assay over another; 2) allow one to use a DNA binding dye (e.g., SYBR Green I) in order to minimize reagent costs; 3) allow one to eliminate the use of a defined copy number reference assay; and 4) allow one to minimize aberrant copy number observed from a minority of copy number assays by diluting the aberrant average CT values in a “double” average CT value. Karyotyping assay formats 5, 6, 7, 8, 9, 10, 11, 12, or 13, or modified versions thereof, can be used to interrogate a single chromosome or multiple chromosomes.
In other embodiments, the copy number assays can serve as a copy number reference assay herein denoted as a “virtual copy number reference assay”. In these embodiments, a multiplicity of copy number assays are performed and each individual copy number assay is compared to one or more of the remaining copy number assays, which then serve as the copy number reference assays. For example, if ten copy number assays are performed, copy number assay no. 1 is compared to each of copy number assays nos. 2-9, copy number assay no. 2 is compared to each of copy number assays nos. 1 and 3-10, and so forth, for a possible total number of 90 different copy number calculations. For each copy number determination, any two or more of the other nine virtual copy number reference assays can be combined (e.g., via average, median, etc.) toward a total number of ten different analyses. This combining takes place at the virtual copy number reference assay CT level or ΔCT level with one copy number assay and up to nine virtual copy number reference assays. These karyotyping assay formats provide the following advantages: 1) the use of incompatible assays that cannot be multiplexed due to assay interaction issues or PCR bias of one assay over another; 2) the use of a DNA binding dye in order to minimize reagent costs; 3) elimination of a properly defined copy number reference assay; 4) confirmation of copy number using each copy number assay as a virtual copy number reference assay; and 5) for troubleshooting of a copy number assay by the use of each copy number assay as a virtual copy number reference assay.
It is also possible to further modify these karyotyping assay formats. For example, in some embodiments, at least ten different copy number analyses can be used in karyotyping assay formats 14, 15, or 16, and the ten different average ΔCT values derived from the ten different copy number analyses are then averaged toward a “double” average ΔCT value. Subsequent copy number analysis by the 2−ΔΔCT method is performed as previously described herein. This copy number analysis mirrors that of the other assays with the exception that average ΔCT values, not average CT values, are being averaged. In other embodiments, assay formats 14, 15 or 16 can be combined with assay formats 5, 6, 7, 11, 12 or 13 by averaging the average CT values from multiple copy number assays prior to calculating ΔCT values. For example, a karyotyping assay with 10 copy number assays in which five target chromosome 1 and five target chromosome 2, the replicate CT values are averaged prior to averaging the five assays targeting each chromosome. Accordingly, only two “double” average CT values exist, one for each chromosome. A variety of subsequent calculations of ΔCT are possible. These karyotyping assay formats with modified copy number analyses provide the following advantages: 1) the use of previously incompatible assays that can not be multiplexed due to assay interaction issues or PCR bias of one assay over another; 2) the use of a DNA binding dye in order to minimize reagent costs; 3) do not rely on the use of a properly defined copy number reference assay; and 4) minimize aberrant copy number results observed in some copy number assays by diluting the aberrant average ΔCT values in a “double” average ΔCT value. Such assay formats can be used to interrogate a single chromosome or multiple chromosomes. It should also be noted that the real-time PCR reactions can utilize increased multiplexing.
Karyotyping methods that use a virtual copy number reference assay (i.e., those that do not utilize either a reference assay or calibrator sample) are also provided. Situations in which a reference assay is not used include the lack of a suitable reference genomic sequence, non-identification of a control assay locus from a reference genomic sequence, assay unavailability, cost reduction, or to facilitate ease of use. Situations in which a calibrator sample is not used include unavailability of such a sample or a reference sequence, to reduce costs, or facilitate ease of use. Such assays are described below.
In the following two embodiments, the copy number assay CT values and copy number reference assay CT values for the virtual calibrator sample are calculated. The copy number assay CT values can be calculated by averaging the replicate CT for each copy number assay and then averaging these values across the multiple test samples. Alternatively, the copy number assay CT values can be calculated by averaging replicate CT values of each copy number assay, averaging these values across multiple copy number assays (those targeting the same chromosome), and averaging these values across the multiple test samples. These methods can also be applied to the copy number reference assay CT values. Depending upon the number of copy number assays and copy number reference assays used, up to four formats can be used to perform this calculation. Furthermore, the copy number assays and the copy number reference assays for the test samples can be treated separately or in combination as outlined in Table 1.
In the first of these embodiments, a “virtual calibrator sample” can be used in lieu of an actual calibrator sample (i.e., virtual calibrator assay no. 1). Typically, in these karyotyping assays, one or more copy number assays targeting one or more test chromosomes, one or more copy number reference assays targeting one or more reference chromosomes, and two or more test samples (e.g., gDNA) are utilized. In these assays, the copy number of the reference chromosome may or may not be previously known. An average CT value is calculated from the replicate CT values of copy number assay. If more than one copy number assay is used, the average CT values of each copy number assay are averaged. This calculation is performed for each of the two or more test samples, and for the one or more reference assay(s). The average CT values or “double” average CT values can then be further averaged across the two or more test samples. This “double” or “triple” average CT value is representative of a “virtual” calibrator sample. If appropriate, other calculations, such as the median, can be used in lieu of the average. The calculated CT values from the copy number assay(s) and copy number reference assay(s) can then be used to calculate the ΔCT value for each of the test samples. Subsequent copy number analysis by the 2−ΔΔCT method (see below) is similar to that as previously described herein. More accurate results will typically result when the test samples are normalized to the same concentration prior to chromosomal copy number analysis and when additional reference assays are utilized. Even in the presence of equal numbers of target loci, different reference assays may provide different CT values due to differences in assay performance. Therefore, the use of additional reference assays can provide a more accurate CT representation of each test gDNA sample. Likewise, the use of additional test samples can provide a more accurate CT representation of the “virtual” calibrator sample.
In another of these embodiments, a modified “virtual calibrator sample” (i.e., virtual calibrator assay no. 2) can be used in lieu of an actual calibrator sample as described in the immediately preceding paragraph. However, after averaging replicate CT values, the average CT values from one copy number reference assay can be averaged across the multiple test samples. This “double” average CT value represents the CT value from the virtual calibrator sample. If multiple copy number reference assays are used, average CT values from replicates can be averaged across multiple copy number reference assays and then averaged across multiple test samples toward a “triple” average CT value representing the virtual calibrator sample. Multiple copy number assays can be treated separately or combined.
In yet other embodiments, multiple copy number assays with virtual calibrator assays are used, for example, multiple copy number assays targeting one or more test loci on a chromosome (virtual calibrator assay no. 1 targeting the test sample); multiple copy number assays targeting one or more test locus on a reference chromosome (virtual calibrator assay no. 2 targeting the virtual calibrator sample); multiple copy number reference assays targeting a reference chromosome (virtual calibrator assay no. 3, copy number reference assay targeting the test sample and the virtual calibrator sample). In these embodiments, the copy number of the reference chromosome may or may not be previously known. The copy number of the targeted chromosome(s) in virtual calibrator assays no. 2 and no. 3 need not be the same. In these embodiments, a calibrator sample is not required. “Double” average CT values are calculated from each of the three assay groups. The “double” average CT values from virtual calibrator assays no. 1 and no. 3 are then used to calculate the ΔCT value for the test sample. If appropriate, other calculations, such as the median, can be used in lieu of the average. The “double” average CT values from virtual calibrator assays no. 2 and no. 3 are used to calculate the ΔCT value for the “virtual” calibrator sample. Subsequent copy number analysis by the 2−ΔΔCT method is similar to that previously described herein. It is noted that the use of additional test and reference assays can provide a more accurate CT representation of the test gDNA sample and “virtual” calibrator gDNA sample.
Yet another embodiment is similar to that described in the immediately preceding paragraph except that virtual calibrator assay #3 is omitted. Because the test sample and “virtual” calibrator sample are derived from the same biological sample, there is no need to normalize for input variation by using one or more reference assays. Therefore, the reference assays can be removed from the workflow. Upon calculating the “double” average CT values from virtual calibrator assays no. 1 and no. 2, the “double” average CT value from virtual calibrator assay no. 2, representing the “virtual” calibrator sample, is subtracted from the “double” average CT value from assay group no. 1, representing the test sample, to calculate a ΔCT value. If appropriate, other calculations, such as the median, can be used in lieu of the average. Because the reference assays are not used, calculation of relative quantity via the equation of 2−ΔΔCT is replaced with calculation of relative quantity via the equation of the 2−ΔCT. Subsequent copy number analysis is similar to that previously described herein.
As described above, data generated using the karyotyping methods described herein (with or without reference or calibrator samples) can be analyzed using any of several well-known methods or algorithms. As an example, relative copy number can be calculated using assay format 1 (of Table 1) as shown in Table 2:
As understood by those of skill in the art, a one cycle difference is theoretically a two-fold difference in quantity of the target sample. The skilled artisan can derive from step 2 that the Test or Unknown Sample (CT value=23) can have eight-fold fewer copies than the Calibrator Sample (CT value=20). However, because the copy number reference assay targets the same number of copies in all samples, the test sample (CT value=23) actually contains four-fold less sample input the calibrator sample (CT value=21). Because of the normalization of sample input variation via the copy number reference assay, it can be determined that the test sample has two-fold fewer copies than the calibrator sample (i.e., 4-fold variation is accounted for as sample input variation from the initial 8-fold variation) indicating that the actual copy number variation between the test and calibrator samples is two-fold.
The data generated using the methods described herein can also be analyzed using CopyCaller™ software available from Applied Biosystems. CopyCaller™ software is based on the 2−ΔΔCT method described by Livak, et al. (supra). In one method, the user is required to select one sample as the calibrator sample and to identify the copy number thereof. For example, if the target loci in the test sample (such as gDNA) is detected using VIC™ dye and the target loci in the calibrator sample is detected using FAM™ dye, the ΔCT value is calculated for each reaction well by subtracting the VIC™ CT value from the FAM™ CT value. The calculation of the ΔCT value normalizes for sample input variation between wells. Because multiple replicates are typically assayed, the ΔCT values from each replicate for each sample are averaged toward an average ΔCT value for each sample. The ΔΔCT value is calculated for each sample by subtracting the average ΔCT value of the calibrator sample from the average ΔCT value of each test sample. The calculation of the ΔΔCT value provides for determining the copy number of a test sample relative to the calibrator sample of known copy number. The ΔΔCT values are converted to a relative quantity (RQ) by the equation of RQ=2−ΔΔCT. The RQ value of each sample is then multiplied by the copy number of the calibrator sample to obtain the copy number determination of each sample.
Another method of copy number analysis by CopyCaller™ utilizes proprietary algorithms for determining copy number without selection of a calibrator sample. The user is only required to identify the most common copy number to be provided by the samples. Both methods of copy number analysis by CopyCaller™ employ outlier removal strategies for all CT and ΔCT values. CopyCaller™ can also provide confidence and absolute Z-score values for the predicted copy number values. It is also noted that copy number analysis methods other than the 2−ΔΔCT method can be utilized (e.g., Cikos, et al. Anal Biochem. 384:1-10 (2009)). Other methods of analysis can also be used, as would be understood by one of skill in the art.
While many of the embodiments described herein relate to PCR, other amplification or product detection methods or reagents can also be utilized. For example, such methods or reagents may include any of several methods that can be used to amplify the target nucleic acid from the sample. The term “amplifying” which typically refers to an “exponential” increase in target nucleic acid is used herein to describe both linear and exponential increases in the numbers of a select target sequence of nucleic acid. The term “amplification reaction mixture” refers to an aqueous solution comprising the various reagents used to amplify a target nucleic acid. These include enzymes, including, but not limited to, polymerases (such as DNA polymerase, RNA polymerase and reverse transcriptase), aqueous buffers, salts, amplification primers, target nucleic acid, and nucleoside triphosphates. Depending upon the context, the mixture can be either a complete or incomplete amplification reaction mixture. The method used to amplify the target nucleic acid can be any method available to one of skill in the art. Any in vitro means for multiplying or amplifying the copies of a target sequence of nucleic acid can be utilized. These include linear, logarithmic, or any other amplification method. Exemplary methods include polymerase chain reaction (PCR; see, e.g., U.S. Pat. Nos. 4,683,202; 4,683,195; 4,965,188; and 5,035,996), isothermal procedures (using one or more RNA polymerases (see, e.g., WO 2006/081222), strand displacement (see, e.g., U.S. Pat. No. RE39,007), partial destruction of primer molecules (see, e.g., WO 2006/087574), ligase chain reaction (LCR) (see, e.g., Wu, et al. Genomics 4:560-569 (1989) and Barany, Proc. Natl. Acad. Sci. USA 88:189-193 (1991)), Qβ RNA replicase systems, RNA transcription-based systems (e.g., TAS, 3SR), rolling circle amplification (RCA) (see, e.g., U.S. Pat. No. 5,854,033; U.S. Pub. No. 2004/0265897; Lizardi, et al. Nat. Genet. 19: 225-232 (1998); and Barter, et al. Nucleic Acid Res. 26:5073-5078 (1998)), and strand displacement amplification (SDA) (Little, et al. Clin. Chem. 45:777-784 (1999)), among others. Many systems are suitable for use in amplifying target nucleic acids and are contemplated herein as would be understood by one of skill in the art.
Any of several methods can be used to detect amplified target nucleic acids using primers or probes. Many different reagents, systems, or detectable labels can be used in the methods described herein. These include, for example, TaqMan® systems, detectable nucleic acid intercalating agents, such as ethidium bromide and SYBR dyes, detectable label-quencher systems (e.g., FRET, salicylate/DTPA ligand systems (see, e.g., Oser et al. Angew. Chem. Int. Ed. Engl. 29:1167-1169 (1990)), displacement hybridization, homologous probes, assays described in EP 070685), molecular beacons (e.g., NASBA), Scorpion, locked nucleic acid (LNA) bases (Singh, et al. Chem. Commun. 4:455-456 (1998)), peptide nucleic acid (PNA) probes (Pellestor, et al. Eur. J. Hum. Gen. 12:694-700 (2004)), Eclipse probes (Afonina, et al. Biotechniques 32:940-949 (2002)), light-up probes (Svanvik, et al. Anal. Biochem. 281:26-35 (2000)), molecular beacons (Tyagi, et al. Nat. Biotechnol. 14:303-308 (1996)), tripartite molecular beacons (Nutiu, Nucleic Acids Res. 30:E94 (2002)), QuantiProbes (www.qiagen.com), HyBeacons (French, et al. Mol. Cell. Probes 15:363-374 (2001)), displacement probes (Li, et al. Nucleic Acids Res. 30:E5 (2002)), HybProbes (Cardullo, et al. Proc. Natl. Acad. Sci. USA 85:8790-8794 (1988)), MGB Alert® probes (www.nanogen.com), Q-PNA (Fiandaca, et al. Genome Res. 11:609-613 (2001)), Plexor (www.Promega.com), LUX primers (Nazarenko, et al. Nucleic Acids Res. 30:E37 (2002)), Scorpion primers (Whitcombe, et al. Nat. Biotechnol. 17:804-807 (1999)), AmpliFluor® (Sunrise) primers (Nazarenko, et al. Nucleic Acids Res. 25:2516-2521 (1997)), DzyNA primers (Todd, et al. Clin. Chem. 46:625-630 (2000)), and the like. In each of these assays, the generation of amplification products can be monitored while the reaction is in progress. An apparatus for detecting the signal generated by the detectable label can be used to detect, measure, and quantify the signal before, during, or after amplification. The particular type of signal may dictate the choice of detection method. For example, in some embodiments, fluorescent dyes are used to label probes or amplified products. The probes bind to single-stranded or double-stranded amplified products, or the dyes intercalate into the double-stranded amplified products, and consequently, the resulting fluorescence increases as the amount of amplified product increases. In some embodiments, the Tm is ascertained by observing a fluorescence decrease as the double-stranded amplified product dissociates and the intercalating dye is released therefrom. The amount of fluorescence can be quantitated using standard equipment such as a spectra-fluorometer, for example. The use of other methods or reagents is also contemplated herein as would be understood by one of skill in the art.
One exemplary method for amplifying and detecting target nucleic acids is the TaqMan® system described above (see, e.g., U.S. Pat. Nos. 4,889,818; 5,079,352; 5,210,015; 5,436,134; 5,487,972; 5,658,751; 5,538,848; 5,618,711; 5,677,152; 5,723,591; 5,773,258; 5,789,224; 5,801,155; 5,804,375; 5,876,930; 5,994,056; 6,030,787; 6,084,102; 6,127,155; 6,171,785; 6,214,979; 6,258,569; 6,814,934; 6,821,727; 7,141,377; and 7,445,900). As described herein and elsewhere, TaqMan® assays are typically carried out by performing nucleic acid amplification on a target polynucleotide using a nucleic acid polymerase having 5′-to-3′ nuclease activity, a primer capable of hybridizing to said target polynucleotide, and an oligonucleotide probe capable of hybridizing to said target polynucleotide 3′ relative to said primer. The oligonucleotide probe typically includes a detectable label (e.g., a fluorescent reporter molecule) and a quencher molecule capable of quenching the fluorescence of said reporter molecule. The oligonucleotide probe typically exists in at least one single-stranded conformation when unhybridized where said quencher molecule quenches the fluorescence of said reporter molecule. When hybridized to a target nucleic acid, the probe exhibits at least one other conformation in which the fluorescence of the reporter molecule is unquenched. Typically, the detectable label and quencher molecule are part of a single probe. As amplification proceeds, the polymerase digests the probe to separate the detectable label from the quencher molecule. The detectable label (e.g., fluorescence) is monitored during the reaction, where detection of the label corresponds to the occurrence of nucleic acid amplification (e.g., the higher the signal the greater the amount of amplification). Variations of TaqMan® assays (e.g., LNA™ spiked TaqMan® assay) are known in the art and would be suitable for use in the methods described herein.
Another exemplary system utilizes double-stranded probes in displacement hybridization methods (see, e.g., Morrison, et al. Anal. Biochem. 183:231-244 (1989); and Li, et al., (supra)). In such methods, the probe typically includes two complementary oligonucleotides of different lengths where one includes a detectable label and the other includes a quencher molecule. When not bound to a target nucleic acid, the quencher suppresses the signal from the detectable label. The probe becomes detectable upon displacement hybridization with a target nucleic acid. Multiple probes can be used, each containing different detectable labels, such that multiple target nucleic acids may be queried in a single reaction.
Additional exemplary methods for amplifying and detecting target nucleic acids involve “molecular beacons”, which are single-stranded hairpin shaped oligonucleotide probes. In the presence of the target sequence, the probe unfolds, binds and emits a signal (e.g., fluoresces). A molecular beacon typically includes at least four components: 1) the “loop”, an 18-30 nucleotide region which is complementary to the target sequence; 2) two 5-7 nucleotide “stems” found on either end of the loop and being complementary to one another; 3) at the 5′ end, a detectable label; and 4) at the 3′ end, a quencher dye that prevents the detectable label from emitting a signal when the probe is in the closed loop shape (i.e., not bound to a target nucleic acid). Thus, in the presence of a complementary target, the “stem” portion of the beacon separates out resulting in the probe hybridizing to the target. Other types of molecular beacons are also known and may be suitable for use in the methods described herein. Molecular beacons can be used in a variety of assay systems. One such system is nucleic acid sequence-based amplification (NASBA®), a single step isothermal process for amplifying RNA to double stranded DNA without temperature cycling. A NASBA® reaction typically requires avian myeloblastosis virus (AMV), reverse transcriptase (RT), T7 RNA polymerase, RNase H, and two oligonucleotide primers. After amplification, the amplified target nucleic acid can be detected using a molecular beacon. Other uses for molecular beacons are known in the art and would be suitable for use in the methods described herein.
The Scorpion system is another exemplary assay format that can be used in the methods described herein. Scorpion primers are bi-functional molecules in which a primer is covalently linked to the probe, along with a detectable label (e.g., a fluorophore) and a quencher. In the presence of a target nucleic acid, the detectable label and the quencher separate which leads to an increase in signal emitted from the detectable label. Typically, a primer used in the amplification reaction includes a probe element at the 5′ end along with a “PCR blocker” element (e.g., HEG monomer) at the start of the hairpin loop. The probe typically includes a self-complementary stem sequence with a detectable label at one end and a quencher at the other. In the initial amplification cycles (e.g., PCR), the primer hybridizes to the target and extension occurs due to the action of the polymerase. The Scorpion system can be used to examine and identify point mutations using multiple probes that have different tags to distinguish between the probes. Using PCR as an example, after one extension cycle is complete, the newly synthesized target region will be attached to the same strand as the probe. Following the second cycle of denaturation and annealing, the probe and the target hybridize. The hairpin sequence then hybridizes to a part of the newly produced PCR product. This results in the separation of the detectable label from the quencher and causes emission of the signal. Other uses for molecular beacons are known in the art and would be suitable for use in the methods described herein.
One or more detectable labels or quenching agents are typically attached to a primer or probe. The detectable label can emit a signal when free or when bound to the target nucleic acid. The detectable label can also emit a signal when in proximity to another detectable label. Detectable labels can also be used with quencher molecules such that the signal is only detectable when not in sufficiently close proximity to the quencher molecule. For instance, in some embodiments, the assay system can cause the detectable label to be liberated from the quenching molecule. Any of several detectable labels can be used to label the primers and probes used in the methods described herein.
As mentioned above, in some embodiments the detectable label can be attached to a probe, which can be incorporated into a primer, or may otherwise bind to amplified target nucleic acid (e.g., a detectable nucleic acid binding agent such as an intercalating or non-intercalating dye). When using more than one detectable label, each should differ in their spectral properties such that the labels can be distinguished from each other, or such that together the detectable labels emit a signal that is not emitted by either detectable label alone. Exemplary detectable labels include, for instance, a fluorescent dye or fluorphore (e.g., a chemical group that can be excited by light to emit fluorescence or phosphorescence), “acceptor dyes” capable of quenching a fluorescent signal from a fluorescent donor dye, and the like. Suitable detectable labels may include, for example, fluorosceins (e.g., 5-carboxy-2,7-dichlorofluorescein; 5-Carboxyfluorescein (5-FAM); 5-HAT (Hydroxy Tryptamine); 6-HAT; 6-JOE; 6-carboxyfluorescein (6-FAM); FITC); Alexa fluors (e.g., 350, 405, 430, 488, 500, 514, 532, 546, 555, 568, 594, 610, 633, 635, 647, 660, 680, 700, 750); BODIPY fluorophores (e.g., 492/515, 493/503, 500/510, 505/515, 530/550, 542/563, 558/568, 564/570, 576/589, 581/591, 630/650-X, 650/665-X, 665/676, FL, FL ATP, FI-Ceramide, R6G SE, TMR, TMR-X conjugate, TMR-X, SE, TR, TR ATP, TR-X SE), coumarins (e.g., 7-amino-4-methylcoumarin, AMC, AMCA, AMCA-S, AMCA-X, ABQ, CPM methylcoumarin, coumarin phalloidin, hydroxycoumarin, CMFDA, methoxycoumarin), calcein, calcein AM, calcein blue, calcium dyes (e.g., calcium crimson, calcium green, calcium orange, calcofluor white), Cascade Blue, Cascade Yellow; Cy™ dyes (e.g., 3, 3.18, 3.5, 5, 5.18, 5.5, 7), cyan GFP, cyclic AMP Fluorosensor (FiCRhR), fluorescent proteins (e.g., green fluorescent protein (e.g., GFP. EGFP), blue fluorescent protein (e.g., BFP, EBFP, EBFP2, Azurite, mKalama1), cyan fluorescent protein (e.g., ECFP, Cerulean, CyPet), yellow fluorescent protein (e.g., YFP, Citrine, Venus, YPet), FRET donor/acceptor pairs (e.g., fluorescein/tetramethylrhodamine, IAEDANS/fluorescein, EDANS/dabcyl, fluorescein/fluorescein, BODIPY FL/BODIPY FL, Fluorescein/QSY7 and QSY9), LysoTracker and LysoSensor (e.g., LysoTracker Blue DND-22, LysoTracker Blue-White DPX, LysoTracker Yellow HCK-123, LysoTracker Green DND-26, LysoTracker Red DND-99, LysoSensor Blue DND-167, LysoSensor Green DND-189, LysoSensor Green DND-153, LysoSensor Yellow/Blue DND-160, LysoSensor Yellow/Blue 10,000 MW dextran), Oregon Green (e.g., 488, 488-X, 500, 514); rhodamines (e.g., 110, 123, B, B 200, BB, BG, B extra, 5-carboxytetramethylrhodamine (5-TAMRA), 5 GLD, 6-Carboxyrhodamine 6G, Lissamine, Lissamine Rhodamine B, Phallicidine, Phalloidine, Red, Rhod-2,5-ROX (carboxy-X-rhodamine), Sulphorhodamine B can C, Sulphorhodamine G Extra, Tetramethylrhodamine (TRITC), WT), Texas Red, Texas Red-X, VIC and other labels described in, e.g., US Pub. No. 2009/0197254), among others as would be known to those of skill in the art.
The amplified target nucleic acid can also be detected using a detectable nucleic acid binding agent (see, e.g., Table 1) which can be, for example, an intercalating agent or a non-intercalating agent. As used herein, an intercalating agent is an agent or moiety capable of non-covalent insertion between stacked base pairs of a double-stranded nucleic acid molecule. A non-intercalating agent is one that does not insert into the double-stranded nucleic acid molecule. The nucleic acid binding agent can produce a detectable signal directly or indirectly. The signal can be detectable directly using, for example, fluorescence or absorbance, or indirectly using, for example, any moiety or ligand that is detectably affected by its proximity to double-stranded nucleic acid is suitable, for example a substituted label moiety or binding ligand attached to the nucleic acid binding agent. It is typically necessary for the nucleic acid binding agent to produce a detectable signal when bound to a double-stranded nucleic acid that is distinguishable from the signal produced when that same agent is in solution or bound to a single-stranded nucleic acid. For example, intercalating agents, such as ethidium bromide and SYBR dyes, fluoresce more intensely when intercalated into double-stranded DNA than when bound to single-stranded DNA, RNA, or in solution (see, e.g., U.S. Pat. Nos. 5,994,056; 6,171,785; and 6,814,934). Similarly, actinomycin D fluoresces red when bound to single-stranded nucleic acids, and green when bound to double-stranded nucleic acids. And in another example, the photoreactive psoralen 4-aminomethyle-4-5′,8-trimethylpsoralen (AMT) has been reported to exhibit decreased absorption at long wavelengths and fluorescence upon intercalation into double-stranded DNA (Johnston et al. Photochem. Photobiol. 33:785-791 (1981). For example, U.S. Pat. No. 4,257,774 describes the direct binding of fluorescent intercalators to DNA (e.g., ethidium salts, daunomycin, mepacrine and acridine orange, 4′,6-diamidino-α-phenylindole). Non-intercalating agents (e.g., minor groove binders such as Hoechst 33258, distamycin, netropsin) may also be suitable for use. For example, Hoechst 33258 (Searle, et al. Nucleic Acids Res. 18:3753-3762 (1990)) exhibits altered fluorescence with an increasing amount of target. Exemplary detectable DNA binding agents may include, for example, acridine derivatives (e.g., acridine homodimer, acridine orange, acridine yellow, 9-amino-6-chloro-2-methoxyacridine (ACMA), proflavin), actinomycins (e.g., actinomycin D (Jain, et al. J. Mol. Biol. 68:1-10 (1972), 7-amino-actinomycin D (7-AAD)), anthramycin, auramine, azure B, BOBO™-1, BOBO™-3, BO-PRO™-1, BO-PRO™-3, chromomycin (e.g., A3), crystal violet, cyanine dyes, DAPI (Kapuściński, et al. Nucleic Acids Res. 6:3519-3534 (1979)), 4′,6-diamidino-2-phenylindole (DAPI), daunomycin, distamycin (e.g., distamycin D), dyes described in U.S. Pat. No. 7,387,887, ellipticine, ethidium salts (e.g., ethidium bromide, ethidium homdimer-1, ethidium homdimer-2, dihydroethidium (also known as hydroethidine), ethidium monoazide), fluorcoumanin, fluorescent intercalators as described in U.S. Pat. No. 4,257,774, GelStar®(Cambrex Bio Science Rockland Inc., Rockland, Me.), hexidium iodide, Hoechst 33258 (Searle, et al., (supra)), Hoechst 33342, Hoechst 34580, homidium, hydroxystilbamidine, JO-JO-1, JO-PRO™-1, LDS 751, LOLO-1, LO-PRO™-1, malachite green, mepacrine (e.g., orange), mithramycin, netropsin, the Nissl substance, 4′,6-diamidino-α-phenylindole, proflavine, POPO™-1, POPO™-3, PO-PRO™-1, propidium iodide, ruthenium polypyridyls, Sevron dyes (e.g., Brilliant Red 2B, Brilliant Red 4G, Brilliant Red B, Orange, Yellow L), SYBR 101, SYBR 102, SYBER 103, SYBR® Gold, SYBR® Green I (U.S. Pat. Nos. 5,436,134 and 5,658,751), SYBR® Green II, SYTOX® Blue, SYTOX® Green, SYTOX® Orange, SYTO® 1, SYTO® 11, SYTO® 13, SYTO® 14, SYTO® 15, SYTO® 16, SYTO® 17, SYTO® 18, SYTO® 20, SYTO® 21, SYTO® 22, SYTO® 23, SYTO® 24, SYTO® 25, SYTO® 40, SYTO® 43, SYTO® 44, SYTO® 45, SYTO® 59, SYTO® 60, SYTO® 61, SYTO® 62, SYTO® 63, SYTO® 64, SYTO® 80, SYTO® 81, SYTO® 82, SYTO® 83, SYTO® 84, SYTO® 85, thiazole orange (Aldrich Chemical Co., Milwaukee, Wis.), TO-PRO-1, TO-PRO-3, TO-PRO-5, TOTO-1, TOTO-2, TOTO™-3, YO-PRO®-1, YO-PRO®-3, YOYO-1, and YOYO®-3 (Molecular Probes, Inc., Eugene, Oreg.), among others. SYBR® Green I (see, e.g., U.S. Pat. Nos. 5,436,134; 5,658,751; and/or 6,569,927), for example, has been used to monitor a PCR reaction by amplifying the target sequence in the presence of the dye, exciting the biological sample with light at a wavelength absorbed by the dye and detecting the emission therefrom; and, determining a melting profile of the amplified target sequence. The presence of amplified products and, therefore, the target sequence in the sample, can thereafter be determined by, for example, performing a melting curve analysis (i.e., non-linear least squares regression of the sum of multiple gaussians). It is to be understood that the use of the SYBR® Green dye is presented as an example and that many such dyes may be used in the methods described herein. Other nucleic acid binding agents can also be suitable as would be understood by one of skill in the art.
The nature of the test and control amplified products can vary depending on the type of assay system that is utilized. For example, the product of PCR amplification of DNA is typically referred to as an “amplicon”. In contrast, the Ligation Chain Reaction product is instead referred to as “LCR product” or “ligation product”. In PCR, primers are utilized but in other methods, such as LCR ligation probes, can be utilized. It is known that both PCR and LCR functions through exponential amplification or linear amplification. It is important to note, however, that regardless of the kind of method used to perform amplification or detection of target loci, any target loci should exhibit similar “amplification efficiency” to one another (e.g., within about 10 percent). This means that within a particular reaction, the target loci being compared should amplify to approximately the same extent (e.g., within about 10 percent) under similar reaction conditions. It is understood in the art that amplification efficiency may be affected by, for example, sample quality (e.g., purity, presence of reaction inhibitors, and the like) or sequence (e.g., G/C content, mismatches, and the like). As used herein, amplification efficiency refers to any product that may be quantified to determine copy number (e.g., a PCR amplicon, an LCR ligation product, or similar product). Variations in amplification due to factors other than copy number may thereby be limited.
For example, when using PCR, it may be optimal to amplify amplicons of a similar size (e.g., within about ten percent of the length of any other amplicon) and amplify with similar amplification efficiency (e.g., within about ten percent of one another). PCR-related amplification efficiencies can be tested as described by Livak, et al. (supra), or as is otherwise known in the art. Briefly, Livak, et al. teach that, as the ΔCT varies with template dilution, amplifications can be performed over a particular range (e.g., a 100-fold range) and the ΔCT calculated at each point. A log of dilution versus ΔCT is made and if the absolute value of the slope is close to zero, the efficiencies are similar. For PCR-based assays, the resultant amplicons are typically about 70 to 125 nucleotides in length, and can be less than or equal to about 110 nucleotides in length. The amplification efficiency is typically controlled by amplifying amplicons of similar G/C content, length, or melting temperature. The primers used in the test and control assays are also typically of similar character (e.g., G/C content, length, melting temperature). Variations of these parameters can also be used, as would be understood by the skilled artisan.
The methods described herein typically require a target locus to be present on the one or more chromosomes being studied (exceptions include, for example, the absence of a target locus such as when the Y chromosome is assayed in a normal (i.e., XX) female sample). Pre-existing bioinformatics algorithms can be used to generate “qualified” pre-characterized genome targets and to design assays against those targets. Qualified genome targets may include, for example, genes, exons, introns, intergenic regions, junctions thereof, evolutionarily conserved regions (“cold” spots not prone to change), or other regions within the genome. In certain embodiments, multiple target loci can be interrogated on a chromosome. This can be used to control for unanticipated abnormalities in regions containing target loci. It is possible that such regions can include a duplication, deletion, or other abnormality that may result in an inaccurate karyotype determination. By targeting multiple target loci of a chromosome, the user can control for such unanticipated abnormalities. Thus, in certain embodiments, one, more than one, at least two, or two or more loci are targeted from at least one chromosome or arm thereof, but optionally include more than one arm, to be interrogated. In certain embodiments, three, four, five, six, seven, eight, nine, 10 or more loci may be interrogated. In one embodiment, four loci are interrogated on each chromosome. In some embodiments, it may also be important to separate the target loci by a particular approximate number of nucleotides which can, in some embodiments, be calculated as is known in the art to further insure against interrogating a CNVR. It may also be beneficial to separate each target loci from any other target loci by approximately the same number of nucleotides (e.g., along the arm of a chromosome). Many suitable targets exist on each chromosome.
Exemplary targets (which can serve as qualified targets) and their corresponding PCR primers that can be used to interrogate human chromosomes are shown in Table 3. The probes that can be used in combination with particular primers or primer pairs, and the sequences amplified using particular primer pairs are shown in Table 4. One or more of such target loci can be amplified from a biological sample (e.g., gDNA). It should be understood that non-CNVR targets and corresponding probes or primers other than those listed in Tables 3 and 4 can be used in the methods disclosed herein. It is noted that chromosomes of other organisms may be similarly interrogated but that different target loci are used.
Any of the exemplary loci shown in Table 3, or any other suitable target loci that are outside of a CNVR, can be used in the karyotyping methods described herein. Such target loci can also be combined with target loci that fall within a CNVR, if desired. Particular loci can be targeted on each chromosome and any such target loci can be combined with any other target loci to form a subset, or screening panel. For instance, any of SEQ ID NOS. 1-192 can be used to interrogate chromosomes 1-22, X or Y as shown below: SEQ ID NOS. 1, 2, 3, 4, 5, 6, 7, 8 (chromosome 1); SEQ ID NOS. 9, 10, 11, 12, 13, 14, 15, 16 (chromosome 2); SEQ ID NOS.: 17, 18, 19, 20, 21, 22, 23, 24 (chromosome 3); SEQ ID NOS. 25, 26, 27, 28, 29, 30, 31, 32 (chromosome 4), 33, 34, 35, 36, 37, 38, 39, 40 (chromosome 5); SEQ ID NOS. 41, 42, 43, 44, 45, 46, 47, 48 (chromosome 6); SEQ ID NOS. 49, 50, 51, 52, 53, 54, 55, 56 (chromosome 7); SEQ ID NOS. 57, 58, 59, 60, 61, 62, 63, 64 (chromosome 8); SEQ ID NOS. 65, 66, 67, 68, 69, 70, 71, 72 (chromosome 9); SEQ ID NOS. 73, 74, 75, 76, 77, 78, 79, 80 (chromosome 10); SEQ ID NOS. 81, 82, 83, 84, 85, 86, 87, 88 (chromosome 11); SEQ ID NOS. 89, 90, 91, 92, 93, 94, 95, 96 (chromosome 12); SEQ ID NOS. 97, 98, 99, 100, 101, 102, 103, 104 (chromosome 13); SEQ ID NOS. 105, 106, 107, 108, 109, 110, 111, 112 (chromosome 14); SEQ ID NOS. 113, 114, 115, 116, 117, 118, 119, 120 (chromosome 15); SEQ ID NOS. 121, 122, 123, 124, 125, 126, 127, 128 (chromosome 16); SEQ ID NOS. 129, 130, 131, 132, 133, 134, 135, 136 (chromosome 17); SEQ ID NOS. 137, 138, 139, 140, 141, 142, 143, 144 (chromosome 18); SEQ ID NOS. 145, 146, 147, 148, 149, 150, 151, 152 (chromosome 19); SEQ ID NOS. 153, 154, 155, 156, 157, 158, 159, 160 (chromosome 20); SEQ ID NOS. 161, 162, 163, 164, 165, 166, 167, 168 (chromosome 21); SEQ ID NOS. 169, 170, 171, 172, 173, 174, 175, 176 (chromosome 22); SEQ ID NOS. 177, 178, 179, 180, 181, 182, 183, 184 (X chromosome); or SEQ ID NOS. 185, 186, 187, 188, 189, 190, 191, 192 (Y chromosome). Typically, the primers are used in pairs as set forth in Tables 3 and 4. It is understood that other primers and primer pairs not listed in Tables 3 and 4 can be used to interrogate chromosomes or parts of chromosomes 1-22, X and Y. The probes (e.g., SEQ ID NOS. 193-288) are typically used to detect their corresponding amplified sequence (e.g., SEQ ID NOS. 289-384) (Table 4). When multiple chromosomes are interrogated, any of SEQ ID NOS. 1-192 can be used to interrogate their respective chromosomes. Any one or more of SEQ ID NOS. 1-192 can be used with any other one or more of SEQ ID NOS. 1-192 to perform such assays. For instance, any one or more of SEQ ID NOS. 1-8 can be used with any one or more of SEQ ID NOS. 9-16 to interrogate chromosomes 1 and 2, respectively. Other similar subsets of SEQ ID NOS. 1-192 can also be suitable for use, as would be contemplated by one of skill in the art. At least one additional locus can also be interrogated as a control (e.g., “control gene”). Such additional loci can be, for example, the aforementioned RNase P (chromosome 14, cytoband 14q11.2) or TERT (chromosome 5, cytoband 5p15.33).
In certain embodiments, the karyotyping methods described herein are used to simultaneously interrogate multiple chromosomes. For example, using a human gDNA sample, more than one chromosome (i.e., chromosomes 1-22, X and Y) can be simultaneously interrogated in the same or separate assays. Furthermore, a subset of the total complement of chromosomes can be interrogated in the same or separate assays. Individual chromosomes can be targeted by amplifying loci specific to such chromosomes that are outside of CNVRs. For example, chromosome 1 can be interrogated along with any one or more of chromosomes 2-22, X or Y; chromosome 2 can be interrogated along with any one or more of chromosomes 1 and 3-22, X or Y; chromosome 3 can be interrogated along with any one or more of chromosomes 1-2, 3-22, X or Y; chromosome 4 can be interrogated along with any one or more of chromosomes 1-3, 5-22, X or Y; chromosome 5 can be interrogated along with any one or more of chromosomes 1-4, 6-22, X or Y; chromosome 6 can be interrogated along with any one or more of chromosomes 1-5, 7-22, X or Y; chromosome 7 can be interrogated along with any one or more of chromosomes 1-6, 8-22, X or Y; chromosome 8 can be interrogated along with any one or more of chromosomes 1-7, 9-22, X or Y; chromosome 9 can be interrogated along with any one or more of chromosomes 1-8, 10-22, X or Y; chromosome 10 can be interrogated along with any one or more of chromosomes 1-9, 11-22, X or Y; chromosome 11 can be interrogated along with any one or more of chromosomes 1-10, 12-22, X or Y; chromosome 12 can be interrogated along with any one or more of chromosomes 1-11, 13-22, X or Y; chromosome 13 can be interrogated along with any one or more of chromosomes 1-12, 14-22, X or Y; chromosome 14 can be interrogated along with any one or more of chromosomes 1-13, 15-22, X or Y; chromosome 15 can be interrogated along with any one or more of chromosomes 1-14, 16-22, X or Y; chromosome 16 can be interrogated along with any one or more of chromosomes 1-15, 17-22, X or Y; chromosome 17 can be interrogated along with any one or more of chromosomes 1-16, 18-22, X or Y; chromosome 18 can be interrogated along with any one or more of chromosomes 1-17, 19-22, X or Y; chromosome 19 can be interrogated along with any one or more of chromosomes 1-18, 20-22, X or Y; chromosome 20 can be interrogated along with any one or more of chromosomes 1-19, 21, 22, X or Y; chromosome 21 can be interrogated along with any one or more of chromosomes 1-20, 22, X or Y; chromosome 22 can be interrogated along with any one or more of chromosomes 1-21, X or Y; X chromosome can be interrogated along with any one or more of chromosomes 1-22 or Y; Y chromosome can be interrogated along with any one or more of chromosomes 1-22 or X. Other subsets of chromosomes can also be interrogated as would be understood by one of skill in the art. For PCR-based assays, any of the primer sets shown in Tables 3 or 4 can be used as may any other suitable primer set.
Described herein are karyotyping methods for determining the copy number of one or more chromosomes in a test biological sample (such as a test gDNA). In certain embodiments, the methods can be used to determine the karyotype of a test genome. The karyotype information typically relates to one or more chromosomes of a genome.
Such methods described herein can be used in prenatal diagnostic assays to screen for and detect chromosomal abnormalities in a fetus or embryo, such as Down's Syndrome (Trisomy 21), Edward's Syndrome (Trisomy 18) and Patau Syndrome (Trisomy 13). These assays can be performed on biological samples such as fetal cells or cell-free fetal DNA in maternal blood, amniotic fluid, trophoblast cells, chorionic villus samples and percutaneous umbilical cord blood. Furthermore, the methods disclosed herein can also be performed on embryos used in in vitro fertilization (IVF) wherein blastocyst cells or cells biopsied from the trophectoderm are analyzed from embryos that are 3 to 6 days post-fertilization, in order to determine if they should be implanted in utero.
The methods described herein can be useful in, for example, stem cell analysis, quality control assays of cell cultures, analysis of samples of limited quantity (e.g. single cell, formalin-fixed paraffin-embedded (FFPE)), comparisons between cell or tissue types, comparisons between diseased or non-diseased tissues, detection of chromosomal polymorphisms, and the like. The methods disclosed herein can be used to compare chromosome copy number between, for example, cell types (e.g., lymphocyte, epithelial cell), tissue types (e.g., neural tissue, skeletal muscle tissue), disease states (e.g., cancerous, non-cancerous), or types of organisms (e.g. human, mouse, plant, fruit fly). Other uses for the methods described herein will be apparent to one of skill in the art.
The methods described herein can be used to detect chromosomal abnormalities. Exemplary abnormalities include, but are not limited to, deletion, duplication, translocation, mosaicism, aneuploidy (e.g., nullisomy, disomy, trisomy and tetrasomy), inversion, ring formation, isochromosome formation, or chromosomal instability syndromes. Pathologies have been associated with abnormalities in particular human chromosomes including, for example, chromosome 1 (e.g., acute lymphoblastic leukemia, acute megakaryoblastic leukemia, alveolar rhabdomyosarcoma), chromosome 2 (e.g., anaplastic large cell lymphoma, alveolar rhabdomyosarcoma), chromosome 3, chromosome 4 (e.g., Wolf-Hirschhorn syndrome, acute lymphoblastic leukemia), chromosome 5 (e.g., Cri du chat, also known as Chromosome 5q deletion syndrome, anaplastic large cell lymphoma), chromosome 6, chromosome 7 (e.g., Williams syndrome), chromosome 8 (e.g., Burkitt's lymphoma, acute lymphoblastic leukemia, acute myeloblastic leukemia with maturation), chromosome 9 (e.g., Trisomy 9, Warkany syndrome 2, acute lymphoblastic leukemia, Philadelphia chromosome), chromosome 10, chromosome 11 (e.g., Jacobsen syndrome, mantle cell lymphoma, multiple myeloma, acute lymphoblastic leukemia, Ewing's sarcoma, desmoplastic small round cell tumor), chromosome 12 (e.g., acute lymphoblastic leukemia, myxoid liposarcoma), chromosome 13 (e.g., Patau syndrome, alveolar rhabdomyosarcoma, breast and ovarian cancers, deafness, Wilson's disease), chromosome 14 (e.g., Burkitt's lymphoma, follicular lymphoma, acute lymphoblastic leukemia), chromosome 15 (e.g., Angelman syndrome, Prader-Willi syndrome, acute promyelocytic leukemia, Marfan syndrome, Tay-Sach's disease), chromosome 16 (e.g., Trisomy 16, myxoid liposarcoma, polycystic kidney disease, α-thalassemia), chromosome 17 (e.g., Miller-Dieker syndrome, Smith-Magenis syndrome, acute promyelocytic leukemia, dermatofibrosarcoma protuberans, Charcot-Marie-Tooth disease), chromosome 18 (e.g., Edwards syndrome, Burkitt's lymphoma, follicular lymphoma, mantle cell lymphoma, multiple myeloma, synovial sarcoma, Nieman-Pick disease, pancreatic cancer), chromosome 19 (e.g., acute lymphoblastic leukemia), chromosome 20, chromosome 21 (e.g., Down syndrome, acute lymphoblastic leukemia, acute myeloblastic leukemia with maturation), chromosome 22 (e.g., Di George's syndrome, Trisomy 22, acute lymphoblastic leukemia, Philadelphia chromosome, acute megakaryoblastic leukemia, Ewing's sarcoma, dermatofibrosarcoma protuberans, desmoplastic small round cell tumor), X chromosome (e.g., Fragile X syndrome, Turner syndrome, Triple X syndrome, Klinefelter's syndrome, synovial sarcoma, mixed gonadal dysgenesis, XX gonadal dysgenesis, uniparental disomy, Duchenne muscular dystrophy, X-linked diseases, hemophilia, adrenoleukodystrophy, Hunter's disease), or Y chromosome (e.g., Klinefelter's syndrome, uniparental disomy, acute myeloid leukemia). Any of these disorders, or any other disorders that may be detected by karyotyping as would be known by the skilled artisan, can be studied using the methods described herein.
Kits for determining the copy number of multiple chromosomes are also provided. For use with a PCR-based assay, the kit can include at least a set of primers for amplifying at least one test locus on a chromosome (e.g., any one or more of SEQ ID NOS. 1-192) or corresponding probes (e.g., any one or more of SEQ ID NOS. 193-288). The kit can also include samples of pre-determined amplified sequences such as those listed in Table 4 (e.g., any one or more of SEQ ID NOS. 289-384), which can optionally be affixed to a solid support (e.g. a membrane, glass slide, multi-well plate) for use as a control reaction. The kit can also optionally include stock solutions, buffers, enzymes, detectable labels or reagents required for detection, tubes, membranes, and the like that can be used to complete the amplification reaction. In some embodiments, multiple primer sets are included. The kits can also comprise the reagents required to use a reference (e.g., primers and probes for amplifying or detecting a reference sequence such as RNase P or TERT) or calibrator sample (i.e., a control genome of known karyotype, a chromosome sample, or the like). In some embodiments, the kits can contain multiple sets of primers for amplifying multiple target loci on a single chromosome or among two or more chromosomes (e.g., any combination of primers sets and probes shown in Tables 3 and 4, or any primer/probe sets that target loci outside of CNVRs). In one embodiment, a multi-well plate contains within its wells an array of replicate wells for use in carrying out assays. For example, a 384-well array containing four replicate wells for each of chromosomes 1-22, X and Y comprising primer/probe sets targeting loci outside of CNVRs for each of the chromosomes, is provided, optionally including reagents such as stock solutions, buffers, enzymes such as polymerases, and detectable labels. Other embodiments of particular systems and kits are also contemplated which would be understood by one of skill in the art. In addition, the arrays can be customized depending on the target samples being analyzed, for example, arrays comprising one, two, three or more chromosomes comprising one, two, or more than two loci to be analyzed per chromosome. Furthermore, the number of replicate wells or samples per chromosome per array can vary (e.g., one, two, three, four, or more than four).
All references cited within this disclosure are hereby incorporated by reference in their entirety. Certain embodiments are further described in the following examples. These embodiments are provided as examples only and are not intended to limit the scope of the claims in any way.
The primers and probes used to interrogate the targets of interest described below were designed using masked genomic DNA sequences (the human genome assembly B36) by an Applied Biosytems proprietary TaqMan® Copy Number Assay Design Pipeline. The existing pre-designed TaqMan® copy number assay protocols were utilized. TaqMan® copy number assays targeting chromosomal regions located outside of copy number variable regions (CNVR) on each of the 24 chromosomes were designed (the data shown below, however, relates only to chromosomes 2, 3, 4, 5, 6, 20 and 22.) As described above, to cover all 24 chromosomes on a 384 well plate or TaqMan® array card (TLDA), four assays can be selected for each chromosome, using two assays on each arm of each chromosome whenever possible. The assays were run as duplex TaqMan® real-time PCR reactions. A FAM™ dye-based assay was used to detect test target loci and the VIC® dye-based assay was used for the reference locus RNase P (PN 4316844 from Applied Biosystems). Each assay utilized 10 ng gDNA, 1×TaqMan® probe/primer mix (900 nM of each primer, 250 nM of each probe) in 1×TaqMan® Genotyping Master Mix in a 10 μl reaction (run in quadruplicate). PCR reactions were incubated in an Applied Biosystems 7900HT SDS instrument for 10 minutes at 95° C., followed by 40 cycles of 15 seconds at 95° C. and 60 seconds at 60° C. Real-time data was collected by the SDS 2.3 software. The relative quantification analysis was performed to calculate estimated copy number of each sample for the gene of interest by CopyCaller™ software.
These methods were used to interrogate loci on chromosomes 2, 3, 4, 5, 6, 20 and 22 of human gDNA samples (a set of 92 African American and Caucasian gDNA samples from Coriell Institute) using the following primer sets and probes:
The resultant data is shown in
Many forms of cancer have been associated with various chromosomal abnormalities (see discussion hereinabove). Detection of diseases or susceptibility to diseases is critical for proper diagnosis or treatment. For example, certain breast and ovarian cancers are associated with abnormalities in chromosome 13. Primers and probes useful for interrogation of non-CNVR loci of chromosome 13 can be used to determine whether a patient has or is susceptible to having a breast or ovarian cancer.
The primers and probes used to interrogate the chromosome 13 target loci are designed using masked genomic DNA sequences (the human genome assembly B36) by an Applied Biosystems proprietary TaqMan® Copy Number Assay Design Pipeline and interrogated using copy number assays of format 1 and using virtual calibrator assay no. 1. Test samples are obtained from normal and tumor cell and tissue. As described above, four assays are selected for chromosome 13, using two assays on each arm of the chromosome. The assays are run as duplex TaqMan® real-time PCR reactions. A FAM™ dye-based assay is used to detect test target loci and the VIC® dye-based assay is used for the reference locus RNase P (PN 4316844 from Applied Biosystems) in assay format 1. Each assay uses 10 ng gDNA, 1×TaqMan® probe/primer mix (900 nM of each primer, 250 nM of each probe) in 1×TaqMan® Genotyping Master Mix in a 10 μl reaction (run in quadruplicate). PCR reactions are incubated in an Applied Biosystems 7900HT SDS instrument for 10 minutes at 95° C., followed by 40 cycles of 15 seconds at 95° C. and 60 seconds at 60° C. Real-time data is collected by the SDS 2.3 software. The relative quantification analysis is performed as described herein above to calculate estimated copy number of each sample for the gene of interest by CopyCaller™ software.
Identification of abnormal chromosomal copy number and karyotype of chromosome 13 indicates the presence or susceptibility to the breast or ovarian cancer associated with chromosome 13.
While certain embodiments have been described in terms of the preferred embodiments, it is understood that variations and modifications will occur to those skilled in the art. Therefore, it is intended that the appended claims cover all such equivalent variations that come within the scope of the following claims.
This application claims the benefit of priority under 35 U.S.C. §119(e) to U.S. Provisional Patent Application No. 61/334,658, filed May 14, 2010, which is herein incorporated by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
61334658 | May 2010 | US |