We developed a massive-scale RNA sequencing protocol, short quantitative random RNA libraries or SQRL, to survey the complexity, dynamics and sequence content of transcriptomes in a near-complete fashion. This method generates directional, random-primed, linear cDNA libraries that are optimized for next-generation short-tag sequencing. We surveyed the poly(A)(+) transcriptomes of undifferentiated mouse embryonic stem cells (ESCs) and embryoid bodies (EBs) at an unprecedented depth (10 Gb), using the Applied Biosystems SOLiD technology. These libraries capture the genomic landscape of expression, state-specific expression, single-nucleotide polymorphisms (SNPs), the transcriptional activity of repeat elements, and both known and new alternative splicing events. We investigated the impact of transcriptional complexity on current models of key signaling pathways controlling ESC pluripotency and differentiation, highlighting how SQRL can be used to characterize transcriptome content and dynamics in a quantitative and reproducible manner, and suggesting that our understanding of transcriptional complexity is far from complete.

Using the massively parallel technique of sequencing by oligonucleotide ligation and detection (SOLiD; Applied Biosystems), we have assessed the in vivo positions of more than 44 million putative nucleosome cores in the multicellular genetic model organism Caenorhabditis elegans. These analyses provide a global view of the chromatin architecture of a multicellular animal at extremely high density and resolution. While we observe some degree of reproducible positioning throughout the genome in our mixed stage population of animals, we note that the major chromatin feature in the worm is a diversity of allowed nucleosome positions at the vast majority of individual loci. While absolute positioning of nucleosomes can vary substantially, relative positioning of nucleosomes (in a repeated array structure likely to be maintained at least in part by steric constraints) appears to be a significant property of chromatin structure. The high density of nucleosomal reads enabled a substantial extension of previous analysis describing the usage of individual oligonucleotide sequences along the span of the nucleosome core and linker. We release this data set, via the UCSC Genome Browser, as a resource for the high-resolution analysis of chromatin conformation and DNA accessibility at individual loci within the C. elegans genome.[Supplemental material is available online at www.genome.org. SOLiD raw sequencing data from this study have been submitted to the Short Read Archive at NCBI under accession no. SRA001023.]The regulation of genetic information within eukaryotic cells involves a high degree of specificity both in the availability of individual DNA-binding factors in individual cells, and in availability in the genome of DNA sequences that are their potential targets. Modulating the accessibility of individual DNA sequences are many complex interactions, the most prevalent of which are the interactions between histone octamers and DNA in compacted chromosomes. Each histone core interacts with 147 bp of DNA, which coil 1.7 times around the histone octamer (Luger et al. 1997;Davey et al. 2002) to form the basic unit of chromatin structure, the nucleosome. Since the first description over three decades ago (Kornberg 1974), the nucleosome and its role in gene regulation has been the subject of intensive study and speculation.As we have an increasingly detailed functional view of the genome, the tools of high-throughput molecular characterization have been used to begin obtaining a genome-wide description of nucleosome positions (Satchwell et al. 1986;Yuan et al. 2005;Johnson et al. 2006;Albert et al. 2007;Lee et al. 2007;Peckham et al. 2007;Schones et al. 2008;Shivaswamy et al. 2008). These data have in turn been used in attempts to reveal nucleosome positioning signals in DNA sequence (Satchwell et al. 1986;Ioshikhes et al. 1996;Segal et al. 2006;Yuan and Liu 2008). Although of great interest and value, sequence-based predictions of nucleosome position have been limited to date in their accuracy and resolution.In the litera...

SUMMARY To define the mutation spectrum in non-Down syndrome acute megkaryoblastic leukemia (non-DS-AMKL), we performed transcriptome sequencing on diagnostic blasts from 14 pediatric patients and validated our findings in a recurrency/validation cohort consisting of 34 pediatric and 28 adult AMKL leukemia samples. Our analysis identified a cryptic chromosome 16 inversion [inv(16)(p13.3q24.3)] in 27% of pediatric cases, which encodes a CBFA2T3-GLIS2 fusion protein. Expression of CBFA2T3-GLIS2 in Drosophila and murine hematopoietic cells induced bone morphogenic protein (BMP) signaling, and resulted in a marked increase in the self-renewal capacity of hematopoietic progenitors. These data suggest that expression of CBFA2T3-GLIS2 directly contributes to leukemogenesis.

We describe the genome sequencing of an anonymous individual of African origin using a novel ligation-based sequencing assay that enables a unique form of error correction that improves the raw accuracy of the aligned reads to >99.9%, allowing us to accurately call SNPs with as few as two reads per allele. We collected several billion mate-paired reads yielding ;183 haploid coverage of aligned sequence and close to 3003 clone coverage. Over 98% of the reference genome is covered with at least one uniquely placed read, and 99.65% is spanned by at least one uniquely placed matepaired clone. We identify over 3.8 million SNPs, 19% of which are novel. Mate-paired data are used to physically resolve haplotype phases of nearly two-thirds of the genotypes obtained and produce phased segments of up to 215 kb. We detect 226,529 intra-read indels, 5590 indels between mate-paired reads, 91 inversions, and four gene fusions. We use a novel approach for detecting indels between mate-paired reads that are smaller than the standard deviation of the insert size of the library and discover deletions in common with those detected with our intra-read approach. Dozens of mutations previously described in OMIM and hundreds of nonsynonymous single-nucleotide and structural variants in genes previously implicated in disease are identified in this individual. There is more genetic variation in the human genome still to be uncovered, and we provide guidance for future surveys in populations and cancer biopsies.[Supplemental material is available online at

The Decapentaplegic and Notch signaling pathways are thought to direct regional specification in the Drosophila eye-antennal epithelium by controlling the expression of selector genes for the eye (Eyeless/Pax6, Eyes absent) and/or antenna (Distal-less). Here, we investigate the function of these signaling pathways in this process. We find that organ primordia formation is indeed controlled at the level of Decapentaplegic expression but critical steps in regional specification occur earlier than previously proposed. Contrary to previous findings, Notch does not specify eye field identity by promoting Eyeless expression but it influences eye primordium formation through its control of proliferation. Our analysis of Notch function reveals an important connection between proliferation, field size, and regional specification. We propose that field size modulates the interaction between the Decapentaplegic and Wingless pathways, thereby linking proliferation and patterning in eye primordium development.

Obtaining high-quality sequence continuity of complex regions of recent segmental duplication remains one of the major challenges of finishing genome assemblies. In the human and mouse genomes, this was achieved by targeting large-insert clones using costly and laborious capillary-based sequencing approaches. Sanger shotgun sequencing of clone inserts, however, has now been largely abandoned, leaving most of these regions unresolved in newer genome assemblies generated primarily by next-generation sequencing hybrid approaches. Here we show that it is possible to resolve regions that are complex in a genome-wide context but simple in isolation for a fraction of the time and cost of traditional methods using long-read single molecule, real-time (SMRT) sequencing and assembly technology from Pacific Biosciences (PacBio). We sequenced and assembled BAC clones corresponding to a 1.3-Mbp complex region of chromosome 17q21.31, demonstrating 99.994% identity to Sanger assemblies of the same clones. We targeted 44 differences using Illumina sequencing and find that PacBio and Sanger assemblies share a comparable number of validated variants, albeit with different sequence context biases. Finally, we targeted a poorly assembled 766-kbp duplicated region of the chimpanzee genome and resolved the structure and organization for a fraction of the cost and time of traditional finishing approaches. Our data suggest a straightforward path for upgrading genomes to a higher quality finished state.

Forward genetic mutational studies, adaptive evolution, and phenotypic screening are powerful tools for creating new variant organisms with desirable traits. However, mutations generated in the process cannot be easily identified with traditional genetic tools. We show that new high-throughput, massively parallel sequencing technologies can completely and accurately characterize a mutant genome relative to a previously sequenced parental (reference) strain. We studied a mutant strain of Pichia stipitis, a yeast capable of converting xylose to ethanol. This unusually efficient mutant strain was developed through repeated rounds of chemical mutagenesis, strain selection, transformation, and genetic manipulation over a period of seven years. We resequenced this strain on three different sequencing platforms. Surprisingly, we found fewer than a dozen mutations in open reading frames. All three sequencing technologies were able to identify each single nucleotide mutation given at least 10-15-fold nominal sequence coverage. Our results show that detecting mutations in evolved and engineered organisms is rapid and cost-effective at the whole-genome level using new sequencing technologies. Identification of specific mutations in strains with altered phenotypes will add insight into specific gene functions and guide further metabolic engineering efforts.[Supplemental material is available online at www.genome.org. Complete data sets are available at the NCBI Short Read Archive under accession no. SRA 001158 (ftp://ftp.ncbi.nih.gov/pub/TraceDB/ShortRead).]Pichia stipitis (Pignal) is a haploid yeast related to endosymbionts of beetles that degrade rotting wood (Suh et al. 2003). It is an important organism for bioenergy production from lignocellulosic materials because of its high capacity to ferment xylose and cellobiose to ethanol (Parekh et al. 1988). We previously sequenced the reference strain, Pichia stipitis CBS-6054, resulting in a completely characterized genome of eight chromosomes totaling 15.4 Mb of sequence (Jeffries et al. 2007). This strain has been subjected to chemical mutagenesis, phenotypic selection, genetic engineering, and adaptive evolution in order to develop strains improved for ethanol production. Chemical mutagenesis and selection resulted in small improvements in ethanol production attributable in part to carbon catabolite derepression (Supplemental Fig. 1; Methods). Disruption of CYC1 (cyctochrome c, isoform 1) to create strain Shi21 increased the specific ethanol production rate by 50% and the ethanol yield by 10%; however, the nature of additional mutational events leading to this phenotype was uncharacterized.Traditional methods for identifying mutations are laborand time-intensive, so we tested the ability of next-generation sequencing technologies to determine the differences in this improved strain's entire genome relative to the reference strain. We generated high-coverage, whole-genome data sets using single fragment end reads from three next-generation sequencing platforms: 454 Life Sc...

During eye development, the selector factors of the Eyeless/Pax6 or Retinal Determination (RD) network control specification of organ-type whereas the bHLH-type proneural factor Atonal drives neurogenesis. Although significant progress has been made in dissecting the acquisition of 'eye identity' at the transcriptional level, the molecular mechanisms underlying the progression from neuronal progenitor to differentiating neuron remain unclear. A recently proposed model for the integration of organ specification and neurogenesis hypothesizes that atonal expression in the eye is RD-network-independent and that Eyeless works in parallel or downstream of atonal to modify the neurogenetic program. We show here that distinct cis-regulatory elements control atonal expression specifically in the eye and that the RD factors Eyeless and Sine oculis function as direct regulators. We find that these transcription factors interact in vitro and provide indirect evidence that this interaction may be required in vivo. The subordination of neurogenesis to the RD pathway in the eye provides a direct mechanism for the coordination of neurogenesis and tissue specification during sensory organ formation.