Middha, Sumit

Abstract Several platforms for genome sequencing and methods for exon enrichment are now available and continue to evolve quickly. In this study, we compared the performance of whole exome sequencing using two instruments, Illumina's GAIIx and HiSeq 2000, comparing single versus duplicate runs. Eight runs of eight samples each were performed with genomic DNA derived from whole blood and sample libraries prepared with the Agilent SureSelect capture array. One run was conducted in duplicate on the GAIIx, and, for some analyses, the combined data were used. We used TREAT (Targeted RE-sequencing and Annotation Tool), developed in-house for data analysis, including sequence alignment (MAQ, BWA), local re-alignment (GATK), variant calling (utilizing MAQ), annotation (SIFT and Seattle Seq), and visualization. Overall, duplicate runs on the GAIIx platform (compared to a single run) increased the number of total reads by approximately two-fold: from ∼66 M to ∼130 M. For both single and duplicate runs, ∼90% of reads mapped to the reference sequence, and ∼55% of reads mapped on-target. Of importance, the percent coverage of the target region increased substantially with duplicate runs; ∼72% of the target region was covered at 40-fold or more when run in duplicate compared to only ∼52% when run once. At 30x, 20x, and 10x, there was a 1.3- (62% v. 78%), 1.2- (72% v. 84%), and 1.1-fold (84% v. 91%) increase respectively, in the percent coverage of the target region for samples run twice. A similar number of filtered, on-target SNPs per sample (∼24 K) was found for both single and duplicate run analyses. However, a 1.3-fold increase in the number of on-target indels was seen in the duplicate run (∼1,200) compared to the single run (∼900). When the combined data for samples run twice on a GAIIx were compared to the samples run using a HiSeq, the total number of reads were similar (∼130 M), although several of the samples had substantially more reads on the HiSeq platform. Samples run on a HiSeq had an increase in the percent coverage of the target region at 40x (72% v. 81%), 30x (78% v. 86%), 20x (84% v. 91%) and 10x (91% v. 95%) compared to the two-run approach on a GAIIx. A similar number of filtered on-target SNPs (∼24 K/sample) and indels (∼1200) was found with the HiSeq compared to two runs on a GAIIx. In summary, our results demonstrate increases in the total number of reads and in overall coverage with the HiSeq 2000 instrument. The total number of SNPs and indel that mapped on target for the version 1 Agilent SureSelect capture array is ∼25 K and 1,200 per sample, respectively. As capture is of variable efficiency for individual runs, higher average coverage of the target region is necessary for sufficient coverage of poorly captured regions. Additional analyses comparing Agilent SureSelect versions 1 and 2 arrays are underway. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 2742. doi:10.1158/1538-7445.AM2011-2742

AbstractThe biological role of vitamin D receptors (VDR), which are abundantly expressed in developing zebrafish (Danio rerio) as early as 48 h after fertilization, and before the development of a mineralized skeleton and mature intestine and kidney, is unknown. We probed the role of VDR in developing zebrafish biology by examining changes in expression of RNA by whole transcriptome shotgun sequencing (RNA-seq) in fish treated with picomolar concentrations of the VDR ligand and hormonal form of vitamin D3, 1α,25-dihydroxyvitamin D3 [1α,25(OH)2D3)].We observed significant changes in RNAs of transcription factors, leptin, peptide hormones, and RNAs encoding proteins of fatty acid, amino acid, xenobiotic metabolism, receptor-activator of NFκB ligand (RANKL), and calcitonin-like ligand receptor pathways. Early highly restricted, and subsequent massive changes in more than 10% of expressed cellular RNA were observed. At days post fertilization (dpf) 2 [24 h 1α,25(OH)2D3-treatment], only four RNAs were differentially expressed (hormone vs. vehicle). On dpf 4 (72 h treatment), 77 RNAs; on dpf 6 (120 h treatment) 1039 RNAs; and on dpf 7 (144 h treatment), 2407 RNAs were differentially expressed in response to 1α,25(OH)2D3. Fewer RNAs (n = 481) were altered in dpf 7 larvae treated for 24 h with 1α,25(OH)2D3vs. those treated with hormone for 144 h. At dpf 7, in 1α,25(OH)2D3-treated larvae, pharyngeal cartilage was larger and mineralization was greater. Changes in expression of RNAs for transcription factors, peptide hormones, and RNAs encoding proteins integral to fatty acid, amino acid, leptin, calcitonin-like ligand receptor, RANKL, and xenobiotic metabolism pathways, demonstrate heretofore unrecognized mechanisms by which 1α,25(OH)2D3 functions in vivo in developing eukaryotes.

Abstract Summary: Reduced representation bisulfite sequencing (RRBS) is a cost-effective approach for genome-wide methylation pattern profiling. Analyzing RRBS sequencing data is challenging and specialized alignment/mapping programs are needed. Although such programs have been developed, a comprehensive solution that provides researchers with good quality and analyzable data is still lacking. To address this need, we have developed a Streamlined Analysis and Annotation Pipeline for RRBS data (SAAP-RRBS) that integrates read quality assessment/clean-up, alignment, methylation data extraction, annotation, reporting and visualization. This package facilitates a rapid transition from sequencing reads to a fully annotated CpG methylation report to biological interpretation. Availability and implementation: SAAP-RRBS is freely available to non-commercial users at the web site http://ndc.mayo.edu/mayo/research/biostat/stand-alone-packages.cfm. Contact: baheti.saurabh@mayo.edu or sun.zhifu@mayo.edu Supplementary Information: Supplementary data are available at Bioinformatics online.

Abstract The longitudinal evolution of a myeloma genome from diagnosis to plasma cell leukemia has not previously been reported. We used whole-genome sequencing (WGS) on 4 purified tumor samples and patient germline DNA drawn over a 5-year period in a t(4;14) multiple myeloma patient. Tumor samples were acquired at diagnosis, first relapse, second relapse, and end-stage secondary plasma cell leukemia (sPCL). In addition to the t(4;14), all tumor time points also shared 10 common single-nucleotide variants (SNVs) on WGS comprising shared initiating events. Interestingly, we observed genomic sequence variants that waxed and waned with time in progressive tumors, suggesting the presence of multiple independent, yet related, clones at diagnosis that rose and fell in dominance. Five newly acquired SNVs, including truncating mutations of RB1 and ZKSCAN3, were observed only in the final sPCL sample suggesting leukemic transformation events. This longitudinal WGS characterization of the natural history of a high-risk myeloma patient demonstrated tumor heterogeneity at diagnosis with shifting dominance of tumor clones over time and has also identified potential mutations contributing to myelomagenesis as well as transformation from myeloma to overt extramedullary disease such as sPCL.

BACKGROUND Next-generation DNA sequencing (NGS) techniques have the potential to revolutionize molecular diagnostics; however, a thorough evaluation of these technologies is necessary to ensure their performance meets or exceeds that of current clinical sequencing methods. METHODS We evaluated the NimbleGen Sequence Capture 385K Human Custom Arrays for enrichment of 22 genes. We sequenced each sample on both the Roche 454 Genome Sequencer FLX (GS-FLX) and the Illumina Genome Analyzer II (GAII) to compare platform performance. RESULTS Although the sequence capture method allowed us to rapidly develop a large number of sequencing assays, we encountered difficulty enriching G+C–rich regions. Although a high proportion of reads consistently mapped outside of the targeted regions, >80% of targeted bases for the GAII and >30% of bases for the GS-FLX were covered by a read depth of ≥20, and >90% of bases for the GAII and >80% of bases for the GS-FLX were covered by a read depth of ≥5. We observed discrepancies among sequence variants identified by the different platforms. CONCLUSIONS Although oligonucleotide arrays are quick and easy to develop, some problematic regions may evade capture, necessitating sequential redesigning for complete optimization. Neither sequencing technology was able to detect every variant identified by Sanger sequencing because of well-known drawbacks of the NGS technologies. The rapidly decreasing error rates and costs of these technologies, however, coupled with advancing bioinformatic capabilities, make them an attractive option for molecular diagnostics in the very near future.