Accurate whole-genome sequencing and haplotyping from 10 to 20 human cells.

Recent advances in whole-genome sequencing have brought the vision of personal genomics and genomic medicine closer to reality. However, current methods lack clinical accuracy and the ability to describe the context (haplotypes) in which genome variants co-occur in a cost-effective manner. Here we describe a low-cost DNA sequencing and haplotyping process, long fragment read (LFR) technology, which is similar to sequencing long single DNA molecules without cloning or separation of metaphase chromosomes. In this study, ten LFR libraries were made using only ∼100 picograms of human DNA per sample. Up to 97% of the heterozygous single nucleotide variants were assembled into long haplotype contigs. Removal of false positive single nucleotide variants not phased by multiple LFR haplotypes resulted in a final genome error rate of 1 in 10 megabases. Cost-effective and accurate genome sequencing and haplotyping from 10-20 human cells, as demonstrated here, will enable comprehensive genetic studies and diverse clinical applications.

The mutation spectrum revealed by paired genome sequences from a lung cancer patient.

Lung cancer is the leading cause of cancer-related mortality worldwide, with non-small-cell lung carcinomas in smokers being the predominant form of the disease. Although previous studies have identified important common somatic mutations in lung cancers, they have primarily focused on a limited set of genes and have thus provided a constrained view of the mutational spectrum. Recent cancer sequencing efforts have used next-generation sequencing technologies to provide a genome-wide view of mutations in leukaemia, breast cancer and cancer cell lines. Here we present the complete sequences of a primary lung tumour (60x coverage) and adjacent normal tissue (46x). Comparing the two genomes, we identify a wide variety of somatic variations, including >50,000 high-confidence single nucleotide variants. We validated 530 somatic single nucleotide variants in this tumour, including one in the KRAS proto-oncogene and 391 others in coding regions, as well as 43 large-scale structural variations. These constitute a large set of new somatic mutations and yield an estimated 17.7 per megabase genome-wide somatic mutation rate. Notably, we observe a distinct pattern of selection against mutations within expressed genes compared to non-expressed genes and in promoter regions up to 5 kilobases upstream of all protein-coding genes. Furthermore, we observe a higher rate of amino acid-changing mutations in kinase genes. We present a comprehensive view of somatic alterations in a single lung tumour, and provide the first evidence, to our knowledge, of distinct selective pressures present within the tumour environment.

Noninvasive prenatal detection and selective analysis of cell-free DNA obtained from maternal blood: evaluation for trisomy 21 and trisomy 18.

OBJECTIVE: We sought to develop a novel biochemical assay and algorithm for the prenatal evaluation of risk for fetal trisomy 21 (T21) and trisomy 18 (T18) using cell-free DNA obtained from maternal blood. STUDY DESIGN: We assayed cell-free DNA from a training set and a blinded validation set of pregnant women, comprising 250 disomy, 72 T21, and 16 T18 pregnancies. We used digital analysis of selected regions in combination with a novel algorithm, fetal-fraction optimized risk of trisomy evaluation (FORTE), to determine trisomy risk for each subject. RESULTS: In all, 163/171 subjects in the training set passed quality control criteria. Using a Z statistic, 35/35 T21 cases and 7/7 T18 cases had Z statistic >3 and 120/121 disomic cases had Z statistic <3. FORTE produced an individualized trisomy risk score for each subject, and correctly discriminated all T21 and T18 cases from disomic cases. All 167 subjects in the blinded validation set passed quality control and FORTE performance matched that observed in the training set correctly discriminating 36/36 T21 cases and 8/8 T18 cases from 123/123 disomic cases. CONCLUSION: Digital analysis of selected regions and FORTE enable accurate, scalable noninvasive fetal aneuploidy detection.

Non-Invasive Chromosomal Evaluation (NICE) Study: results of a multicenter prospective cohort study for detection of fetal trisomy 21 and trisomy 18.

OBJECTIVE: We sought to evaluate performance of a noninvasive prenatal test for fetal trisomy 21 (T21) and trisomy 18 (T18). STUDY DESIGN: A multicenter cohort study was performed whereby cell-free DNA from maternal plasma was analyzed. Chromosome-selective sequencing on chromosomes 21 and 18 was performed with reporting of an aneuploidy risk (High Risk or Low Risk) for each subject. RESULTS: Of the 81 T21 cases, all were classified as High Risk for T21 and there was 1 false-positive result among the 2888 normal cases, for a sensitivity of 100% (95% confidence interval [CI], 95.5-100%) and a false-positive rate of 0.03% (95% CI, 0.002-0.20%). Of the 38 T18 cases, 37 were classified as High Risk and there were 2 false-positive results among the 2888 normal cases, for a sensitivity of 97.4% (95% CI, 86.5-99.9%) and a false-positive rate of 0.07% (95% CI, 0.02-0.25%). CONCLUSION: Chromosome-selective sequencing of cell-free DNA and application of an individualized risk algorithm is effective in the detection of fetal T21 and T18.

Selective analysis of cell-free DNA in maternal blood for evaluation of fetal trisomy.

OBJECTIVE: To develop a novel prenatal assay based on selective analysis of cell-free DNA in maternal blood for evaluation of fetal Trisomy 21 (T21) and Trisomy 18 (T18). METHODS: Two hundred ninety-eight pregnancies, including 39 T21 and seven T18 confirmed fetal aneuploidies, were analyzed using a novel, highly multiplexed assay, termed digital analysis of selected regions (DANSR™). Cell-free DNA from maternal blood samples was analyzed using DANSR assays for loci on chromosomes 21 and 18. Products from 96 separate patients were pooled and sequenced together. A standard Z-test of chromosomal proportions was used to distinguish aneuploid samples from average-risk pregnancy samples. DANSR aneuploidy discrimination was evaluated at various sequence depths. RESULTS: At the lowest sequencing depth, corresponding to 204,000 sequencing counts per sample, average-risk cases where distinguished from T21 and T18 cases, with Z statistics for all cases exceeding 3.6. Increasing the sequencing depth to 410,000 counts per sample substantially improved separation of aneuploid and average-risk cases. A further increase to 620,000 counts per sample resulted in only marginal improvement. This depth of sequencing represents less than 5% of that required by massively parallel shotgun sequencing approaches. CONCLUSION: Digital analysis of selected regions enables highly accurate, cost efficient, and scalable noninvasive fetal aneuploidy assessment.

Using the transcriptome to annotate the genome.

A remaining challenge for the human genome project involves the identification and annotation of expressed genes. The public and private sequencing efforts have identified approximately 15,000 sequences that meet stringent criteria for genes, such as correspondence with known genes from humans or other species, and have made another approximately 10,000-20,000 gene predictions of lower confidence, supported by various types of in silico evidence, including homology studies, domain searches, and ab initio gene predictions. These computational methods have limitations, both because they are unable to identify a significant fraction of genes and exons and because they are unable to provide definitive evidence about whether a hypothetical gene is actually expressed. As the in silico approaches identified a smaller number of genes than anticipated, we wondered whether high-throughput experimental analyses could be used to provide evidence for the expression of hypothetical genes and to reveal previously undiscovered genes. We describe here the development of such a method--called long serial analysis of gene expression (LongSAGE), an adaption of the original SAGE approach--that can be used to rapidly identify novel genes and exons.

Human genome sequencing using unchained base reads on self-assembling DNA nanoarrays.

Genome sequencing of large numbers of individuals promises to advance the understanding, treatment, and prevention of human diseases, among other applications. We describe a genome sequencing platform that achieves efficient imaging and low reagent consumption with combinatorial probe anchor ligation chemistry to independently assay each base from patterned nanoarrays of self-assembling DNA nanoballs. We sequenced three human genomes with this platform, generating an average of 45- to 87-fold coverage per genome and identifying 3.2 to 4.5 million sequence variants per genome. Validation of one genome data set demonstrates a sequence accuracy of about 1 false variant per 100 kilobases. The high accuracy, affordable cost of $4400 for sequencing consumables, and scalability of this platform enable complete human genome sequencing for the detection of rare variants in large-scale genetic studies.

The genomic landscapes of human breast and colorectal cancers.

Human cancer is caused by the accumulation of mutations in oncogenes and tumor suppressor genes. To catalog the genetic changes that occur during tumorigenesis, we isolated DNA from 11 breast and 11 colorectal tumors and determined the sequences of the genes in the Reference Sequence database in these samples. Based on analysis of exons representing 20,857 transcripts from 18,191 genes, we conclude that the genomic landscapes of breast and colorectal cancers are composed of a handful of commonly mutated gene "mountains" and a much larger number of gene "hills" that are mutated at low frequency. We describe statistical and bioinformatic tools that may help identify mutations with a role in tumorigenesis. These results have implications for understanding the nature and heterogeneity of human cancers and for using personal genomics for tumor diagnosis and therapy.