Diagnosing fetal chromosomal aneuploidy using genomic sequencing

Abstract

To provide methods, systems, and apparatuses for determining whether a fetal chromosomal aneuploidy exists in a biological sample obtained from a pregnant female.SOLUTION: A method for determining whether a fetal chromosomal aneuploidy exists comprises: sequencing nucleic acid molecules in a biological sample obtained from a pregnant female subject; determining respective amounts of a clinically relevant chromosome and of background chromosomes on the basis of the sequencing; comparing a parameter calculated from these amounts to one or more cutoff values; and, on the basis of the comparison, determining a classification of whether a fetal chromosomal aneuploidy exists.SELECTED DRAWING: Figure 1

Description

[Technical field]

[0001] Claiming priority This application claims priority from U.S. Provisional Application No. 60 / 951,438, entitled "DETERMINATION OF NUCLEIC ACID SEQUENCE IMBALANCES," filed July 23, 2007 (Attorney Docket No. 016285-005200US), the entire contents of which are incorporated herein by reference for all purposes.

[0002] CROSS-REFERENCE TO RELATED APPLICATIONS This application is also related to concurrently filed regular application (Attorney Docket No. 016285-005210US), entitled "Determining Nucleic Acid Sequence Imbalances," the entire contents of which are incorporated herein by reference for all purposes.

[0003] TECHNICAL FIELD OF THEINVENTION The present invention relates generally to diagnostic testing for fetal chromosomal aneuploidies by determining imbalances between different nucleic acid sequences, and more particularly to identifying Trisomy 21 (Down's Syndrome) and other chromosomal aneuploidies by testing maternal samples (e.g., blood). [Background technology]

[0004] background Fetal chromosomal aneuploidy results from the presence of an abnormal amount (or amounts) of chromosomes or chromosomal regions. The abnormal amount (or amounts) can be abnormally large, such as the presence of an extra chromosome 21, or the presence of a chromosomal region in trisomy 21; or abnormally small, such as the absence of a copy of the X chromosome in Turner syndrome.

[0005] Conventional methods for prenatal diagnosis of fetal chromosomal aneuploidies, e.g., trisomy 21, involve the collection of fetal samples by invasive procedures, e.g., amniocentesis or chorionic villus sampling, which carry the risk of fetal death. Non-invasive procedures, e.g., ultrasound and biochemical marker screening, have been used in risk strata of pregnant women prior to the final invasive diagnostic procedure. However, these screening methods usually measure epiphenomena associated with chromosomal aneuploidies, e.g., trisomy 21, instead of core chromosomal abnormalities, and therefore have suboptimal diagnostic accuracy and are significantly affected by other disadvantages, e.g., gestational age.

[0006] The discovery of circulating cell-free fetal DNA in maternal plasma in 1997 offered new possibilities for non-invasive diagnosis (Lo, YMD and Chiu, RWK 2007 Nat Rev Genet 8, 71-77). This method was quickly applied to the prenatal diagnosis of sex-linked genetic disorders (Costa, JM et al. 2002 N Engl J Med 346, 1502) and certain single-gene disorders (Lo, YMD et al. 1998 N Engl J Med 339, 1734-1738), but its application to the prenatal detection of fetal chromosomal aneuploidies presented considerable challenges (Lo, YMD and Chiu, RWK 2007, see above). First, fetal nucleic acids are usually present in maternal plasma with a high background of maternal nucleic acids that can interfere with the analysis of fetal nucleic acids (Lo, YMD et al. 1998 Am J Hum Genet 62, 768-775). Second, fetal nucleic acids circulate in maternal plasma mainly in acellular form, which makes it difficult to obtain information on gene dosage or chromosome dosage in the fetal genome.

[0007] Significant progress has been made in recent years to overcome these challenges (Benachi, A & Costa, JM 2007 Lancet 369, 440-442). One method is to detect fetal-specific nucleic acids in maternal plasma, thus overcoming the problem of maternal background disturbances (Lo, YMD and Chiu, RWK 2007, see above). The dosage of chromosome 21 was estimated from the ratio of polymorphic alleles in DNA / RNA molecules from the placenta. However, this method is less accurate when the sample contains a smaller amount of target nucleic acid and can only be applied to fetuses that are heterozygous for the targeted polymorphism, and only a portion of the population if one polymorphism is used.

[0008] Dhallan et al (Dhallan, R, et al. 2007, supra; Dhallan, R, et al. 2007 Lancet 369, 474-481) describe an alternative method to increase the proportion of circulating fetal DNA by adding formaldehyde to maternal plasma. The proportion of fetal chromosome 21 sequences in maternal plasma was determined by assessing the ratio of paternally inherited fetal-specific to non-fetal specific alleles for single nucleotide polymorphisms (SNPs) on chromosome 21. SNP ratios were calculated similarly for the reference chromosome. Fetal chromosome 21 imbalances were then inferred by detecting statistically significant differences between the SNP ratios of chromosome 21 and the reference chromosome, where significance is defined using a fixed p-value of 0.05 or less. To ensure high population coverage, more than 500 SNPs per chromosome were targeted. However, there is debate about the effectiveness of formaldehyde to increase fetal DNA to a high ratio (Chung, GTY, et al. 2005 Clin Chem 51, 655-658), and therefore the reproducibility of the method needs to be further evaluated. Similarly, since fetuses and mothers each provide information on a different number of SNPs for each chromosome, the statistical power of comparing SNP ratios will vary from case to case (Lo, YMD & Chiu, RWK. 2007 Lancet 369, 1997). Furthermore, since these methods rely on the detection of gene polymorphisms, they are limited to fetuses that are heterozygous for these polymorphisms.

[0009] Using polymerase chain reaction (PCR) and DNA quantification of chromosome 21 and reference loci in amniotic cell cultures derived from trisomy 21 and euploid fetuses, Zimmermann et al. (2002 Clin Chem 48,362-363) were able to distinguish between the two groups of fetuses based on a 1.5-fold increase in the former chromosome 21 DNA sequence. Discrimination of a 1.5-fold difference was the limit of conventional real-time PCR, since a 2-fold difference in DNA template concentration constitutes only one threshold cycle (Ct) difference. Alternative methods are needed to achieve a finer degree of quantitative discrimination.

[0010] Digital PCR was developed for the detection of skewed allele ratios in nucleic acid samples (Chang, HW et al. 2002 J Natl Cancer Inst 94, 1697-1703). Digital PCR is an amplification-based nucleic acid analysis technique that requires the partitioning of a nucleic acid-containing specimen into many individual samples that contain, on average, about one or less target sequence per sample. Specific nucleic acid targets are amplified by digital PCR using sequence-specific primer species or groups to generate specific amplicons. The nucleic acid loci to be targeted and the sequence-specific primer species or groups to be included in the reaction are determined or selected prior to nucleic acid analysis.

[0011] Clinically, it has been shown to be useful in detecting loss of heterozygosity (LOH) in tumor DNA samples (Zhou,W.et al.2002 Lancet 359,219-225).A previous study employed sequential probability ratio test (SPRT) for analyzing the results of digital PCR to classify whether the experimental results indicate the presence of LOH in the sample (El Karoui at al.2006 Stat Med 25,3124-3133).

[0012] In the methods used in previous studies, the amount of data collected from digital PCR is very low, and therefore precision may be compromised by the small number of data points and typical statistical fluctuations.

[0013] Therefore, it is desirable that non-invasive test has high sensitivity and specificity to minimize false negative and false positive respectively.However, fetal DNA exists in low absolute concentration and represents a small portion of the total DNA sequence in maternal plasma and serum.Therefore, it is also desirable to have a method that can perform non-invasive detection of fetal chromosomal aneuploidy by maximizing the amount of genetic information that can be inferred from the limited amount of fetal nucleic acid that exists as a small group in biological samples that contain maternal background nucleic acid. Summary of the Invention [Means for solving the problem]

[0014] overview The present invention provides a method, system, and device for determining whether a nucleic acid sequence imbalance (e.g., chromosomal imbalance) exists in a biological sample obtained from a pregnant woman. This determination can be made by using a parameter of the amount of a clinically significant chromosomal region in relation to other clinically insignificant chromosomal regions (background regions) in the biological sample. In one embodiment, the amount of a chromosome is determined from sequencing a nucleic acid molecule in a maternal sample, such as urine, plasma, serum, and other suitable biological samples. The nucleic acid molecules in the biological sample are sequenced, resulting in sequencing of genomic fragments. One or more cutoff values ​​are selected to determine whether a change (i.e., an imbalance) exists compared to a reference amount, for example, with respect to the ratio of the amount of two chromosomal regions (or sets of regions).

[0015] According to one exemplary embodiment, a biological sample obtained from a pregnant woman is analyzed to perform prenatal diagnosis of fetal chromosomal aneuploidy. The biological sample contains nucleic acid molecules. A portion of the nucleic acid molecules contained in the biological sample is sequenced. In one embodiment, the amount of genetic information obtained is sufficient for accurate diagnosis without being excessive overall, including the cost and amount of input biological sample required.

[0016] Based on the sequencing, a first amount of a first chromosome is determined from a sequence identified as originating from the first chromosome. A second amount of one or more second chromosomes is determined from a sequence identified as originating from one of the second chromosomes. A parameter from the first amount and the second amount is then compared to one or more cutoff values. Based on the comparison, a classification of whether or not a fetal chromosomal aneuploidy is present in the first chromosome is determined. The sequencing significantly maximizes the amount of genetic information that can be inferred from the limited amount of fetal nucleic acid present as a minority group in a biological sample containing maternal background nucleic acid.

[0017] According to one exemplary embodiment, a biological sample obtained from a pregnant woman is analyzed for prenatal diagnosis of fetal chromosomal aneuploidy. The biological sample comprises nucleic acid molecules. A percentage of fetal DNA in the biological sample is identified. The number N of sequences to be analyzed based on a desired accuracy is calculated based on the percentage. At least N of the nucleic acid molecules contained in the biological sample are randomly sequenced.

[0018] Based on the random sequencing, a first amount of a first chromosome is determined from a sequence identified as originating from the first chromosome. A second amount of one or more second chromosomes is determined from a sequence identified as originating from one of the second chromosomes. A parameter from the first amount and the second amount is then compared to one or more cutoff values. Based on the comparison, a classification of whether or not a fetal chromosomal aneuploidy is present in the first chromosome is determined. The random sequencing significantly maximizes the amount of genetic information that can be inferred from the limited amount of fetal nucleic acid present as a minority group in a biological sample containing background nucleic acid from the mother.

[0019] Other embodiments of the invention relate to systems and computer readable media associated with the methods described herein.

[0020] A further understanding of the nature and advantages of the present invention will be realized by reference to the following detailed description and the accompanying drawings. [Brief description of the drawings]

[0021] [Figure 1] FIG. 1 is a flow chart of a method 100 for performing prenatal diagnosis of fetal chromosomal aneuploidies in a biological sample obtained from a pregnant female subject, according to an embodiment of the present invention. [Diagram 2] FIG. 2 is a flow chart of a method 200 for performing prenatal diagnosis of fetal chromosomal aneuploidies using random sequencing, according to an embodiment of the present invention. [Figure 3A] FIG. 3A shows a plot of the percentage representation of chromosome 21 sequences in maternal plasma samples associated with trisomy 21 or euploid fetuses according to an embodiment of the invention. [Figure 3B] FIG. 3B shows the correlation between fractional fetal DNA concentration in maternal plasma as determined by massively parallel sequencing and microfluidics digital PCR according to an embodiment of the present invention. [Figure 4A]FIG. 4A shows a plot of the percentage of aligned sequences per chromosome according to an embodiment of the present invention. [Figure 4B] FIG. 4B shows a plot of the % difference in percentage expression per chromosome between the trisomy 21 and euploid cases shown in FIG. 4A. [Diagram 5] FIG. 5 shows the correlation between the degree of over-representation of chromosome 21 sequences in maternal plasma associated with trisomy 21 fetuses and partial fetal DNA concentration, in accordance with an embodiment of the present invention. [Figure 6] 6 shows a table of portions of the human genome that have been analyzed according to an embodiment of the present invention, where T21 indicates a sample obtained from a pregnant woman carrying a trisomy 21 fetus. [Figure 7] FIG. 7 shows a table of the number of sequences required to distinguish between euploid and trisomy 21 fetuses in accordance with an embodiment of the present invention. [Figure 8A] FIG. 8A shows a table of the top 10 start positions of sequenced tags aligned to chromosome 21 according to an embodiment of the present invention. [Figure 8B] FIG. 8B shows a table of the top 10 start positions of sequenced tags aligned to chromosome 22 according to an embodiment of the present invention. [Figure 9] FIG. 9 illustrates a block diagram of an exemplary computing device that can be used with the systems and methods of embodiments of the present invention. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0022] definition As used herein, the term "biological sample" refers to any sample taken from a subject (e.g., a human, such as a pregnant female) and contains one or more nucleic acid molecule(s) of interest.

[0023] The term "nucleic acid" or "polynucleotide" refers to deoxyribonucleic acid (DNA) or ribonucleic acid (RNA) and polymers thereof in either single- or double-stranded form. Unless specifically limited, the term includes nucleic acids containing known natural nucleotide analogs that have similar binding properties as the reference nucleic acid and are metabolized in a manner similar to naturally occurring nucleotides. Unless otherwise indicated, a particular nucleic acid sequence also implicitly includes its conservatively modified variants (e.g., degenerate codon substitutions), alleles, orthologs, SNPs, and complementary sequences, as well as the sequence explicitly indicated. Specifically, degenerate codon substitutions can be achieved by generating sequences in which the third position of one or more selected (or all) codons is substituted with mixed-base and / or deoxyinosine residues (Batzer et al., Nucleic Acid Res. 19:5081 (1991); Ohtsuka et al, J. Biol. Chem. 260:2605-2608 (1985); and Rossolini et al, Mol. Cell. Probes 8:91-98 (1994)). The term nucleic acid is used interchangeably with gene, cDNA, mRNA, small non-coding RNA, microRNA (miRNA), Piwi-interacting RNA, and short hairpin RNA (shRNA) encoded by a gene or locus.

[0024] The term "gene" refers to a DNA segment involved in producing a polypeptide chain. It can include regions preceding and following the coding region (leader and trailer), as well as intervening sequences between individual coding segments (exons).

[0025] The term "reaction" as used herein refers to any process involving chemical, enzymatic, or physical action that indicates the presence or absence of a particular polynucleotide sequence of interest. An example of a "reaction" is an amplification reaction, such as the polymerase chain reaction (PCR). Another example of a "reaction" is a sequencing reaction, either by synthesis or ligation. An "informative reaction" is one that indicates the presence of one or more particular polynucleotide sequences of interest in the presence of only one sequence of interest. The term "well" as used herein refers to a reaction in a predetermined location within a closed structure, such as a well-shaped vial, cell, or chamber in a PCR array.

[0026] The term "clinically significant nucleic acid sequence" as used herein can refer to a polynucleotide sequence that corresponds to a larger gene sequence fragment whose potential imbalance is being tested, or to a polynucleotide sequence that corresponds to the larger gene sequence itself. One example is the sequence of chromosome 21. Other examples include chromosomes 18, 13, X, and Y. Still other examples include mutated gene sequences that a fetus may inherit from one or both of its parents, or genetic polymorphisms, or copy number variations. Still other examples include sequences that are mutated, deleted, or amplified in malignant tumors, such as sequences that undergo loss of heterozygosity or gene duplication. In some embodiments, multiple clinically significant nucleic acid sequences, or even multiple manufacturers of clinically significant nucleic acid sequences, can be used to provide data for detecting the imbalance. For example, data from five discontinuous sequences on chromosome 21 can be used in an additive method for determining possible chromosome 21 imbalances, effectively reducing the volume of the sample by 5-fold.

[0027] The term "background nucleic acid sequence" as used herein refers to a nucleic acid sequence whose normal ratio to a clinically significant nucleic acid sequence is known, for example, a 1:1 ratio.In one example, the background nucleic acid sequence and the clinically significant nucleic acid sequence are two alleles from the same chromosome that are separate due to heterozygosity.In another example, the background nucleic acid sequence is one allele that is heterozygous to another allele that is a clinically significant nucleic acid sequence.Furthermore, each of the background nucleic acid sequence and some of the clinically significant nucleic acid sequences may be from different individuals.

[0028] As used herein, the term "reference nucleic acid sequence" refers to a nucleic acid sequence whose average concentration per reaction is known, or equivalently measured.

[0029] As used herein, the term "overrepresented nucleic acid sequence" refers to a nucleic acid sequence in two sequences of interest (e.g., a clinically important sequence and a background sequence) that is more abundant than other sequences in a biological sample.

[0030] The term "based on" as used herein means "based at least in part on" and refers to a value (or result) used in the determination of one value that occurs in the relationship between a method's inputs and the output of that method. The term "derive" as used herein also refers to the relationship between a method's inputs and the output of that method, such as occurs when the derivation is a formula calculation.

[0031] As used herein, the term "quantitative data" refers to data obtained from one or more reactions and providing one or more numerical values. For example, the number of wells that display a fluorescent marker for a particular sequence would be quantitative data.

[0032] The term "parameter" as used herein refers to a numerical value that characterizes a set of quantitative data and / or a numerical relationship between quantitative data sets.For example, the ratio (or a function of the ratio) between a first amount of a first nucleic acid sequence and a second amount of a second nucleic acid sequence is a parameter.

[0033] As used herein, the term "cutoff value" refers to a numerical value used to mediate between two or more states of classification of a biological sample (e.g., diseased state and non-diseased state). For example, if the parameter is greater than the cutoff value, a first classification of the quantitative data is made (e.g., diseased state); or if the parameter is less than the cutoff value, a different classification of the quantitative data is made (e.g., non-diseased state).

[0034] The term "imbalance" as used herein refers to any significant deviation from a reference amount, defined by at least one cutoff value in the amount of a clinically significant nucleic acid sequence.For example, the reference amount may be a ratio of 3 / 5, and thus an imbalance would occur if the measured ratio were 1:1.

[0035] As used herein, the term "chromosomal aneuploidy" refers to a variation in chromosomal dosage from that of a diploid genome. The variation can be an increase or a deletion. It can involve an entire chromosome or a region of a chromosome.

[0036] The term "random sequencing" as used herein refers to sequencing where the sequenced nucleic acid fragments were not specifically identified or targeted prior to the sequencing procedure. No sequence-specific primers were required to target specific loci. The pool of sequenced nucleic acids varies from sample to sample, and even from analysis to analysis on the same sample. The identity of the sequenced nucleic acid is only revealed from the resulting sequencing results. In some embodiments of the present invention, random sequencing can be performed by a procedure that enriches biological samples with a specific group of nucleic acid molecules that share a common feature. In one embodiment, each fragment in a biological sample has an equal chance of being sequenced.

[0037] The term "human genome fragment" or "part of human genome" as used herein refers to less than 100% of the nucleotide sequence in the human genome, which contains 3 billion base pairs of nucleotides. In sequencing, it refers to less than 1-fold coverage of the nucleotide sequence in the human genome. The term can be expressed as a percentage or absolute number of nucleotides / base pairs. As an example of use, the term can be used to refer to the actual amount of sequencing that has been performed. An embodiment can determine the minimum amount of sequenced fragments of the human genome that are required to perform an accurate diagnosis. As another example of use, the term can refer to the amount of sequenced data that is used to derive a parameter or amount for disease classification.

[0038] The term "sequenced tag" as used herein refers to a nucleotide sequence that is sequenced from any part or all of a nucleic acid molecule.For example, a sequenced tag can be a short nucleotide sequence that is sequenced from a nucleic acid fragment, a short nucleotide sequence at both ends of a nucleic acid fragment, or the sequencing of the complete nucleic acid fragment present in a biological sample.A nucleic acid fragment is any part of a larger nucleic acid molecule.A fragment (e.g., a gene) can be present separately (i.e., unlinked) from other parts of a larger nucleic acid molecule.

[0039] Detailed Description The present invention provides a method, system, and device for determining whether there is an increase or decrease in a clinically significant chromosomal region (disease state) compared to a non-disease state. This determination can be made by using a parameter of the amount of the clinically significant chromosomal region in relation to other clinically insignificant chromosomal regions (background regions) in the biological sample. The nucleic acid molecules in the biological sample can be sequenced, resulting in sequencing of genomic fragments, and the amount can be determined from the results of the sequencing. For example, one or more cutoff values ​​are selected to determine whether there is an alteration (i.e., imbalance) in the ratio of the amounts of two chromosomal regions (or sets of regions) compared to a reference amount.

[0040] The change detected in the reference amount can be any deviation (upward or downward) in the relationship between the clinically significant nucleic acid sequence and other sequences that are not clinically significant. Thus, the reference state can be any ratio or other amount (e.g., other than a 1-1 correspondence), and the measured state indicating the change can be any ratio or other amount that differs from the reference amount as determined by one or more cutoff values.

[0041] The clinically significant chromosomal regions (also referred to as clinically significant nucleic acid sequences) and background nucleic acid sequences can be from a first type of cell and can be from one or more second type of cells. For example, fetal nucleic acid sequences from fetal / placental cells are present in a biological sample, such as maternal plasma, that contains a background of maternal nucleic acid sequences from maternal cells. In one embodiment, the cutoff value is determined at least in part based on the percentage of cells of the first type in the biological sample. It should be noted that the percentage of fetal sequences in the sample can be determined by any fetal-derived locus and is not limited to measuring clinically significant nucleic acid sequences. In another embodiment, the cutoff value is determined at least in part based on the percentage of tumor sequences in a biological sample, such as plasma, serum, saliva, or urine, that contains a background of nucleic acid sequences from non-malignant cells in the body.

[0042] I. General method FIG. 1 is a flow chart of a method 100 for performing prenatal diagnosis of fetal chromosomal aneuploidies in a biological sample obtained from a pregnant female subject, according to an embodiment of the present invention.

[0043] In step 110, a biological sample is obtained from the pregnant female. The biological sample may be plasma, urine, serum, or any other suitable sample. The sample includes nucleic acid molecules from the fetus and the pregnant female. For example, the nucleic acid molecules may be fragments of chromosomes.

[0044] In step 120, at least a portion of the plurality of nucleic acid molecules contained in the biological sample is sequenced. The sequenced portions represent human genome fragments. In one embodiment, the nucleic acid molecules are fragments of respective chromosomes. One end (e.g., 35 base pairs (bp)), both ends, or the complete fragment may be sequenced. All of the nucleic acid molecules in the sample may be sequenced, or only a subset may be sequenced. The subset may be randomly selected, as described in more detail below.

[0045] In one embodiment, the sequencing is performed using massively parallel sequencing. Massive sequencing, such as 454 platform (Roche) (Margulies, M. et al. 2005 Nature 437, 376-380), Illumina Genome Analyzer (or Solexa platform), or SOLiD system (Applied Biosystems), or Helicos True Single Molecule DNA sequencing technology (Harris TD et al. 2008 Science, 320, 106-109), Pacific Biosciences' single molecule, real-time (SMRT®) technology, and nanopore sequencing (Soni GV and Meller A. 2007 Clin Chem 53: 1996-2001), allows the sequencing of many nucleic acid molecules isolated from a sample in a highly multiplexed parallel manner (Dear Brief Funct Genomic Proteomic 2003; 1: 397-416). Each of these platforms sequences single molecules of nucleic acid fragments that have not been clonally expanded or even amplified.

[0046] Since a large number of sequencing reads are generated from each sample in each run, on the order of 100,000 to 1 million or even millions or billions, the resulting sequenced reads form a representative signature of the mixture of nucleic acid species in the original specimen. For example, the haplotype, transcriptome, and methylation profile of the sequenced reads are similar to those of the original specimen (Brenner et al Nat Biotech 2000;18:630-634; Taylor et al Cancer Res 2007;67:8511-8518). Due to the large sampling of sequences from each specimen, the number of identical sequences, such as those resulting from sequencing a nucleic acid pool with several-fold coverage or high redundancy, is also a good quantitative indication of the number of specific nucleic acid species or loci in the original sample.

[0047] In step 130, a first amount of a first chromosome (e.g., a clinically important chromosome) is determined based on sequencing (e.g., data from sequencing). The first amount is determined from sequences identified as originating from the first chromosome. For example, bioinformatics techniques can then be used to place each of these DNA sequences into the human genome. Such sequences may be present in repetitive regions of the human genome, or in regions subject to inter-individual variation, e.g., copy number variation, and therefore such sequence ratios may be excluded from subsequent analysis. The amount of the chromosome of interest and the amount of one or more other chromosomes may be determined in this manner.

[0048] In step 140, based on the sequencing, a second amount of one or more second chromosomes is determined from the sequence identified as coming from one of the second chromosomes. In one embodiment, the second chromosome is all other chromosomes other than the first chromosome (i.e., the chromosome being tested). In another embodiment, the second chromosome is just a single other chromosome.

[0049] There are many methods for determining chromosome quantity, including but not limited to the number of tags sequenced, the number of nucleotides (base pairs) sequenced, or a tally of the cumulative length of sequenced nucleotides (base pairs) from a particular chromosome(s) or chromosomal region.

[0050] In another embodiment, rules can be imposed on the results of sequencing to determine what is counted.In one embodiment, a certain amount can be obtained based on the ratio of the results of sequencing.For example, the results of sequencing corresponding to the nucleic acid fragments of a certain size range are selected after bioinformatics analysis.Examples of size ranges are less than about 300bp, less than 200bp, or less than 100bp.

[0051] In step 150, a parameter is determined from the first amount and the second amount. For example, the parameter can be a simple ratio of the first amount to the second amount, or a simple ratio of the first amount to the second amount and the first amount. In one embodiment, each amount can be an independent variable of one or another function, where a ratio can be obtained from these other functions. Those skilled in the art will recognize a number of different suitable parameters.

[0052] In one embodiment, a parameter (e.g., fractional representation) of a chromosome potentially involved in chromosome aneuploidy, such as chromosome 21, chromosome 18, or chromosome 13, can thus be calculated from the results of the bioinformatics procedure. The fractional representation can be obtained based on the abundance of the entire chromosome sequence (e.g., some measurement of all chromosomes, including the clinically important chromosomes), or of a specific chromosome subset (e.g., just one other than the one being tested).

[0053] In step 150, the parameters are compared to one or more cutoff values. The cutoff values ​​may be determined from a variety of suitable methods. Such methods include Bayesian likelihood methods, sequential probability ratio tests (SPRT), false discovery, confidence intervals, and receiver operating characteristics (ROC). Examples of these methods and the application of sample-specific methods are described in co-filed application "DETERMINING A NUCLEIC ACID SEQUENCE IMBALANCE" (Attorney Docket No. 016285-005210US), which is incorporated by reference.

[0054] In one embodiment, the parameters (e.g., fractional representation of clinically significant chromosomes) are then compared to reference ranges established in pregnancies involving normal (i.e., euploid) fetuses. In some variations of the procedure, the reference ranges (i.e., cutoff values) may be adjusted according to the fractional concentration (f) of fetal DNA in a particular maternal plasma sample. The value of f may be determined from a sequencing dataset, for example, using sequences mappable to the Y chromosome if the fetus is male. The value of f may also be determined from separate analyses, for example, using fetal epigenetic markers (Chan KCA et al 2006 Clin Chem 52,2211-8), or from analysis of single nucleotide polymorphisms.

[0055] At step 160, a classification is determined of whether or not fetal chromosomal aneuploidy is present at the first chromosome based on the comparison. In one embodiment, the classification is a definite yes or no. In another embodiment, the classification may be non-deterministic or unclear. In yet another embodiment, the classification may be a score that can be interpreted at a later date, for example, by a physician.

[0056] II. Sequencing, Alignment, and Quantitation As mentioned above, only genome fragments are sequenced.In one embodiment, even when the pool of nucleic acids in a sample is sequenced with less than 100% gene coverage instead of several-fold coverage, in the proportion of captured nucleic acid molecules, most of each nucleic acid species is sequenced only once.Similarly, the imbalance of the amount of a particular chromosome or chromosome region can be quantitatively determined.That is, the imbalance of the amount of a chromosome or chromosome region is inferred from the percentage representation of said locus in other sequenced tags that can be mapped in the sample.

[0057] This is in contrast to situations where the same pool of nucleic acids are sequenced multiple times to achieve large excess, or several-fold coverage, such that individual nucleic acid species are sequenced multiple times. In such situations, the number of times a particular nucleic acid species is sequenced correlates with their relative concentration in the original sample compared to that of another nucleic acid species. The cost of sequencing increases with the number of fold coverages required to achieve accurate representation of the nucleic acid species.

[0058] In one embodiment, the ratio of such sequences will be from the chromosome involved in the aneuploidy, for example chromosome 21 in this example. Further sequences from such sequencing runs will be from other chromosomes. By taking into account the relative size of chromosome 21 compared to other chromosomes, a normalized frequency of chromosome 21-specific sequences from such sequencing runs within the reference range can be obtained. If the fetus has trisomy 21, the normalized frequency of chromosome 21-derived sequences from such sequencing runs will increase, thus allowing the detection of trisomy 21. The extent of the change in normalized frequency will depend on the partial concentration of fetal nucleic acids in the analyzed sample.

[0059] In one embodiment, we used an Illumina Genome Analyzer for single-end sequencing of human genome and human plasma DNA samples. The Illumina Genome Analyzer sequences clonally expanded single DNA molecules captured on a solid surface called a flow cell. Each flow cell has eight lanes for sequencing eight individual samples or pools of samples. Each lane can yield approximately 200 Mb of sequence, which is only a portion of the 3 billion base pairs of sequence in the human genome. Each genomic DNA or plasma DNA sample was sequenced using one lane of a flow cell. The resulting short sequence tags were aligned to the human reference genome sequence, and the origins of replication of the chromosomes were indicated. The total number of individual sequenced tags aligned to each chromosome was tabulated and compared to the relative size of each chromosome expected from the reference human genome or a non-disease representative sample. Chromosome gains or losses were then identified.

[0060] The described approach is only one example of the gene / chromosome dosage method described herein. Alternatively, paired end sequencing can be performed. Instead of comparing the length of the sequenced fragments with that expected in the reference genome described by Campbell et al. (Nat Genet 2008;40:722-729), the aligned and sequenced tags are counted and sorted according to chromosomal location. The gain or loss of chromosomal regions, or entire chromosomes, was determined by comparing the number of tags with the expected chromosomal size in the reference genome or that of a representative non-disease specimen. Since paired end sequencing allows for an estimate of the size of the original nucleic acid fragment, one example is to focus on counting paired sequenced tags that correspond to nucleic acid fragments of a certain size, for example, less than 300bp, less than 200bp, or less than 100bp.

[0061] In another embodiment, the nucleic acid pool fractions to be sequenced in one run are further selected before sequencing. For example, hybridization-type techniques, such as oligonucleotide arrays, can be used to first select nucleic acid sequences from specific chromosomes, such as potentially aneuploid chromosomes and other chromosomes that are not involved in the tested aneuploidy. Another example is that a specific subgroup of nucleic acid sequences from the sample pool is selected or enriched before sequencing. For example, as discussed above, it has been reported that fetal DNA molecules in maternal plasma contain shorter fragments than maternal background DNA molecules (Chan et al Clin Chem 2004;50:88-92). Therefore, one or more methods known to those skilled in the art can be used to fractionate nucleic acid sequences in a sample according to molecular size, for example by gel electrophoresis, or size exclusion column, or by microfluidics-type techniques. Further alternatively, in the example of analyzing cell-free fetal DNA in maternal plasma, the fetal nucleic acid portion can be amplified by methods that suppress maternal background, such as by the addition of formaldehyde (Dhallan et al JAMA 2004;291:1114-9). In one embodiment, a portion or subset of the pool of auxiliary selected nucleic acids is randomly sequenced.

[0062] Other single molecule sequencing methods, such as the Roche454 platform, the Applied Biosystems SOLiD platform, the Helicos True Single Molecule DNA sequencing technology, Pacific Biosciences' Single Molecule Real Time (SMRT®) technology, and nanopore sequencing methods may be used in this application as well.

[0063] III. Determining chromosome dosage from sequencing results After massively parallel sequencing, bioinformatics analysis was performed to locate the chromosomal origin of replication of the sequenced tags. After this procedure, the tags identified as potentially aneuploid chromosomes, i.e., chromosome 21 in this study, are quantitatively compared with all sequenced tags, or with tags from one or more chromosomes not associated with aneuploidy. The relationship between the sequencing results from chromosome 21 and the sequencing results of other non-chromosomes 21 for the test sample is compared with the cutoff value obtained from the method described in the above section to determine whether the sample is obtained from a pregnant woman involved in aneuploidy or trisomy 21 fetus.

[0064] Many different quantities can be obtained from the sequenced tags, including but not limited to the following. For example, the number, i.e., absolute number, of sequenced tags aligned to a particular chromosome can be compared to the absolute number of sequenced tags aligned to other chromosomes. Alternatively, the fractional count of the amount of sequenced tags from chromosome 21 relative to all or some other sequenced tags can be compared to that of other non-aneuploid chromosomes. In this experiment, 36bp was sequenced from each DNA fragment, so the number of nucleotides sequenced from a particular chromosome was easily obtained by multiplying 36bp by the number of tags sequenced.

[0065] Furthermore, because each maternal plasma specimen was sequenced statistically using one flow cell capable of sequencing human genome fragments, most maternal plasma DNA fragment species would each be sequenced to produce one sequenced tag number. That is, the nucleic acid fragments present in the maternal plasma specimen were sequenced to a range of less than 1x. Therefore, the total number of sequenced nucleotides for any particular chromosome will largely correspond to the amount, proportion or length of the portion of said chromosome sequenced. Thus, the quantitative determination of the representation of a chromosome that is potentially aneuploid is obtained from the portion of the number of nucleotides sequenced from that chromosome or a portion of equivalent length relative to the amount for other chromosomes similarly obtained.

[0066] IV. Increasing the Pool of Nucleic Acids for Sequencing As mentioned above and described in the following examples, only a portion of the human genome needs to be sequenced to distinguish between euploid cases and trisomy 21. Therefore, it would be possible and cost-effective to increase the pool of nucleic acids to be sequenced prior to random sequencing of a fraction of the increased pool. For example, fetal DNA molecules in maternal plasma contain shorter fragments than maternal background DNA molecules (Chan et al Clin Chem 2004;50:88-92). Therefore, one or more methods known to those skilled in the art can be used to fractionate nucleic acid sequences in a sample according to molecular size, for example by gel electrophoresis, or size exclusion columns, or by microfluidics-type techniques.

[0067] Furthermore, alternatively, in the example of analyzing cell-free fetal DNA in maternal plasma, the fetal nucleic acid portion can be increased by methods that suppress maternal background, such as by the addition of formaldehyde (Dhallan et al JAMA 2004;291:1114-9). The fetal sequence ratio will be increased in nucleic acid pools containing shorter fragments. According to Figure 7, the number of sequenced tags required to distinguish euploid from trisomy 21 will be reduced by increasing the partial fetal DNA concentration.

[0068] Alternatively, sequences from the potentially aneuploid chromosome and one or more chromosomes not involved in the aneuploidy can be enriched, for example by hybridization techniques on oligonucleotide microarrays. The enriched pool of nucleic acids would then be subjected to random screening. This could reduce the cost of sequencing.

[0069] V. Random Screening Figure 2 is a flow chart of a method 200 for performing prenatal diagnosis of fetal chromosomal aneuploidy using random screening of an embodiment of the present invention. In one embodiment for massively parallel sequencing, representative data from all chromosomes can be generated simultaneously. The origin of replication of a particular fragment is not selected in advance of time. The sequencing is done randomly, and then a database search can be performed to see where a particular fragment originates from. This is in contrast to the situation where a particular fragment from chromosome 21 and another fragment from chromosome 1 are amplified.

[0070] In step 210, a biological sample from a pregnant woman is obtained. In step 220, the N number of sequences to be analyzed is calculated with the desired accuracy. In one embodiment, the percentage of fetal DNA in the biological sample is first identified. This can be done by any suitable means known to those skilled in the art. The identification can simply read a value measured by another. In this embodiment, the calculation of the N number of sequences to be analyzed is based on a percentage. For example, the number of sequences required for analysis can be increased when the fetal DNA percentage falls and decreased when fetal DNA increases. The N number can be a fixed number or a relative number, such as a percentage. In another embodiment, an N number known to be sufficient for accurate disease diagnosis can be sequenced. The N number can be sufficient even in pregnant women with fetal DNA concentrations that are at the lower end of the normal range.

[0071] In step 230, at least N of the plurality of nucleic acid molecules contained in the biological sample are randomly sequenced. A feature of this described approach is that the nucleic acids to be sequenced are not specifically identified or targeted prior to sample analysis, i.e., sequencing. No sequence-specific primers targeting specific loci are required in the sequencing. The pool of nucleic acids to be sequenced varies from sample to sample, and even from analysis to analysis on the same sample. Furthermore, as will be seen from the following description (FIG. 6), the amount of sequencing results required in a diagnostic case may vary between the amount of specimen tested and the amount of the reference group. These aspects are in stark contrast to most molecular diagnostic approaches, such as those based on fluorescent in situ hybridization, quantitative fluorescent PCR, quantitative real-time PCR, digital PCR, comparative fluorescent PCR, macroarray comparative gene hybridization, etc., in which the targeted loci are required in advance prior to their determination, thus necessitating the use of locus-specific primers or probe sets or groups thereof.

[0072] In one embodiment, random screening is performed on DNA fragments present in the plasma of pregnant women, and obtains gene sequences that may originate from either the fetus or its mother. Random screening involves sampling (sequencing) a random portion of nucleic acid molecules present in a biological sample. Because the sequencing is random, a different subset (fraction) of nucleic acid molecules (and their genome) may be sequenced in each analysis. Embodiments may be performed when this subset varies from sample to sample, and even from analysis to analysis that may use the same sample. Examples of such fractions are about 0.1%, 0.5%, 1%, 5%, 10%, 20%, or 30% of the genome. In other embodiments, the fraction is at least any of these values.

[0073] The remaining steps 240 - 270 may proceed in a similar manner to method 100 .

[0074] VI. Post-sequencing selection of the pool of sequenced tags As described in the following Examples II and III, a subset of sequenced data is sufficient to distinguish between trisomy 21 and euploid cases. A subset of sequenced data can be the proportion of sequenced tags that pass a certain quality parameter. For example, in Example II, sequenced tags that are individually aligned to a repeat-masked reference human genome are used. Alternatively, a representative pool of nucleic acid fragments from all chromosomes can be sequenced, but with a focus on the comparison between data associated with potentially aneuploid chromosomes and data associated with many non-aneuploid chromosomes.

[0075] Alternatively, a subset of sequencing results including sequenced tags resulting from nucleic acid fragments corresponding to a particular size window in the original specimen can be selected during post-sequencing analysis. For example, paired-end sequencing can be used using an Illumina genome analyzer, which involves sequencing two ends of a nucleic acid fragment. The sequenced data from each paired end is then aligned to a reference human genome sequence. The distance or number of nucleotides spanning between the two ends is then estimated. The full length of the original nucleic acid fragment can also be estimated. Alternatively, sequencing platforms such as the 454 platform, and possibly some single molecule sequencing techniques, can sequence the full length of short nucleic acid sequences, for example 200bp. In this way, the actual length of the nucleic acid fragment will be known quickly from the sequenced data.

[0076] Such paired-end analysis is also possible using other sequencing platforms, such as the Applied Biosystems SOLiD system. For the Roche454 platform, due to the increased read length compared to other massively parallel sequencing systems, it is also possible to determine fragment lengths from the complete sequence.

[0077] The advantage of focusing data analysis on the subset of sequenced tags that represent short nucleic acid fragments in the original maternal plasma sample is that the data set is effectively enriched for fetal DNA sequences, since fetal DNA molecules in maternal plasma contain more short fragments than maternal background DNA molecules (Chan et al Clin Chem 2004;50:88-92). According to Figure 7, the number of sequenced tags required to distinguish between euploid and trisomy 21 cases will decrease due to the increased fractional fetal DNA concentration.

[0078] Post-sequencing selection of a subset of the nucleic acid pool differs from other nucleic acid enrichment methods performed prior to sample analysis, such as gel electrophoresis or the use of size exclusion columns to select nucleic acids of a particular size, which require physical separation of the enriched pool from the background pool of nucleic acids. Physical procedures introduce more experimental steps and can create problems such as contamination. Post-sequencing in silico selection of a subset of the sequencing results may also vary depending on the sensitivity and specificity required for disease determination.

[0079] The bioinformatics, computer and statistical methods used to determine whether a maternal plasma sample is obtained from a pregnant woman carrying a trisomy 21 or euploid fetus can be compiled into a computer program product used to determine parameters from the results of sequencing. The operation of the computer program will relate to the determination of a quantitative amount from the essentially aneuploid chromosome and an amount(s) from one or more other chromosomes. The parameters will be determined and compared using suitable cut-off values ​​to determine whether fetal chromosomal aneuploidy is present in the essentially aneuploid chromosome. EXAMPLES

[0080] The following examples are offered to illustrate, but not to limit, the claimed invention. I. Prenatal diagnosis of fetal trisomy 21 Eight pregnant women were recruited for the study. All pregnant women were in the first or second trimester and had singleton pregnancies. Four of them were carrying fetuses with trisomy 21, and the other four were carrying euploid fetuses. Twenty milliliters of peripheral venous blood were collected from each subject. Maternal plasma was collected after centrifugation at 1600xg for 10 min and at 16000xg for another 10 min. DNA was then extracted from 5-10 mL of each plasma sample. Maternal plasma DNA was then used for massively parallel sequencing by an Illumina genome analyzer, according to the manufacturer's instructions. The sequencing technicians were temporarily blinded to the fetal diagnosis during sequencing and sequence data analysis.

[0081] Briefly, approximately 50 ng of maternal plasma DNA was used for DNA library preparation. It is possible to start with smaller amounts, e.g., 15 ng or 10 ng of maternal plasma DNA. Maternal plasma DNA fragments were blunt-end ligated with Solexa adapters, and fragments of 150–300 bp were selected by gel purification. Alternatively, blunt-end and adapter-ligated maternal plasma DNA fragments can be passed through a column (e.g., AMPure, Agencourt) to remove unligated adapters without size selection prior to cluster generation. Adapter-ligated DNA was hybridized to the flow cell surface, and DNA clusters were generated using an Illumina cluster station, followed by 36 cycles of sequencing using an Illumina genome analyzer. DNA from each maternal plasma specimen was sequenced by one flow cell. Sequenced reads were compiled using the Solexa analysis pipeline. All reads were then aligned to the repeat-masked reference human genome sequence, the NCBI36 assembly (GenBank accession numbers: NC_000001 to NC_000024), using the Eland application.

[0082] In this study, to reduce the complexity of data analysis, only sequences that were located at unique positions in the repeat-masked human genome reference were further considered. Other subsets or complete sets of sequenced data could alternatively be used. The total number of individually mappable sequences for each specimen was counted. The number of sequences that were individually aligned to chromosome 21 was expressed as a ratio to the total number of aligned sequences for each specimen. Since maternal plasma contains fetal DNA in the background of DNA from its mother, a fetus with trisomy 21 would contribute extra sequenced tags derived from chromosome 21 due to the presence of extra copies of chromosome 21 in the fetal genome. Thus, the percentage of chromosome 21 sequences in maternal plasma from pregnant women carrying a fetus with trisomy 21 would be higher than that from pregnant women carrying a euploid fetus. The analysis does not require targeting of fetal-specific sequences. It also does not require prior physical separation of fetal nucleic acid from maternal nucleic acid. It is also not necessary to distinguish or differentiate between fetal and maternal sequences after sequencing.

[0083] FIG. 3A shows the percentage of sequences mapped to chromosome 21 (percentage representation of chromosome 21) for each of the eight maternal plasma DNA samples. The percentage representation of chromosome 21 was significantly higher in maternal plasma of trisomy 21 pregnancies than that of euploid pregnancies. These data suggest that non-invasive prenatal diagnosis of fetal aneuploidy can be achieved by determining the percentage of aneuploid chromosomes compared to that of a reference group. Alternatively, over-representation of chromosome 21 can be detected by comparing the percentage representation of chromosome 21 obtained experimentally with the percentage representation of chromosome 21 expected for a euploid human genome. This can be done by masking or not masking repetitive regions in the human genome.

[0084] Of the eight pregnant women, five were carrying a male fetus each. Sequences mapping to the Y chromosome would be fetal-specific. The percentage of sequences mapping to the Y chromosome was used to calculate the partial fetal DNA concentration in the original maternal plasma samples. In addition, partial fetal DNA concentrations were also determined by using microfluidics digital PCR involving X-linked (ZFX) zinc finger proteins, Y-linked (ZFY) paralogous genes.

[0085] Figure 3B shows the correlation between the partial fetal DNA concentration estimated by the percentage representation of the Y chromosome by sequencing and that determined by ZFY / ZFX microfluidics digital PCR. There was a positive correlation between the partial fetal DNA concentrations in maternal plasma determined by these two methods. The correlation coefficient (r) was 0.917 in the Pearson correlation analysis.

[0086] The percentage of maternal plasma DNA sequences aligned to each of the 24 chromosomes (22 autosomes, plus chromosomes X and Y) for two representative cases is shown in Figure 4A. One pregnant woman was carrying a trisomy 21 fetus and the other was carrying a euploid fetus. The percentage representation of sequences mapping to chromosome 21 is higher in pregnant women carrying trisomy 21 fetuses compared to pregnant women carrying normal fetuses.

[0087] The percentage difference (%) per chromosome between the maternal plasma DNA samples in the above two cases is shown in Figure 4B. The percentage difference for a particular chromosome was calculated using the following formula: Percentage difference (%) = (P 21 -P E ) / P E x100% Here, P 21 is the percentage of plasma DNA sequences that aligned to a particular chromosome in pregnant women carrying a trisomy 21 fetus; and P Eis the percentage of plasma DNA sequences that aligned to a particular chromosome in pregnant women carrying euploid fetuses.

[0088] As shown in FIG. 4B, in the plasma of pregnant women carrying trisomy 21 sequences, there is an 11% over-representation of chromosome 21 when compared to pregnant women carrying euploid fetuses. For sequences aligned to other chromosomes, the difference between the two cases was within 5%. Because the percentage representation for chromosome 21 is increased in trisomy 21 samples compared to euploid maternal plasma samples, the difference (%) can be alternatively referred to as the degree of over-representation between chromosome 21 sequences. In addition to the difference (%) and absolute difference between the percentage representation of chromosome 21, the ratio of numbers from the test sample and the reference sample can also be calculated, and the reference can be calculated similarly, and will indicate the degree of over-representation of chromosome 21 in trisomy 21 compared to euploid samples.

[0089] For four pregnant women, each carrying a euploid fetus, an average of 1.3455% of their plasma DNA sequences were aligned to chromosome 21. For the four women carrying a trisomy 21 fetus, three of their fetuses were female. The percentage representation of chromosome 21 was calculated for each of these three cases. The difference (%) in the percentage representation of chromosome 21 for these three trisomy 21 cases from the average percentage representation of chromosome 21 obtained from the values ​​for the four euploid cases was determined as described above; i.e., the average of the four cases carrying euploid fetuses was used as the basis in this calculation. The fractional fetal DNA concentrations for these three male trisomy 21 cases were inferred from the percentage representation of each of the Y chromosome sequences.

[0090] The correlation between the degree of over-representation of chromosome 21 and the partial fetal DNA concentration is shown in FIG. 5. There was a significant positive correlation between the two parameters. The correlation coefficient (r) was 0.898 in the Pearson correlation analysis. These results indicate that the degree of over-representation of chromosome 21 in maternal plasma is associated with the partial concentration of fetal DNA in maternal plasma samples. Therefore, a cut-off value in the degree of over-representation of chromosome 21 sequences in relation to the partial fetal DNA concentration can be determined to identify pregnancies involving trisomy 21 fetuses.

[0091] The determination of the fractional concentration of fetal DNA in maternal plasma can also be made separately from the sequencing run. For example, the Y chromosome DNA concentration can be determined in advance using real-time PCR, microfluidic PCR, or mass spectrometry. For example, we have demonstrated in FIG. 3 that there is an excellent correlation between the fetal DNA concentration estimated based on the Y chromosome number generated during the sequencing run and the ZFY / ZFX ratio generated separately from the sequencing run. In fact, fetal DNA concentration can be determined using loci applicable in female fetuses other than the Y chromosome. For example, Chan et al. have shown that fetal-derived methylated RASSF1A sequences are detected in the plasma of pregnant women in the background of maternal-derived unmethylated RASSF1A sequences (Chan et al, Clin Chem 2006;52:2211-8). Thus, the fractional fetal DNA concentration can be determined by dividing the amount of methylated RASSF1A sequences by the total RASSF1A (methylated and unmethylated) sequences.

[0092] It is expected that in carrying out the present invention, maternal plasma is more preferable than maternal serum, because DNA is released from maternal blood cells during blood clotting.Therefore, if serum is used, it is expected that the partial concentration of fetal DNA will be lower in maternal plasma than in maternal serum.That is, if maternal serum is used, it is expected that more sequences need to be generated to diagnose fetal chromosomal aneuploidy compared with the plasma sample obtained from the same pregnant woman at the same time.

[0093] Yet another alternative method for determining the fractional concentration of fetal DNA would be through quantification of polymorphic differences between the pregnant woman and the fetus (Dhallan R, et al.2007 Lancet,369,474-481). An example of this method would be to target polymorphic sites where the pregnant woman is homozygous and the fetus is heterozygous. The amount of fetal-specific alleles can be compared to the amount of common alleles to determine the fractional concentration of fetal DNA.

[0094] In contrast to existing techniques for detecting chromosomal abnormalities, including comparative genomic hybridization, microarray comparative genomic hybridization, and quantitative real-time polymerase chain reaction, which detect and quantify one or more specific sequence(s), massively parallel sequencing does not rely on the detection or analysis of a set of pre-determined or defined DNA sequences. Random representative fragments of DNA molecules from a pool of samples were sequenced. The number of different sequenced tags aligned to various chromosomal regions was compared between samples that contained or did not contain the DNA species of interest. Chromosomal abnormalities would be revealed by differences in the number (or percentage) of sequences aligned to any given chromosomal region in the samples.

[0095] In another example, sequencing technology in cell-free DNA of plasma may be used to detect chromosomal abnormalities in plasma DNA for the detection of specific cancer. Different cancers have a set of typical chromosomal abnormalities. Alterations (amplification and deletion) in multiple chromosomal regions may be used. Thus, there will be amplified sequence ratios aligned to amplified regions and reduced sequence ratios aligned to reduced regions. The percentage representation per chromosome may be compared with the size for each corresponding chromosome in a reference genome, expressed as a percentage of the genome representation of any given chromosome relative to the whole genome. Direct comparison or comparison with a reference chromosome may be used as well.

[0096] II. Sequencing only a portion of the human genome In the experiment described in Example I above, maternal plasma DNA from each individual specimen was sequenced using only one flow cell. The number of sequenced tags generated from each of the specimens tested by the sequencing run is shown in Figure 6. T21 indicates a sample obtained from a pregnant woman associated with a trisomy 21 fetus.

[0097] From each sequenced maternal plasma DNA fragment, 36 bp were sequenced, so the number of nucleotides / base pairs sequenced from each specimen can be determined by multiplying 36 bp by the number of tags sequenced, and is similarly shown in Figure 6. With approximately 3 billion base pairs in the human genome, the amount of sequencing data generated from individual maternal plasma specimens, of which only a portion is shown, ranged from approximately 10% to 13%.

[0098] Furthermore, in this study, only uniquely mappable sequenced tags, designated U0, a nomenclature from the Eland software, were used to demonstrate the presence of over-representation in the amount of chromosome 21 sequences in maternal plasma samples from pregnant women each carrying a fetus with trisomy 21 as described in Example 1 above. As shown in Figure 6, the U0 sequences represent only a subset of all sequenced tags generated from each specimen, and further represent a smaller proportion of the human genome, approximately 2%. These data indicate that sequencing only a portion of the human genome sequences present in the tested samples is sufficient to achieve a diagnosis of fetal aneuploidy.

[0099] III. Determining the Number of Sequences Needed The results of sequencing plasma DNA from a pregnant woman carrying a euploid male fetus were used for this analysis. The number of sequenced tags located without mismatches in the reference human genome sequence was 1,990,000. Subsets of sequences were randomly selected from these 1,990,000 tags, and the percentage of sequences aligned to chromosome 21 was calculated within each subset. The number of sequences in the subsets was varied from 60,000 to 540,000 sequences. For the number of individual subsets, multiple subsets of the same number of sequenced tags were compiled by random selection of sequenced tags from the total pool until no other combinations were possible. The average percentage of sequences aligned to chromosome 21 and its standard deviation (SD) were then calculated from multiple subsets within each subset size. These data were compared between different subset sizes to determine the effect of subset size on the percentage distribution of sequences aligned to chromosome 21. The 5th and 95th percentiles of the percentages were then calculated according to the mean and SD.

[0100] When a pregnant woman is carrying a trisomy 21 fetus, sequenced tags aligned to chromosome 21 should be over-represented in maternal plasma due to an overdose of chromosome 21 from the fetus. The degree of over-representation depends on the fetal DNA percentage in the maternal plasma DNA sample according to the following equation: Per T21 =Per Eu ×(1+f / 2) Here, Per T21 shows the percentage of sequences that aligned to chromosome 21 in women with trisomy 21 fetuses; and Per Eu indicates the percentage of sequences that aligned to chromosome 21 in women with euploid fetuses; and f indicates the fetal DNA percentage in maternal plasma DNA.

[0101] As shown in Figure 7, the SD for the percentage of sequences aligned to chromosome 21 decreases with increasing number of sequences in each subset. Thus, the interval between the 5th and 95th percentiles decreases as the number of sequences in each subset increases. When the 5% and 95% intervals for euploid and trisomy 21 do not overlap, discrimination between the two case groups will be possible with greater than 95% accuracy.

[0102] As shown in Figure 7, the minimum subset size for distinguishing between trisomy 21 and euploid cases depends on the fetal DNA percentage. The minimum subset sizes for distinguishing between trisomy 21 and euploid cases were 120,000, 180,000, and 540,000 sequences for fetal DNA percentages of 20%, 10%, and 5%, respectively. That is, when the maternal plasma DNA sample contains 20% fetal DNA, the number of sequences required to be analyzed to determine whether the fetus has trisomy 21 would be 120,000. When the fetal DNA percentage falls to 5%, the number of sequences required to be analyzed would increase to 540,000.

[0103] Since data is generated using 36 base pair sequencing, 120,000, 180,000, and 540,000 sequences correspond to 0.14%, 0.22%, and 0.65% of the human genome, respectively.Since the low range fetal DNA concentration in maternal plasma obtained from early pregnancy is reported to be about 5% (Lo, YMD et al.1998 Am J Hum Genet 62,768-775), sequencing about 0.6% of the human genome may represent the minimum amount of sequencing required for the diagnosis of fetal chromosomal aneuploidy for any pregnant woman with at least 95% accuracy.

[0104] IV. Random Screening To demonstrate that the sequenced DNA fragments are randomly selected during the sequencing run, we obtained sequenced tags resulting from the eight maternal plasma samples analyzed in Example I. For each maternal plasma specimen, we determined the start position of each of the 36 bp sequenced tags that were uniquely aligned to chromosome 21 without mismatches relative to the reference human genome sequence, NCBI assembly 36. We then ranked the start position numbers for the pool of aligned sequenced tags from each specimen in ascending order. We performed a similar analysis for chromosome 22. For illustration, the top 10 start positions for chromosomes 21 and 22 for each maternal plasma specimen are shown in Figures 8A and 8B, respectively. As can be seen from these tables, the sequenced pools of DNA fragments were not identical between samples.

[0105] Any software components or functions described in this application may be implemented using a suitable computer language, such as Java, C++, or Perl, for example using conventional or object-oriented techniques. The software code may be stored as a series of instructions or commands in a suitable medium, including a computer-readable medium for storage and / or transmission, a random access memory (RAM), a read-only memory (ROM), a magnetic medium, such as a hard drive or a floppy disk, or an optical medium, such as a compact disk (CD) or a DVD (digital versatile disk), a flash memory, etc. The computer-readable medium may be any combination of such storage or transmission devices.

[0106] Such programs may also be encoded and transmitted using carrier signals adapted for transmission over wired, optical, and / or wireless networks conforming to various protocols, including the Internet. Thus, computer readable media of the present invention may be created using data signals encoded with such programs. Such program-encoded computer readable media may be packaged with a compatible device or provided separately from other devices (e.g., via Internet download). Any such computer readable media may be provided on or within a single program product (e.g., a hard drive or an entire computer system) and may reside on or within different computer program products within a system or network. The computer system may include a monitor, printer, or other suitable display to provide a user with any of the results mentioned herein.

[0107] An example of a computer system is shown in FIG. 9. The subsystems shown in FIG. 9 are interconnected through a system bus 975. Additional subsystems are shown, such as a printer 974, a keyboard 978, a fixed disk 979, a monitor 976 coupled to a display adapter 982, and others. Peripherals and input / output (I / O) devices coupled to an I / O controller 971 may be connected to the computer system by any number of means known in the art, such as a serial port 977. For example, the serial port 977 or an external interface 981 may be used to connect the computer system to a wide area network such as the Internet, a mouse input device, or a scanner. The interconnection via the system bus allows the central processor 973 to communicate with the individual subsystems and to execute instructions from the system memory 972 or fixed disk 979, and to exchange information between the subsystems. The system memory 972 and / or the fixed disk 979 may embody a computer-readable medium.

[0108] The above description of exemplary embodiments of the present invention is presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise form described, and many modifications and variations are possible in light of the above teachings. The embodiments have been chosen and described to best explain the principles of the invention and its practical application, with various modifications suitable for the particular use intended, so that others skilled in the art may best utilize the invention in various embodiments.

[0109] All applications, patents, and patent applications cited herein are hereby incorporated by reference in their entirety for all purposes.

Claims

1. 1. A method for determining whether a fetal chromosomal aneuploidy is present in a biological sample obtained from a female subject carrying a fetus, comprising: wherein the biological sample is maternal plasma and comprises cell-free nucleic acid molecules derived from the genome of the female subject and cell-free nucleic acid molecules derived from the genome of one of the fetuses; Steps below: obtaining said biological sample; randomly sequencing the cell-free nucleic acid molecules derived from the genome of the female subject and the cell-free nucleic acid molecules derived from the genome of the single fetus contained in the biological sample to obtain sequenced tags, wherein the sequenced tags include sequenced tags corresponding to the cell-free nucleic acid molecules derived from the genome of the female subject and sequenced tags corresponding to the cell-free nucleic acid molecules derived from the genome of the single fetus; aligning at least a portion of the sequenced tags to the human genome; determining a first amount of sequenced tags identified as aligning to a first chromosome of the human genome and not aligning to one or more second chromosomes of the human genome; determining a second amount of sequenced tags identified as aligning to the one or more second chromosomes and not aligning to the first chromosome; using said first amount and said second amount to determine a percentage representation of sequenced tags identified as aligning to said first chromosome, wherein said percentage representation is a measure of the proportion of cell free nucleic acid molecules in the biological sample that are derived from said first chromosome; comparing the percentage representation to one or more cutoff values; and determining whether a fetal chromosomal aneuploidy is present in said first chromosome; Including, wherein the first amount is the total number or total length of sequenced tags identified as aligning to the first chromosome and not aligning to the one or more second chromosomes; The second amount is the total number or total length of sequenced tags identified as aligning to the one or more second chromosomes and not aligning to the first chromosome; and The method, wherein the first chromosome is chromosome 21, chromosome 18, chromosome 13, an X chromosome, or a Y chromosome.

2. 2. The method of claim 1, wherein the sequenced tags identified as aligning to the first chromosome and not aligning to the one or more second chromosomes are less than 300 bp.

3. 2. The method of claim 1, wherein the sequenced tags identified as aligning to the first chromosome and not aligning to the one or more second chromosomes are less than 200 bp.

4. 2. The method of claim 1, wherein the sequenced tags identified as aligning to the first chromosome and not aligning to the one or more second chromosomes are less than 100 bp.

5. 10. The method of claim 1, wherein the cell-free nucleic acid molecules of the biological sample are enriched for sequences that are derived from at least one specific chromosome.

6. 10. The method of claim 1, wherein the cell-free nucleic acid molecules of the biological sample are enriched for sequences that are less than 300 bp.

7. 10. The method of claim 1, wherein the cell-free nucleic acid molecules of the biological sample are enriched for sequences that are less than 200 bp.

8. The method of claim 1, wherein the nucleic acid molecules sequenced by random sequencing are at least 0.1% of the human genome.

9. The method of claim 1, wherein the nucleic acid molecules sequenced by random sequencing represent at least 0.5% of the human genome.

10. The method of claim 1 , wherein the cutoff value is a reference value established in normal biological samples.

11. 1. A computer program for controlling a computing system to perform operations for performing prenatal diagnosis of fetal chromosomal aneuploidies in a biological sample obtained from a female subject carrying a fetus, the computing system comprising an input device, one or more processors and an output device, the biological sample being maternal plasma and comprising cell-free nucleic acid molecules derived from the genome of the female subject and cell-free nucleic acid molecules derived from the genome of the fetus, the operations comprising the steps of: obtaining, by the input device, sequenced tag data resulting from random sequencing of a portion of the cell-free nucleic acid molecules contained in the biological sample of the female subject; aligning, by the processor, at least a portion of the sequenced tags to the human genome; determining a first amount of the first chromosome from sequences identified by the processor as aligning to a first chromosome of the human genome and not aligning to one or more second chromosomes of the human genome; determining a second amount of the one or more second chromosomes from sequences identified by the processor as aligning to the one or more second chromosomes and not aligning to the first chromosome; using the first amount and the second amount to determine a percentage representation of sequences identified by the processor as aligning to the first chromosome, wherein the percentage representation is a measure of the proportion of cell free nucleic acid molecules in the biological sample that are derived from the first chromosome; comparing, by the processor, the percentage representation to one or more cutoff values; and outputting a result of the comparison by the output device; wherein based on a result of the comparison, it is determined whether a fetal chromosomal aneuploidy is present in the first chromosome; wherein the first amount is the total number or total length of sequenced tags identified as aligning to the first chromosome and not aligning to the one or more second chromosomes; The second amount is the total number or total length of sequenced tags identified as aligning to the one or more second chromosomes and not aligning to the first chromosome; and The first chromosome is chromosome 21, chromosome 18, chromosome 13, chromosome X, or chromosome Y. The computer program.

Images